title: A projection method for approximating fixed points of quasi nonexpansive mappings without the usual demiclosedness condition

abstract: We introduce and analyze an abstract algorithm that aims to find the projection onto a closed convex subset of a Hilbert space. When specialized to the fixed point set of a quasi nonexpansive mapping, the required sufficient condition (termed "fixed-point closed") is less restrictive than the usual conditions based on the demiclosedness principle. A concrete example of a subgradient projector is presented which illustrates the applicability of this generalization.


title: Tensor products of tautological bundles under the Bridgeland-King-Reid-Haiman equivalence

abstract: We give a description of the image of tensor products of tautological bundles on Hilbert schemes of points on surfaces under the Bridgeland-King-Reid-Haiman equivalence. Using this, some new formulas for cohomological invariants of these bundles are obtained. In particular, we give formulas for the Euler characteristic of arbitrary tensor products on the Hilbert scheme of two points and of triple tensor products in general.


title: Probing W+W-gamma Production and Anomalous Quartic Gauge Boson Couplings at the CERN LHC

abstract: Triple gauge boson associated production at the LHC serves as an interesting channel to test the robustness of the Standard Model. Any deviation from its SM prediction may indicate possible existence of relevant new physics, e.g., anomalous quartic gauge boson couplings. In this paper, a Monte-Carlo feasibility study of measuring WWA production with pure leptonic decays and probing anomalous quartic gauge-boson (e.g., WWAA) couplings, is presented in detail for the first time, with parton shower and detector simulation effects taken into account. Our results show that at the sqrt{s} = 14 TeV LHC with an integrated luminosity of 100 (30) fb-1, one can reach a significance of 9 (5) sigma to observe the SM WWA production, and can constrain at the 95% CL the anomalous WWAA coupling parameters, e.g., a_{0,c}^W/\Lambda^2 (see Ref.[15] for their definitions), at 10^{-5} GeV^{-2}, respectively.


title: Randomized Dimension Reduction on Massive Data

abstract: Scalability of statistical estimators is of increasing importance in modern applications and dimension reduction is often used to extract relevant information from data. A variety of popular dimension reduction approaches can be framed as symmetric generalized eigendecomposition problems. In this paper we outline how taking into account the low rank structure assumption implicit in these dimension reduction approaches provides both computational and statistical advantages. We adapt recent randomized low-rank approximation algorithms to provide efficient solutions to three dimension reduction methods: Principal Component Analysis (PCA), Sliced Inverse Regression (SIR), and Localized Sliced Inverse Regression (LSIR). A key observation in this paper is that randomization serves a dual role, improving both computational and statistical performance. This point is highlighted in our experiments on real and simulated data.


title: Hybrid Behaviour of Markov Population Models

abstract: We investigate the behaviour of population models written in Stochastic Concurrent Constraint Programming (sCCP), a stochastic extension of Concurrent Constraint Programming. In particular, we focus on models from which we can define a semantics of sCCP both in terms of Continuous Time Markov Chains (CTMC) and in terms of Stochastic Hybrid Systems, in which some populations are approximated continuously, while others are kept discrete. We will prove the correctness of the hybrid semantics from the point of view of the limiting behaviour of a sequence of models for increasing population size. More specifically, we prove that, under suitable regularity conditions, the sequence of CTMC constructed from sCCP programs for increasing population size converges to the hybrid system constructed by means of the hybrid semantics. We investigate in particular what happens for sCCP models in which some transitions are guarded by boolean predicates or in the presence of instantaneous transitions.


title: Heat equation approach to geometric changes of the torus Laughlin-state

abstract: We study the second quantized -or guiding center- description of the torus Laughlin state. Our main focus is the change of the guiding center degrees of freedom with the torus geometry, which we show to be generated by a two-body operator. We demonstrate that this operator can be used to evolve the full torus Laughlin state at given modular parameter \tau\ from its simple (Slater-determinant) thin torus limit, thus giving rise to a new presentation of the torus Laughlin state in terms of its "root partition" and an exponential of a two-body operator. This operator therefore generates in particular the adiabatic evolution between Laughlin states on regular tori and the quasi-one-dimensional thin torus limit. We make contact with the recently introduced notion of a "Hall viscosity" for fractional quantum Hall states, to which our two-body operator is naturally related, and which serves as a demonstration of our method to generate the Laughlin state on the torus.


title: Hard-disk equation of state: First-order liquid-hexatic transition in two dimensions with three simulation methods

abstract: We report large-scale computer simulations of the hard-disk system at high densities in the region of the melting transition. Our simulations reproduce the equation of state, previously obtained using the event-chain Monte Carlo algorithm, with a massively parallel implementation of the local Monte Carlo method and with event-driven molecular dynamics. We analyze the relative performance of these simulation methods to sample configuration space and approach equilibrium. Our results confirm the first-order nature of the melting phase transition in hard disks. Phase coexistence is visualized for individual configurations via the orientational order parameter field. The analysis of positional order confirms the existence of the hexatic phase.


title: Massively parallel Monte Carlo for many-particle simulations on GPUs

abstract: Current trends in parallel processors call for the design of efficient massively parallel algorithms for scientific computing. Parallel algorithms for Monte Carlo simulations of thermodynamic ensembles of particles have received little attention because of the inherent serial nature of the statistical sampling. In this paper, we present a massively parallel method that obeys detailed balance and implement it for a system of hard disks on the GPU. We reproduce results of serial high-precision Monte Carlo runs to verify the method. This is a good test case because the hard disk equation of state over the range where the liquid transforms into the solid is particularly sensitive to small deviations away from the balance conditions. On a Tesla K20, our GPU implementation executes over one billion trial moves per second, which is 148 times faster than on a single Intel Xeon E5540 CPU core, enables 27 times better performance per dollar, and cuts energy usage by a factor of 13. With this improved performance we are able to calculate the equation of state for systems of up to one million hard disks. These large system sizes are required in order to probe the nature of the melting transition, which has been debated for the last forty years. In this paper we present the details of our computational method, and discuss the thermodynamics of hard disks separately in a companion paper.

categories: physics.comp-ph cond-mat.mtrl-sci cond-mat.soft cond-mat.stat-mech


title: Deformation theory and rational homotopy type

abstract: We regard the classification of rational homotopy types as a problem in algebraic deformation theory: any space with given cohomology is a perturbation, or deformation, of the "formal" space with that cohomology. The classifying space is then a "moduli" space --- a certain quotient of an algebraic variety of perturbations. The description we give of this moduli space links it with corresponding structures in homotopy theory, especially the classification of fibres spaces with fixed fibre F in terms of homotopy classes of maps of the base B into a classifying space constructed from the monoid of homotopy equivalences of F to itself. We adopt the philosophy, later promoted by Deligne in response to Goldman and Millson, that any problem in deformation theory is "controlled" by a differential graded Lie algebra, unique up to homology equivalence (quasi-isomorphism) of dg Lie algebras. Here we extend this philosophy further to control by sh-Lie-algebras.


title: Syzygies and singularities of tensor product surfaces of bidegree (2,1)

abstract: Let U be a basepoint free four-dimensional subspace of the space of sections of O(2,1) on P^1 x P^1. The sections corresponding to U determine a regular map p_U: P^1 x P^1 --> P^3. We study the associated bigraded ideal I_U in k[s,t;u,v] from the standpoint of commutative algebra, proving that there are exactly six numerical types of possible bigraded minimal free resolution. These resolutions play a key role in determining the implicit equation of the image p_U(P^1 x P^1), via work of Buse-Jouanolou, Buse-Chardin, Botbol and Botbol-Dickenstein-Dohm on the approximation complex. In four of the six cases I_U has a linear first syzygy; remarkably from this we obtain all differentials in the minimal free resolution. In particular this allows us to describe the implicit equation and singular locus of the image.


title: Information content of the low-energy electric dipole strength: correlation analysis

abstract: We study the information content carried by the electric dipole strength with respect to isovector and isoscalar indicators characterizing bulk nuclear matter and finite nuclei. To separate isoscalar and isovector modes, and low-energy strength and giant resonances, we analyze the E1 strength as a function of excitation energy $E$ and momentum transfer $q$. We use the self-consistent nuclear density functional theory with Skyrme energy density functionals, augmented by the random phase, to compute the E1 strength, and covariance analysis to assess correlations between observables. Calculations are performed for spherical, doubly-magic nuclei $^{208}$Pb and $^{132}$Sn. We demonstrate that E1 transition densities in the low-energy region below the giant dipole resonance exhibit appreciable state dependence and multi-nodal structures, which are fingerprints of weak collectivity. The correlation between the accumulated low-energy strength and symmetry energy is weak, and dramatically depends on the energy cutoff assumed. On the other hand, a strong correlation is predicted between isovector indicators and the accumulated isovector strength at $E$ around 20 MeV and momentum transfer $q\sim 0.65$ fm$^{-1}$. Momentum- and coordinate-space pattern of the low-energy dipole modes indicate a strong fragmentation into individual particle-hole excitations. The global measure of low-energy dipole strength poorly correlates with the nuclear symmetry energy and other isovector characteristics. Consequently, our results do not support the suggestion that there exists a collective "pygmy dipole resonance," which is a strong indicator of nuclear isovector properties.


title: Different Operating Systems Compatible for Image Prepress Process in Color Management: Analysis and Performance Testing

abstract: Image computing has become a real catchphrase over the past few years and the interpretations of the meaning of the term vary greatly. The Imagecomputing market is currently rapidly evolving with high growth prospects and almost daily announcements of new devices and application platforms, which results in an increasing diversification of devices, operating system and development platforms. Compared to more traditional information technology markets like the one of desktop computing, mobile computing is much less consolidated and neither standards nor even industry standards have yet been established. There are various platforms and interfaces which may be used to perform the desired tasks through the device. We have tried to compare the various mobile operating systems and their trade-offs.


title: Unstable analogues of the Lichtenbaum-Quillen conjecture

abstract: This survey is mostly concerned with unstable analogues of the Lichtenbaum-Quillen Conjecture. The Lichtenbaum-Quillen Conjecture (now implied by the Voevodsky-Rost Theorem) attempts to describe the algebraic K-theory of rings of integers in number fields in terms of much more accessible "etale models". Suitable versions of the conjecture predict the cohomology of infinite general linear groups of rings of S-integers over suitable number fields; our survey focuses on an unstable version of this form of the conjecture.


