title: Improved Measures of Integrated Information

abstract: Although there is growing interest in measuring integrated information in computational and cognitive systems, current methods for doing so in practice are computationally unfeasible. Existing and novel integration measures are investigated and classified by various desirable properties. A simple taxonomy of Phi-measures is presented where they are each characterized by their choice of factorization method (5 options), choice of probability distributions to compare (3x4 options) and choice of measure for comparing probability distributions (7 options). When requiring the Phi-measures to satisfy a minimum of attractive properties, these hundreds of options reduce to a mere handful, some of which turn out to be identical. Useful exact and approximate formulas are derived that can be applied to real-world data from laboratory experiments without posing unreasonable computational demands.


title: Certification of Boson Sampling Devices with Coarse-Grained Measurements

abstract: A boson sampling device could efficiently sample from the output probability distribution of noninteracting bosons undergoing many-body interference. This problem is not only classically intractable, but its solution is also believed to be classically unverifiable. Hence, a major difficulty in experiment is to ensure a boson sampling device performs correctly. We present an experimental friendly scheme to extract useful and robust information from the quantum boson samplers based on coarse-grained measurements. The procedure can be applied to certify the equivalence of boson sampling devices while ruling out alternative fraudulent devices. We perform numerical simulations to demonstrate the feasibility of the method and consider the effects of realistic noise. Our approach is expected to be generally applicable to other many-body certification tasks beyond the boson sampling problem.


title: Evidence for High-Frequency QPOs with a 3:2 Frequency Ratio from a 5000 Solar Mass Black Hole

abstract: Following the discovery of 3:2 resonance quasi-periodic oscillations (QPOs) in M82X-1 (Pasham et al. 2014), we have constructed power density spectra (PDS) of all 15 (sufficiently long) {\it XMM-Newton} observations of the ultraluminous X-ray source NGC1313X-1 ($L_{X}$ $\approx$ 2$\times$10$^{40}$ erg/sec). We detect a strong QPO at a frequency of 0.29$\pm$0.01 Hz in data obtained on 2012 December 16. Subsequent searching of all the remaining observations for a 3:2/2:3 frequency pair revealed a feature at 0.46$\pm$0.02 Hz on 2003 Dec 13 (frequency ratio of 1.59$\pm$0.09). The global significance of the 0.29 Hz feature considering all frequencies between 0.1 and 4 Hz is $>$ 3.5 $\sigma$. The significance of the 0.46$\pm$0.02 Hz QPO is $>$ 3.5$\sigma$ for a search at 2/3 and 3/2 of 0.29 Hz. We also detect lower frequency QPOs (32.9$\pm$2.6 and 79.7$\pm$1.2 mHz). All the QPOs are super-imposed on a continuum consisting of flat-topped, band-limited noise, breaking into a power-law at a frequency of 16$\pm$3 mHz and white noise at $\gtrsim$ 0.1 Hz. NGC1313X-1's PDS is analogous to stellar-mass black holes' (StMBHs) PDS in the so-called steep power-law state, but with the respective frequencies (both QPOs and break frequencies) scaled down by a factor of $\sim$ 1000. Using the inverse mass-to-high-frequency QPO scaling of StMBHs, we estimate NGC1313X-1's black hole mass to be 5000$\pm$1300 $M_{\odot}$, consistent with an inference from the scaling of the break frequency. However, the implied Eddington ratio, L$_{Edd}$ $>$ 0.03$\pm$0.01, is significantly lower compared to StMBHs in the steep power-law state (L$_{Edd}$ $\gtrsim$ 0.2).


title: The Relation between Galaxy Structure and Spectral Type: Implications for the Buildup of the Quiescent Galaxy Population at 0.5<z<2.0

abstract: We present the relation between galaxy structure and spectral type, using a K-selected galaxy sample at 0.5<z<2.0. Based on similarities between the UV-to-NIR spectral energy distributions, we classify galaxies into 32 spectral types. The different types span a wide range in evolutionary phases, and thus -- in combination with available CANDELS/F160W imaging -- are ideal to study the structural evolution of galaxies. Effective radii (R_e) and Sersic parameters (n) have been measured for 572 individual galaxies, and for each type, we determine R_e at fixed stellar mass by correcting for the mass-size relation. We use the rest-frame U-V vs. V-J diagram to investigate evolutionary trends. When moving into the direction perpendicular to the star-forming sequence, in which we see the Halpha equivalent width and the specific star formation rate (sSFR) decrease, we find a decrease in R_e and an increase in n. On the quiescent sequence we find an opposite trend, with older redder galaxies being larger. When splitting the sample into redshift bins, we find that young post-starburst galaxies are most prevalent at z>1.5 and significantly smaller than all other galaxy types at the same redshift. This result suggests that the suppression of star formation may be associated with significant structural evolution at z>1.5. At z<1, galaxy types with intermediate sSFRs (10^{-11.5}-10^{-10.5} yr^-1) do not have post-starburst SED shapes. These galaxies have similar sizes as older quiescent galaxies, implying that they can passively evolve onto the quiescent sequence, without increasing the average size of the quiescent galaxy population.


title: The SCUBA-2 Cosmology Legacy Survey: Multi-wavelengths counterparts to 10$^3$ submillimeter galaxies in the UKIDSS-UDS field

abstract: We present multiwavelength identifications for the counterparts of 1088 submillimeter sources detected at 850$\mu$m in the SCUBA-2 Cosmology Legacy Survey study of the UKIDSS-UDS field. By utilising an ALMA pilot study on a subset of our bright SCUBA-2 sample as a training set, along with the deep optical-near-infrared data available in this field, we develop a novel technique, Optical-IR Triple Color (OIRTC), using $z-K$, $K-[3.6]$, $[3.6]-[4.5]$ colors to select the candidate submillimeter galaxy (SMG) counterparts. By combining radio identification and the OIRTC technique, we find counterpart candidates for 80% of the Class = 1 $\geq4\,\sigma$ SCUBA-2 sample, defined as those that are covered by both radio and OIR imaging and the base sample for our scientific analyses. Based on the ALMA training set, we expect the accuracy of these identifications to be $82\pm20$%, with a completeness of $69\pm16$%, essentially as accurate as the traditional $p$-value technique but with higher completeness. We find that the fraction of SCUBA-2 sources having candidate counterparts is lower for fainter 850$\mu$m sources, and we argue that for follow-up observations sensitive to SMGs with $S_{850}\gtrsim 1$ mJy across the whole ALMA beam, the fraction with multiple counterparts is likely to be $>40$% for SCUBA-2 sources at $S_{850} \gtrsim 4$ mJy. We find that the photometric redshift distribution for the SMGs is well fit by a lognormal distribution, with a median redshift of $z=2.3\pm0.1$. After accounting for the sources without any radio and/or OIRTC counterpart, we estimate the median redshift to be $z=2.6\pm0.1$ for SMGs with $S_{850} >1$ mJy. We also use this new large sample to study the clustering of SMGs and the the far-infrared properties of the unidentified submillimeter sources by stacking their Herschel SPIRE far-infrared emission.


title: Weakened magnetic braking as the origin of anomalously rapid rotation in old field stars

abstract: A knowledge of stellar ages is crucial for our understanding of many astrophysical phenomena, and yet ages can be difficult to determine. As they become older, stars lose mass and angular momentum, resulting in an observed slowdown in surface rotation. The technique of 'gyrochronology' uses the rotation period of a star to calculate its age. However, stars of known age must be used for calibration, and, until recently, the approach was untested for old stars (older than 1 gigayear, Gyr). Rotation periods are now known for stars in an open cluster of intermediate age (NGC 6819; 2.5 Gyr old), and for old field stars whose ages have been determined with asteroseismology. The data for the cluster agree with previous period-age relations, but these relations fail to describe the asteroseismic sample. Here we report stellar evolutionary modelling, and confirm the presence of unexpectedly rapid rotation in stars that are more evolved than the Sun. We demonstrate that models that incorporate dramatically weakened magnetic braking for old stars can---unlike existing models---reproduce both the asteroseismic and the cluster data. Our findings might suggest a fundamental change in the nature of ageing stellar dynamos, with the Sun being close to the critical transition to much weaker magnetized winds. This weakened braking limits the diagnostic power of gyrochronology for those stars that are more than halfway through their main-sequence lifetimes.


title: A Bosonic Analog of a Topological Dirac Semi-Metal: Effective Theory, Neighboring Phases, and Wire Construction

abstract: We construct a bosonic analog of a two-dimensional topological Dirac Semi-Metal (DSM). The low-energy description of the most basic 2D DSM model consists of two Dirac cones at positions $\pm\mathbf{k}_0$ in momentum space. The local stability of the Dirac cones is guaranteed by a composite symmetry $Z_2^{\mathcal{TI}}$, where $\mathcal{T}$ is time-reversal and $\mathcal{I}$ is inversion. This model also exhibits interesting time-reversal and inversion symmetry breaking electromagnetic responses. In this work we construct a bosonic version by replacing each Dirac cone with a copy of the $O(4)$ Nonlinear Sigma Model (NLSM) with topological theta term and theta angle $\theta=\pm \pi$. One copy of this NLSM also describes the gapless surface termination of the 3D Bosonic Topological Insulator (BTI). We compute the time-reversal and inversion symmetry breaking electromagnetic responses for our model and show that they are twice the value one gets in the DSM case matching what one might expect from, for example, a bosonic Chern insulator. We also investigate the stability of the BSM model and find that the composite $Z_2^{\mathcal{TI}}$ symmetry again plays an important role. Along the way we clarify many aspects of the surface theory of the BTI including the electromagnetic response, the charges and statistics of vortex excitations, and the stability to symmetry-allowed perturbations. We briefly comment on the relation between the various descriptions of the $O(4)$ NLSM with $\theta=\pi$ used in this paper (a dual vortex description and a description in terms of four massless fermions) and the recently proposed dual description of the BTI surface in terms of $2+1$ dimensional Quantum Electrodynamics with two flavors of fermion ($N=2$ QED$_3$). In a set of four Appendixes we review some of the tools used in the paper, and also derive some of the more technical results.


title: ADE Spectral Networks

abstract: We introduce a new perspective and a generalization of spectral networks for 4d $\mathcal{N}=2$ theories of class $\mathcal{S}$ associated to Lie algebras $\mathfrak{g} = \textrm{A}_n$, $\textrm{D}_n$, $\textrm{E}_{6}$, and $\textrm{E}_{7}$. Spectral networks directly compute the BPS spectra of 2d theories on surface defects coupled to the 4d theories. A Lie algebraic interpretation of these spectra emerges naturally from our construction, leading to a new description of 2d-4d wall-crossing phenomena. Our construction also provides an efficient framework for the study of BPS spectra of the 4d theories. In addition, we consider novel types of surface defects associated with minuscule representations of $\mathfrak{g}$.


title: Generalized Gravitational Entropy for Warped Anti-de Sitter Space

abstract: For spacetimes that are not asymptotic to anti-de Sitter Space (non AAdS), we adapt the Lewkowycz-Maldacena procedure to find the holographic entanglement entropy. The key observation, which to our knowledge is not very well appreciated, is that asymptotic boundary conditions play an essential role on extending the replica trick to the bulk. For non AAdS, we expect the following three main modifications: (1) the expansion near the special surface has to be compatible with the asymptotic expansion; (2) periodic conditions are imposed to coordinates on the phase space with diagonalized symplectic structure, not to all fields appearing in the action; (3) evaluating the entanglement functional using the boundary term method amounts to evaluating the presymplectic structure at the special surface, where some additional exact form may contribute. An explicit calculation is carried out for three-dimensional warped anti-de Sitter spacetime (WAdS3) in a consistent truncation of string theory, the so-called S-dual dipole theory. It turns out that the generalized gravitational entropy in WAdS3 is captured by the least action of a charged particle in WAdS3 space, or equivalently, by the geodesic length in an auxiliary AdS3. Consequently, the bulk calculation agrees with the CFT results, providing another piece of evidence for the WAdS3/CFT2 correspondence.


title: Modelling the circumstellar medium in RS Ophiuchi and its link to Type Ia supernovae

abstract: Recent interpretations of narrow, variable absorption lines detected in some Type Ia supernovae suggest that their progenitors are surrounded by dense, circumstellar material. Similar variations detected in the symbiotic recurrent nova system RS Oph, which undergoes thermonuclear outbursts every ~20 years, making it an ideal candidate to investigate the origin of these lines. To this end, we present simulations of multiple mass transfer-nova cycles in RS Oph. We find that the quiescent mass transfer produces a dense, equatorial outflow, i.e., concentrated towards the binary orbital plane, and an accretion disc forms around the white dwarf. The interaction of a spherical nova outburst with these aspherical circumstellar structures produces a bipolar outflow, similar to that seen in HST imaging of the 2006 outburst. In order to produce an ionization structure that is consistent with observations, a mass-loss rate of $5 \times 10^{-7}\,\mathrm{M}_{\odot}\,\mathrm{yr}^{-1}$ from the red giant is required. The simulations also produce a polar accretion flow, which may explain the broad wings of the quiescent H {\alpha} line and hard X-rays. By comparing simulated absorption line profiles to observations of the 2006 outburst, we are able to determine which components arise in the wind and which are due to the novae. We explore the possible behaviour of absorption line profiles as they may appear should a supernova occur in a system like RS Oph. Our models show similarities to supernovae like SN 2006X, but require a high mass-loss rate, $\dot{M} \sim 10^{-6}$ to $10^{-5}\,\mathrm{M_\odot}\,\mathrm{yr}^{-1}$, to explain the variability in SN 2006X.


title: Light Amplification in Semiconductor-Superconductor Structures

abstract: We study a new effect of Cooper-pair-based two-photon gain in semiconductor-superconductor structures, showing broadband enhancement of ultrafast two-photon amplification. We further show that with the superconducting enhancement, at moderately high seed intensities the two-photon gain contribution approaches that of the one-photon gain. A full quantum-optical model of singly- and fully-stimulated two-photon emission is developed. Our results provide new insights on nonlinear light-matter interaction in the superconducting state, including the possibility of coherent control in two-photon semiconductor-superconductor devices. The theoretically demonstrated effects can have important implications in optoelectronics and in coherent control applications.


title: Number of trials required to estimate a free-energy difference, using fluctuation relations

abstract: The difference Delta F between free energies has applications in biology, chemistry, and pharmacology. The value of Delta F can be estimated from experiments or simulations, via fluctuation theorems developed in statistical mechanics. Calculating the error in a Delta F estimate is difficult. Worse, atypical trials dominate estimates. How many trials one should perform was estimated roughly in [Jarzynski, Phys. Rev. E 73, 046105 (2006)]. We enhance the approximation with information-theoretic strategies: We quantify "dominance" with a tolerance parameter chosen by the experimenter or simulator. We bound the number of trials one should expect to perform, using the order-infinity Renyi entropy. The bound can be estimated if one implements the "good practice" of bidirectionality, known to improve estimates of Delta F. Estimating Delta F from this number of trials leads to an error that we bound approximately. Numerical experiments on a weakly interacting dilute classical gas support our analytical calculations.


title: First Search for an X-ray -- Optical Reverberation Signal in an Ultraluminous X-ray Source

abstract: Using simultaneous optical (VLT/FORS2) and X-ray (XMM-Newton) data of NGC 5408, we present the first ever attempt to search for a reverberation signal in an ultraluminous X-ray source (NGC 5408 X-1). The idea is similar to AGN broad line reverberation mapping where a lag measurement between the X-ray and the optical flux combined with a Keplerian velocity estimate should enable us to weigh the central compact object. We find that although NGC 5408 X-1's X-rays are variable on a timescale of a few hundred seconds (RMS of 9.0$\pm$0.5%), the optical emission does not show any statistically significant variations. We set a 3$\sigma$ upper limit on the RMS optical variability of 3.3%. The ratio of the X-ray to the optical variability is an indicator of X-ray reprocessing efficiency. In X-ray binaries, this ratio is roughly 5. Assuming a similar ratio for NGC 5408 X-1, the expected RMS optical variability is $\approx$2% which is still a factor of roughly two lower than what was possible with the VLT observations in this study. We find marginal evidence (3$\sigma$) for optical variability on a $\sim$ 24 hour timescale. Our results demonstrate that such measurements can be made, but photometric conditions, low sky background levels and longer simultaneous observations will be required to reach the optical variability levels similar to X-ray binaries.


title: Positivity Sectors and the Amplituhedron

abstract: We initiate a detailed investigation into the assembly of simple amplituhedron-like building blocks to obtain spaces of physical interest. In particular, we describe the geometric process through which the building blocks, which we call positivity sectors, glue together to form the desired geometries. Positivity sectors are seen to naturally segment the space describing the Lth power of the one-loop amplitude. In this way, we obtain a good understanding of how the geometric complexity of the building blocks can be washed out in the formation of larger spaces. Conversely, the tools we develop allow us to form spaces of ever greater complexity, a process which is crucial to the construction of the amplituhedron from its triangulations, which remains an important open question. We present the full boundary structure of all positivity sectors related to the three-loop amplituhedron. We also construct a practical algorithm that achieves the desired geometric assembly of positivity sectors, and make available supporting Mathematica files containing the full boundary structure of all positivity sectors at three loops.


title: High-Dimensional Stochastic Design Optimization by Adaptive-Sparse Polynomial Dimensional Decomposition

abstract: This paper presents a novel adaptive-sparse polynomial dimensional decomposition (PDD) method for stochastic design optimization of complex systems. The method entails an adaptive-sparse PDD approximation of a high-dimensional stochastic response for statistical moment and reliability analyses; a novel integration of the adaptive-sparse PDD approximation and score functions for estimating the first-order design sensitivities of the statistical moments and failure probability; and standard gradient-based optimization algorithms. New analytical formulae are presented for the design sensitivities that are simultaneously determined along with the moments or the failure probability. Numerical results stemming from mathematical functions indicate that the new method provides more computationally efficient design solutions than the existing methods. Finally, stochastic shape optimization of a jet engine bracket with 79 variables was performed, demonstrating the power of the new method to tackle practical engineering problems.


title: Expansion of a Valley-Polarized Exciton Cloud in a 2D Heterostructure

abstract: Heterostructures comprising different monolayer semiconductors provide a new system for fundamental science and device technologies, such as in the emerging field of valleytronics. Here, we realize valley-specific interlayer excitons in monolayer WSe2-MoSe2 vertical heterostructures. We create interlayer exciton spin-valley polarization by circularly polarized optical pumping and determine a valley lifetime of 40 nanoseconds. This long-lived polarization enables the visualization of the expansion of a valley-polarized exciton cloud over several microns. The spatial pattern of the polarization evolves into a ring with increasing exciton density, a manifestation of valley exciton exchange interactions. Our work lays a foundation for quantum optoelectronics and valleytronics based on interlayer excitons in van der Waals heterostructures.


title: Resolved Gas Interior to the Dust Rings of the HD 141569 Disk

abstract: The disk around HD 141569 is one of a handful of systems whose weak infrared emission is consistent with a debris disk, but still has a significant reservoir of gas. Here we report spatially resolved mm observations of the CO(3-2) and CO(1-0) emission as seen with the SMA and CARMA. We find that the excitation temperature for CO is lower than expected from cospatial blackbody grains, similar to previous observations of analogous systems, and derive a gas mass that lies between that of gas-rich primordial disks and gas-poor debris disks. The data also indicate a large inner hole in the CO gas distribution and an outer radius that lies interior to the outer scattered light rings. This spatial distribution, with the dust rings just outside the gaseous disk, is consistent with the expected interactions between gas and dust in an optically thin disk. This indicates that gas can have a significant effect on the location of the dust within debris disks.


title: Deterministic phase slips in mesoscopic superconducting rings

abstract: The properties of one-dimensional superconductors are strongly influenced by topological fluctuations of the order parameter, known as phase slips, which cause the decay of persistent current in superconducting rings and the appearance of resistance in superconducting wires. Despite extensive work, quantitative studies of phase slips have been limited by uncertainty regarding the order parameter's free energy landscape. Here we show detailed agreement between measurements of the persistent current in isolated flux-biased rings and Ginzburg-Landau theory over a wide range of temperature, magnetic field, and ring size; this agreement provides a quantitative picture of the free energy landscape. We also demonstrate that phase slips occur deterministically as the barrier separating two competing order parameter configurations vanishes. These results will enable studies of quantum and thermal phase slips in a well-characterized system and will provide access to outstanding questions regarding the nature of one-dimensional superconductivity.