title: MKID development for SuperSpec: an on-chip, mm-wave, filter-bank spectrometer

abstract: SuperSpec is an ultra-compact spectrometer-on-a-chip for millimeter and submillimeter wavelength astronomy. Its very small size, wide spectral bandwidth, and highly multiplexed readout will enable construction of powerful multibeam spectrometers for high-redshift observations. The spectrometer consists of a horn-coupled microstrip feedline, a bank of narrow-band superconducting resonator filters that provide spectral selectivity, and Kinetic Inductance Detectors (KIDs) that detect the power admitted by each filter resonator. The design is realized using thin-film lithographic structures on a silicon wafer. The mm-wave microstrip feedline and spectral filters of the first prototype are designed to operate in the band from 195-310 GHz and are fabricated from niobium with at Tc of 9.2K. The KIDs are designed to operate at hundreds of MHz and are fabricated from titanium nitride with a Tc of 2K. Radiation incident on the horn travels along the mm-wave microstrip, passes through the frequency-selective filter, and is finally absorbed by the corresponding KID where it causes a measurable shift in the resonant frequency. In this proceedings, we present the design of the KIDs employed in SuperSpec and the results of initial laboratory testing of a prototype device. We will also briefly describe the ongoing development of a demonstration instrument that will consist of two 500-channel, R=700 spectrometers, one operating in the 1-mm atmospheric window and the other covering the 650 and 850 micron bands.


title: Accuracy of a new hybrid finite element method for solving a scattering Schroedinger equation

abstract: This hybrid method (FE-DVR), introduced by Resigno and McCurdy, Phys. Rev. A 62, 032706 (2000), uses Lagrange polynomials in each partition, rather than "hat" functions or Gaussian functions. These polynomials are discrete variable representation functions, and are.orthogonal under the Gauss-Lobatto quadrature discretization approximation. Accuracy analyses of this method are performed for the case of a one dimensional Schroedinger equation with various types of local and nonlocal potentials for scattering boundary conditions. The accuracy is ascertained by comparison with a spectral Chebyshev integral equation method, accurate to 1:10^(-11). For an accuracy of the phase shift of 1:10^(-8) The FE-DVR method is found to be 100 times faster than a sixth order finite difference method (Numerov), is easy to program, and can routinely achieve an accuracy of better than 1:10^(-8) for the numerical examples studied.


title: A New Randomness Evaluation Method with Applications to Image Shuffling and Encryption

abstract: This letter discusses the problem of testing the degree of randomness within an image, particularly for a shuffled or encrypted image. Its key contributions are: 1) a mathematical model of perfectly shuffled images; 2) the derivation of the theoretical distribution of pixel differences; 3) a new $Z$-test based approach to differentiate whether or not a test image is perfectly shuffled; and 4) a randomized algorithm to unbiasedly evaluate the degree of randomness within a given image. Simulation results show that the proposed method is robust and effective in evaluating the degree of randomness within an image, and may often be more suitable for image applications than commonly used testing schemes designed for binary data like NIST 800-22. The developed method may be also useful as a first step in determining whether or not a shuffling or encryption scheme is suitable for a particular cryptographic application.


title: Jet spectral breaks in black hole X-ray binaries

abstract: In X-ray binaries, compact jets are known to commonly radiate at radio to infrared frequencies, whereas at optical to gamma-ray energies, the contribution of the jet is debated. The total luminosity, and hence power of the jet is critically dependent on the position of the break in its spectrum, between optically thick (self-absorbed) and optically thin synchrotron emission. This break, or turnover, has been reported in just one black hole X-ray binary (BHXB) thus far, GX 339-4, and inferred via spectral fitting in two others, A0620-00 and Cyg X-1. Here, we collect a wealth of multiwavelength data from the outbursts of BHXBs during hard X-ray states, in order to search for jet breaks as yet unidentified in their spectral energy distributions. In particular, we report the direct detection of the jet break in the spectrum of V404 Cyg during its 1989 outburst, at nu_b = (1.8 +- 0.3) x 10^14 Hz (1.7 +- 0.2 microns). We increase the number of BHXBs with measured jet breaks from three to eight. Jet breaks are found at frequencies spanning more than two orders of magnitude, from nu_b = (4.5 +- 0.8) x 10^12 Hz for XTE J1118+480 during its 2005 outburst, to nu_b > 4.7 x 10^14 Hz for V4641 Sgr in outburst. A positive correlation between jet break frequency and luminosity is expected theoretically; nu_b \propto L_nu,jet^~0.5 if other parameters are constant. With constraints on the jet break in a total of 12 BHXBs including two quiescent systems, we find a large range of jet break frequencies at similar luminosities and no obvious global relation (but such a relation cannot be ruled out for individual sources). We speculate that different magnetic field strengths and/or different radii of the acceleration zone in the inner regions of the jet are likely to be responsible for the observed scatter between sources. (abridged)


title: James-Stein Type Center Pixel Weights for Non-Local Means Image Denoising

abstract: Non-Local Means (NLM) and variants have been proven to be effective and robust in many image denoising tasks. In this letter, we study the parameter selection problem of center pixel weights (CPW) in NLM. Our key contributions are: 1) we give a novel formulation of the CPW problem from the statistical shrinkage perspective; 2) we introduce the James-Stein type CPWs for NLM; and 3) we propose a new adaptive CPW that is locally tuned for each image pixel. Our experimental results showed that compared to existing CPW solutions, the new proposed CPWs are more robust and effective under various noise levels. In particular, the NLM with the James-Stein type CPWs attain higher means with smaller variances in terms of the peak signal and noise ratio, implying they improve the NLM robustness and make it less sensitive to parameter selection.


title: New 1/N expansions in random tensor models