title: 3D Gaze Estimation from 2D Pupil Positions on Monocular Head-Mounted Eye Trackers

abstract: 3D gaze information is important for scene-centric attention analysis but accurate estimation and analysis of 3D gaze in real-world environments remains challenging. We present a novel 3D gaze estimation method for monocular head-mounted eye trackers. In contrast to previous work, our method does not aim to infer 3D eyeball poses but directly maps 2D pupil positions to 3D gaze directions in scene camera coordinate space. We first provide a detailed discussion of the 3D gaze estimation task and summarize different methods, including our own. We then evaluate the performance of different 3D gaze estimation approaches using both simulated and real data. Through experimental validation, we demonstrate the effectiveness of our method in reducing parallax error, and we identify research challenges for the design of 3D calibration procedures.


title: Modulating Functions-Based Method for Parameters and Source Estimation in One-Dimensional Partial Differential Equations

abstract: In this paper, modulating functions-based method is proposed for estimating space-time dependent unknowns in one-dimensional partial differential equations. The proposed method simplified the problem into a system of algebraic equations linear in unknown parameters. The well-posedness of modulating functions-based solution is proven. The wave and the fifth order KdV equations are used as examples to show the effectiveness of the proposed method in both noise-free and noisy cases.


title: Compact Groups of Galaxies with Complete Spectroscopic Redshifts in the Local Universe

abstract: Dynamical analysis of compact groups provides important tests of models of compact group formation and evolution. By compiling 2066 redshifts from FLWO/FAST, from the literature, and from SDSS DR12 in the fields of compact groups in \citet{McC09}, we construct the largest sample of compact groups with complete spectroscopic redshifts in the redshift range $0.01 < z < 0.22$. This large redshift sample shows that the interloper fraction in the \citet{McC09} compact group candidates is $\sim 42\%$. A secure sample of 332 compact groups includes 192 groups with four or more member galaxies and 140 groups with three members. The fraction of early-type galaxies in these compact groups is 62\%, slightly higher than for the original Hickson compact groups. The velocity dispersions of early- and late-type galaxies in compact groups change little with groupcentric radius; the radii sampled are less than $100 ~h^{-1}$ kpc, smaller than the radii typically sampled by members of massive clusters of galaxies. The physical properties of our sample compact groups include size, number density, velocity dispersion, and local environment; these properties slightly differ from those derived for the original Hickson compact groups and for the DPOSS II compact groups. Differences result from subtle differences in the way the group candidates were originally selected. The space density of the compact groups changes little with redshift over the range covered by this sample. The approximate constancy of the space density for this sample is a potential constraint on the evolution of compact groups on a few Gigayear timescale.


title: Charged Perfect Fluid Sphere in Higher Dimensional Spacetime

abstract: Present paper provides a new model for perfect fluid sphere filled with charge in higher dimensional spacetime admitting conformal symmetry. We consider a linear equation of state with coefficients fixed by the matching conditions at the boundary of the source corresponding to the exterior Reissner-Nordstr\"{o}m higher dimensional spacetime. Several physical features for different dimensions, starting from four up to eleven, are briefly discussed. It is shown that all the features as obtained from the present model are physically desirable and valid as far as the observed data set for the compact star $SAX~J1808.4-3658~(SS2)$ is concerned.


title: Bianchi Type-III Dark Energy Model in f(R,T) gravity with a variable Deceleration parameter

abstract: In this paper, we have investigated Bianchi Type-III Dark Energy Model in f(R,T) theory of gravity proposed by Harko et al. (Phys. Rev. D84 : 024020, 2011). To find a determinate solution of the field equations we have used (i) the fact that scalar expansion of the space-time is proportional to the shear scalar and (ii) a variable deceleration parameter given by Mishra et al. (Int. J. Theor. Phys. 52: 2546, 2013). Physical and kinematical properties of the model are also discussed. It is observed that our model is in agreement with the recent cosmological observations.


title: On the reduction of generalized polylogarithms to $\text{Li}_n$ and $\text{Li}_{2,2}$ and on the evaluation thereof

abstract: We give expressions for all generalized polylogarithms up to weight four in terms of the functions log, $\text{Li}_n$, and $\text{Li}_{2,2}$, valid for arbitrary complex variables. Furthermore we provide algorithms for manipulation and numerical evaluation of $\text{Li}_n$ and $\text{Li}_{2,2}$, and add codes in Mathematica and C++ implementing the results. With these results we calculate a number of previously unknown integrals, which we add in App. C.


title: Inference rules for RDF(S) and OWL in N3Logic

abstract: This paper presents inference rules for Resource Description Framework (RDF), RDF Schema (RDFS) and Web Ontology Language (OWL). Our formalization is based on Notation 3 Logic, which extended RDF by logical symbols and created Semantic Web logic for deductive RDF graph stores. We also propose OWL-P that is a lightweight formalism of OWL and supports soft inferences by omitting complex language constructs.


title: Relaxation dynamics of the asymmetric simple exclusion process with Langmuir kinetics on a ring

abstract: We consider the asymmetric simple exclusion process with Langmuir kinetics on a periodic lattice. We analytically obtain the exact time evolution of correlation functions with arbitrary length starting from the initial state with no particle in the system. The exact stationary state of this model has been known for the totally asymmetric case. We propose a basis transformation which simplifies the proof of the stationarity of this state and enables the generalization to the partially asymmetric case. Moreover, we construct low-energy excitations and obtain the exact relaxation time.

categories: cond-mat.stat-mech math-ph math.MP nlin.SI


title: The Poisson Boltzmann equation and the charge separation phenomenon at the silica-water interface: A holistic approach

abstract: The Poisson Boltzmann equation is known for its success in describing the Debye layer that arises from the charge separation phenomenon at the silica/water interface. However, by treating only the mobile ionic charges in the liquid, the Poisson Boltzmann equation accounts for only half of the electrical double layer, with the other half, the surface charge layer, being beyond its computational domain. In this work, we take a holistic approach to the charge separation phenomenon at the silica/water interface by treating, within a single computational domain, the electrical double layer that comprises both the mobile ions in the liquid and the surface charge density. The Poisson Nernst Planck equations are used as the rigorous basis for our methodology. The holistic approach has the advantage of being able to predict surface charge variations that arise either from the addition of salt and acid to the liquid, or from the decrease of the liquid channel width to below twice the Debye length. As the electrical double layer must be overall neutral, we use this constraint to derive both the form of the static limit of the Poisson Nernst Planck equations, as well as a global chemical potential that replaces the classical zeta potential as the boundary value for the PB equation, which can be re-derived from our formalism. We present several predictions of our theory that are beyond the framework of the PB equation alone, e.g., the surface capacitance and the so-called pK and pL values, the isoelectronic point at which the surface charge layer is neutralized, and the appearance of a Donnan potential that arises from the formation of an electrical double layer at the inlet regions of a nano-channel connected to the bulk reservoir. All theory predictions are shown to be in good agreement with the experimental observations.


title: The nuclear dust lane of Circinus: collimation without a torus

abstract: In some AGN, nuclear dust lanes connected to kpc-scale dust structures provide all the extinction required to obscure the nucleus, challenging the role of the dusty torus proposed by the Unified Model. In this letter we show the pc-scale dust and ionized gas maps of Circinus constructed using sub-arcsec-accuracy registration of infrared VLT AO images with optical \textit{Hubble Space Telescope} images. We find that the collimation of the ionized gas does not require a torus but is caused by the distribution of dust lanes of the host galaxy on $\sim$10 pc scales. This finding questions the presumed torus morphology and its role at parsec scales, as one of its main attributes is to collimate the nuclear radiation, and is in line with interferometric observations which show that most of the pc-scale dust is in the polar direction. We estimate that the nuclear dust lane in Circinus provides $1/3$ of the extinction required to obscure the nucleus. This constitutes a conservative lower limit to the obscuration at the central parsecs, where the dust filaments might get optically thicker if they are the channels that transport material from $\sim$100 pc scales to the centre.


title: Using young massive star clusters to understand star formation and feedback in high-redshift-like environments

abstract: The formation environment of stars in massive stellar clusters is similar to the environment of stars forming in galaxies at a redshift of 1 - 3, at the peak star formation rate density of the Universe. As massive clusters are still forming at the present day at a fraction of the distance to high-redshift galaxies they offer an opportunity to understand the processes controlling star formation and feedback in conditions similar to those in which most stars in the Universe formed. Here we describe a system of massive clusters and their progenitor gas clouds in the centre of the Milky Way, and outline how detailed observations of this system may be able to: (i) help answer some of the fundamental open questions in star formation and (ii) quantify how stellar feedback couples to the surrounding interstellar medium in this high-pressure, high-redshift analogue environment.


title: Proper image subtraction - optimal transient detection, photometry and hypothesis testing

abstract: Transient detection and flux measurement via image subtraction stand at the base of time domain astronomy. Due to the varying seeing conditions, the image subtraction process is non-trivial, and existing solutions suffer from a variety of problems. Starting from basic statistical principles, we develop the optimal statistic for transient detection, flux measurement and any image-difference hypothesis testing. We derive a closed-form statistic that: (i) Is mathematically proven to be the optimal transient detection statistic in the limit of background-dominated noise; (ii) Is numerically stable; (iii) For accurately registered, adequately sampled images, does not leave subtraction or deconvolution artifacts; (iv) Allows automatic transient detection to the theoretical sensitivity limit by providing credible detection significance; (v) Has uncorrelated white noise; (vi) Is a sufficient statistic for any further statistical test on the difference image, and in particular, allows to distinguish particle hits and other image artifacts from real transients; (vii) Is symmetric to the exchange of the new and reference images; (viii) Is at least an order of magnitude faster to compute than some popular methods; and (ix) Is straightforward to implement. Furthermore, we present extensions of this method that make it resilient to registration errors, color-refraction errors, and any noise source that can be modelled. In addition, we show that the optimal way to prepare a reference image is the proper image coaddition presented in Zackay \& Ofek (2015b). We demonstrate this method on simulated data and real observations from the Palomar Transient Factory data release 2. We provide an implementation of this algorithm in MATLAB and Python.


title: Localizing energy in granular materials

abstract: A device for absorbing and storing short duration impulses in an initially uncompressed one-dimensional granular chain is presented. Simply stated, short regions of sufficiently soft grains are embedded in a hard granular chain. These grains exhibit long-lived standing waves of predictable frequencies regardless of the timing of the arrival of solitary waves from the larger matrix. We explore the origins, symmetry, and energy content of the soft region and its intrinsic modes.


title: Accurate Tight-Binding Hamiltonians for 2D and Layered Materials

abstract: We present a scheme to controllably improve the accuracy of tight-binding Hamiltonian matrices derived by projecting the solutions of plane-wave ab initio calculations on atomic orbital basis sets. By systematically increasing the completeness of the basis set of atomic orbitals, we are able to optimize the quality of the band structure interpolation over wide energy ranges including unoccupied states. This methodology is applied to the case of interlayer and image states, which appear several eV above the Fermi level in materials with large interstitial regions or surfaces such as graphite and graphene. Due to their spatial localization in the empty regions inside or outside of the system, these states have been inaccessible to traditional tight-binding models and even to ab initio calculations with atom-centered basis functions.


title: Information-theoretic thresholds for community detection in sparse networks

abstract: We give upper and lower bounds on the information-theoretic threshold for community detection in the stochastic block model. Specifically, let $k$ be the number of groups, $d$ be the average degree, the probability of edges between vertices within and between groups be $c_\mathrm{in}/n$ and $c_\mathrm{out}/n$ respectively, and let $\lambda = (c_\mathrm{in}-c_\mathrm{out})/(kd)$. We show that, when $k$ is large, and $\lambda = O(1/k)$, the critical value of $d$ at which community detection becomes possible -- in physical terms, the condensation threshold -- is \[ d_c = \Theta\!\left( \frac{\log k}{k \lambda^2} \right) \, , \] with tighter results in certain regimes. Above this threshold, we show that the only partitions of the nodes into $k$ groups are correlated with the ground truth, giving an exponential-time algorithm that performs better than chance -- in particular, detection is possible for $k \ge 5$ in the disassortative case $\lambda < 0$ and for $k \ge 11$ in the assortative case $\lambda > 0$. (Similar upper bounds were obtained independently by Abbe and Sandon.) Below this threshold, we use recent results of Neeman and Netrapalli (who generalized arguments of Mossel, Neeman, and Sly) to show that no algorithm can label the vertices better than chance, or even distinguish the block model from an Erd\H{o}s-R\'enyi random graph with high probability. We also rely on bounds on certain functions of doubly stochastic matrices due to Achlioptas and Naor; indeed, our lower bound on $d_c$ is the second moment lower bound on the $k$-colorability threshold for random graphs with a certain effective degree.

categories: math.PR cond-mat.stat-mech cs.CC cs.IT cs.SI math.IT


title: Two Higgs doublet models with an $S_3$ symmetry

abstract: We study all implementations of the $S_3$ symmetry in the two Higgs doublet model with quarks, consistent with non-zero quark masses and a Cabibbo-Kobayashi-Maskawa (CKM) matrix which is not block diagonal. We study the impact of the various soft-breaking terms and vacuum expectation values, and find an interesting relation between $\alpha$ and $\beta$. We also show that, in this minimal setting, only two types of assignments are possible: either all field sectors are in singlets or all field sectors have a doublet.


title: Three-Point Disc Amplitudes in the RNS Formalism

abstract: We calculate all tree level string theory vacuum to Dp-brane disc amplitudes involving an arbitrary RR-state and two NS-NS vertex operators. This computation was earlier performed by K. Becker, Guo, and Robbins for the simplest case of a RR-state of type C_{p-3}. Here we use the aid of a computer to calculate all possible three-point amplitudes involving a RR-vertex operator of type C_{p+1+2k}.


title: Systematic errors in low latency gravitational wave parameter estimation impact electromagnetic follow-up observations

abstract: Among the most eagerly anticipated opportunities made possible by Advanced LIGO/Virgo are multimessenger observations of compact mergers. Optical counterparts may be short-lived so rapid characterization of gravitational wave (GW) events is paramount for discovering electromagnetic signatures. One way to meet the demand for rapid GW parameter estimation is to trade off accuracy for speed, using waveform models with simplified treatment of the compact objects' spin. We report on the systematic errors in GW parameter estimation suffered when using different spin approximations to recover generic signals. Component mass measurements can be biased by $>5\sigma$ using simple-precession waveforms and in excess of $20\sigma$ when non-spinning templates are employed. This suggests that electromagnetic observing campaigns should not take a strict approach to selecting which LIGO/Virgo candidates warrant follow-up observations based on low-latency mass estimates. For sky localization, we find searched areas are up to a factor of ${\sim}$2 larger for non-spinning analyses, and are systematically larger for any of the simplified waveforms considered in our analysis. Distance biases for the non-precessing waveforms can be in excess of 100\% and are largest when the spin angular momenta are in the orbital plane of the binary. We confirm that spin-aligned waveforms should be used for low-latency parameter estimation at the minimum. Including simple precession, though more computationally costly, mitigates biases except for signals with extreme precession effects. Our results shine a spotlight on the critical need for development of computationally inexpensive precessing waveforms and/or massively parallel algorithms for parameter estimation.


title: On Disjunctive and Conjunctive Set-Labelings of Graphs

abstract: Let $X$ be a non-empty set and $\sP(X)$ be its power set. A set-valuation or a set-labeling of a given graph $G$ is an injective function $f:V(G) \to \sP(X)$ such that the induced function $f^{\ast}:E(G) \to \sP(X)$ defined by $f^{\ast} (uv) = f(u)\ast f(v)$, where $\ast$ is a binary operation on sets. A set-indexer of a graph $G$ is an injective set-valued function $f:V(G) \to \sP(X)$ such that the induced function $f^{\ast}:E(G) \to \sP(X)$ is also injective. In this paper, two types of set-labelings, called conjunctive set-labeling and disjunctive set-labeling, of graphs are introduced and some properties and characteristics of these types of set-labelings of graphs are studied.


title: Damping of Landau levels in neutral graphene at low magnetic fields: A phonon Raman scattering study

abstract: Landau level broadening mechanisms in electrically neutral and quasineutral graphene were investigated through micro-magneto-Raman experiments in three different samples, namely, a natural single-layer graphene flake and a back-gated single-layer device, both deposited over Si/SiO2 substrates, and a multilayer epitaxial graphene employed as a reference sample. Interband Landau level transition widths were estimated through a quantitative analysis of the magnetophonon resonances associated with optically active Landau level transitions crossing the energy of the E_2g Raman-active phonon. Contrary to multilayer graphene, the single-layer graphene samples show a strong damping of the low-field resonances, consistent with an additional broadening contribution of the Landau level energies arising from a random strain field. This extra contribution is properly quantified in terms of a pseudomagnetic field distribution Delta_B = 1.0-1.7 T in our single-layer samples.


title: A Comparative Study of Knots of Star Formation in Interacting vs. Spiral Galaxies

abstract: Interacting galaxies are known to have higher global rates of star formation on average than normal galaxies, relative to their stellar masses. Using UV and IR photometry combined with new and published H-alpha images, we have compared the star formation rates of ~700 star forming complexes in 46 nearby interacting galaxy pairs with those of regions in 39 normal spiral galaxies. The interacting galaxies have proportionally more regions with high star formation rates than the spirals. The most extreme regions in the interacting systems lie at the intersections of spiral/tidal structures, where gas is expected to pile up and trigger star formation. Published Hubble Telescope images show unusually large and luminous star clusters in the highest luminosity regions. The star formation rates of the clumps correlate with measures of the dust attenuation, consistent with the idea that regions with more interstellar gas have more star formation. For the clumps with the highest star formation rates, the apparent dust attenuation is consistent with the Calzetti starburst dust attenuation law. This suggests that the high luminosity regions are dominated by a central group of young stars surrounded by a shell of clumpy interstellar gas. In contrast, the lower luminosity clumps are bright in the UV relative to H-alpha, suggesting either a high differential attenuation between the ionized gas and the stars, or a post-starburst population bright in the UV but faded in H-alpha. The fraction of the global light of the galaxies in the clumps is higher on average for the interacting galaxies than for the spirals. Thus the star forming regions in interacting galaxies are more luminous, dustier, or younger on average.


title: HerMES: A search for high-redshift dusty galaxies in the HerMES Large Mode Survey - Catalogue, number counts and early results

abstract: Selecting sources with rising flux densities towards longer wavelengths from Herschel/SPIRE maps is an efficient way to produce a catalogue rich in high-redshift (z > 4) dusty star-forming galaxies. The effectiveness of this approach has already been confirmed by spectroscopic follow-up observations, but the previously available catalogues made this way are limited by small survey areas. Here we apply a map-based search method to 274 deg$^2$ of the HerMES Large Mode Survey (HeLMS) and create a catalogue of 477 objects with SPIRE flux densities $S_{500} > S_{350} >S_{250}$ and a $5 \sigma$ cut-off $S_{500} > $ 52 mJy. From this catalogue we determine that the total number of these "red" sources is at least an order of magnitude higher than predicted by galaxy evolution models. These results are in agreement with previous findings in smaller HerMES fields; however, due to our significantly larger sample size we are also able to investigate the shape of the red source counts for the first time. We have obtained spectroscopic redshift measurements for two of our sources using the Atacama Large Millimeter/submillimeter Array (ALMA). The redshifts z = 5.1 and z = 3.8 confirm that with our selection method we can indeed find high-redshift dusty star-forming galaxies.


title: Experimentally accessible witnesses of many-body localisation

abstract: The phenomenon of many-body localised (MBL) systems has attracted significant interest in recent years, for its intriguing implications from a perspective of both condensed-matter and statistical physics: they are insulators even at non-zero temperature and fail to thermalise, violating expectations from quantum statistical mechanics. What is more, recent seminal experimental developments with ultra-cold atoms in optical lattices constituting analog quantum simulators have pushed many-body localised systems into the realm of physical systems that can be measured with high accuracy. In this work, we introduce experimentally accessible witnesses that directly probe distinct features of MBL, distinguishing it from its Anderson counterpart. We insist on building our toolbox from techniques available in the laboratory, including on-site addressing, super-lattices, and time-of-flight measurements, identifying witnesses based on fluctuations, density-density correlators, densities, and entanglement. We build upon the theory of out of equilibrium quantum systems, in conjunction with tensor network and exact simulations, showing the effectiveness of the tools for realistic models.