abstract: Although random tensor models were introduced twenty years ago, it is only in 2011 that Gurau proved the existence of a 1/N expansion. Here we show that there actually is more than a single 1/N expansion, depending on the dimension. These new expansions can be used to define tensor models for `rectangular' tensors (whose indices have different sizes). In the large N limit, they retain more than the melonic graphs. Still, in most cases, the large N limit is found to be Gaussian, and therefore extends the scope of the universality theorem for large random tensors. Nevertheless, a scaling which leads to non-Gaussian large N limits, in even dimensions, is identified for the first time.


title: Extending Task Parallelism for Frequent Pattern Mining

abstract: Algorithms for frequent pattern mining, a popular informatics application, have unique requirements that are not met by any of the existing parallel tools. In particular, such applications operate on extremely large data sets and have irregular memory access patterns. For efficient parallelization of such applications, it is necessary to support dynamic load balancing along with scheduling mechanisms that allow users to exploit data locality. Given these requirements, task parallelism is the most promising of the available parallel programming models. However, existing solutions for task parallelism schedule tasks implicitly and hence, custom scheduling policies that can exploit data locality cannot be easily employed. In this paper we demonstrate and characterize the speedup obtained in a frequent pattern mining application using a custom clustered scheduling policy in place of the popular Cilk-style policy. We present PFunc, a novel task parallel library whose customizable task scheduling and task priorities facilitated the implementation of our clustered scheduling policy.


title: Measurement of the Coulomb Logarithm in a Radio-Frequency Paul Trap

abstract: Samples of ultracold 174 Yb+ ions, confined in a linear radio-frequency Paul trap, are heated via self-induced micromotion interruption, while their temperature, density, and therefore structural phase are monitored and simulated. The observed time evolution of the ion temperature is compared to a theoretical model for ion-ion heating allowing a direct measurement of the Coulomb logarithm in a linear Paul trap. This result permits a simple, yet accurate, analytical description of ion cloud thermodynamic properties, e.g. density, temperature, and structural phase, as well as suggests limits to and improvements for on-going trapped-ion quantum information efforts.


title: Secret-Key Agreement Capacity over Reciprocal Fading Channels: A Separation Approach

abstract: Fundamental limits of secret-key agreement over reciprocal wireless channels are investigated. We consider a two-way block-fading channel where the channel gains in the forward and reverse links between the legitimate terminals are correlated. The channel gains between the legitimate terminals are not revealed to any terminal, whereas the channel gains of the eavesdropper are revealed to the eavesdropper. We propose a two-phase transmission scheme, that reserves a certain portion of each coherence block for channel estimation, and the remainder of the coherence block for correlated source generation. The resulting secret-key involves contributions of both channel sequences and source sequences, with the contribution of the latter becoming dominant as the coherence period increases. We also establish an upper bound on the secret-key capacity, which has a form structurally similar to the lower bound. Our upper and lower bounds coincide in the limit of high signal-to-noise-ratio (SNR) and large coherence period, thus establishing the secret-key agreement capacity in this asymptotic regime. Numerical results indicate that the proposed scheme achieves significant gains over training-only schemes, even for moderate SNR and small coherence periods, thus implying the necessity of randomness-sharing in practical secret-key generation systems.


title: A One-Vertex Decomposition Algorithm for Generating Algebraic Expressions of Square Rhomboids

abstract: The paper investigates relationship between algebraic expressions and graphs. We consider a digraph called a square rhomboid that is an example of non-series-parallel graphs. Our intention is to simplify the expressions of square rhomboids and eventually find their shortest representations. With that end in view, we describe the new algorithm for generating square rhomboid expressions, which improves on our previous algorithms.


title: Non-abelian vortices on CP^1 and Grassmannians

abstract: Many properties of the moduli space of abelian vortices on a compact Riemann surface are known. For non-abelian vortices the moduli space is less well understood. Here we consider non-abelian vortices on the Riemann sphere CP^1, and we study their moduli spaces near the Bradlow limit. We give an explicit description of the moduli space as a Kahler quotient of a finite-dimensional linear space. The dimensions of some of these moduli spaces are derived. Strikingly, there exist non-abelian vortex configurations on CP^1, with non-trivial vortex number, for which the moduli space is a point. This is in stark contrast to the moduli space of abelian vortices. For a special class of non-abelian vortices the moduli space is a Grassmannian, and the metric near the Bradlow limit is a natural generalization of the Fubini--Study metric on complex projective space. We use this metric to investigate the statistical mechanics of non-abelian vortices. The partition function is found to be analogous to the one for abelian vortices.


title: Exchange Field-Mediated Magnetoresistance in the Correlated Insulator Phase of Be Films

abstract: We present a study of the proximity effect between a ferromagnet and a paramagnetic metal of varying disorder. Thin beryllium films are deposited onto a 5 nm-thick layer of the ferromagnetic insulator EuS. This bilayer arrangement induces an exchange field, $H_{ex}$, of a few tesla in low resistance Be films with sheet resistance $R\ll R_Q$, where $R_Q=h/e^2$ is the quantum resistance. We show that $H_{ex}$ survives in very high resistance films and, in fact, appears to be relatively insensitive to the Be disorder. We exploit this fact to produce a giant low-field magnetoresistance in the correlated insulator phase of Be films with $R\gg R_Q$.


title: The Faint End of the Cluster Galaxy Luminosity Function at High Redshift

abstract: We measure the faint end slope of the galaxy luminosity function (LF) for cluster galaxies at 1<z<1.5 using Spitzer IRAC data. We investigate whether this slope, alpha, differs from that of the field LF at these redshifts, and with the cluster LF at low redshifts. The latter is of particular interest as low-luminosity galaxies are expected to undergo significant evolution. We use seven high-redshift spectroscopically confirmed galaxy clusters drawn from the IRAC Shallow Cluster Survey to measure the cluster galaxy LF down to depths of M* + 3 (3.6 microns) and M* + 2.5 (4.5 microns). The summed LF at our median cluster redshift (z=1.35) is well fit by a Schechter distribution with alpha[3.6] = -0.97 +/- 0.14 and alpha[4.5] = -0.91 +/- 0.28, consistent with a flat faint end slope and is in agreement with measurements of the field LF in similar bands at these redshifts. A comparison to alpha in low-redshift clusters finds no statistically significant evidence of evolution. Combined with past studies which show that M* is passively evolving out to z~1.3, this means that the shape of the cluster LF is largely in place by z~1.3. This suggests that the processes that govern the build up of the mass of low-mass cluster galaxies have no net effect on the faint end slope of the cluster LF at z<1.3.


title: Spin torque switching in perpendicular films at finite temperature

abstract: We show how the phase diagram for spin torque switching in the case of perpendicular anisotropy is altered at nonzero temperature. The hysteresis region in which the parallel and antiparallel states coexist shrinks, and a new region of telegraph noise appears. In a small sample, the region of coexistence of a precessional and parallel state can disappear entirely. We show that the phase diagram for both zero and nonzero temperatures can be understood and calculated by plotting an effective energy as a function of angle. A combinatorial analysis is useful for systematically describing the phase diagram.


title: Pairings on Generalized Huff Curves

abstract: This paper presents the Tate pairing computation on generalized Huff curves proposed by Wu and Feng. In fact, we extend the results of the Tate pairing computation on the standard Huff elliptic curves done previously by Joye, Tibouchi and Vergnaud. We show that the addition step of the Miller loop can be performed in $1\mathbf{M}+(k+15)\mathbf{m}+2\mathbf{c}$ and the doubling one in $1\mathbf{M} + 1\mathbf{S} + (k + 12) \mathbf{m} + 5\mathbf{s} + 2\mathbf{c}$ on the generalized Huff curve.


title: Giant atomic displacement induced by built-in strain in metastable Mn$_3$O$_4$

abstract: We present x-ray, neutron scattering and heat capacity data that reveal a coupled first-order magnetic and structural phase transition of the metastable mixed-valence post-spinel compound Mn$_3$O$_4$ at 210 K. Powder neutron diffraction measurements reveal a magnetic structure in which Mn$^{3+}$ spins align antiferromagnetically along the edge-sharing \emph{a}-axis, with a magnetic propagation vector k = [1/2, 0, 0]. In contrast, the Mn$^{2+}$ spins, which are geometrically frustrated, do not order until a much lower temperature. Although the Mn$^{2+}$ spins do not directly participate in the magnetic phase transition at 210 K, structural refinements reveal a large atomic shift at this phase transition, corresponding to a physical motion of approximately 0.25 {\AA} even though the crystal symmetry remains unchanged. This "giant" response is due to the coupled effect of built-in strain in the metastable post-spinel structure with the orbital realignment of the Mn$^{3+}$ ion.


title: Anisotropic CR diffusion and gamma-ray production close to supernova remnants, with an application to W28

abstract: Cosmic rays that escape their acceleration site interact with the ambient medium and produce gamma rays as the result of inelastic proton-proton collisions. The detection of such diffuse emission may reveal the presence of an accelerator of cosmic rays, and also constrain the cosmic ray diffusion coefficient in its vicinity. Preliminary results in this direction have been obtained in the last years from studies of the gamma-ray emission from molecular clouds located in the vicinity of supernova remnants, which are the prime candidate for cosmic ray production. Hints have been found for a significant suppression of the diffusion coefficient with respect to the average one in the Galaxy. However, most of these studies rely on the assumption of isotropic diffusion, which may not be very well justified. Here, we extend this study to the case in which cosmic rays that escape an accelerator diffuse preferentially along the magnetic field lines. As a first approximation, we further assume that particles are strongly magnetized and that their transport perpendicular to the magnetic field is mainly due to the wandering of the field lines. The resulting spatial distribution of runaway cosmic rays around the accelerator is, in this case, strongly anisotropic. An application of the model to the case of the supernova remnant W28 demonstrates how the estimates of the diffusion coefficient from gamma-ray observations strongly depend on the assumptions made on the isotropy (or anisotropy) of diffusion. For higher levels of anisotropy of the diffusion, larger values of the diffusion coefficient are found to provide a good fit to data. Thus, detailed models for the propagation of cosmic rays are needed in order to interpret in a correct way the gamma-ray observations.


title: A Deep X-ray View of the Hot Halo in the Edge-on Spiral Galaxy NGC 891

abstract: NGC 891 is a nearby edge-on galaxy that is similar to the Milky Way and has a hot X-ray emitting halo that could arise from accretion, a galactic fountain, or a combination of the two. The metallicity of the gas can help distinguish between these models, and here we report on results that use 138 ks of archival Chandra data and 92 ks of new XMM-Newton data to measure the temperature and metallicity of the hot halo of the galaxy. We find good fits for a thermal model with kT ~ 0.2 keV and Z ~ 0.1 solar, and rule out solar metallicity to more than 99% confidence. This result suggests accretion from the intergalactic medium as the origin for the hot halo. However, it is also possible to fit a two-temperature thermal model with solar metallicity where kT_1 = 0.1 keV and kT_2 = 0.25 keV. A consideration of the cooling rate and scale height prefers the single-temperature model. We also find that the cooling rate in the hot gas cannot explain the massive HI halo in the steady state. In addition, a galactic fountain model cannot eject enough mass to account for the HI halo, and we speculate that the neutral halo may be gas from a prior outflow that has since cooled.


title: The First Billion Years project - IV: Proto-galaxies reionising the Universe

abstract: The contribution of stars in galaxies to cosmic reionisation depends on the star formation history in the Universe, the abundance of galaxies during reionisation, the escape fraction of ionising photons and the clumping factor of the inter-galactic medium (IGM). We compute the star formation rate and clumping factor during reionisation in a cosmological volume using a high-resolution hydrodynamical simulation. We post-process the output with detailed radiative transfer simulations to compute the escape fraction of ionising photons. Together, this gives us the opportunity to assess the contribution of galaxies to reionisation self-consistently. The strong mass and redshift dependence of the escape fraction indicates that reionisation occurred between z=15 and z=10 and was mainly driven by proto-galaxies forming in dark-matter haloes with masses between 1e7 and 1e8 solar mass. More massive galaxies that are rare at these redshifts and have significantly lower escape fractions contribute less photons to the reionisation process than the more-abundant low-mass galaxies. Star formation in the low-mass haloes is suppressed by radiative feedback from reionisation, therefore these proto-galaxies only contribute when the part of the Universe they live in is still neutral. After z~10, massive galaxies become more abundant and provide most of the ionising photons. In addition, we find that Population (Pop) III stars are too short-lived and not frequent enough to have a major contribution to reionisation. Although the stellar component of the proto-galaxies that produce the bulk of ionising photons during reionisation is too faint to be detected by the James Webb Space Telescope (JWST), these sources are brightest in the H-alpha and Ly-alpha recombination lines, which will likely be detected by JWST in deep surveys.


title: The Imprint of The Extragalactic Background Light in the Gamma-Ray Spectra of Blazars

abstract: The light emitted by stars and accreting compact objects through the history of the Universe is encoded in the intensity of the extragalactic background light (EBL). Knowledge of the EBL is important to understand the nature of star formation and galaxy evolution, but direct measurements of the EBL are lim- ited by Galactic and other foreground emissions. Here we report an absorption feature seen in the combined spectra of a sample of gamma-ray blazars out to a redshift of z$\sim$1.6. This feature is caused by attenuation of gamma rays by the EBL at optical to UV frequencies, and allowed us to measure the EBL flux density in this frequency band.


title: The Formation and Evolution of Wind-Capture Disks In Binary Systems

abstract: We study the formation, evolution and physical properties of accretion disks formed via wind capture in binary systems. Using the AMR code AstroBEAR, we have carried out high resolution 3D simulations that follow a stellar mass secondary in the co-rotating frame as it orbits a wind producing AGB primary. We first derive a resolution criteria, based on considerations of Bondi-Hoyle flows, that must be met in order to properly resolve the formation of accretion disks around the secondary. We then compare simulations of binaries with three different orbital radii (10, 15, 20 AU). Disks are formed in all three cases, however the size of the disk and, most importantly, its accretion rate decreases with orbital radii. In addition, the shape of the orbital motions of material within the disk becomes increasingly elliptical with increasing binary separation. The flow is mildly unsteady with "fluttering" around the bow shock observed. The disks are generally well aligned with the orbital plane after a few binary orbits. We do not observe the presence of any large scale, violent instabilities (such as the flip-flop mode). For the first time, moreover, it is observed that the wind component that is accreted towards the secondary has a vortex tube-like structure, rather than a column-like one as it was previously thought. In the context of AGB binary systems that might be precursors to Pre-Planetary and Planetary Nebula, we find that the wind accretion rates at the chosen orbital separations are generally too small to produce the most powerful outflows observed in these systems if the companions are main sequence stars but marginally capable if the companions are white dwarfs. It is likely that many of the more powerful PPN and PN involve closer binaries than the ones considered here. The results also demonstrate principles of broad relevance to all wind-capture binary systems.


title: The Minimum-Mass Extrasolar Nebula: In-Situ Formation of Close-In Super-Earths

abstract: Close-in super-Earths, with radii R = 2-5 R_Earth and orbital periods P < 100 days, orbit more than half, and perhaps nearly all Sun-like stars in the universe. We use this omnipresent population to construct the minimum-mass extrasolar nebula (MMEN), the circumstellar disk of solar-composition solids and gas from which such planets formed, if they formed near their current locations and did not migrate. In a series of back-of-the-envelope calculations, we demonstrate how in-situ formation in the MMEN is fast, efficient, and can reproduce many of the observed properties of close-in super-Earths, including their gas-to-rock fractions. Testable predictions are discussed.


title: XMM-Newton Observations of Three Interacting Luminous Infrared Galaxies

abstract: We investigate the X-ray properties of three interacting luminous infrared galaxy systems. In one of these systems, IRAS 18329+5950, we resolve two separate sources. A second, IRAS 20550+1656, and third, IRAS 19354+4559, have only a single X-ray source detected. We compare the observed emission to PSF profiles and determine that three are extended in emission. One is compact, which is suggestive of an AGN, although all of our profiles have large uncertainties. We then model the spectra to determine soft (0.5--2 keV) and hard (2--10 keV) luminosities for the resolved sources and then compare these to relationships found in the literature between infrared and X-ray luminosities for starburst galaxies. We obtain luminosities of $\log(L_{\textrm{soft}}/\textrm{L}_{\odot}) = 7.32,\:7.06,\:7.68$ and $\log(L_{\textrm{hard}}/\textrm{L}_{\odot}) = 7.33,\: 7.07,\: 7.88$ for IRAS 18329+5950, IRAS 19354+4559, and IRAS 20550+1656, respectively. These are intermediate to two separate predictions in the literature for star-formation-dominated sources. Our highest quality spectrum of IRAS 20550+1656 suggests super-solar abundance of alpha elements at $2\sigma$ significance, with $\log(\frac{\alpha}{\alpha_{\odot}}) = [\alpha] = 0.4\pm0.2$. This is suggestive of recent enrichment with Type II supernovae, consistent with a starburst environment. The X-ray properties of the target galaxies are most likely due to starbursts, but we cannot conclusively rule out AGN.


title: Stability of edge states in strained graphene

abstract: Spatially inhomogeneous strains in graphene can simulate the effects of valley-dependent magnetic fields. As demonstrated in recent experiments, the realizable magnetic fields are large enough to give rise to well-defined flat pseudo-Landau levels, potentially having counter-propagating edge modes. In the present work we address the conditions under which such edge modes are visible. We find that, whereas armchair edges do not support counter-propagating edge modes, zigzag edges do so, through a novel selective-hybridization mechanism. We then discuss effects of interactions on the stability of counter-propagating edge modes, and find that, for the experimentally relevant case of Coulomb interactions, interactions typically decrease the stability of the edge modes. Finally, we generalize our analysis to address the case of spontaneous valley polarization, which is expected to occur in charge-neutral strained graphene.


title: Ground States of Spin-1/2 Triangular Antiferromagnets in a Magnetic Field

abstract: We use a combination of density matrix renormalization group calculations and analytical approaches to study a simplified model for a spatially anisotropic spin-1/2 triangular lattice Heisenberg antiferromagnet: the three-leg triangular spin tube (TST). The model is described by three Heisenberg chains, with exchange constant J, coupled antiferromagnetically with exchange constant J' along the diagonals of the ladder system, with periodic boundary conditions in the shorter direction. We determine the full phase diagram of this model as a function of spatial anisotropy, J'/J, and magnetic field. We find a rich phase diagram, which is dominated by quantum states - phases corresponding to the classical ground state appears only in a small region. Among the dominant phases generated by quantum effects are commensurate and incommensurate coplanar quasi-ordered states, which appear in the vicinity of the isotropic region for most fields, and in the high field region for most anisotropies. The coplanar states, while not classical ground states, can be understood semiclassically. Even more strikingly, the largest region of phase space is occupied by a spin density wave phase, which has incommensurate collinear correlations along the field. This phase has no semiclassical analog, and may be ascribed to enhanced one-dimensional (1d) fluctuations due to frustration. Cutting across the phase diagram is a magnetization plateau, with a gap to all excitations and up up down spin order, with a quantized magnetization equal to 1/3 of the saturation value. In the TST, this plateau extends almost but not quite to the decoupled chains limit. Most of the above features are expected to carry over to the two dimensional system, which we also discuss. At low field, a dimerized phase appears, which is particular to the 1d nature of the TST, and which can be understood from quantum Berry phase arguments.


title: The Co-evolution of Cosmic Entropy and Structures in the Universe

abstract: According to the second law of thermodynamics, the arrow of time points to an ever increasing entropy of the Universe. However, exactly how the entropy evolves with time and what drives the growth remain largely unknown. Here, for the first time, we quantify the evolving entropy of cosmic structures using a large-scale cosmological hydrodynamical simulation. Our simulation starts from initial conditions predicted by the leading LambdaCDM cosmology, self-consistently evolves the dynamics of both dark and baryonic matter, star formation, black hole growth and feedback processes, from the cosmic dawn to the present day. Tracing the entropy contributions of these distinct components in the simulation, we find a strong link between entropy growth and structure formation. The entropy is dominated by that of the black holes in all epochs, and its evolution follows the same path as that of galaxies: it increases rapidly from a low-entropy state at high redshift until z~2, then transits to a slower growth. Our results suggest that cosmic entropy may co-evolve with cosmic structure, and that its growth may be driven mainly by the formation of black holes in galaxies. We predict that the entropy will continue to increase in the near future, but likely at a constant rate.


title: On trace anomaly in 2+1 flavor QCD

abstract: We report on recent progress by the HotQCD collaboration in studying the trace anomaly at non-zero temperature in 2+1 flavor QCD on lattices with the temporal extent N_tau=4, 6, 8, 10 and 12 using the highly improved staggered quark (HISQ) action as well as the asqtad action. We discuss the dependence of our lattice results on the scale setting procedure and compare them with hadron resonance gas (HRG) calculations as well as to resummed perturbative results.


title: Dynamics around black holes: Radiation Emission and Tidal Effects

abstract: In this thesis we study several dynamical processes involving black holes in four and higher dimensions. First, using perturbative techniques, we compare the massless and massive scalar radiation emitted by a particle radially infalling into a Schwarzchild black hole. We show that the late-time waveform of massive scalar perturbations is dominated by a universal oscillatory decaying tail, which appears due to curvature effects. We also show that the energy spectrum is in perfect agreement with a ZFL calculation once no-hair properties of black holes are taken into account. In the second part, we study the phenomenon of superradiance in higher dimensions and conjecture that the maximum energy extracted from a rotating black hole can be understood in terms of the ergoregion proper volume. We then study some consequences of superradiance in the dynamics of moons orbiting around higher-dimensional rotating black holes. In four-dimensional spacetime, moons around black holes generate low-amplitude tides, and the energy extracted from the hole's rotation is always smaller than the gravitational radiation lost to infinity. We show that in dimensions larger than five the energy extracted from the black hole through superradiance is larger than the energy carried out to infinity. Our results lend strong support to the conjecture that tidal acceleration is the rule, rather than the exception, in higher dimensions. Superradiance dominates the energy budget and moons "outspiral"; for some particular orbital frequency, the energy extracted at the horizon equals the energy emitted to infinity and "floating orbits" generically occur. We give an interpretation of this phenomenon in terms of the membrane paradigm and of tidal acceleration due to energy dissipation across the horizon.


title: S2LET: A code to perform fast wavelet analysis on the sphere

abstract: We describe S2LET, a fast and robust implementation of the scale-discretised wavelet transform on the sphere. Wavelets are constructed through a tiling of the harmonic line and can be used to probe spatially localised, scale-depended features of signals on the sphere. The scale-discretised wavelet transform was developed previously and reduces to the needlet transform in the axisymmetric case. The reconstruction of a signal from its wavelets coefficients is made exact here through the use of a sampling theorem on the sphere. Moreover, a multiresolution algorithm is presented to capture all information of each wavelet scale in the minimal number of samples on the sphere. In addition S2LET supports the HEALPix pixelisation scheme, in which case the transform is not exact but nevertheless achieves good numerical accuracy. The core routines of S2LET are written in C and have interfaces in Matlab, IDL and Java. Real signals can be written to and read from FITS files and plotted as Mollweide projections. The S2LET code is made publicly available, is extensively documented, and ships with several examples in the four languages supported. At present the code is restricted to axisymmetric wavelets but will be extended to directional, steerable wavelets in a future release.


title: Kiloparsec-Scale Simulations of Star Formation in Disk Galaxies. I. The unmagnetized and zero-feedback limit

abstract: We present hydrodynamic simulations of the evolution of self-gravitating dense gas on scales of 1 kiloparsec down to < parsec in a galactic disk, designed to study dense clump formation from giant molecular clouds (GMCs). These structures are expected to be the precursors to star clusters and this process may be the rate limiting step controling star formation rates in galactic systems as described by the Kennicutt-Schmidt relation. We follow the thermal evolution of the gas down to ~5K using extinction-dependent heating and cooling functions. We do not yet include magnetic fields or localized stellar feedback, so the evolution of the GMCs and clumps is determined solely by self-gravity balanced by thermal and turbulent pressure support and the large scale galactic shear. While cloud structures and densities change significantly during the simulation, GMC virial parameters remain mostly above unity for time scales exceeding the free-fall time of GMCs indicating that energy from galactic shear and large-scale cloud motions continuously cascades down to and within the GMCs. We implement star formation at a slow, inefficient rate of 2% per local free-fall time, but even this yields global star formation rates that are about two orders of magnitude larger than the observed Kennicutt-Schmidt relation due to over-production of dense gas clumps. We expect a combination of magnetic support and localized stellar feedback is required to inhibit dense clump formation to ~1% of the rate that results from the nonmagnetic, zero-feedback limit.


title: Solution of the Anderson impurity model via the functional renormalization group

abstract: We show that the functional renormalization group is a numerically cheap method to obtain the low-energy behavior of the Anderson impurity model describing a localized interacting electron coupled to a bath of conduction electrons. Our approach uses an external magnetic field as flow parameter, partial bosonization of the transverse spin fluctuations, and frequency-independent interaction vertices which are fixed by Ward identities. We calculate the quasi-particle residue and the spin susceptibility in the particle-hole symmetric case and obtain excellent agreement with the Bethe ansatz results for arbitrary strengths of the interaction.


title: Electrodynamic Response and Stability of Molecular Crystals

abstract: We show that electrodynamic dipolar interactions, responsible for long-range fluctuations in matter, play a significant role in the stability of molecular crystals. Density functional theory calculations with van der Waals interactions determined from a semilocal "atom-in-a-molecule" model result in a large overestimation of the dielectric constants and sublimation enthalpies for polyacene crystals from naphthalene to pentacene, whereas an accurate treatment of non-local electrodynamic response leads to an agreement with the measured values for both quantities. Our findings suggest that collective response effects play a substantial role not only for optical excitations, but also for cohesive properties of non-covalently bound molecular crystals.


title: Complex-mass scheme and perturbative unitarity

abstract: We derive cutting rules for loop integrals containing propagators with complex masses. Using a field-theoretical model of a heavy vector boson interacting with a light fermion, we demonstrate that the complex-mass scheme respects unitarity order by order in a perturbative expansion provided that the renormalized coupling constant remains real.


title: Domain wall description of superconductivity

abstract: In the present work we shall address the issue of electrical conductivity in superconductors in the perspective of superconducting domain wall solutions in the realm of field theory. We take our set up made out of a dynamical complex scalar field coupled to gauge field to be responsible for superconductivity and an extra scalar real field that plays the role of superconducting domain walls. The temperature of the system is interpreted through the fact that the soliton following accelerating orbits is a Rindler observer experiencing a thermal bath.


title: Tidal heating and torquing of a Kerr black hole to next-to-leading order in the tidal coupling

abstract: We calculate the linear vacuum perturbations of a Kerr black hole surrounded by a slowly-varying external spacetime to third order in the ratio of the black-hole mass to the radius of curvature of the external spacetime. This expansion applies to two relevant physical scenarios: (i) a small Kerr black hole immersed in the gravitational field of a much larger external black hole, and (ii) a Kerr black hole moving slowly around another external black hole of comparable mass. This small-hole/slow-motion approximation allows us to parametrize the perturbation through slowly-varying, time-dependent electric and magnetic tidal tensors, which then enables us to separate the Teukolsky equation and compute the Newman-Penrose scalar analytically to third order in our expansion parameter. We obtain generic expressions for the mass and angular momentum flux through the perturbed black hole horizon, as well as the rate of change of the horizon surface area, in terms of certain invariants constructed from the electric and magnetic tidal tensors. We conclude by applying these results to the two scenarios described above.


title: Robustness of a topological phase: Topological color code in parallel magnetic field

abstract: The robustness of the topological color code, which is a class of error correcting quantum codes, is investigated under the influence of an uniform magnetic field on the honeycomb lattice. Our study relies on two high-order series expansions using perturbative continuous unitary transformations in the limit of low and high fields, exact diagonalization and a classical approximation. We show that the topological color code in a single parallel field is isospectral to the Baxter-Wu model in a transverse field on the triangular lattice. It is found that the topological phase is stable up to a critical field beyond which it breaks down to the polarized phase by a first-order phase transition. The results also suggest that the topological color code is more robust than the toric code, in the parallel magnetic field.


title: Star Formation Histories of Globular Clusters with Multiple Populations. I. \omega\ Cen, M22, and NGC 1851

abstract: There is increasing evidence that some massive globular clusters (GCs) host multiple stellar populations having different heavy element abundances enriched by supernovae. They usually accompany multiple red giant branches (RGBs) in the color-magnitude diagrams (CMDs), and are distinguished from most of the other GCs which display variations only in light element abundances. In order to investigate the star formation histories of these peculiar GCs, we have constructed synthetic CMDs for \omega\ Cen, M22, and NGC 1851. Our models are based on the updated versions of Yonsei-Yale (Y^2) isochrones and horizontal branch (HB) evolutionary tracks which include the cases of enhancements in both helium and the total CNO abundances. To estimate ages and helium abundances of subpopulations in each GC, we have compared our models with the observations on the Hess diagram by employing a \chi^2 minimization technique. We find that metal-rich subpopulations in each of these GCs are also enhanced in helium abundance, and the age differences between the metal-rich and metal-poor subpopulations are fairly small (~0.3-1.7 Gyr), even in the models with the observed variations in the total CNO content. These are required to simultaneously reproduce the observed extended HB and the splits on the main sequence, subgiant branch, and RGB. Our results are consistent with the hypothesis that these GCs are the relics of more massive primeval dwarf galaxies that merged and disrupted to form the proto-Galaxy.


title: Spectrum of hyperplane arrangements in four variables

abstract: One of the most important invariants in singularity theory is the Hodge spectrum. Calculating the Hodge spectrum is a difficult task and formulas exist for only a few cases. In this article the main result is the formula for reduced hyperplane arrangements in four variables.


title: Learning Monocular Reactive UAV Control in Cluttered Natural Environments

abstract: Autonomous navigation for large Unmanned Aerial Vehicles (UAVs) is fairly straight-forward, as expensive sensors and monitoring devices can be employed. In contrast, obstacle avoidance remains a challenging task for Micro Aerial Vehicles (MAVs) which operate at low altitude in cluttered environments. Unlike large vehicles, MAVs can only carry very light sensors, such as cameras, making autonomous navigation through obstacles much more challenging. In this paper, we describe a system that navigates a small quadrotor helicopter autonomously at low altitude through natural forest environments. Using only a single cheap camera to perceive the environment, we are able to maintain a constant velocity of up to 1.5m/s. Given a small set of human pilot demonstrations, we use recent state-of-the-art imitation learning techniques to train a controller that can avoid trees by adapting the MAVs heading. We demonstrate the performance of our system in a more controlled environment indoors, and in real natural forest environments outdoors.


title: Spin Orbit Magnetism and Unconventional Superconductivity

abstract: We find an exotic spin excitation in a magnetically ordered system with spin orbit magnetism in 2D, where the order parameter has a net spin current and no net magnetization. Starting from a Fermi liquid theory, similar to that for a weak ferromagnet, we show that this excitation emerges from an exotic magnetic Fermi liquid state that is protected by a generalized Pomeranchuck condition. We derive the propagating mode using the Landau kinetic equation, and find that the dispersion of the mode has a $\sqrt q$ behavior in leading order in 2D. We find an instability toward superconductivity induced by this exotic mode, and a further analysis based on the forward scattering sum rule strongly suggests that this superconductivity has p-wave pairing symmetry. We perform similar studies in the 3D case, with a slightly different magnetic system and find that the mode leads to a Lifshitz-like instability most likely toward an inhomogeneous magnetic state in one of the phases.


title: Divisorial models of normal varieties

abstract: We prove that the canonical ring of a canonical variety in the sense of de Fernex and Hacon is finitely generated. We prove that canonical varieties are klt if and only if R(-K_X) is finitely generated. We introduce a notion of nefness for non-Q-Gorenstein varieties and study some of its properties. We then focus on these properties for non-Q-Gorenstein toric varieties.


title: The Constrained NMSSM with a 126 GeV Higgs boson: A global analysis

abstract: We present the first global analysis of the Constrained NMSSM that investigates the impact of the recent discovery of a 126 GeV Higgs-like boson, of the observation of a signal for BR(B_s->\mu^+\mu^-), and of constraints on supersymmetry from ~5/fb of data accumulated at the LHC, as well as of other relevant constraints from colliders, flavor physics and dark matter. We consider three possible cases, assuming in turn that the discovered Higgs boson is: i) the lightest Higgs boson of the model; ii) the next-to-lightest Higgs boson; and iii) a combination of both roughly degenerate in mass. The likelihood function for the Higgs signal uses signal rates in the \gamma\gamma\ and ZZ --> 4l channels, while that for the Higgs exclusion limits assumes decay through the \gamma\gamma, \tau\tau, ZZ and W^+W^- channels. In all cases considered we identify the 68% and 95% credible posterior probability regions in a Bayesian approach. We find that, when the constraints are applied with their respective uncertainties, the first case shows strong CMSSM-like behavior, with the stau coannihilation region featuring highest posterior probability, the best-fit point, a correct mass of the lightest Higgs boson and the lighter stop mass in the ballpark of 1 TeV. The second and the third cases, on the other hand, while allowed are disfavored by the constraints from direct detection of dark matter and from BR(B_s->\mu^+\mu^-). Without the anomalous magnetic moment of the muon the fit improves considerably, especially for negative effective \mu\ parameter. We discuss how the considered scenarios could be tested furtherly at the LHC and in dark matter searches.


title: Does a Daily Deal Promotion Signal a Distressed Business? An Empirical Investigation of Small Business Survival

abstract: In the last four years, daily deals have emerged from nowhere to become a multi-billion dollar industry world-wide. Daily deal sites such as Groupon and Livingsocial offer products and services at deep discounts to consumers via email and social networks. As the industry matures, there are many questions regarding the impact of daily deals on the marketplace. Important questions in this regard concern the reasons why businesses decide to offer daily deals and their longer-term impact on businesses. In the present paper, we investigate whether the unobserved factors that make marketers run daily deals are correlated with the unobserved factors that influence the business, In particular, we employ the framework of seemingly unrelated regression to model the correlation between the errors in predicting whether a business uses a daily deal and the errors in predicting the business' survival. Our analysis consists of the survival of 985 small businesses that offered daily deals between January and July 2011 in the city of Chicago. Our results indicate that there is a statistically significant correlation between the unobserved factors that influence the business' decision to offer a daily deal and the unobserved factors that impact its survival. Furthermore, our results indicate that the correlation coefficient is significant in certain business categories (e.g. restaurants).


title: Properties of gas clumps and gas clumping factor in the intra cluster medium

abstract: The spatial distribution of gas matter inside galaxy clusters is not completely smooth, but may host gas clumps associated with substructures. These overdense gas substructures are generally a source of unresolved bias of X-ray observations towards high density gas, but their bright luminosity peaks may be resolved sources within the ICM, that deep X-ray exposures may be (already) capable to detect. In this paper we aim at investigating both features, using a set of high-resolution cosmological simulations with ENZO. First, we monitor how the bias by unresolved gas clumping may yield incorrect estimates of global cluster parameters and affects the measurements of baryon fractions by X-ray observations. We find that based on X-ray observations of narrow radial strips, it is difficult to recover the real baryon fraction to better than 10 - 20 percent uncertainty. Second, we investigated the possibility of observing bright X-ray clumps in the nearby Universe (z<=0.3). We produced simple mock X-ray observations for several instruments (XMM, Suzaku and ROSAT) and extracted the statistics of potentially detectable bright clumps. Some of the brightest clumps predicted by simulations may already have been already detected in X- ray images with a large field of view. However, their small projected size makes it difficult to prove their existence based on X-ray morphology only. Preheating, AGN feedback and cosmic rays are found to have little impact on the statistical properties of gas clumps.


title: On the Economic Value and Price-Responsiveness of Ramp-Constrained Storage

abstract: The primary concerns of this paper are twofold: to understand the economic value of storage in the presence of ramp constraints and exogenous electricity prices, and to understand the implications of the associated optimal storage management policy on qualitative and quantitative characteristics of storage response to real-time prices. We present an analytic characterization of the optimal policy, along with the associated finite-horizon time-averaged value of storage. We also derive an analytical upperbound on the infinite-horizon time-averaged value of storage. This bound is valid for any achievable realization of prices when the support of the distribution is fixed, and highlights the dependence of the value of storage on ramp constraints and storage capacity. While the value of storage is a non-decreasing function of price volatility, due to the finite ramp rate, the value of storage saturates quickly as the capacity increases, regardless of volatility. To study the implications of the optimal policy, we first present computational experiments that suggest that optimal utilization of storage can, in expectation, induce a considerable amount of price elasticity near the average price, but little or no elasticity far from it. We then present a computational framework for understanding the behavior of storage as a function of price and the amount of stored energy, and for characterization of the buy/sell phase transition region in the price-state plane. Finally, we study the impact of market-based operation of storage on the required reserves, and show that the reserves may need to be expanded to accommodate market-based storage.


title: Modeling of octave-spanning Kerr frequency combs using a generalized mean-field Lugiato-Lefever model

abstract: A generalized Lugiato-Lefever equation is numerically solved with a Newton-Raphson method to model Kerr frequency combs. We obtain excellent agreement with past experiments, even for an octave-spanning comb. Simulations are much faster than with any other technique despite including more modes than ever before. Our study reveals that Kerr combs are associated with temporal cavity solitons and dispersive waves, and opens up new avenues for the understanding of Kerr comb formation.


title: Synthesis of ReN3 thin films by magnetron sputtering

abstract: Recently was reported a novel compound between rhenium and nitrogen, announced with ReN2 composition. This compound was synthesized by the high temperature and high pressure method. We found that the diffraction peaks of this compound are in agreement with the x-ray pattern of a rhenium-nitrogen film, under the assumption that the film is oriented on the substrate. The film was prepared by reactive magnetron sputtering, at room temperature, and deposited on a silicon wafer. From the analysis of the diffractograms it could be concluded that both materials share the same structure. By density functional calculation was found that the composition could be ReN3, instead of ReN2, as stated before. The ReN3 fits in the Ama2 (40) orthorhombic space group, and by the existence of N3 anions it should be categorized as an azide; that is, a nitrogen-rich compound. To reach high nitrogen concentrations by sputtering a crucial step is the target-poisoning. Under this regime of deposition is ensured that the compound is formed simultaneously on the substrate and the target. The poisoned target is rarely used because of a reduced sputtering yield, but as shall see, it can be used as a novel synthetic technique.


title: Optimal Auctions via the Multiplicative Weight Method

abstract: We show that the multiplicative weight update method provides a simple recipe for designing and analyzing optimal Bayesian Incentive Compatible (BIC) auctions, and reduces the time complexity of the problem to pseudo-polynomial in parameters that depend on single agent instead of depending on the size of the joint type space. We use this framework to design computationally efficient optimal auctions that satisfy ex-post Individual Rationality in the presence of constraints such as (hard, private) budgets and envy-freeness. We also design optimal auctions when buyers and a seller's utility functions are non-linear. Scenarios with such functions include (a) auctions with "quitting rights", (b) cost to borrow money beyond budget, (c) a seller's and buyers' risk aversion. Finally, we show how our framework also yields optimal auctions for variety of auction settings considered in Cai et al, Alaei et al, albeit with pseudo-polynomial running times.


title: Finite-Momentum Dimer Bound State in Spin-Orbit Coupled Fermi Gas

abstract: We investigate the two-body properties of a spin-1/2 Fermi gas subject to a spin-orbit coupling induced by laser fields. When an attractive s-wave interaction between unlike spins is present, the system may form a dimer bound state. Surprisingly, in the presence of a Zeeman field along the direction of the spin-orbit coupling, the bound state obtains finite center-of-mass mechanical momentum, whereas under the same condition but in the absence of the two-body interaction, the system has zero total momentum. This unusual result can be regarded as a consequence of the broken Galilean invariance by the spin-orbit coupling. Such a finite-momentum bound state will have profound effects on the many-body properties of the system.


title: Gaussification through decoherence

abstract: We investigate the loss of nonclassicality and non-Gaussianity of a single-mode state of the radiation field in contact with a thermal reservoir. The damped density matrix for a Fock-diagonal input is written using the Weyl expansion of the density operator. Analysis of the evolution of the quasiprobability densities reveals the existence of two successive characteristic times of the reservoir which are sufficient to assure the positivity of the Wigner function and, respectively, of the $P$ representation. We examine the time evolution of non-Gaussianity using three recently introduced distance-type measures. They are based on the Hilbert-Schmidt metric, the relative entropy, and the Bures metric. Specifically, for an $M$-photon-added thermal state, we obtain a compact analytic formula of the time-dependent density matrix that is used to evaluate and compare the three non-Gaussianity measures. We find a good consistency of these measures on the sets of damped states. The explicit damped quasiprobability densities are shown to support our general findings regarding the loss of negativities of Wigner and $P$ functions during decoherence. Finally, we point out that Gaussification of the attenuated field mode is accompanied by a nonmonotonic evolution of the von Neumann entropy of its state conditioned by the initial value of the mean photon number.


title: c-axis resistivity, pseudogap, superconductivity and Widom line in doped Mott insulators

abstract: Layered doped Mott insulators, such as the cuprates, show unusual temperature dependence of the resistivity. Intriguingly, the resistivity perpendicular to the CuO$_2$ planes, $\rho_c(T)$, shows both metallic ($d\rho_c/dT > 0$) and semi-conducting ($d\rho_c/dT<0$) behavior. We shed light on this puzzle by calculating $\rho_c$ for the two-dimensional Hubbard model within plaquette cellular dynamical mean-field theory and strong-coupling continuous-time quantum Monte Carlo as the impurity solver. The temperature, $T$, and doping, $\delta$, dependencies of $\rho_c$ are controlled by the first-order transition between pseudogap and correlated metal phases from which superconductivity can emerge. On the large doping side of the transition $\rho_c(T)$ is metallic, while on the low-doping side $\rho_c(T)$ changes from metallic to semi-conducting behavior with decreasing $T$. As a function of doping, the jump in $\rho_c$ across the first-order transition evolves into a sharp crossover at higher temperatures. This crossover coincides with the pseudogap temperature $T^*$ in the single-particle density of states, the spin susceptibility and other observables. Such coincidence in crossovers is expected along the continuation of the first-order transition into the super-critical regime, called the Widom line. This implies that not only the dynamic and the thermodynamic properties but also the DC transport in the normal state are governed by the hidden first-order transition. $\rho_c(T)$ has a high-temperature quasi-linear regime where it can exceed the Mott-Ioffe-Regel limit and when it has a minimum it is nearly parallel to the Widom line.


title: The Complexity of Optimal Mechanism Design

abstract: Myerson's seminal work provides a computationally efficient revenue-optimal auction for selling one item to multiple bidders. Generalizing this work to selling multiple items at once has been a central question in economics and algorithmic game theory, but its complexity has remained poorly understood. We answer this question by showing that a revenue-optimal auction in multi-item settings cannot be found and implemented computationally efficiently, unless ZPP contains P^#P. This is true even for a single additive bidder whose values for the items are independently distributed on two rational numbers with rational probabilities. Our result is very general: we show that it is hard to compute any encoding of an optimal auction of any format (direct or indirect, truthful or non-truthful) that can be implemented in expected polynomial time. In particular, under well-believed complexity-theoretic assumptions, revenue-optimization in very simple multi-item settings can only be tractably approximated. We note that our hardness result applies to randomized mechanisms in a very simple setting, and is not an artifact of introducing combinatorial structure to the problem by allowing correlation among item values, introducing combinatorial valuations, or requiring the mechanism to be deterministic (whose structure is readily combinatorial). Our proof is enabled by a flow-interpretation of the solutions of an exponential-size linear program for revenue maximization with an additional supermodularity constraint.


title: Study on the formation and the decomposition of AgN3 and a hypothetical compound ReN3 by using density functional calculations

abstract: We present a comparative study between ReN3 and AgN3 by using density functional theory. The ReN3 is a hypothetical compound proposed by us to interpret the Re to Re interplanar spacing of thin films grown by sputtering. Both, the AgN3 as the ReN3, are calculated as positive enthalpy compounds. The enthalpy might give a clue about the spontaneous decomposition of the solid form, but it cannot be interpreted as a restriction of its synthesizability. As from the calculated total-energy, we discuss the route for the formation of AgN3 starting from atomic species in aqueous solution. We propose that their synthesizability is conditioned by the energy of free nitrogen atoms, and the kinetics of reaction. We conclude that the intrinsic stability of a certain atomic arrangement depends only of the equilibrium of atomic forces, and not from the energy value associated with that structure.


title: Implementing Quantum Walks Using Orbital Angular Momentum of Classical Light

abstract: We present an implementation scheme for a quantum walk in the orbital angular momentum space of a laser beam. The scheme makes use of a ring interferometer, containing a quarter-wave plate and a q plate. This setup enables one to perform an arbitrary number of quantum walk steps. In addition, the classical nature of the implementation scheme makes it possible to observe the quantum walk evolution in real time. We use nonquantum entanglement of the laser beam's polarization with its orbital angular momentum to implement the quantum walk.


title: Convergence analysis for a primal-dual monotone + skew splitting algorithm with applications to total variation minimization

abstract: In this paper we investigate the convergence behavior of a primal-dual splitting method for solving monotone inclusions involving mixtures of composite, Lipschitzian and parallel sum type operators proposed by Combettes and Pesquet in [7]. Firstly, in the particular case of convex minimization problems, we derive convergence rates for the sequence of objective function values by making use of conjugate duality techniques. Secondly, we propose for the general monotone inclusion problem two new schemes which accelerate the sequences of primal and/or dual iterates, provided strong monotonicity assumptions for some of the involved operators are fulfilled. Finally, we apply the theoretical achievements in the context of different types of image restoration problems solved via total variation regularization.


title: Inflection Points and the Power Spectrum

abstract: Inflection point inflation generically includes a deviation from slow-roll when the inflaton approaches the inflection point. Such deviations are shown to be generated by transitions between singular trajectories. The effects on the power spectrum are studied within the context of universality classes for small-field models. These effects are shown to scale with universality parameters, and can explain the anomalously low power at large scales observed in the CMB. The reduction of power is related to the inflection point's basin of attraction. Implications for the likelihood of inflation are discussed.


title: A Nonlinear Force-Free Magnetic Field Approximation Suitable for Fast Forward-Fitting to Coronal Loops. III. The Free Energy

abstract: An analytical approximation of a nonlinear force-free magnetic field (NLFFF) solution was developed in Paper I, while a numerical code that performs fast forward-fitting of this NLFFF approximation to a line-of-sight magnetogram and coronal 3D loops has been described and tested in Paper II. Here we calculate the free magnetic energy $E_{\rm free}=E_{\rm N}-E_{\rm P}$, i.e. the difference of the magnetic energies between the nonpotential field and the potential field. A second method to estimate the free energy is obtained from the mean misalignment angle change $\Delta\mu=\mu_{\rm P}-\mu_{\rm N}$ between the potential and nonpotential field, which scales as $E_{\rm free}/E_{\rm P} \approx \tan^2{(\Delta\mu)}$. For four active regions observed with STEREO in 2007 we find free energies in the range of $q_{\rm free}=(E_{\rm free}/E_{\rm P}) \approx 1%-10%$, with an uncertainty of less than $\pm 2%$ between the two methods, while the free energies obtained from 11 other NLFFF codes exhibit a larger scatter of order $\approx\pm10%$. We find also a correlation between the free magnetic energy and the GOES flux of the largest flare that occurred during the observing period, which can be quantified by an exponential relationship, $F_{\rm GOES} \propto \exp{(q_{\rm free}/0.015)}$, implying an exponentiation of the dissipated currents.


title: Three Dimensional Radiative Hydrodynamical Simulations of the Highly Irradiated Short Period Exoplanet HD189733b

abstract: We present a detailed three-dimensional radiative-hydrodynamical simulation of the well known irradiated exoplanet HD189733b. Our model solves the fully compressible Navier-Stokes equations coupled to wavelength-dependent radiative transfer throughout the entire planetary envelope. We provide detailed comparisons between the extensive observations of this system and predictions calculated directly from the numerical models. The atmospheric dynamics is characterized by supersonic winds that fairly efficiently advect energy from the dayside to the nightside. Super-rotating equatorial jets form for a wide range of pressures from 10^-5 to 10 bars while counter rotating jets form at higher latitudes. Calculated transit spectrum agree well with the data from the infrared to the UV including the strong Rayleigh scattering seen at short wavelength, though we slightly under-predict the observations at wavelengths shorter then ~0.6 microns. Our predicted emission spectrum agrees remarkably well at 5.8 and 8 microns, but slightly over-predicts the emission at 3.6 and 4.5 microns when compared to the latest analysis by Knutson et. al (2012). Our simulated IRAC phasecurves agree fairly well with the amplitudes of variations, shape, and phases of minimum and maximum flux. However, we over-predict the peak amplitude at 3.6 and 4.5 microns, and slightly under-predict the location of the phasecurve maximum and minimum. These simulations include, for the first time in a multi-dimensional simulation, a strong Rayleigh scattering component to the absorption opacity, necessary to explain observations in the optical and UV. The agreement between our models and observations suggest that including the effects of condensates in simulations as the dominant form of opacity will be very important in future models.


title: X-Ray Determination of the Variable Rate of Mass Accretion onto TW Hydrae

abstract: Diagnostics of electron temperature (T_e), electron density (n_e), and hydrogen column density (N_H) from the Chandra High Energy Transmission Grating spectrum of He-like Ne IX in TW Hydrae (TW Hya), in conjunction with a classical accretion model, allow us to infer the accretion rate onto the star directly from measurements of the accreting material. The new method introduces the use of the absorption of Ne IX lines as a measure of the column density of the intervening, accreting material. On average, the derived mass accretion rate for TW Hya is 1.5 x 10^{-9} M_{\odot} yr^{-1}, for a stellar magnetic field strength of 600 Gauss and a filling factor of 3.5%. Three individual Chandra exposures show statistically significant differences in the Ne IX line ratios, indicating changes in N_H, T_e, and n_e by factors of 0.28, 1.6, and 1.3, respectively. In exposures separated by 2.7 days, the observations reported here suggest a five-fold reduction in the accretion rate. This powerful new technique promises to substantially improve our understanding of the accretion process in young stars.


title: Waveguide-QED-Based Photonic Quantum Computation

abstract: We propose a new scheme for quantum computation using flying qubits--propagating photons in a one-dimensional waveguide--interacting with matter qubits. Photon-photon interactions are mediated by the coupling to a three- or four-level system, based on which photon-photon \pi-phase gates (Controlled-NOT) can be implemented for universal quantum computation. We show that high gate fidelity is possible given recent dramatic experimental progress in superconducting circuits and photonic-crystal waveguides. The proposed system can be an important building block for future on-chip quantum networks.


title: Network models of dissolution of porous media

abstract: We investigate the chemical dissolution of porous media using a network model in which the system is represented as a series of interconnected pipes with the diameter of each segment increasing in proportion to the local reactant consumption. Moreover, the topology of the network is allowed to change dynamically during the simulation: as the diameters of the eroding pores become comparable with the interpore distances, the pores are joined together thus changing the interconnections within the network. With this model, we investigate different growth regimes in an evolving porous medium, identifying the mechanisms responsible for the emergence of specific patterns. We consider both the random and regular network and study the effect of the network geometry on the patterns. Finally, we consider practically important problem of finding an optimum flow rate that gives a maximum increase in permeability for a given amount of reactant.


title: Pseudobulge Formation as a Dynamical Rather than a Secular Process

abstract: We investigate the formation and evolution of the pseudobulge in "Eris", a high-resolution N-body + smoothed particle hydrodynamics (SPH) cosmological simulation that successfully reproduces a Milky Way-like massive late-type spiral in a cold dark matter (LCDM) universe. At the present epoch, Eris has a virial mass Mvir=8x10^11 Msun, a photometric stellar mass M*=3.2x10^10 Msun, a bulge-to-total ratio B/T = 0.26, and a weak nuclear bar. We find that the bulk of the pseudobulge forms quickly at high redshift via a combination of non-axisymmetric disk instabilities and tidal interactions or mergers both occurring on dynamical timescales, not through slow secular processes at lower redshift. Its subsequent evolution is not strictly secular either, and is closely intertwined with the evolution of the stellar bar. In fact, the structure that we recognize as a pseudobulge today evolves from a stellar bar that formed at high redshift, was destroyed by minor mergers at z~3, reformed shortly after, and weakened again following a steady gas inflow at z<1. The gradual dissolution of the bar ensues at z=1 and continues until the present without increasing the stellar velocity dispersion in the inner regions. In this scenario the pseudobulge is not a separate component from the inner disk in terms of formation path, rather it is the first step in the inside-out formation of the baryonic disk, in agreement with the fact that pseudobulges of massive spiral galaxies have typically a dominant old stellar population. If our simulations do indeed reproduce the formation mechanisms of massive spirals, then the progenitors of late-type galaxies should have strong bars and small photometric pseudobulges at high redshift.


title: Even Subdivision-Factors of Cubic Graphs

abstract: We call a set $\mathcal S$ of graphs an "even subdivison-factor" of a cubic graph $G$ if $G$ contains a spanning subgraph $H$ such that every component of $H$ has an even number of vertices and is a subdivision of an element of $\mathcal S$. We show that any set of 2-connected graphs which is an even subdivison-factor of every 3-connected cubic graph, satisfies certain properties. As a consequence, we disprove a conjecture which was stated in an attempt to solve the circuit double cover conjecture.


title: Exact Lagrangian immersions with one double point revisited

abstract: We study exact Lagrangian immersions with one double point of a closed orientable manifold K into n-complex-dimensional Euclidean space. Our main result is that if the Maslov grading of the double point does not equal 1 then K is homotopy equivalent to the sphere, and if, in addition, the Lagrangian Gauss map of the immersion is stably homotopic to that of the Whitney immersion, then K bounds a parallelizable (n+1)-manifold. The hypothesis on the Gauss map always holds when n=2k or when n=8k-1. The argument studies a filling of K obtained from solutions to perturbed Cauchy-Riemann equations with boundary on the image f(K) of the immersion. This leads to a new and simplified proof of some of the main results of arXiv:1111.5932, which treated Lagrangian immersions in the case n=2k by applying similar techniques to a Lagrange surgery of the immersion, as well as to an extension of these results to the odd-dimensional case.


title: Blind Signal Separation in the Presence of Gaussian Noise

abstract: A prototypical blind signal separation problem is the so-called cocktail party problem, with n people talking simultaneously and n different microphones within a room. The goal is to recover each speech signal from the microphone inputs. Mathematically this can be modeled by assuming that we are given samples from an n-dimensional random variable X=AS, where S is a vector whose coordinates are independent random variables corresponding to each speaker. The objective is to recover the matrix A^{-1} given random samples from X. A range of techniques collectively known as Independent Component Analysis (ICA) have been proposed to address this problem in the signal processing and machine learning literature. Many of these techniques are based on using the kurtosis or other cumulants to recover the components. In this paper we propose a new algorithm for solving the blind signal separation problem in the presence of additive Gaussian noise, when we are given samples from X=AS+\eta, where \eta is drawn from an unknown, not necessarily spherical n-dimensional Gaussian distribution. Our approach is based on a method for decorrelating a sample with additive Gaussian noise under the assumption that the underlying distribution is a linear transformation of a distribution with independent components. Our decorrelation routine is based on the properties of cumulant tensors and can be combined with any standard cumulant-based method for ICA to get an algorithm that is provably robust in the presence of Gaussian noise. We derive polynomial bounds for the sample complexity and error propagation of our method.


title: Bayesian Learning and Predictability in a Stochastic Nonlinear Dynamical Model

abstract: Bayesian inference methods are applied within a Bayesian hierarchical modelling framework to the problems of joint state and parameter estimation, and of state forecasting. We explore and demonstrate the ideas in the context of a simple nonlinear marine biogeochemical model. A novel approach is proposed to the formulation of the stochastic process model, in which ecophysiological properties of plankton communities are represented by autoregressive stochastic processes. This approach captures the effects of changes in plankton communities over time, and it allows the incorporation of literature metadata on individual species into prior distributions for process model parameters. The approach is applied to a case study at Ocean Station Papa, using Particle Markov chain Monte Carlo computational techniques. The results suggest that, by drawing on objective prior information, it is possible to extract useful information about model state and a subset of parameters, and even to make useful long-term forecasts, based on sparse and noisy observations.


title: Space-time diagrammatics

abstract: We introduce a new class of two-dimensional diagrams, the \emph{projection diagrams}, as a tool to visualize the global structure of space-times. We construct the diagrams for several metrics of interest, including the Kerr-Newman - (anti) de Sitter family, with or without cosmological constant, and the Emparan-Reall black rings.


title: Indestructibility of compact spaces

abstract: In this article we investigate which compact spaces remain compact under countably closed forcing. We prove that, assuming the Continuum Hypothesis, the natural generalizations to $\omega_1$-sequences of the selection principle and topological game versions of the Rothberger property are not equivalent, even for compact spaces. We also show that Tall and Usuba's "$\aleph_1$-Borel Conjecture" is equiconsistent with the existence of an inaccessible cardinal.


title: A Novel Model Order Reduction Approach for Navier-Stokes Equations at High Reynolds Number

abstract: A new approach to model order reduction of the Navier-Stokes equations at high Reynolds number is proposed. Unlike traditional approaches, this method does not rely on empirical turbulence modeling or modification of the Navier-Stokes equations. It provides spatial basis functions different from the usual proper orthogonal decomposition basis function in that, in addition to optimally representing the training data set, the new basis functions also provide stable and accurate reduced-order models. The proposed approach is illustrated with two test cases: two-dimensional flow inside a square lid-driven cavity and a two-dimensional mixing layer.


title: Review of measurements of Kerr nonlinearities in lithium niobate: the role of the delayed Raman response

abstract: We review measurements in the literature regarding Kerr nonlinearity of lithium niobate. Of particular interest is to separate the value of the electronic Kerr nonlinearity and the strength of the delayed Raman nonlinearity. The spectral shape of the Raman response in LN is discussed and we restrict ourselves to the case where the pump is polarized along the optical axis. This configuration is important for so-called type 0 noncritical interaction, which is used in QPM devices and in recent supercontinuum generation and few-cycle soliton compression experiments.


title: Inverse problems in approximate uniform generation

abstract: We initiate the study of \emph{inverse} problems in approximate uniform generation, focusing on uniform generation of satisfying assignments of various types of Boolean functions. In such an inverse problem, the algorithm is given uniform random satisfying assignments of an unknown function $f$ belonging to a class $\C$ of Boolean functions, and the goal is to output a probability distribution $D$ which is $\epsilon$-close, in total variation distance, to the uniform distribution over $f^{-1}(1)$. Positive results: We prove a general positive result establishing sufficient conditions for efficient inverse approximate uniform generation for a class $\C$. We define a new type of algorithm called a \emph{densifier} for $\C$, and show (roughly speaking) how to combine (i) a densifier, (ii) an approximate counting / uniform generation algorithm, and (iii) a Statistical Query learning algorithm, to obtain an inverse approximate uniform generation algorithm. We apply this general result to obtain a poly$(n,1/\eps)$-time algorithm for the class of halfspaces; and a quasipoly$(n,1/\eps)$-time algorithm for the class of $\poly(n)$-size DNF formulas. Negative results: We prove a general negative result establishing that the existence of certain types of signature schemes in cryptography implies the hardness of certain inverse approximate uniform generation problems. This implies that there are no {subexponential}-time inverse approximate uniform generation algorithms for 3-CNF formulas; for intersections of two halfspaces; for degree-2 polynomial threshold functions; and for monotone 2-CNF formulas. Finally, we show that there is no general relationship between the complexity of the "forward" approximate uniform generation problem and the complexity of the inverse problem for a class $\C$ -- it is possible for either one to be easy while the other is hard.


title: Dispersion law for a one-dimensional weakly interacting Bose gas with zero boundary conditions

abstract: From the time-dependent Gross equation, we find the quasiparticle dispersion law for a one-dimensional weakly interacting Bose gas with a non-point interatomic potential and zero boundary conditions (BCs). The result coincides with the dispersion law for periodic BCs, i.e. the Bogolyubov law $E_{B}(k) = \sqrt{\left (\frac{\hbar^{2} k^2}{2m}\right )^{2} + n_{0}\nu(k)\frac{\hbar^2 k^2}{m}}$. In the case of periodic BCs, the dispersion law can be easily derived from Gross' equation. However, for zero BCs, the analysis is not so simple.


title: Coherent feedback that beats all measurement-based feedback protocols

abstract: We show that when the speed of control is bounded, there is a widely applicable minimal-time control problem for which a coherent feedback protocol is optimal, and is faster than all measurement-based feedback protocols, where the latter are defined in a strict sense. The superiority of the coherent protocol is due to the fact that it can exploit a geodesic path in Hilbert space, a path that measurement-based protocols cannot follow.


title: Bahadur efficiency of nonparametric test for independence based on $L_1$-error

abstract: We introduce new test statistic to test the independence of two multi-dimensional random variables. Based on the $L_1$-distance and the historgram density estimation method, the test is compared via Bahadur relative efficiency to several tests available in the literature. It arises that our test reaches better performances than a number of usual tests among whom we cite the Kolmogorov-Smirnov test. Beforehand, large deviation result is stated for the associated statistic. The local asymptotic optimality relative to the test is also studied.


title: Stringy Models of Modified Gravity: Space-time defects and Structure Formation

abstract: Starting from microscopic models of space-time foam, based on brane universes propagating in bulk space-times populated by D0-brane defects ("D-particles"), we arrive at effective actions used by a low-energy observer on the brane world to describe his/her observations of the Universe. These actions include, apart from the metric tensor field, also scalar (dilaton) and vector fields, the latter describing the interactions of low-energy matter on the brane world with the recoiling point-like space-time defect (D-particle). The vector field is proportional to the recoil velocity of the D-particle and as such it satisfies a certain constraint. The vector breaks locally Lorentz invariance, which however is assumed to be conserved on average in a space-time foam situation, involving the interaction of matter with populations of D-particle defects. In this paper we demonstrate that, already at the end of the radiation era, the (constrained) vector field associated with the recoil of the defects provides the seeds for a growing mode in the evolution of the Universe. Such a growing mode survives during the matter dominated era, provided the variance of the D-particle recoil velocities on the brane is larger than a critical value.


title: An Alternative Proof and Generalization of Ferrero's Computations of Iwasawa \lambda-Invariants

abstract: We prove a slight generalization of Iwasawa's `Riemann-Hurwitz' formula for number fields and use it to generalize Ferrero's and Kida's well-known computations of Iwasawa \lambda-invariants for the cyclotomic Z_2-extensions of imaginary quadratic number fields. In particular, we show that if p is a Fermat prime, then similar computations of Iwasawa \lambda-invariants hold for certain imaginary quadratic extensions of the unique subfield k of Q(\zeta_{p^2}) such that [k:Q] = p. In fact, we actually prove more by explicitly computing cohomology groups of principal ideals. The computation of lambda invariants obtained is a special case of a much more general result concerning relative lambda invariants for cyclotomic Z_2-extensions of CM number fields due to Yuji Kida. However, the approach used here significantly differs from that of Kida, and the intermediate computations of cohomology groups found here do not hold in Kida's more general setting.


title: Maximum Distance Separable Codes for Symbol-Pair Read Channels

abstract: We study (symbol-pair) codes for symbol-pair read channels introduced recently by Cassuto and Blaum (2010). A Singleton-type bound on symbol-pair codes is established and infinite families of optimal symbol-pair codes are constructed. These codes are maximum distance separable (MDS) in the sense that they meet the Singleton-type bound. In contrast to classical codes, where all known q-ary MDS codes have length O(q), we show that q-ary MDS symbol-pair codes can have length \Omega(q^2). In addition, we completely determine the existence of MDS symbol-pair codes for certain parameters.