title: Kirchhoff migration without phases

abstract: We present a simple, frequency domain, preprocessing step to Kirchhoff migration that allows the method to image scatterers when the wave field phase information is lost at the receivers, and only intensities are measured. The resulting imaging method does not require knowing the phases of the probing field or manipulating the phase of the wave field at the receivers. In a regime where the scattered field is small compared to the probing field, the problem of recovering the full-waveform scattered field from intensity data can be formulated as an embarrassingly simple least-squares problem. Although this only recovers the projection (on a known subspace) of the full-waveform scattered field, we show that, for high frequencies, this projection gives Kirchhoff images asymptotically identical to the images obtained with full waveform data. Our method can also be used when the source is modulated by a Gaussian process and autocorrelations are measured at an array of receivers.


title: Compositions in YBa2Cu3Oy for 3-D bond order superconductors

abstract: Self-doping in YBa2Cu3Oy starts to switch from a charge balancing of chain metal reduction and plane oxidation to one of chain ligand oxidation (O-) and plane reduction on quenching from >600K. Amongst the responses are plane expansions and types of unusually strong superconductivity such as elevated temperature superconductivity (ETS), observed through laser pulsing (Tc=552K*) and upon shot quenching (Tc=200K*). We ascribe ETS to limited 3-D superconductivity due to a correlated system of bond ordering within chain-plane sandwiches and propose how to stabilize it. Accordingly, plane expanding n-doping arises from self-doping charge equilibration with local chain Cu of two-fold O coordination (2). The dumbbell type bonds of both apical O- are a result of high-energy environment and comparable in their metric with electron pairs on the plane. Following empirical Tc=%(2)x11 we suggest increasing %(2) to cause the observed retrograde rise in Tc (photo-induced reflectivity edge) to a theoretical limit of Tc=1100K* at 100%(2). We propose compositions and heat treatments based on charge-lattice commensurability. Paired charge concentrations on planes, determined by bond ordering based on magic number counts, suggest several promising candidates such as c=0.22=2/3a0x3b0, c=0.17=2/3x4 or c=0.080=2/5x5. The latter could be achieved with 32%(2) and display Tc=352K* not only in laser pulsed but in shot quenched materials.


title: Microfabrication of large area high-stress silicon nitride membranes for optomechanical devices

abstract: In view of the integration of membrane resonators with more complex MEMS structures, we developed a general fabrication procedure for circular shape SiN$_x$ membranes using Deep Reactive Ion Etching (DRIE). Large area and high-stress SiN$_x$ membranes were fabricated and used as optomechanical resonators in a Michelson interferometer and in a Fabry-P\'erot cavity. The measurements show that the fabrication process preserves both the optical quality and the mechanical quality factor of the membrane.


title: On an extension of the Iwatsuka model

abstract: We prove absolute continuity for an extended class of two-dimensional magnetic Hamiltonians that were initially studied by A. Iwatsuka. In particular, we add an electric field that is translation invariant in the same direction as the magnetic field is. As an example, we study the effective Hamiltonian for a thin quantum layer in a homogeneous magnetic field.


title: Internal wave pressure, velocity, and energy flux from density perturbations

abstract: Determination of energy transport is crucial for understanding the energy budget and fluid circulation in density varying fluids such as the ocean and the atmosphere. However, it is rarely possible to determine the energy flux field $\mathbf{J} = p \mathbf{u}$, which requires simultaneous measurements of the pressure and velocity perturbation fields, $p$ and $\mathbf{u}$. We present a method for obtaining the instantaneous $\mathbf{J}(x,z,t)$ from density perturbations alone: a Green's function-based calculation yields $p$, and $\mathbf{u}$ is obtained by integrating the continuity equation and the incompressibility condition. We validate our method with results from Navier-Stokes simulations: the Green's function method is applied to the density perturbation field from the simulations, and the result for $\mathbf{J}$ is found to agree typically to within $1\%$ with $\mathbf{J}$ computed directly using $p$ and $ \mathbf{u}$ from the Navier-Stokes simulation. We also apply the Green's function method to density perturbation data from laboratory schlieren measurements of internal waves in a stratified fluid, and the result for $\mathbf{J}$ agrees to within $6\%$ with results from Navier-Stokes simulations. Our method for determining the instantaneous velocity, pressure, and energy flux fields applies to any system described by a linear approximation of the density perturbation field, e.g., to small amplitude lee waves and propagating vertical modes. The method can be applied using our Matlab graphical user interface EnergyFlux.


title: Extreme residues of Dedekind zeta functions

abstract: In a family of $S_{d+1}$-fields ($d=2,3,4$), we obtain the true upper and lower bound of the residues of Dedekind zeta functions except for a density zero set. For $S_5$-fields, we need to assume the strong Artin conjecture. We also show that there exists an infinite family of number fields with the upper and lower bound, resp.


title: Chemodynamics of the Milky Way and disc formation history: insight from the RAVE and Gaia-ESO surveys

abstract: Multi-object spectrographs have opened a new window on the analyses of the chemo-dynamical properties of old Milky Way stars. These analyses allow us to trace back the internal mechanisms and the external factors that have influenced the evolution of our Galaxy, and therefore understand fundamental aspects of galaxy evolution in general. Here, we present recent results from the RAdial Velocity Experiment (RAVE) and the Gaia-ESO survey. These surveys explore the Milky Way properties in different ways, in terms of sample size and selection, magnitude range, and spectral resolution. We focus here on (i) the first direct detection of evidence for radial migration within the thin disc, providing insight into the history of spiral structure of the Milky Way, and (ii) the chemo-dynamical characterisation of the metal-weak thick and thin discs, for which chemo-dynamical models still have difficulties in reproducing.


title: An analogue of the Aluffi algebra for modules

abstract: P. Aluffi introduced in [1] a new graded algebra in order to conveniently express characteristic cycles in the theory of singular varieties. This algebra is attached to a surjective ring homomorphism $A\surjects B$ by taking a suitable inverse limit of graded algebras, one for each representation of $A$ as a residue ring of a given "ambient" ring $R.$ Since giving a ring surjection $A\surjects B$ is tantamount to giving an ideal $I \subset A$, it would seem natural to ask for an analogous notion for $A$-modules. This is the central purpose of this work. Since a given module may not admit any embedding into a free module, a preparatory toil includes dealing with this technical point at the outset. On the bright side, the intrusion of modules raises a few algebraic questions interesting on their own. It is to expect that this extension to modules may be transcribed in terms of coherent sheaves, thus possibly providing an answer to a question by Aluffi in this regard. Two main bodies of examples are treated in detail to illustrate how the theory works and to show the relation to finer properties of other algebras.


title: Statistical Decoupling of Lagrangian Fluid Parcel in Newtonian Cosmology

abstract: The Lagrangian dynamics of a single fluid element within a self-gravitational matter field is intrinsically non-local due to the presence of the tidal force. This complicates the theoretical investigation of the non-linear evolution of various cosmic objects, e.g. dark matter halos, in the context of Lagrangian fluid dynamics, since a fluid parcel with given initial density and shape may evolve differently depending on their environments. In this paper, we provide a statistical solution that could decouple this environmental dependence. After deriving the probability distribution evolution equation of the matter field, our method produces a set of closed ordinary differential equations whose solution is uniquely determined by the initial condition of the fluid element. Mathematically, it corresponds to the projected characteristic curve of the transport equation of the density-weighted probability density function (PDF). Consequently it is guaranteed that the one-point PDF would be preserved by evolving these local, yet non-linear, curves with the same set of initial data as the real system. Physically, these trajectories describe the mean evolution averaged over all environments by substituting the tidal tensor with its conditional average. For Gaussian distributed dynamical variables, this mean tidal tensor is simply proportional to the velocity shear tensor, and the dynamical system would recover the prediction of Zel'dovich approximation (ZA) with the further assumption of the linearized continuity equation. For Weakly non-Gaussian field, the averaged tidal tensor could be expanded perturbatively as a function of all relevant dynamical variables whose coefficients are determined by the statistics of the field.


title: Effect of four-fermion operators on the mass of the composite particles

abstract: We propose a theoretical framework for evaluating the effect of four-fermion operators on the mass of composite particles in confining strongly-coupled gauge theories. The confining sector is modelled by a non-local Nambu-Jona Lasinio action, whereas the four-fermion operators, arising from a different sector, are local. In order to illustrate the method, we investigate a simple toy model with a global $SU(2)_L\times SU(2)_R\to SU(2)_V$ symmetry breaking, and a four-fermion operator breaking $SU(2)_L\times SU(2)_R$ but preserving $SU(2)_V$. In the particle spectrum we only include the pseudoscalar isospin triplet, that is the pseudo-Nambu-Goldstone bosons associated with chiral symmetry breaking, and the lightest scalar singlet. After checking that the nonlocal model successfully accounts for the experimental results in two-flavour QCD, we investigate the mass spectrum as a function of the four-fermion coupling. For our specific choice of four-fermion operator, we find that the mass of the pseudoscalar triplet grows, whereas the mass of the lightest scalar singlet is approximately unaffected, as the four-fermion coupling grows. We argue that these results can be directly tested on the lattice, and briefly discuss possible applications of this technique to models of dynamical electroweak symmetry breaking.


title: Dependence of technological improvement on artifact interactions

abstract: Empirical research has shown performance improvement of many different technological domains occurs exponentially but with widely varying improvement rates. What causes some technologies to improve faster than others do? Previous quantitative modeling research has identified artifact interactions, where a design change in one component influences others, as an important determinant of improvement rates. The models predict that improvement rate for a domain is proportional to the inverse of the domain interaction parameter. However, no empirical research has previously studied and tested the dependence of improvement rates on artifact interactions. A challenge to testing the dependence is that any method for measuring interactions has to be applicable to a wide variety of technologies. Here we propose a patent-based method that is both technology domain-agnostic and less costly than alternative methods. We use textual content from patent sets in 27 domains to find the influence of interactions on improvement rates. Qualitative analysis identified six specific keywords that signal artifact interactions. Patent sets from each domain were then examined to determine the total count of these 6 keywords in each domain, giving an estimate of artifact interactions in each domain. It is found that improvement rates are positively correlated with the inverse of the total count of keywords with correlation coefficient of +0.56 with a p-value of 0.002. The empirical results agree with model predictions and support the suggestion that domains with higher number of artifacts interactions (higher complexity) will improve at a slower pace.


title: Radiative decays $V\rightarrow P\gamma^{*}$ in the instant form of relativistic quantum mechanics

abstract: Calculations of form factor for the radiative decays $V\rightarrow P\gamma^{*}$ process are performed in the framework of an instant form of relativistic quantum mechanics. The electromagnetic current operator for this decay is constructed. The transition form factor is obtained in the so called relativistic modified impulse approximation (MIA). The current operator satisfies the conditions of Lorentz-covariance and current conservation in MIA. The results of the calculations are compared with the analogous results in the light-front dynamics and in the model of vector meson dominance


title: On Calculating the Current-Voltage Characteristic of Multi-Diode Models for Organic Solar Cells

abstract: We provide an alternative formulation of the exact calculation of the current-voltage characteristic of solar cells which have been modeled with a lumped parameters equivalent circuit with one or two diodes. Such models, for instance, are suitable for describing organic solar cells whose current-voltage characteristic curve has an inflection point, also known as an S-shaped anomaly. Our formulation avoids the risk of numerical overflow in the calculation. It is suitable for implementation in Fortran, C or on micro-controllers.


title: Using SVM to pre-classify government purchases

abstract: The Brazilian government often misclassifies the goods it buys. That makes it hard to audit government expenditures. We cannot know whether the price paid for a ballpoint pen (code #7510) was reasonable if the pen was misclassified as a technical drawing pen (code #6675) or as any other good. This paper shows how we can use machine learning to reduce misclassification. I trained a support vector machine (SVM) classifier that takes a product description as input and returns the most likely category codes as output. I trained the classifier using 20 million goods purchased by the Brazilian government between 1999-04-01 and 2015-04-02. In 83.3% of the cases the correct category code was one of the three most likely category codes identified by the classifier. I used the trained classifier to develop a web app that might help the government reduce misclassification. I open sourced the code on GitHub; anyone can use and modify it.


title: Probability density function of in-plane permeability of fibrous media: Constant Kozeny coefficient

abstract: This study addresses the issue of whether or not a normal distribution appropriately represents permeability variations. To do so, (i) the distribution of local fibre volume fraction for each tow is experimentaly determined by estimation of each pair of local areal density and thickness, (ii) the Kozeny-Carmen equation together with the change of variable technique are used to compute the PDF of permeability, (iii) using the local values of fibre volume fraction, the distribution of local average permeability is computed and subsequently the goodness of fit of the computed PDF is compared with the distribution of the permeability at microscale level. Finally variability of local permeability at the microscale level is determined. The first set of results reveals that (1) the relationship between the local areal density and local thickness in a woven carbon-epoxy composite is modelled by a bivariate normal distribution, (2) while fibre volume fraction follows a normal distribution, permeability follows a gamma distribution, (3) this work also shows that there is significant agreement between the analytical approach and the simulation results. The second set of results shows that the coefficient of variation of permeability is one order of magnitude larger than that of fibre volume fraction. Future work will consider other variables, such as type of fabrics, the degree of fibre preform compaction to determine whether or not the bivariate normal model is applicable for a broad range of fabrics.


title: Parametric Number Covariance in Quantum Chaotic Spectra

abstract: We study spectral parametric correlations in quantum chaotic systems and introduce the number covariance as a measure of such correlations. We derive analytic results for the classical random matrix ensembles using the binary correlation method and obtain compact expressions for the co- variance. We illustrate the universality of this measure by presenting the spectral analysis of the quantum kicked rotors for the time-reversal invariant and time-reversal non-invariant cases. A local version of the parametric number variance introduced earlier is also investigated.


title: Software for enumerative and analytic combinatorics

abstract: We survey some general-purpose symbolic software packages that implement algorithms from enumerative and analytic combinatorics. Software for the following areas is covered: basic combinatorial objects, symbolic combinatorics, P\'olya theory, combinatorial species, and asymptotics. We describe the capabilities that the packages offer as well as some of the algorithms used, and provide links to original documentation. Most of the packages are freely downloadable from the web.


title: Wish-list for Run II top quark measurements (for BSM)

abstract: In this contribution I will highlight the new challenges for top quark physics at LHC Run II, focusing in particular on the interplay between precision studies on the top quark and searches for new physics. A new strategy to search for subtle scenarios of new physics is envisaged. The ability of compute and measure very accurately top quark properties such as its production rate, decay properties, e.g rates and distributions, and variables sensitive to the top quark mass is put at the center of this strategy to probe new physics.


title: Coherent Controllable Transport of a Surface Plasmon Coupled to Plasmonic Waveguide with a Metal Nano Particle-Semiconductor Quantum Dot Hybrid System

abstract: By using the real-space method, switching of a single plasmon interacting with a hybrid nanosystem composed of a semiconductor quantum dot (SQD) and a metallic nanoparticle (MNP) coupled to one-dimensional surface plasmonic waveguide is investigated theoretically. We discussed that the dipole coupling between an exciton and a localized surface plasmon results in the formation of a hybrid exciton and the transmission and reflection of the propagating single plasmon could be controlled by changing the interparticle distance between the SQD and the MNP and the size of the nanoparticles. The controllable transport of the propagating single surface plasmon by such a nanosystem discussed here could find the significant potential in the design of next-generation quantum devices such as plasmonic switch, single photon transistor and nanolaser and quantum information.


title: High-temperature superconductivity from fine-tuning of Fermi-surface singularities in iron oxypnictides

abstract: In the family of the iron-based superconductors, the $RE$FeAsO-type compounds (with $RE$ being a rare-earth metal) exhibit the highest bulk superconducting transition temperatures ($T_{\mathrm{c}}$) up to $55\ \textrm{K}$ and thus hold the key to the elusive pairing mechanism. Recently, it has been demonstrated that the intrinsic electronic structure of SmFe$_{0.92}$Co$_{0.08}$AsO ($T_{\mathrm{c}}=18\ \textrm{K}$) is highly nontrivial and consists of multiple band-edge singularities in close proximity to the Fermi level. However, it remains unclear whether these singularities are generic to the $RE$FeAsO-type materials and if so, whether their exact topology is responsible for the aforementioned record $T_{\mathrm{c}}$. In this work, we use angle-resolved photoemission spectroscopy (ARPES) to investigate the inherent electronic structure of the NdFeAsO$_{0.6}$F$_{0.4}$ compound with a twice higher $T_{\mathrm{c}}=38\ \textrm{K}$. We find a similarly singular Fermi surface and further demonstrate that the dramatic enhancement of superconductivity in this compound correlates closely with the fine-tuning of one of the band-edge singularities to within a fraction of the superconducting energy gap $\Delta$ below the Fermi level. Our results provide compelling evidence that the band-structure singularities near the Fermi level in the iron-based superconductors must be explicitly accounted for in any attempt to understand the mechanism of superconducting pairing in these materials.


title: Compressing molecular dynamics trajectories: breaking the one-bit-per-sample barrier

abstract: Molecular dynamics simulations yield large amounts of trajectory data. For their durable storage and accessibility an efficient compression algorithm is paramount. State of the art domain-specific algorithms combine quantization, Huffman encoding and occasionally domain knowledge. We propose the high resolution trajectory compression scheme (HRTC) that relies on piecewise linear functions to approximate quantized trajectories. By splitting the error budget between quantization and approximation, our approach beats the current state of the art by several orders of magnitude given the same error tolerance. It allows storing samples at far less than one bit per sample. It is simple and fast enough to be integrated into the inner simulation loop, store every time step, and become the primary representation of trajectory data.

categories: cs.DC cs.IT math.IT physics.comp-ph q-bio.QM


title: On the irrationality of generalized $q$-logarithm

abstract: For integer $p$, $|p|>1$, and generic rational $x$ and $z$, we establish the irrationality of the series $$\ell_p(x,z)=x\sum_{n=1}^\infty\frac{z^n}{p^n-x}.$$ It is a symmetric ($\ell_p(x,z)=\ell_p(z,x)$) generalization of the $q$-logarithmic function ($x=1$ and $p=1/q$ where $|q|<1$), which in turn generalizes the $q$-harmonic series ($x=z=1$). Our proof makes use of the Hankel determinants built on the Pad\'e approximations to $\ell_p(x,z)$.


title: Observation of Strong Radiation Pressure Forces from Squeezed Light on a Mechanical Oscillator

abstract: Quantum enhanced sensing is a powerful technique in which nonclassical states are used to improve the sensitivity of a measurement. For enhanced mechanical displacement sensing, squeezed states of light have been shown to reduce the photon counting noise that limits the measurement noise floor. It has long been predicted, however, that suppressing the noise floor with squeezed light should produce an unavoidable increase in radiation pressure noise that drives the mechanical system. Such nonclassical radiation pressure forces have thus far been hidden by insufficient measurement strengths and residual thermal mechanical motion. Since the ultimate measurement sensitivity relies on the delicate balance between these two noise sources, the limits of the quantum enhancement have not been observed. Using a microwave cavity optomechanical system, we observe the nonclassical radiation pressure noise that necessarily accompanies any quantum enhancement of the measurement precision. By varying both the magnitude and phase of the squeezing, we optimize the fundamental trade-off between mechanical imprecision and backaction noise in accordance with the Heisenberg uncertainty principle. As the strength of the measurement is further increased, radiation pressure forces eventually dominate the mechanical motion. In this regime, the optomechanical interaction can be exploited as an efficient quantum nondemolition (QND) measurement of the amplitude fluctuations of the light field. By overwhelming mechanical thermal noise with radiation pressure by two orders of magnitude, we demonstrate a mechanically-mediated measurement of the squeezing with an effective homodyne efficiency of 94%. Thus, with strong radiation pressures forces, mechanical motion enhances the measurement of nonclassical light, just as nonclassical light enhances the measurement of the motion.


title: Measurement scheme and analysis for weak ground state hyperfine transition moments through two-pathway coherent control

abstract: We report our detailed analysis of a table-top system for the measurement of the weak-force-induced electric dipole moment of a ground state hyperfine transition carried out in an atomic beam geometry. We describe an experimental configuration of conductors for application of orthogonal r.f. and static electric fields, with cavity enhancement of the r.f. field amplitude, that allows confinement of the r.f. field to a region in which the static fields are uniform and well-characterized. We carry out detailed numerical simulations of the field modes, and analyze the expected magnitude of statistical and systematic limits to the measurement of this transition amplitude in atomic cesium. The combination of an atomic beam with this configuration leads to strong suppression of magnetic dipole contributions to the atomic signal. The application of this technique to the measurement of extremely weak transition amplitudes in other atomic systems, especially alkali metals, seems very feasible.


title: An observation on Classification of Lie point symmetries for quadratic Li\'{e}nard type equation $\ddot{x}+f\left(x\right) \dot{x}^{2}+g\left( x\right) =0$ [JMP 54, 053506 (2013)] and its erratum [JMP 55, 059901 (2014)]

abstract: We demonstrate a simplification of some recent works on the classification of the Lie symmetries for a quadratic equation of Li\'{e}nard type. We observe that the problem could have been resolved more simply.


title: Embryo as an active granular fluid: stress-coordinated cellular constriction chains

abstract: Mechanical stress plays an intricate role in gene expression in individual cells and sculpting of developing tissues. However, systematic methods of studying how mechanical stress and feedback help to harmonize cellular activities within a tissue have yet to be developed. Motivated by our observation of the cellular constriction chains (CCCs) during the initial phase of ventral furrow formation in the Drosophila melanogaster embryo, we propose an active granular fluid (AGF) model that provides valuable insights into cellular coordination in the apical constriction process. In our model, cells are treated as circular particles connected by a predefined force network, and they undergo a random constriction process in which the particle constriction probability P is a function of the stress exerted on the particle by its neighbors. We find that when P favors tensile stress, constricted particles tend to form chain-like structures. In contrast, constricted particles tend to form compact clusters when P favors compression. A remarkable similarity of constricted-particle chains and CCCs observed in vivo provides indirect evidence that tensile-stress feedback coordinates the apical constriction activity. We expect that our particle-based AGF model will be useful in analyzing mechanical feedback effects in a wide variety of morphogenesis and organogenesis phenomena.


title: The effects of assembly bias on cosmological inference from galaxy-galaxy lensing and galaxy clusters

abstract: The combination of galaxy-galaxy lensing (GGL) and galaxy clustering is a promising route to measuring the amplitude of matter clustering and testing modified gravity theories of cosmic acceleration. Halo occupation distribution (HOD) modeling can extend the approach down to nonlinear scales, but galaxy assembly bias could introduce systematic errors by causing the HOD to vary with large scale environment at fixed halo mass. We investigate this problem using the mock galaxy catalogs created by Hearin & Watson (2013, HW13), which exhibit significant assembly bias because galaxy luminosity is tied to halo peak circular velocity and galaxy colour is tied to halo formation time. The preferential placement of galaxies (especially red galaxies) in older halos affects the cutoff of the mean occupation function $\langle N_\text{cen}(M_\text{min}) \rangle$ for central galaxies, with halos in overdense regions more likely to host galaxies. The effect of assembly bias on the satellite galaxy HOD is minimal. We introduce an extended, environment dependent HOD (EDHOD) prescription to describe these results and fit galaxy correlation measurements. Crucially, we find that the galaxy-matter cross-correlation coefficient, $r_{gm} \equiv \xi_{gm} \cdot [ \xi_{mm} \xi_{gg} ]^{-1/2}$, is insensitive to assembly bias on scales $r \gtrsim 1 \; h^{-1}\text{Mpc}$, even though $\xi_{gm} $ and $\xi_{gg} $ are both affected individually. We can therefore recover the correct $\xi_{mm} $ from the HW13 galaxy-galaxy and galaxy-matter correlations using either a standard HOD or EDHOD fitting method. For $M_r \leq -19$ or $M_r \leq -20$ samples the recovery of $\xi_{mm}$ is accurate to 2% or better. For a sample of red $M_r \leq -20$ galaxies we achieve 2% recovery at $r \gtrsim 2\;h^{-1}\text{Mpc}$ with EDHOD modeling but lower accuracy at smaller scales or with a standard HOD fit.


title: Breakdown of the Kondo insulating state in SmB6 by introducing Sm vacancies

abstract: SmB6 is a proposed topological Kondo insulator where the presence of topological nontriviality can be tuned by variations in the Sm valence. Experimentally, Sm valence can be changed by tuning stoichiometry of SmB6. We show that Raman scattering can detect vacancies lower than 1% of Sm sites in SmB6 crystal by probing the intensity of defect-induced scattering of the acoustic phonon branch at 10~meV. In the electronic Raman spectra of SmB6 at temperatures below 130~K, we observe features developing in A$_{1g}$ and E$_g$ symmetries at 100 and 41~meV which we assign to excitations between hybridized bands, and depressed spectral weight below 20~meV associated with the hybridization gap. With the increased number of Sm vacancies up to 1% we observe an increase of spectral weight below 20~meV showing that the gap is filling in with electronic states. For the sample with the lowest number of vacancies the in-gap exciton excitations with long lifetimes protected by hybridization gap are observed at 16-18~meV in E$_g$ and T$_{2g}$ symmetries. These excitations broaden as a decrease in the lifetime with increasing number of vacancies and are quenched by the presence of in-gap states at concentration of Sm vacancies of about 1%. Based on this study we suggest that only the most stoichiometric SmB6 samples have a bulk gap necessary for topological Kondo insulators.


title: On Milstein approximations with varying coefficients: the case of super-linear diffusion coefficients

abstract: A new class of explicit Milstein schemes, which approximate stochastic differential equations (SDEs) with superlinearly growing drift and diffusion coefficients, is proposed in this article. It is shown, under very mild conditions, that these explicit schemes converge in $\mathcal L^p$ to the solution of the corresponding SDEs with optimal rate.


title: High-fidelity spatial and polarization addressing of Ca-43 qubits using near-field microwave control

abstract: Individual addressing of qubits is essential for scalable quantum computation. Spatial addressing allows unlimited numbers of qubits to share the same frequency, whilst enabling arbitrary parallel operations. We demonstrate addressing of long-lived $^{43}\text{Ca}^+$ "atomic clock" qubits held in separate zones ($960\mu$m apart) of a microfabricated surface trap with integrated microwave electrodes. Such zones could form part of a "quantum CCD" architecture for a large-scale quantum information processor. By coherently cancelling the microwave field in one zone we measure a ratio of Rabi frequencies between addressed and non-addressed qubits of up to 1400, from which we calculate a spin-flip probability on the qubit transition of the non-addressed ion of $1.3\times 10^{-6}$. Off-resonant excitation then becomes the dominant error process, at around $5 \times 10^{-3}$. It can be prevented either by working at higher magnetic field, or by polarization control of the microwave field. We implement polarization control with error $2 \times 10^{-5}$, which would suffice to suppress off-resonant excitation to the $\sim 10^{-9}$ level if combined with spatial addressing. Such polarization control could also enable fast microwave operations.


title: Minimum Average Delay of Routing Trees

abstract: The general communication tree embedding problem is the problem of mapping a set of communicating terminals, represented by a graph G, into the set of vertices of some physical network represented by a tree T. In the case where the vertices of G are mapped into the leaves of the host tree T the underlying tree is called a routing tree and if the internal vertices of T are forced to have degree 3, the host tree is known as layout tree. Different optimization problems have been studied in the class of communication tree problems such as well-known minimum edge dilation and minimum edge congestion problems. In this report we study the less investigate measure i.e. tree length, which is a representative for average edge dilation (communication delay) measure and also for average edge congestion measure. We show that finding a routing tree T for an arbitrary graph G with minimum tree length is an NP-Hard problem.


title: Efficient Markov Chain Monte Carlo Sampling for Hierarchical Hidden Markov Models

abstract: Traditional Markov chain Monte Carlo (MCMC) sampling of hidden Markov models (HMMs) involves latent states underlying an imperfect observation process, and generates posterior samples for top-level parameters concurrently with nuisance latent variables. When potentially many HMMs are embedded within a hierarchical model, this can result in prohibitively long MCMC runtimes. We study combinations of existing methods, which are shown to vastly improve computational efficiency for these hierarchical models while maintaining the modeling flexibility provided by embedded HMMs. The methods include discrete filtering of the HMM likelihood to remove latent states, reduced data representations, and a novel procedure for dynamic block sampling of posterior dimensions. The first two methods have been used in isolation in existing application-specific software, but are not generally available for incorporation in arbitrary model structures. Using the NIMBLE package for R, we develop and test combined computational approaches using three examples from ecological capture-recapture, although our methods are generally applicable to any embedded discrete HMMs. These combinations provide several orders of magnitude improvement in MCMC sampling efficiency, defined as the rate of generating effectively independent posterior samples. In addition to being computationally significant for this class of hierarchical models, this result underscores the potential for vast improvements to MCMC sampling efficiency which can result from combinations of known algorithms.


title: Group Communication Over LTE : A Radio Access Perspective

abstract: Long Term Evolution (LTE), which has its root on commercial mobile communications, recently becomes an influential solution to future public safety communications. To verify the feasibility of LTE for public safety, it is essential to investigate whether an LTE system optimized for one-to-one communications is capable of providing group communication, which is one of the most important service concepts in public safety. In general, a number of first responders for public safety need to form a group for communicating with each other or sharing the common data for collaboration on their mission. In this article, we analyze how the current LTE system can support group communication in a radio access aspect. Based on the requirements for group communication, we validate whether each LTE-enabled radio access method can efficiently support group communication. In addition, we propose a new multicast transmission scheme, termed index-coded Hybrid Automatic Retransmission reQuest (HARQ). By applying the index coding concept to HARQ operations, we show that the LTE system can provide group communication more sophisticatedly in terms of radio resource efficiency and scalability. We finally evaluate the performance of LTE-enabled group communication using several radio access methods and show how the proposed transmission scheme brings the performance enhancement via system level simulations.


title: Vibrational Signatures in the THz Spectrum of 1,3-DNB: A First-Principles and Experimental Study

abstract: Understanding the fundamental processes of light-matter interaction is important for detection of explosives and other energetic materials, which are active in the infrared and terahertz (THz) region. We report a comprehensive study on electronic and vibrational lattice properties of structurally similar 1,3-dinitrobenzene (1,3- DNB) crystals through first-principles electronic structure calculations and THz spectroscopy measurements on polycrystalline samples. Starting from reported x-ray crystal structures, we use density-functional theory (DFT) with periodic boundary conditions to optimize the structures and perform linear response calculations of the vibrational properties at zero phonon momentum. The theoretically identified normal modes agree qualitatively with those obtained experimentally in a frequency range up to 2.5 THz and quantitatively at much higher frequencies. The latter frequencies are set by intra-molecular forces. Our results suggest that van der Waals dispersion forces need to be included to improve the agreement between theory and experiment in the THz region, which is dominated by intermolecular modes and sensitive to details in the DFT calculation. An improved comparison is needed to assess and distinguish between intra- and intermolecular vibrational modes characteristic of energetic materials.


title: Imaging polarimetry and spectropolarimetry of comet C/2013 R1 (Lovejoy)

abstract: We have obtained imaging polarimetry of the comet C/2013 R1 (Lovejoy) with 2-Channel-Focal-Reducer Rozhen instrument at 2m Ritchey-Chr\'etien-Coud\'e telescope of the Bulgarian National Astronomical Observatory Rozhen in two dust continuum filters covering wavelength intervals clear from molecular emissions and centred at $4430\AA$ in (blue filter) and at $6840\AA$ in (red filter). In imaging mode we measured the degree of linear polarisation $17.01\pm0.09$% in the blue and $18.81\pm0.02$% in the red, which is in a very good agreement with measurements of other comets at the similar phase angle. We have also obtained polarisation maps in both filters. We found a strong correlation between the spatial distribution of the polarisation and the dust colour. Spectropolarimetry of the nucleus region shows an increase of the polarisation with wavelength, and a depolarisation in the spectral regions with gas emission lines, most noticeable in C$_2$ emission band, which shows a polarisation of $6.0\pm1.1$%.


title: Suzaku Confirms NGC~3660 is an Unabsorbed Seyfert 2

abstract: An enigmatic group of objects, unabsorbed Seyfert 2s may have intrinsically weak broad line regions, obscuration in the line of sight to the BLR but not to the X-ray corona, or so much obscuration that the X-ray continuum is completely suppressed and the observed spectrum is actually scattered into the line of sight from nearby material. NGC 3660 has been shown to have weak broad optical/near infrared lines, no obscuration in the soft X-ray band, and no indication of "changing look" behavior. The only previous hard X-ray detection of this source by Beppo-SAX seemed to indicate that the source might harbor a heavily obscured nucleus. However, our analysis of a long-look Suzaku observation of this source shows that this is not the case, and that this source has a typical power law X-ray continuum with normal reflection and no obscuration. We conclude that NGC 3660 is confirmed to have no unidentified obscuration and that the anomolously high Beppo-SAX measurement must be due to source confusion or similar, being inconsistent with our Suzaku measurements as well as non-detections from Swift-BAT and RXTE.


title: Steepest-entropy-ascent quantum thermodynamic modeling of the far-from-equilibrium interactions between nonequilibrium systems of indistinguishable particle ensembles

abstract: This paper presents a nonequilibrium, first-principles, thermodynamic-ensemble based model for the relaxation process of interacting non-equilibrium systems. This model is formulated using steepest-entropy-ascent quantum thermodynamics (SEAQT) and its equation of motion for a grand canonical ensemble and is applied to a many particle system of classical or indistinguishable particles. Two kinds of interactions are discussed, including pure heat diffusion and heat and mass diffusion together. Since no local equilibrium assumption is made, the conjugate fluxes and forces are intrinsic to the subspaces of the state space of one system and/or of the state space of the two interacting systems. They are derived via the concepts of hypoequilibrium state and nonequilibrium intensive properties, which describe the nonmutual equilibrium status between subspaces of the thermodynamic state space of a single system and/or of the state space of the two interacting systems. The Onsager relations are shown to be thermodynamic kinematic features of the system and are found without knowledge of the detailed mechanics of the dynamic process. A fundamental thermodynamic explanation for the measurement of each intensive property of a system in a nonequilibrium state is given. The fundamental thermodynamic definition of reservoir is also discussed. Finally, the equation of motion for a system undergoing multiple interactions is provided, which permits the modeling of a network of local systems in nonequilibrium at any spatial and temporal scale.


title: Global Classical Solutions to the Relativistic Boltzmann Equation without Angular Cut-off

abstract: We prove the unique existence and exponential decay of global in time classical solutions to the special relativistic Boltzmann equation without any angular cut-off assumptions with initial perturbations in some weighted Sobolev spaces. We consider perturbations of the relativistic Maxwellian equilibrium states. We work in the case of a spatially periodic box. We consider the general conditions on the collision kernel from Dudy\'nski and Ekiel-Je\'zewska (Commun Math Phys \textbf{115}(4):607--629, 1985). Additionally, we prove sharp constructive upper and coercive lower bounds for the linearized relativistic Boltzmann collision operator in terms of a geometric fractional Sobolev norm; this shows that a spectral gap exists and that this behavior is similar to that of the non-relativistic case as shown by Gressman and Strain (Journal of AMS \textbf{24}(3), 771--847, 2011). Lastly, we derive the relativistic analogue of Carleman dual representation of Boltzmann collision operator. This is the first global existence and stability result for relativistic Boltzmann equation without angular cutoff and this resolves the open question of perturbative global existence for the relativistic kinetic theory without the Grad's angular cut-off assumption.


title: Robobarista: Learning to Manipulate Novel Objects via Deep Multimodal Embedding

abstract: There is a large variety of objects and appliances in human environments, such as stoves, coffee dispensers, juice extractors, and so on. It is challenging for a roboticist to program a robot for each of these object types and for each of their instantiations. In this work, we present a novel approach to manipulation planning based on the idea that many household objects share similarly-operated object parts. We formulate the manipulation planning as a structured prediction problem and learn to transfer manipulation strategy across different objects by embedding point-cloud, natural language, and manipulation trajectory data into a shared embedding space using a deep neural network. In order to learn semantically meaningful spaces throughout our network, we introduce a method for pre-training its lower layers for multimodal feature embedding and a method for fine-tuning this embedding space using a loss-based margin. In order to collect a large number of manipulation demonstrations for different objects, we develop a new crowd-sourcing platform called Robobarista. We test our model on our dataset consisting of 116 objects and appliances with 249 parts along with 250 language instructions, for which there are 1225 crowd-sourced manipulation demonstrations. We further show that our robot with our model can even prepare a cup of a latte with appliances it has never seen before.


title: Two Small Temperate Planets Transiting Nearby M Dwarfs in K2 Campaigns 0 and 1

abstract: The prime Kepler mission revealed that small planets (<4 R_earth) are common, especially around low-mass M dwarfs. K2, the re-purposed Kepler mission, continues this exploration of small planets around small stars. Here we combine K2 photometry with spectroscopy, adaptive optics imaging, and archival survey images to analyze two small planets orbiting the nearby, field age, M dwarfs K2-26 (EPIC 202083828) and K2-9. K2-26 is an M1.0 +/- 0.5 dwarf at 93 +/- 7 pc from K2 Campaign 0. We validate its 14.5665 d period planet and estimate a radius of 2.67^+0.46_-0.42 R_earth. K2-9 is an M2.5 +/- 0.5 dwarf at 110 +/- 12 pc from K2 Campaign 1. K2-9b was first identified by Montet et al. 2015; here we present spectra and adaptive optics imaging of the host star and independently validate and characterize the planet. Our analyses indicate K2-9b is a 2.25^+0.53_-0.96 R_earth planet with a 18.4498 d period. K2-26b exhibits a transit duration that is too long to be consistent with a circular orbit given the measured stellar radius. Thus, the long transits are likely due to the photoeccentric effect and our transit fits hint at an eccentric orbit. Both planets receive low incident flux from their host stars and have estimated equilibrium temperatures <500 K. K2-9b may receive approximately Earth-like insolation. However, its host star exhibits strong GALEX UV emission which could affect any atmosphere it harbors. K2-26b and K2-9b are representatives of a poorly studied class of small planets with cool temperatures that have radii intermediate to Earth and Neptune. Future study of these systems can provide key insight into trends in bulk composition and atmospheric properties at the transition from silicate dominated to volatile rich bodies.


title: Noncommutative Solenoids and the Gromov-Hausdorff Propinquity

abstract: We prove that noncommutative solenoids are limits, in the sense of the Gromov-Hausdorff propinquity, of quantum tori. From this observation, we prove that noncommutative solenoids can be approximated by finite dimensional quantum compact metric spaces, and that they form a continuous family of quantum compact metric spaces over the space of multipliers of the solenoid, properly metrized.


title: A hybrid multi-time-step framework for pore-scale and continuum-scale modeling of solute transport in porous media

abstract: In this paper, we propose a computational framework,which is based on a domain decomposition technique, to employ both finite element method (which is a popular continuum modeling approach) and lattice Boltzmann method (which is a popular pore-scale modeling approach) in the same computational domain. To bridge the gap across the disparate length and time-scales, we first propose a new method to enforce continuum-scale boundary conditions (i.e., Dirichlet and Neumann boundary conditions) onto the numerical solution from the lattice Boltzmann method. This method are based on maximization of entropy and preserve the non-negativity of discrete distributions under the lattice Boltzmann method. The proposed computational framework allows different grid sizes, orders of interpolation, and time-steps in different subdomains. This allows for different desired resolutions in the numerical solution in different subdomains. Through numerical experiments, the effect of grid and time-step refinement, disparity of time-steps in different subdomains, domain partitioning, and the number of iteration steps on the accuracy and rate of convergence of the proposed methodology are studied. Finally, to showcase the performance of this framework in porous media applications, we use it to simulate the dissolution of calcium carbonate in a porous structure.