title: The SIMPlest Miracle

abstract: It has recently been proposed that dark matter could be a thermal relic of 3-to-2 scatterings in a strongly coupled hidden sector. We present explicit classes of strongly coupled gauge theories that admit this behavior. These are QCD-like theories of dynamical chiral symmetry breaking, where the pions play the role of dark matter. The number-changing 3-to-2 process, which sets the dark matter relic abundance, arises from the Wess-Zumino-Witten term. The theories give an explicit relationship between the 3-to-2 annihilation rate and the 2-to-2 self-scattering rate, which alters predictions for structure formation. This is a simple calculable realization of the strongly-interacting-massive-particle (SIMP) mechanism.


title: Flat space (higher spin) gravity with chemical potentials

abstract: We introduce flat space spin-3 gravity in the presence of chemical potentials and discuss some applications to flat space cosmology solutions, their entropy, free energy and flat space orbifold singularity resolution. Our results include flat space Einstein gravity with chemical potentials as special case. We discover novel types of phase transitions between flat space cosmologies with spin-3 hair and show that the branch that continuously connects to spin-2 gravity becomes thermodynamically unstable for sufficiently large temperature or spin-3 chemical potential.


title: The stability of stellar disks in Milky-Way sized dark matter halos

abstract: We employ an improved methodology to insert live stellar disks into high-resolution dark matter simulations of Milky Way sized halos, allowing us to investigate the fate of thin stellar disks in the tumultuous environment of cold dark matter structures. We study a set of eight different halos, drawn from the Aquarius simulation project, in which stellar disks are adiabatically grown with a prescribed structure, and then allowed to self-consistently evolve. The initial velocity distribution is set-up in very good equilibrium with the help of the GALIC code. We find that the residual triaxiality of the halos leads to significant disk tumbling, qualitatively confirming earlier work. We show that the disk turning motion is unaffected by structural properties of the galaxies such as the presence or absence of a bulge or bar. In typical Milky Way sized dark matter halos, we expect an average turning of the disks by about 40 degrees between z=1 and z=0, over the coarse of 6 Gyr. We also investigate the impact of the disks on substructures, and conversely, the disk heating rate caused by the dark matter halo substructures. The presence of disks reduces the central subhalo abundance by a about a factor of two, due to an increased evaporation rate by gravitational shocks from disk passages. We find that substructures are important for heating the outer parts of stellar disks but do not appear to significantly affect their inner parts.


title: Self-Interacting Dark Matter through the Higgs Portal

abstract: We study self-interacting dark matter coupled to the Standard Model via the Higgs portal. We consider a scenario where dark matter is a thermal relic with strong enough self interactions that can alleviate the problems of collisionless cold dark matter. We study constraints from direct detection searches, the LHC, and Big Bang nucleosynthesis. We show that the tension between these constraints and the need for sufficiently strong self-interactions with light mediators can be alleviated by coupling the mediator to either active or sterile neutrinos. Future direct detection data offers great potential and can be used to find evidence of a light mediator and verify that dark matter scatters via long-range self-interactions.


title: On the two-loop corrections to the Higgs mass in trilinear R-parity violation

abstract: We study the impact of large trilinear R-parity violating couplings on the lightest CP-even Higgs boson mass in supersymmetric models. We use the publicly available computer codes SARAH and SPheno to compute the leading two-loop corrections. We use the effective potential approach. For not too heavy third generation squarks (< 1 TeV) and couplings close to the unitarity bound we find positive corrections up to a few GeV in the Higgs mass.


title: Emission line tomography of the short period cataclysmic variables CC Scl and V2051 Oph

abstract: We present time-series spectroscopy of two short period cataclysmic variables, CC Scl and V2051 Oph, to test the efficiency of Doppler tomography-based methods in constraining orbital parameters of evolved cataclysmic variables. We find that the Ca~II triplet lines offer superior diagnostics, revealing emission components from the mass donors and sharp images of the accretion discs. Furthermore, we use Monte-Carlo methods to estimate the uncertainties from ensembles of Doppler maps. We compare our new methods against traditional radial velocity methods and show that they offer a valid route towards system parameter determination. Our analysis of CC Scl suggests a low mass ratio of $q=0.08\pm0.03$ with a primary velocity of $K_1=37\pm14$ km/s. This mass ratio is in between the pre- and post-period minimum status, however our $K_1$ solution favours a post-period minimum system. Our derived parameters for V2051 Oph ($q= 0.16\pm 0.03$, $K_1=97\pm10$ km/s) are in agreement with the eclipse solution ($q=0.19\pm0.03$), offering a direct validation of our methods.


title: Low-energy phonons and superconductivity in Sn$_{0.8}$In$_{0.2}$Te

abstract: We present neutron scattering measurements on low-energy phonons from a superconducting ($T_{c} = 2.7$~K) Sn$_{0.8}$In$_{0.2}$Te single crystal sample. The longitudinal acoustic phonon mode and one transverse acoustic branch have been mapped out around the (002) Bragg peak for temperatures of 1.7 K and 4.2 K. We observe a substantial energy width of the transverse phonons at energies comparable to twice the superconducting gap; however, there is no change in this width between the superconducting and normal states. We also confirm that the compound is well ordered, with no indications of structural instability.


title: Isomorphism invariants of enveloping algebras

abstract: Let $L$ be a Lie algebra with its enveloping algebra $U(L)$ over a field. In this paper we survey results concerning the isomorphism problem for enveloping algebras: given another Lie algebra $H$ for which $U(L)$ and $U(H)$ are isomorphic as associative algebras, can we deduce that $L$ and $H$ are isomorphic Lie algebras? Over a field of positive characteristic we consider a similar problem for restricted Lie algebras, that is, given restricted Lie algebras $L$ and $H$ for which their restricted enveloping algebras are isomorphic as algebras, can we deduce that $L$ and $H$ are isomorphic?


title: Importance of repulsive interactions for the equation of state and other properties of strongly interacting matter

abstract: We illustrate the role of repulsive interactions in a hadron-resonance gas at freeze-out and in a gas of quark-gluon bags. Taking into account non-zero size of particles in hadron gas leads to a significant decrease and shift of the net-baryon density maximum. The transition point from baryon to meson dominated matter depends on the difference between baryon and meson radii. We also show that depending on the properties of the quark-gluon bags one may obtain any type of the phase transition from hadron gas to quark-gluon plasma: the first or second order, as well as four types of the crossover.


title: Energy-Efficient Routing in Wireless Networks in the Presence of Jamming

abstract: The effectiveness and simple implementation of physical layer jammers make them an essential threat for wireless networks. In a multihop wireless network, where jammers can interfere with the transmission of user messages at intermediate nodes along the path, one can employ jamming oblivious routing and then employ physical-layer techniques (e.g. spread spectrum) to suppress jamming. However, whereas these approaches can provide significant gains, the residual jamming can still severely limit system performance. This motivates the consideration of routing approaches that account for the differences in the jamming environment between different paths. First, we take a straightforward approach where an equal outage probability is allocated to each link along a path and develop a minimum energy routing solution. Next, we demonstrate the shortcomings of this approach and then consider the joint problem of outage allocation and routing by employing an approximation to the link outage probability. This yields an efficient and effective routing algorithm that only requires knowledge of the measured jamming at each node. Numerical results demonstrate that the amount of energy saved by the proposed methods with respect to a standard minimum energy routing algorithm, especially for parameters appropriate for terrestrial wireless networks, is substantial.


title: Holistic Collaborative Privacy Framework for Users' Privacy in Social Recommender Service

abstract: The current business model for existing recommender services is centered around the availability of users' personal data at their side whereas consumers have to trust that the recommender service providers will not use their data in a malicious way. With the increasing number of cases for privacy breaches, different countries and corporations have issued privacy laws and regulations to define the best practices for the protection of personal information. The data protection directive 95/46/EC and the privacy principles established by the Organization for Economic Cooperation and Development (OECD) are examples of such regulation frameworks. In this paper, we assert that utilizing third-party recommender services to generate accurate referrals are feasible, while preserving the privacy of the users' sensitive information which will be residing on a clear form only on his/her own device. As a result, each user who benefits from the third-party recommender service will have absolute control over what to release from his/her own preferences. We proposed a collaborative privacy middleware that executes a two stage concealment process within a distributed data collection protocol in order to attain this claim. Additionally, the proposed solution complies with one of the common privacy regulation frameworks for fair information practice in a natural and functional way -which is OECD privacy principles. The approach presented in this paper is easily integrated into the current business model as it is implemented using a middleware that runs at the end-users side and utilizes the social nature of content distribution services to implement a topological data collection protocol.


title: Field emission in diode and triode vacuum nanostructures

abstract: The paper discusses the nano-diode and nano-triode structures of vacuum electronics. Such structures may have the dielectric film with the thickness of several nanometers wich is located on the cathode. Such film with a large dielectric constant reduces the thickness of the potential barrier by about film thickness and reduces the height of barrier. The grid electrode structure may contain several periods of metallization with a thickness of tens of nanometers, which allows one to obtain the resonant under barrier tunneling. For all structures we obtained electrostatic Green's function, the built profiles potential barriers, the calculated tunneling coefficients and Volt-Ampere Characteristics (VAC) with regard to the distribution of electron energies. The structures under consideration require to use of low voltages on the gate (grid) electrodes. They are promising for the electron gun with a large current controlled by low voltage on the grid, for example, in TWT with subsequent acceleration of the electron beam. Since a major influence on the tunneling provides the grid, the electrons emitted from the grid with approximately the Fermi velocity, and the calculation of the electron gun with multiple electrodes after the grid is greatly simplified.


title: Big Bang Synthesis of Nuclear Dark Matter

abstract: We investigate the physics of dark matter models featuring composite bound states carrying a large conserved dark "nucleon" number. The properties of sufficiently large dark nuclei may obey simple scaling laws, and we find that this scaling can determine the number distribution of nuclei resulting from Big Bang Dark Nucleosynthesis. For plausible models of asymmetric dark matter, dark nuclei of large nucleon number, e.g. > 10^8, may be synthesised, with the number distribution taking one of two characteristic forms. If small-nucleon-number fusions are sufficiently fast, the distribution of dark nuclei takes on a logarithmically-peaked, universal form, independent of many details of the initial conditions and small-number interactions. In the case of a substantial bottleneck to nucleosynthesis for small dark nuclei, we find the surprising result that even larger nuclei, with size >> 10^8, are often finally synthesised, again with a simple number distribution. We briefly discuss the constraints arising from the novel dark sector energetics, and the extended set of (often parametrically light) dark sector states that can occur in complete models of nuclear dark matter. The physics of the coherent enhancement of direct detection signals, the nature of the accompanying dark-sector form factors, and the possible modifications to astrophysical processes are discussed in detail in a companion paper.


title: Multi-brid DBI Inflation

abstract: It is shown that it is possible to apply $\delta \cN$ to some special non-canonic cases. We extended the multi-brid idea to the multi-field separable model with a non-canonical kinetic term, mainly the DBI model. Assuming a specific surface for the end of inflation and introducing new fields, enable us to find an explicit expression for the number of e-folds in terms of the new fields. By using $\delta$N formalism, we arrived at the cosmological parameters. We considered the DBI model for two different limits, viz., speed limit and constant sound speed.


title: Star formation in Chamaeleon I and III: a molecular line study of the starless core population

abstract: The Chamaeleon clouds are excellent targets for low-mass star formation studies. Cha I and II are actively forming stars while Cha III shows no sign of ongoing star formation. We aim to determine the driving factors that have led to the very different levels of star formation activity in Cha I and III and examine the dynamical state and possible evolution of the starless cores within them. Observations were performed in various molecular transitions with APEX and Mopra. Five cores are gravitationally bound in Cha I and one in Cha III. The infall signature is seen toward 8-17 cores in Cha I and 2-5 cores in Cha III, which leads to a range of 13-28% of the cores in Cha I and 10-25% of the cores in Cha III that are contracting and may become prestellar. Future dynamical interactions between the cores will not be dynamically significant in either Cha I or III, but the subregion Cha I North may experience collisions between cores within ~0.7 Myr. Turbulence dissipation in the cores of both clouds is seen in the high-density tracers N2H+ 1-0 and HC3N 10-9. Evidence of depletion in the Cha I core interiors is seen in the abundance distributions of C17O, C18O, and C34S. Both contraction and static chemical models indicate that the HC3N to N2H+ abundance ratio is a good evolutionary indicator in the prestellar phase for both gravitationally bound and unbound cores. In the framework of these models, we find that the cores in Cha III and the southern part of Cha I are in a similar evolutionary stage and are less chemically evolved than the central region of Cha I. The measured HC3N/N2H+ abundance ratio and the evidence for contraction motions seen towards the Cha III starless cores suggest that Cha III is younger than Cha I Centre and that some of its cores may form stars in the future. The cores in Cha I South may on the other hand be transient structures. (abridged)


title: TCP Congestion Control Identification

abstract: Transmission Control Protocol (TCP) carries most of the traffic on the Internet these days. There are several implementations of TCP, and the most important difference among them is their mechanism for controlling congestion. One of the methods for determining type of a TCP is active probing. Active probing considers a TCP implementation as a black box, sends different streams of data to the appropriate host. According to the response received from the host, it figures out the type of TCP version implemented. TCP Behavior Inference Tool (TBIT) is an implemented tool that uses active probing to check the running TCP on web servers. It can check several aspects of the running TCP including initial value of congestion window, congestion control algorithm, conformant congestion control, response to selective acknowledgment, response to Explicit Congestion Notification (ECN) and time wait duration. In this paper we focus on congestion control algorithm aspect of it, explain the mechanism used by TBIT and present the results.


title: Proton Charge Radius and Precision Tests of QED

abstract: The "proton radius puzzle" remains unsolved since it was established in 2010. This paper summarizes the current state and gives an overview over upcoming experiments.


title: Spatially Extended NaI D Resonant Emission and Absorption in the Galactic Wind of the Nearby Infrared-Luminous Quasar F05189-2524

abstract: Emission from metal resonant lines has recently emerged as a potentially powerful probe of the structure of galactic winds at low and high redshift. In this work, we present only the second example of spatially resolved observations of NaI D emission from a galactic wind in a nearby galaxy (and the first 3D observations at any redshift). F05189-2524, a nearby (z=0.043) ultra luminous infrared galaxy powered by a quasar, was observed with the integral field unit on the Gemini Multi-Object Spectrograph (GMOS) at Gemini North. NaI D absorption in the system traces dusty filaments on the near side of an extended, AGN-driven galactic wind (with projected velocities up to 2000 km/s). These filaments (A_V < 4) and N(H) < 10^22 cm^-2) simultaneously obscure the stellar continuum and NaI D emission lines. The NaI D emission lines serve as a complementary probe of the wind; they are strongest in regions of low foreground obscuration and extend up to the limits of the field of view (galactocentric radii of 4 kpc). An azimuthally symmetric Sersic model extincted by the same foreground screen as the stellar continuum reproduces the NaI D emission line surface brightness distribution except in the inner regions of the wind, where some emission-line filling of absorption lines may occur. The presence of detectable NaI D emission in F05189-2524 may be due to its high continuum surface brightness at the rest wavelength of NaI D. These data uniquely constrain current models of cool gas in galactic winds and serve as a benchmark for future observations and models.


title: Invisibility's Flicker: Detecting Thermal Cloaks via Transient Effects

abstract: Recent research on the development of a thermal cloak has concentrated on engineering an inhomogeneous thermal conductivity and homogeneous volumetric heat capacity. While the perfect cloak of inhomogeneous $\kappa$ and $\rho c_p$ is known to be exact (no signals scattering or penetrating to the cloak's interior), no such analysis has been considered for this case. Using analytic, computational, and experimental techniques, we demonstrate that these approximate cloaks are detectable. Although they work as perfect cloaks in the steady-state, their transient (time-dependent) response is imperfect and a detectable amount of heat is scattered. This is sufficient to determine the presence of a cloak and any heat source it contains, but the material composition hidden within the cloak is not detectable in practice.


title: A compact, large-aperture tunable lens with adaptive spherical correction

abstract: In this paper, we present the proof of concept of a very fast adaptive glass membrane lens with a large aperture/diameter ratio, spherical aberration correction and integrated actuation. The membrane is directly deformed using two piezo actuators that can tune the focal length and the conical parameter. This operating principle allows for a usable aperture of the whole membrane diameter. Together with the efficient actuation mechanism, the aperture is around 2/3 of the total system diameter - at a thickness of less than 2mm. The response time is a few milliseconds at 12mm aperture, which is fast compared to similar systems.


title: Cryptographically Blinded Games: Leveraging Players' Limitations for Equilibria and Profit

abstract: In this work we apply methods from cryptography to enable any number of mutually distrusting players to implement broad classes of mediated equilibria of strategic games without the need for trusted mediation. Our implementation makes use of a (standard) pre-play "cheap talk" phase, in which players engage in free and non-binding communication prior to playing in the original game. In our cheap talk phase, the players execute a secure multi-party computation protocol to sample an action profile from an equilibrium of a "cryptographically blinded" version of the original game, in which actions are encrypted. The essence of our approach is to exploit the power of encryption to selectively restrict the information available to players about sampled action profiles, such that these desirable equilibria can be stably achieved. In contrast to previous applications of cryptography to game theory, this work is the first to employ the paradigm of using encryption to allow players to benefit from hiding information \emph{from themselves}, rather than from others; and we stress that rational players would \emph{choose} to hide the information from themselves.


title: A generalised-Lagrangian-mean model of the interactions between near-inertial waves and mean flow

abstract: Wind forcing of the ocean generates a spectrum of inertia-gravity waves that is sharply peaked near the local inertial (or Coriolis) frequency. The corresponding near-inertial waves (NIWs) are highly energetic and play a significant role in the slow, large-scale dynamics of the ocean. To analyse this role, we develop a new model of the nondissipative interactions between NIWs and balanced motion. The model is derived using the generalised-Lagrangian-mean (GLM) framework (specifically, the glm variant of Soward & Roberts (2010)), taking advantage of the time-scale separation between the two types of motion to average over the short NIW period. We combine Salmon's (2013) variational formulation of GLM with Whitham averaging to obtain a system of equations governing the joint evolution of NIWs and mean flow. Assuming that the mean flow is geostrophically balanced reduces this system to a simple model coupling Young & Ben Jelloul's (1997) equation for NIWs with a modified quasi-geostrophic equation. In this coupled model, the mean flow affects the NIWs through advection and refraction; conversely, the NIWs affect the mean flow by modifying the potential-vorticity inversion - the relation between advected potential vorticity and advecting mean velocity - through a quadratic wave term, consistent with the GLM results of Buhler & McIntyre (1998). The coupled model is Hamiltonian and its conservation laws, for wave action and energy in particular, prove illuminating: on their basis, we identify a new interaction mechanism whereby NIWs forced at large scales extract energy from the balanced flow as their horizontal scale is reduced by differential advection and refraction so that their potential energy increases. A rough estimate suggests that this mechanism could provide a significant sink of energy for mesoscale motion and play a part in the global energetics of the ocean.


title: Anomaly Detection in Dynamic Networks of Varying Size

abstract: Dynamic networks, also called network streams, are an important data representation that applies to many real-world domains. Many sets of network data such as e-mail networks, social networks, or internet traffic networks are best represented by a dynamic network due to the temporal component of the data. One important application in the domain of dynamic network analysis is anomaly detection. Here the task is to identify points in time where the network exhibits behavior radically different from a typical time, either due to some event (like the failure of machines in a computer network) or a shift in the network properties. This problem is made more difficult by the fluid nature of what is considered "normal" network behavior. The volume of traffic on a network, for example, can change over the course of a month or even vary based on the time of the day without being considered unusual. Anomaly detection tests using traditional network statistics have difficulty in these scenarios due to their Density Dependence: as the volume of edges changes the value of the statistics changes as well making it difficult to determine if the change in signal is due to the traffic volume or due to some fundamental shift in the behavior of the network. To more accurately detect anomalies in dynamic networks, we introduce the concept of Density-Consistent network statistics. On synthetically generated graphs anomaly detectors using these statistics show a a 20-400% improvement in the recall when distinguishing graphs drawn from different distributions. When applied to several real datasets Density-Consistent statistics recover multiple network events which standard statistics failed to find.


title: Probing GPDs in Ultraperipheral Collisions

abstract: Ultraperipheral collisions in hadron colliders give new opportunites to investigate the hadron stucture through exclusive photoproduction processes. We describe the possibility of measuring the Generalized Parton Distributions in the Timelike Compton Scattering process and in the production of heavy vector meson.


title: Digital Signal Processing using Stream High Performance Computing: A 512-input Broadband Correlator for Radio Astronomy

abstract: A "large-N" correlator that makes use of Field Programmable Gate Arrays and Graphics Processing Units has been deployed as the digital signal processing system for the Long Wavelength Array station at Owens Valley Radio Observatory (LWA-OV), to enable the Large Aperture Experiment to Detect the Dark Ages (LEDA). The system samples a ~100MHz baseband and processes signals from 512 antennas (256 dual polarization) over a ~58MHz instantaneous sub-band, achieving 16.8Tops/s and 0.236 Tbit/s throughput in a 9kW envelope and single rack footprint. The output data rate is 260MB/s for 9 second time averaging of cross-power and 1 second averaging of total-power data. At deployment, the LWA-OV correlator was the largest in production in terms of N and is the third largest in terms of complex multiply accumulations, after the Very Large Array and Atacama Large Millimeter Array. The correlator's comparatively fast development time and low cost establish a practical foundation for the scalability of a modular, heterogeneous, computing architecture.


title: Stability and bifurcation for the Kuramoto model

abstract: We study the mean-field limit of the Kuramoto model of globally coupled oscillators. By studying the evolution in Fourier space and understanding the domain of dependence, we show a global stability result. Moreover, we can identify function norms to show damping of the order parameter for velocity distributions and perturbations in $\mathcal{W}^{n,1}$ for $n > 1$. Finally, for sufficiently regular velocity distributions we can identify exponential decay in the stable case and otherwise identify finitely many eigenmodes. For these eigenmodes we can show a center-unstable manifold reduction, which gives a rigorous tool to obtain the bifurcation behaviour. The damping is similar to Landau damping for the Vlasov equation.


title: On a Phi-Kirchhoff multivalued problem with critical growth in an Orlicz-Sobolev space

abstract: This paper is concerned with the multiplicity of nontrivial solutions in an Orlicz-Sobolev space for a nonlocal problem with critical growth, involving N-functions and theory of locally Lispchitz continuous functionals.


title: Second laws under control restrictions

abstract: The second law of thermodynamics, formulated as an ultimate bound on the maximum extractable work, has been rigorously derived in multiple scenarios. However, the unavoidable limitations that emerge due to the lack of control on small systems are often disregarded when deriving such bounds, which is specifically important in the context of quantum thermodynamics. Here, we study the maximum extractable work with limited control over the working system and its interaction with the heat bath. We derive a general second law when the set of accessible Hamiltonians of the working-system is arbitrarily restricted. We then apply our bound to particular scenarios that are important in realistic implementations: limitations on the maximum energy gap and local control over many-body systems. We hence demonstrate in what precise way the lack of control affects the second law. In particular, contrary to the unrestricted case, we show that the optimal work extraction is not achieved by simple thermal contacts. Our results do not only generalize the second law to scenarios of practical relevance, but also take first steps in the direction of local thermodynamics.


title: Two-neutrino double-beta decay of $^{150}$Nd to excited final states in $^{150}$Sm

abstract: Double-beta decay is a rare nuclear process in which two neutrons in the nucleus are converted to two protons with the emission of two electrons and two electron anti-neutrinos. We measured the half life of the two-neutrino double-beta decay of $^{150}$Nd to excited final states of $^{150}$Sm by detecting the de-excitation gamma rays of the daughter nucleus. This study yields the first detection of the coincidence gamma rays from the 0$^+_1$ excited state of $^{150}$Sm. These gamma rays have energies of 333.97 keV and 406.52 keV, and are emitted in coincidence through a 0$^+_1\rightarrow$2$^+_1\rightarrow$0$^+_{gs}$ transition. The enriched Nd$_2$O$_3$ sample consisted of 40.13 g $^{150}$Nd and was observed for 642.8 days at the Kimballton Underground Research Facility, producing 21.6 net events in the region of interest. This count rate gives a half life of $T_{1/2}=(1.07^{+0.45}_{-0.25}(stat)\pm0.07(syst.))\times 10^{20}$ years. The effective nuclear matrix element was found to be 0.0465$^{+0.0098}_{-0.0054}$. Finally, lower limits were obtained for decays to higher excited final states. Our half-life measurement agrees within uncertainties with another recent measurement in which no coincidence was employed. Our nuclear matrix element calculation may have an impact on a recent neutrinoless double-beta decay nuclear matrix element calculation which implies the decay to the first excited state in $^{150}$Sm is favored over that to the ground state.


title: Goedel's other legacy and the imperative of a self-reflective science

abstract: The Goedelian approach is discussed as a prime example of a science towards the origins. While mere selfreferential objectification locks in to its own byproducts, self-releasing objectification informs the formation of objects at hand and their different levels of interconnection. Guided by the spirit of Goedel's work a self-reflective science can open the road where old tenets see only blocked paths. "This is, as it were, an analysis of the analysis itself, but if that is done it forms the fundamental of human science, as far as this kind of things is concerned." G. Leibniz, ('Methodus Nova ...', 1673)


title: When do wireless network signals appear Poisson?

abstract: We consider the point process of signal strengths from transmitters in a wireless network observed from a fixed position under models with general signal path loss and random propagation effects. We show via coupling arguments that under general conditions this point process of signal strengths can be well-approximated by an inhomogeneous Poisson or a Cox point processes on the positive real line. We also provide some bounds on the total variation distance between the laws of these point processes and both Poisson and Cox point processes. Under appropriate conditions, these results support the use of a spatial Poisson point process for the underlying positioning of transmitters in models of wireless networks, even if in reality the positioning does not appear Poisson. We apply the results to a number of models with popular choices for positioning of transmitters, path loss functions, and distributions of propagation effects.


title: Interactions of a Stabilized Radion and Duality

abstract: We determine the couplings of the graviscalar radion in Randall-Sundrum models to Standard Model fields propagating in the bulk of the space, taking into account effects arising from the dynamics of the Goldberger-Wise scalar that stabilizes the size of the extra dimension. The leading corrections to the radion couplings are shown to arise from direct contact interactions between the Goldberger-Wise scalar and the Standard Model fields. We obtain a detailed interpretation of the results in terms of the holographic dual of the radion, the dilaton. In doing so, we determine how the familiar identification of the parameters on the two sides of the AdS/CFT correspondence is modified in the presence of couplings of the bulk Standard Model fields to the Goldberger-Wise scalar. We find that corrections to the form of the dilaton couplings from effects associated with the stabilization of the extra dimension are suppressed by the square of the ratio of the dilaton mass to the Kaluza-Klein scale, in good agreement with results from the CFT side of the correspondence.


title: Finiteness of unramified deformation rings

abstract: We prove that the universal unramified deformation ring $R^{\mathrm{unr}}$ of a continuous Galois representation $\overline{\rho}: G_{F^{+}} \rightarrow \mathrm{GL}_n(k)$ (for a totally real field $F^{+}$ and finite field $k$) is finite over $\mathcal{O} = W(k)$ in many cases. We also prove (under similar hypotheses) that the universal deformation ring $R^{\mathrm{univ}}$ is finite over the local deformation ring $R^{\mathrm{loc}}$.


title: Simulating tidal streams in a high resolution dark matter halo

abstract: We simulate tidal streams in the presence and absence of substructures inside the zero redshift snapshot of the Via Lactea II (VL-2) simulation. A halo finder is used to remove and isolate the subhalos found inside the high resolution dark matter halo of VL-2, and the potentials for both the main halo and all the subhalos are constructed individually using the self-consistent field (SCF) method. This allows us to make direct comparison of tidal streams between a smooth halo and a lumpy halo without assuming idealized profiles or triaxial fits. We simulate the kinematics of a star cluster starting with the same orbital position but two different velocities. Although these two orbits are only moderately eccentric and have similar apo- and pericentric distances, we find that the two streams have very different morphologies. We conclude that our model of the potential of VL-2 can provide insights about tidal streams that have not been explored by previous studies using idealized or axisymmetric models.


title: A Hybrid Approach to Finding Relevant Social Media Content for Complex Domain Specific Information Needs

abstract: While contemporary semantic search systems offer to improve classical keyword-based search, they are not always adequate for complex domain specific information needs. The domain of prescription drug abuse, for example, requires knowledge of both ontological concepts and 'intelligible constructs' not typically modeled in ontologies. These intelligible constructs convey essential information that include notions of intensity, frequency, interval, dosage and sentiments, which could be important to the holistic needs of the information seeker. We present a hybrid approach to domain specific information retrieval (or knowledge-aware search) that integrates ontology-driven query interpretation with synonym-based query expansion and domain specific rules, to facilitate search in social media. Our framework is based on a context-free grammar (CFG) that defines the query language of constructs interpretable by the search system. The grammar provides two levels of semantic interpretation: 1) a top-level CFG that facilitates retrieval of diverse textual patterns, which belong to broad templates and 2) a low-level CFG that enables interpretation of certain specific expressions that belong to such patterns. These low-level expressions occur as concepts from four different categories of data: 1) ontological concepts, 2) concepts in lexicons (such as emotions and sentiments), 3) concepts in lexicons with only partial ontology representation, called lexico-ontology concepts (such as side effects and routes of administration (ROA)), and 4) domain specific expressions (such as date, time, interval, frequency and dosage) derived solely through rules. Our approach is embodied in a novel Semantic Web platform called PREDOSE developed for prescription drug abuse epidemiology. Keywords: Knowledge-Aware Search, Ontology, Semantic Search, Background Knowledge, Context-Free Grammar


title: Lagrangian Formulation for Particle Interpretations of Quantum Mechanics: Single-Particle Case

abstract: A Lagrangian description is presented which can be used in conjunction with particle interpretations of quantum mechanics. A special example of such an interpretation is the well-known Bohm model. The Lagrangian density introduced here also contains a potential for guiding the particle. The advantages of this description are that the field equations and the particle equations of motion can both be deduced from a single Lagrangian density expression and that conservation of energy and momentum are assured. After being developed in a general form, this Lagrangian formulation is then applied to the special case of the Bohm model as an example. It is thereby demonstrated that such a Lagrangian description is compatible with the predictions of quantum mechanics.


title: The double radiative $\bar{B}\to X_s\gamma \gamma$ decay at $O(\alpha_{s})$ in QCD

abstract: In these proceedings, we briefly review the individual interference contribution of the electromagnetic dipole operator $\mathcal{O}_{7}$ to the double differential decay width $d\Gamma_{77}/(ds_1\, ds_2)$ for the process $\bar{B} \to X_s \gamma \gamma$ at $O(\alpha_s)$ in QCD, which is based on our previous work. We define two kinematical variables $s_1$ and $s_2$ as $s_i=(p_b - q_i)^2/m_b^2$, where $p_b$, $q_1$, $q_2$ are the momenta of b-quark and two photons. While the (renormalized) virtual corrections are worked out exactly for a certain range of $s_1$ and $s_2$, we retained in the gluon bremsstrahlung process only the leading power w.r.t. the (normalized) hadronic mass $s_3=(p_b-q_1-q_2)^2/m_b^2$ in the underlying triple differential decay width $d\Gamma_{77}/(ds_1 ds_2 ds_3)$. We found that the double differential decay width, based on this approximation, is free of infrared- and collinear singularities when summing up the virtual- and real-radiation corrections, while this was not the case when keeping all powers in $s_3$ in the gluon bremsstrahlung process due to the configurations allowing collinear photon emission from the (massless) $s$-quark. Lastly, we compare our analytical results with those obtained in a recently extended work, where a non-zero strange quark mass was introduced to regulate the collinear photon configurations.


title: Mechanical generation of networks with surplus complexity

abstract: In previous work I examined an information based complexity measure of networks with weighted links. The measure was compared with that obtained from by randomly shuffling the original network, forming an Erdos-Renyi random network preserving the original link weight distribution. It was found that real world networks almost invariably had higher complexity than their shuffled counterparts, whereas networks mechanically generated via preferential attachment did not. The same experiment was performed on foodwebs generated by an artificial life system, Tierra, and a couple of evolutionary ecology systems, EcoLab and WebWorld. These latter systems often exhibited the same complexity excess shown by real world networks, suggesting that the complexity surplus indicates the presence of evolutionary dynamics. In this paper, I report on a mechanical network generation system that does produce this complexity surplus. The heart of the idea is to construct the network of state transitions of a chaotic dynamical system, such as the Lorenz equation. This indicates that complexity surplus is a more fundamental trait than that of being an evolutionary system.


title: Temporal and spectral properties of quantum light

abstract: The modes of the electromagnetic field are solutions of Maxwell's equations taking into account the material boundary conditions. The field modes of classical optics - properly normalized - are also the mode functions of quantum optics. Quantum physics adds that the excitation within each mode is quantized in close analogy to the harmonic oscillator. A complete set of mode functions forms a basis with which any new modes can be reconstructed. In full generality each electromagnetic mode function in the four dimensional space-time is mathematically equivalent to a harmonic oscillator. The quantization of the electromagnetic field defines the excitation per mode and the correlation between modes. In classical optics there can be oscillations and stochastic fluctuations of amplitude, phase, polarization et cetera. In quantum optics there are in addition uncertain quantum field components, quantum correlations and quantized energies. Here, we present selected topics from classical to quantum optics. We start in the second chapter with the classical optics description of a light field and its spectral densities, their measurement and their interpretation. In the third chapter the quantum properties of a single light mode are reviewed as well as ways to measure these quantum properties. Gaussian states of a light mode are emphasized, i. e. states for which the Wigner function has a two dimensional Gaussian shape. The fourth chapter will be concerned with more than one mode presenting a unifying approach to quadratic Hamiltonians including phase conjugation which is related to time reversal.


title: Short time-scale frequency and amplitude variations in the pulsations of an roAp star: HD 217522

abstract: Photometric observations of HD~217522 in 1981 revealed only one pulsation frequency $\nu_1=1.21529$~mHz. Subsequent observations in 1989 showed the presence of an additional frequency $\nu_2=2.0174$~mHz. New observations in 2008 confirm the presence of the mode with $\nu_2=2.0174$~mHz. Examination of the 1989 data shows amplitude modulation over a time scale of the order of a day, much shorter than what has been observed in other roAp stars. High spectral and time resolution data obtained using the VLT in 2008 confirm the presence of $\nu_2$ and short term modulations in the radial velocity amplitudes of rare earth elements. This suggests growth and decay times shorter than a day, more typical of solar-like oscillations. The driving mechanism of roAp stars and the Sun are different, and the growth and decay seen in the Sun are due to stochastic nature of the driving mechanism. The driving mechanism in roAp stars usually leads to mode stability on a longer time-scale than in the Sun. We interpret the reported change in $\nu_1$ between the 1982 and 1989 data as part of the general frequency variability observed in this star on many time scales.


title: Broadband resonances in ITO nanorod arrays

abstract: In the nanophotonics community, there is an active discussion regarding the origin of the selective absorption/scattering of light by the resonances with nanorod arrays. Here we report a study of the resonances in ordered indium-tin-oxide (ITO) nanorod arrays resulted from the waveguide modes. We discover that with only 2.4% geometrical coverage, the micron-length nanorod arrays strongly interact with light across an extra-wide band from visible to mid-infrared resulting in less than 10% transmission. Simulations show excellent agreement with our experimental observation. Near-field profile obtained from simulation reveals the electric field is mainly localized on the surfaces of the nanorods at all the resonances. Theoretical analysis is then applied to explain the resonances and it was found that the resonances in the visible are different from those in the infrared. When the light arrives at the array, part of the light wave propagates through the free space in between the nanorods and part of the wave is guided inside the nanorods and the phase difference at the ends of the rod interactions forms the basis of the resonances in the visible region; while the resonances in the infrared are Fabry-Perot resonances of the surface guided waves between the two ends of the nanorods. The simple analytical formulae developed predict the spectral positions of these resonances well. This information can be used to design devices like wavelength-selective photodetector, modulators, and nanorod-based solar cells.


title: Local representations of the loop braid group

abstract: We study representations of the loop braid group $LB_n$ from the perspective of extending representations of the braid group $B_n$. We also pursue a generalization of the braid/Hecke/Temperlely-Lieb paradigm---uniform finite dimensional quotient algebras of the loop braid group algebras.


title: Solving Large-Scale Robust Stability Problems by Exploiting the Parallel Structure of Polya's Theorem

abstract: In this paper, we propose a distributed computing approach to solving large-scale robust stability problems on the simplex. Our approach is to formulate the robust stability problem as an optimization problem with polynomial variables and polynomial inequality constraints. We use Polya's theorem to convert the polynomial optimization problem to a set of highly structured Linear Matrix Inequalities (LMIs). We then use a slight modification of a common interior-point primal-dual algorithm to solve the structured LMI constraints. This yields a set of extremely large yet structured computations. We then map the structure of the computations to a decentralized computing environment consisting of independent processing nodes with a structured adjacency matrix. The result is an algorithm which can solve the robust stability problem with the same per-core complexity as the deterministic stability problem with a conservatism which is only a function of the number of processors available. Numerical tests on cluster computers and supercomputers demonstrate the ability of the algorithm to efficiently utilize hundreds and potentially thousands of processors and analyze systems with 100+ dimensional state-space. The proposed algorithms can be extended to perform stability analysis of nonlinear systems and robust controller synthesis.


title: Highly tunable ultra-narrow-resonances with optical nano-antenna phased arrays in the infrared

abstract: We report our recent development in pursuing high Quality-Factor (high-Q factor) plasmonic resonances, with vertically aligned two dimensional (2-D) periodic nanorod arrays. The 2-D vertically aligned nano-antenna array can have high-Q resonances varying arbitrarily from near infrared to terahertz regime, as the antenna resonances of the nanorod are highly tunable through material properties, the length of the nanorod, and the orthogonal polarization direction with respect to the lattice surface,. The high-Q in combination with the small optical mode volume gives a very high Purcell factor, which could potentially be applied to various enhanced nonlinear photonics or optoelectronic devices. The 'hot spots' around the nanorods can be easily harvested as no index-matching is necessary. The resonances maintain their high-Q factor with the change of the environmental refractive index, which is of great interest for molecular sensing.


title: Pouring Cloud Virtualization Security Inside Out

abstract: In this article, virtualization security concerns in the cloud computing domain are reviewed. The focus is toward virtual machine (VM) security where attacks and vulnerabilities such as VM escape, VM hopping, cross-VM side-channel, VM-based rootkits (VMBRs), VM mobility, and VM remote are mentioned and discussed according to their relevance in the clouds. For each attack we outline how they affect the security of cloud systems. Countermeasures and security measures to detect or prevent them through techniques such as VM detection, GuardHype, VM introspection, VM image scanning, etc. are also discussed. Through the surveyed work we present a classification of VM threats within the clouds. Finally, we include our observations and those of other researchers on this matter of cloud virtualization security.


title: General bounds for quantum discord and discord distance

abstract: For any bipartite state, how strongly can one subsystem be quantum correlated with another? Using the Koashi-Winter relation, we study the upper bound of purified quantum discord, which is given by the sum of the von Neumann entropy of the unmeasured subsystem and the entanglement of formation shared between the unmeasured subsystem with the environment. In particular, we find that the Luo et al.'s conjecture on the quantum correlations and the Lindblad conjecture are all ture, when the entanglement of formation vanishes. Let the difference between the left discord and the right discord be captured by the discord distance. If the Lindblad conjecture is true, we show that the joint entropy is a tight upper bound for the discord distance. Further, we obtain a necessary and sufficient condition for saturating upper bounds of purified quantum discord and discord distance separately with the equality conditions for the Araki-Lieb inequality and the Lindblad conjecture. Furthermore, we show that the subadditive relation holds for any bipartite quantum discord.


title: Inclusive $B$ decays and exclusive radiative decays by Belle

abstract: The $b \rightarrow s \gamma$, $b \rightarrow d \gamma$ and $b \rightarrow s \ell^+ \ell^-$ processes are allowed at higher order via the electroweak loop or box diagrams in the Standard model. It is sensitive probe to search for new physics beyond the Standard model because new particles might enter in the loop. We present preliminary results of branching fraction of the $\bar{B} \rightarrow X_s \gamma$, $CP$ asymmetry in the $\bar{B} \rightarrow X_{s+d} \gamma$, and the forward-backward asymmetry in the $\bar{B} \rightarrow X_s \ell^+ \ell^-$. The results are based on a data sample containing $772 \times 10^6 B\bar{B}$ pairs recorded at the $\Upsilon(4S)$ resonance with the Belle detector at the KEKB $e^+ e^-$ collider.


title: Extraction of the Schottky parameters in metal-semiconductor-metal diodes from a single current-voltage measurement

abstract: In order to develop a method to extract the parameters of the two inherent Schottky contacts from a single current-voltage (I-V) characteristic curve, the I-V characteristics of metal-semiconductor-metal (MSM) diodes with asymmetric Schottky barrier heights are theoretically investigated using the thermionic emission model. The MSM diode structure is commonly used because an additional MS interface is required for the electrical characterization of MS diodes. A finite charge-injection barrier is generally formed at the additional interface. When a local maximum was detected in the first-order derivative of the measured I-V characteristics for a MSM diode, the parameters for the Schottky contacts, the zero-bias barrier heights of both MS interfaces, the series resistance of the MSM diode and the effective ideality factor for the MS diode with a higher barrier could be extracted using this method.


title: Absorption-Ablation-Excitation Mechanisms of Laser-Cluster Interactions in a Nanoaerosol System

abstract: The absorption-ablation-excitation mechanism in laser-cluster interactions is investigated by measuring Rayleigh scattering of aerosol clusters along with atomic emission from phase-selective laser-induced breakdown spectroscopy (PS-LIBS). As the excitation laser intensity is increased beyond 0.16GW/cm2, the scattering cross-section of TiO_2 clusters begins to decrease, concurrent with the onset of atomic emission of Ti, indicating a scattering-to-ablation transition and the formation of nanoplasmas. To better clarify the process, time-resolved measurements of scattering signals are examined for different excitation laser intensities. For increasing laser intensities, the cross-sections of clusters decrease during a single pulse, evincing the shorter ablation delay time and larger ratios of ablation clusters. Assessment of the electron energy distribution during the ablation process is conducted by non-dimensionalizing the Fokker-Planck equation, with analogous Strouhal Sl_E, Peclet Pe_E, and Damkohler Da_E numbers defined to characterize the laser-induced aerothermochemical environment. For conditions of Sl_E>>1, Pe_E>>1, and Da_E<<1, the electrons are excited to the conduction band by two-photon absorption, then relax to bottom of the conduction band by collisional electron energy loss to the lattice, and finally serve as the energy transfer media between laser field and lattice. The relation between delay time and excitation intensity is well predicted by this simplified model with quasi-steady assumption.


title: Measuring Influence in Twitter Ecosystems using a Counting Process Modeling Framework

abstract: Data extracted from social media platforms, such as Twitter, are both large in scale and complex in nature, since they contain both unstructured text, as well as structured data, such as time stamps and interactions between users. A key question for such platforms is to determine influential users, in the sense that they generate interactions between members of the platform. Common measures used both in the academic literature and by companies that provide analytics services are variants of the popular web-search PageRank algorithm applied to networks that capture connections between users. In this work, we develop a modeling framework using multivariate interacting counting processes to capture the detailed actions that users undertake on such platforms, namely posting original content, reposting and/or mentioning other users' postings. Based on the proposed model, we also derive a novel influence measure. We discuss estimation of the model parameters through maximum likelihood and establish their asymptotic properties. The proposed model and the accompanying influence measure are illustrated on a data set covering a five year period of the Twitter actions of the members of the US Senate, as well as mainstream news organizations and media personalities.


title: Glider: A GPU Library Driver for Improved System Security

abstract: Legacy device drivers implement both device resource management and isolation. This results in a large code base with a wide high-level interface making the driver vulnerable to security attacks. This is particularly problematic for increasingly popular accelerators like GPUs that have large, complex drivers. We solve this problem with library drivers, a new driver architecture. A library driver implements resource management as an untrusted library in the application process address space, and implements isolation as a kernel module that is smaller and has a narrower lower-level interface (i.e., closer to hardware) than a legacy driver. We articulate a set of device and platform hardware properties that are required to retrofit a legacy driver into a library driver. To demonstrate the feasibility and superiority of library drivers, we present Glider, a library driver implementation for two GPUs of popular brands, Radeon and Intel. Glider reduces the TCB size and attack surface by about 35% and 84% respectively for a Radeon HD 6450 GPU and by about 38% and 90% respectively for an Intel Ivy Bridge GPU. Moreover, it incurs no performance cost. Indeed, Glider outperforms a legacy driver for applications requiring intensive interactions with the device driver, such as applications using the OpenGL immediate mode API.


title: Constraining the Amount of Circumstellar Matter and Dust around Type Ia Supernovae through Near-Infrared Echoes

abstract: The circumstellar (CS) environment is key to understanding progenitors of type Ia supernovae (SNe Ia), as well as the origin of a peculiar extinction property toward SNe Ia for cosmological application. It has been suggested that multiple scatterings of SN photons by CS dust may explain the non-standard reddening law. In this paper, we examine the effect of re-emission of SN photons by CS dust in the infrared (IR) wavelength regime. This effect allows the observed IR light curves to be used as a constraint on the position/size and the amount of CS dust. The method was applied to observed near-infrared (NIR) SN Ia samples; meaningful upper limits on the CS dust mass were derived even under conservative assumptions. We thereby clarify a difficulty associated with the CS dust scattering model as a general explanation for the peculiar reddening law, while it may still apply to a sub-sample of highly reddened SNe Ia. For SNe Ia in general, the environment at the interstellar scale appears to be responsible for the non-standard extinction law. Furthermore, deeper limits can be obtained using the standard nature of SN Ia NIR light curves. In this application, an upper limit of Mdot ~10^{-8}-10^{-7} Msun/yr (for the wind velocity of ~10 km/s) is obtained for a mass loss rate from a progenitor up to ~0.01 pc, and Mdot ~10^{-7}-10^{-6} Msun/yr up to ~0.1 pc.


title: Two-step pulse observation for Raman-Ramsey coherent population trapping atomic clocks

abstract: We propose a two-step pulse observation method to enhance frequency stability for coherent population trapping (CPT) atomic clocks. The proposed method is a Raman-Ramsey scheme with low light intensity at resonance observation, and provides a Ramsey-CPT resonance with both reduced frequency sensitivity to the light intensity and a high signal-to-noise ratio by reducing the repumping into a steady dark state. The resonance characteristics were calculated based on density matrix analysis of a $\Lambda$-type three level system that was modeled on the $^{133}$Cs-D$_1$ line, and the characteristics were also measured using a vertical-cavity surface-emitting laser and a Cs vapor cell.


title: Computational Techniques in FJRW Theory with Applications to Landau-Ginzburg Mirror Symmetry

abstract: The Landau-Ginzburg A-model, given by FJRW theory, defines a cohomological field theory, but in most examples is very difficult to compute, especially when the symmetry group is not maximal. We give some methods for finding the A-model structure. In many cases our methods completely determine the previously unknown A-model Frobenius manifold structure. In the case where these Frobenius manifolds are semisimple, this can be shown to determine the structure of the higher genus potential as well. We compute the Frobenius manifold structure for 27 of the previously unknown unimodal and bimodal singularities and corresponding groups, including 13 cases using a non-maximal symmetry group.


title: On the Structure of the Lambda 1405

abstract: For almost 50 years the structure of the Lambda 1405 resonance has been a mystery. Recently, a new lattice QCD simulation showing that its strange magnetic form factor vanishes, together with a comprehensive Hamiltonian analysis of the lattice QCD energy levels, has unambiguously established that the structure is dominated by a bound anti-kaon--nucleon component [1]. Here we present supplementary information for Ref. [1] including a presentation of the relevant Hamiltonian effective field theory and an illustration of the volume dependence of the results and their connection to the infinite volume limit of Nature.


title: Quantum discord in the central spin model

abstract: The quantum properties of dynamic correlations in a system of an electron spin surrounded by nuclear spins under the conditions of free induction decay (FID) and spin echo have been studied. Analytical results for the time evolution of mutual information, classical part of correlations, and quantum part characterized by quantum discord have been obtained within the central spin model in the high-temperature approximation. The same formulas describe the quantum discord in both the FID and the spin echo although the forms of the dependences are different because of difference in the parameters entering into the formulas. Discord of the spin echo compared with the FID has a strong dependence on time at short times, and it tends to zero with decreasing of the magnetic field, whereas in the case of the FID it reaches a plateau.


title: Dark Matter Annihilation in the First Galaxy Halos

abstract: We investigate the impact of energy released from self-annihilating dark matter on heating of gas in the small, high-redshift dark matter halos thought to host the first stars. A SUSY neutralino like particle is implemented as our dark matter candidate. The PYTHIA code is used to model the final, stable particle distributions produced during the annihilation process. We use an analytic treatment in conjunction with the code MEDEA2 to find the energy transfer and subsequent partition into heating, ionizing and Lyman alpha photon components. We consider a number of halo density models, dark matter particle masses and annihilation channels. We find that the injected energy from dark matter exceeds the binding energy of the gas within a $10^5$ - $10^6$ M$_\odot$ halo at redshifts above 20, preventing star formation in early halos in which primordial gas would otherwise cool. Thus we find that DM annihilation could delay the formation of the first galaxies.


title: Deep Narrow Boltzmann Machines are Universal Approximators

abstract: We show that deep narrow Boltzmann machines are universal approximators of probability distributions on the activities of their visible units, provided they have sufficiently many hidden layers, each containing the same number of units as the visible layer. We show that, within certain parameter domains, deep Boltzmann machines can be studied as feedforward networks. We provide upper and lower bounds on the sufficient depth and width of universal approximators. These results settle various intuitions regarding undirected networks and, in particular, they show that deep narrow Boltzmann machines are at least as compact universal approximators as narrow sigmoid belief networks and restricted Boltzmann machines, with respect to the currently available bounds for those models.


title: Doyen-Wilson results for odd length cycle systems

abstract: For each odd $m \geq 3$ we completely solve the problem of when an $m$-cycle system of order $u$ can be embedded in an $m$-cycle system of order $v$, barring a finite number of possible exceptions. In cases where $u$ is large compared to $m$, where $m$ is a prime power, or where $m \leq 15$, the problem is completely resolved. In other cases, the only possible exceptions occur when $v-u$ is small compared to $m$. This result is proved as a consequence of a more general result which gives necessary and sufficient conditions for the existence of an $m$-cycle decomposition of a complete graph of order $v$ with a hole of size $u$ in the case where $u \geq m-2$ and $v-u \geq m+1$ both hold.


title: Bayesian Analysis for Stellar Evolution with Nine Parameters (BASE-9): User's Manual

abstract: BASE-9 is a Bayesian software suite that recovers star cluster and stellar parameters from photometry. BASE-9 is useful for analyzing single-age, single-metallicity star clusters, binaries, or single stars, and for simulating such systems. BASE-9 uses Markov chain Monte Carlo and brute-force numerical integration techniques to estimate the posterior probability distributions for the age, metallicity, helium abundance, distance modulus, and line-of-sight absorption for a cluster, and the mass, binary mass ratio, and cluster membership probability for every stellar object. BASE-9 is provided as open source code on a version-controlled web server. The executables are also available as Amazon Elastic Compute Cloud images. This manual provides potential users with an overview of BASE-9, including instructions for installation and use.


title: Asymmetric Minwise Hashing

abstract: Minwise hashing (Minhash) is a widely popular indexing scheme in practice. Minhash is designed for estimating set resemblance and is known to be suboptimal in many applications where the desired measure is set overlap (i.e., inner product between binary vectors) or set containment. Minhash has inherent bias towards smaller sets, which adversely affects its performance in applications where such a penalization is not desirable. In this paper, we propose asymmetric minwise hashing (MH-ALSH), to provide a solution to this problem. The new scheme utilizes asymmetric transformations to cancel the bias of traditional minhash towards smaller sets, making the final "collision probability" monotonic in the inner product. Our theoretical comparisons show that for the task of retrieving with binary inner products asymmetric minhash is provably better than traditional minhash and other recently proposed hashing algorithms for general inner products. Thus, we obtain an algorithmic improvement over existing approaches in the literature. Experimental evaluations on four publicly available high-dimensional datasets validate our claims and the proposed scheme outperforms, often significantly, other hashing algorithms on the task of near neighbor retrieval with set containment. Our proposal is simple and easy to implement in practice.


title: Weight modules for current algebras

abstract: For any finite-dimensional simple Lie algebra $\mathfrak{g}$ and commutative associative algebra $S$ of finite type, we give a complete classification of the simple weight modules of $\mathfrak{g}\otimes S$ with bounded weight multiplicities.


title: Kinematics of classical Cepheids in the Nuclear Stellar Disk

abstract: Classical Cepheids are useful tracers of the Galactic young stellar population because their distances and ages can be determined from their period-luminosity and period-age relations. In addition, the radial velocities and chemical abundance of the Cepheids can be derived from spectroscopic observations, providing further insights into the structure and evolution of the Galaxy. Here, we report the radial velocities of classical Cepheids near the Galactic Center, three of which were reported in 2011, the other reported for the first time. The velocities of these Cepheids suggest that the stars orbit within the Nuclear Stellar Disk, a group of stars and interstellar matter occupying a region of 200 pc around the Center, although the three-dimensional velocities cannot be determined until the proper motions are known. According to our simulation, these four Cepheids formed within the Nuclear Stellar Disk like younger stars and stellar clusters therein.


title: Monotonic Abstraction Techniques: from Parametric to Software Model Checking

abstract: Monotonic abstraction is a technique introduced in model checking parameterized distributed systems in order to cope with transitions containing global conditions within guards. The technique has been re-interpreted in a declarative setting in previous papers of ours and applied to the verification of fault tolerant systems under the so-called "stopping failures" model. The declarative reinterpretation consists in logical techniques (quantifier relativizations and, especially, quantifier instantiations) making sense in a broader context. In fact, we recently showed that such techniques can over-approximate array accelerations, so that they can be employed as a meaningful (and practically effective) component of CEGAR loops in software model checking too.


title: Choreographies and Behavioural Contracts on the Way to Dynamic Updates

abstract: We survey our work on choreographies and behavioural contracts in multiparty interactions. In particular theories of behavioural contracts are presented which enable reasoning about correct service composition (contract compliance) and service substitutability (contract refinement preorder) under different assumptions concerning service communication: synchronous address or name based communication with patient non-preemptable or impatient invocations, or asynchronous communication. Correspondingly relations between behavioural contracts and choreographic descriptions are considered, where a contract for each communicating party is, e.g., derived by projection. The considered relations are induced as the maximal preoders which preserve contract compliance and global traces: we show maximality to hold (permitting services to be discovered/substituted independently for each party) when contract refinement preorders with all the above asymmetric communication means are considered and, instead, not to hold if the standard symmetric CCS/pi-calculus communication is considered (or when directly relating choreographies to behavioral contracts via a preorder, no matter the communication mean). The obtained maximal preorders are then characterized in terms of a new form of testing, called compliance testing, where not only tests must succeed but also the system under test (thus relating to controllability theory), and compared with classical preorders such as may/must testing, trace inclusion, etc. Finally, recent work about adaptable choreographies and behavioural contracts is presented, where the theory above is extended to update mechanisms allowing choreographies/contracts to be modified at run-time by internal (self-adaptation) or external intervention.


title: An Approach to Model Checking of Multi-agent Data Analysis

abstract: The paper presents an approach to verification of a multi-agent data analysis algorithm. We base correct simulation of the multi-agent system by a finite integer model. For verification we use model checking tool SPIN. Protocols of agents are written in Promela language and properties of the multi-agent data analysis system are expressed in logic LTL. We run several experiments with SPIN and the model.


title: A Language Support for Exhaustive Fault-Injection in Message-Passing System Models

abstract: This paper presents an approach towards specifying and verifying adaptive distributed systems. We here take fault-handling as an example of adaptive behavior and propose a modeling language Sandal for describing fault-prone message-passing systems. One of the unique mechanisms of the language is a linguistic support for abstracting typical faults such as unexpected termination of processes and random loss of messages. The Sandal compiler translates a model into a set of NuSMV modules. During the compilation process, faults specified in the model will be woven into the output. One can thus enjoy full-automatic exhaustive fault-injection without writing faulty behaviors explicitly. We demonstrate the advantage of the language by verifying a model of the two-phase commit protocol under faulty environment.


title: 3-D Reciprocal Collision Avoidance on Physical Quadrotor Helicopters with On-Board Sensing for Relative Positioning

abstract: In this paper, we present an implementation of 3-D reciprocal collision avoidance on real quadrotor helicopters where each quadrotor senses the relative position and velocity of other quadrotors using an on-board camera. We show that using our approach, quadrotors are able to successfully avoid pairwise collisions in GPS and motion-capture denied environments, without communication between the quadrotors, and even when human operators deliberately attempt to induce collisions. To our knowledge, this is the first time that reciprocal collision avoidance has been successfully implemented on real robots where each agent independently observes the others using on-board sensors. We theoretically analyze the response of the collision-avoidance algorithm to the violated assumptions by the use of real robots. We quantitatively analyze our experimental results. A particularly striking observation is that at times the quadrotors exhibit "reciprocal dance" behavior, which is also observed when humans move past each other in constrained environments. This seems to be the result of sensing uncertainty, which causes both robots involved to have a different belief about the relative positions and velocities and, as a result, choose the same side on which to pass.


title: The Kohn-Luttinger superconductivity in idealized doped graphene

abstract: Idealized graphene monolayer is considered neglecting the van der Waals potential of the substrate and the role of the nonmagnetic impurities. The effect of the long-range Coulomb repulsion in an ensemble of Dirac fermions on the formation of the superconducting pairing in a monolayer is studied in the framework of the Kohn-Luttinger mechanism. The electronic structure of graphene is described in the strong coupling Wannier representation on the hexagonal lattice. We use the Shubin-Vonsowsky model which takes into account the intra- and intersite Coulomb repulsions of electrons. The Cooper instability is established by solving the Bethe-Salpeter integral equation, in which the role of the effective interaction is played by the renormalized scattering amplitude. The renormalized amplitude contains the Kohn-Luttinger polarization contributions up to and including the second-order terms in the Coulomb repulsion. We construct the superconductive phase diagram for the idealized graphene monolayer and show that the Kohn-Luttinger renormalizations and the intersite Coulomb repulsion significantly affect the interplay between the superconducting phases with $f-$, $d+id-$, and $p+ip-$wave symmetries of the order parameter.


title: Autoregressive Cascades on Random Networks

abstract: This paper considers a model for cascades on random networks in which the cascade propagation at any node depends on the load at the failed neighbor, the degree of the neighbor as well as the load at that node. Each node in the network bears an initial load that is below the capacity of the node. The trigger for the cascade emanates at a single node or a small fraction of the nodes from some external shock. Upon failure, the load at the failed node gets divided randomly and added to the existing load at those neighboring nodes that have not yet failed. Subsequently, a neighboring node fails if its accumulated load exceeds its capacity. The failed node then plays no further part in the process. The cascade process stops as soon as the accumulated load at all nodes that have not yet failed is below their respective capacities. The model is shown to operate in two regimes, one in which the cascade terminates with only a finite number of node failures. In the other regime there is a positive probability that the cascade continues indefinitely. Bounds are obtained on the critical parameter where the phase transition occurs.


title: A direct algorithm of one-dimensional optimal system for the group invariant solutions

abstract: A direct and systematic algorithm is proposed to find one-dimensional optimal system for the group invariant solutions, which is attributed to the classification of its corresponding one-dimensional Lie algebra. Since the method is based on different values of all the invariants, the process itself can both guarantee the comprehensiveness and demonstrate the inequivalence of the optimal system, with no further proof. To illustrate our method more clearly , we give a couple of well-known examples: the Korteweg-de Vries (KdV) equation and the heat equation.


title: The construction of two-dimensional optimal systems for the invariant solutions

abstract: To search for inequivalent group invariant solutions, a general and systematic approach is established to construct two-dimensional optimal systems, which is based on commutator relations, adjoint matrix and the invariants. The details of computing all the invariants for two-dimensional subalgebras is presented and the optimality of twodimensional optimal systems is shown clearly under different values of invariants, with no further proof. Applying the algorithm to (1+1)-dimensional heat equation and (2+1)-dimensional Navier-Stokes (NS) equation, their twodimensional optimal systems are obtained, respectively. For the heat equation, eleven two-parameter elements in the optimal system are found one by one, which are discovered more comprehensive. The two-dimensional optimal system of NS equations is used to generate intrinsically different reduced ordinary differential equations and some interesting explicit solutions are provided.


title: Query complexity of sampling and small geometric partitions

abstract: In this paper we study the following problem: Discrete partitioning problem (DPP): Let $\mathbb{F}_q P^n$ denote the $n$-dimensional finite projective space over $\mathbb{F}_q$. For positive integer $k \leq n$, let $\{ A^i\}_{i=1}^N$ be a partition of $(\mathbb{F}_q P^n)^k$ such that (1) for all $i \leq N$, $A^i = \prod_{j=1}^k A^i_j$ (partition into product sets), (2) for all $i \leq N$, there is a $(k-1)$-dimensional subspace $L^i \subseteq \mathbb{F}_q P^n$ such that $A^i \subseteq (L^i)^k$. What is the minimum value of $N$ as a function of $q,n,k$? We will be mainly interested in the case $k=n$.


title: A Sharp First Order Analysis of Feynman-Kac Particle Models

abstract: This article provides a new theory for the analysis of forward and backward particle approximations of Feynman-Kac models. Such formulae are found in a wide variety of applications and their numerical (particle) approximation are required due to their intractability. Under mild assumptions, we provide sharp and non-asymptotic first order expansions of these particle methods, potentially on path space and for possibly unbounded functions. These expansions allows one to consider upper and lower bound bias type estimates for a given time horizon $n$ and particle number $N$; these non-asymptotic estimates are of order $\mathcal{O}(n/N)$. Our approach is extended to tensor products of particle density profiles, leading to new sharp and non-asymptotic propagation of chaos estimates. The resulting upper and lower bound propagation of chaos estimates seems to be the first result of this kind for mean field particle models. As a by-product of our results, we also provide some analysis of the particle Gibbs sampler, providing first order expansions of the kernel and minorization estimates.


title: Fixation times in differentiation and evolution in the presence of bottlenecks, deserts, and oases

abstract: Cellular differentiation and evolution are stochastic processes that can involve multiple types (or states) of particles moving on a complex, high-dimensional state-space or "fitness" landscape. Cells of each specific type can thus be quantified by their population at a corresponding node within a network of states. Their dynamics across the state-space network involve genotypic or phenotypic transitions that can occur upon cell division, such as during symmetric or asymmetric cell differentiation, or upon spontaneous mutation. Waiting times between transitions can be nonexponentially distributed and reflect e.g., the cell cycle. Here, we use a multi-type branching processes to study first passage time statistics for a single cell to appear in a specific state. We present results for a sequential evolutionary process in which $L$ successive transitions propel a population from a "wild-type" state to a given "terminally differentiated," "resistant," or "cancerous" state. Analytic and numeric results are also found for first passage times across an evolutionary chain containing a node with increased death or proliferation rate, representing a desert/bottleneck or an oasis. Processes involving cell proliferation are shown to be "nonlinear" (even though mean-field equations for the expected particle numbers are linear) resulting in first passage time statistics that depend on the position of the bottleneck or oasis. Our results highlight the sensitivity of stochastic measures to cell division fate and quantify the limitations of using certain approximations and assumptions (such as fixed-population and mean-field assumptions) in evaluating fixation times.


title: Evolution of angular-momentum-losing exoplanetary systems : Revisiting Darwin stability

abstract: We assess the importance of tidal evolution and its interplay with magnetic braking in the population of hot-Jupiter planetary systems. By minimizing the total mechanical energy of a given system under the constraint of stellar angular momentum loss, we rigorously find the conditions for the existence of dynamical equilibrium states. We estimate their duration, in particular when the wind torque spinning down the star is almost compensated by the tidal torque spinning it up. We introduce dimensionless variables to characterize the tidal evolution of observed hot Jupiter systems and discuss their spin and orbital states using generalized Darwin diagrams based on our new approach. We show that their orbital properties are related to the effective temperature of their host stars. The long-term evolution of planets orbiting F- and G-type stars is significantly different owing to the combined effect of magnetic braking and tidal dissipation. The existence of a quasi-stationary state, in the case of short-period planets, can significantly delay their tidal evolution that would otherwise bring the planet to fall into its host star. Most of the planets known to orbit F-type stars are presently found to be near this stationary state, probably in a configuration not too far from that they had when their host star settled on the zero-age main sequence. Considering the importance of angular momentum loss in the early stages of stellar evolution, our results indicate that it has to be taken into account also to properly test the migration scenarios of planetary system formation.


title: Stochastic Compositional Gradient Descent: Algorithms for Minimizing Compositions of Expected-Value Functions

abstract: Classical stochastic gradient methods are well suited for minimizing expected-value objective functions. However, they do not apply to the minimization of a nonlinear function involving expected values or a composition of two expected-value functions, i.e., problems of the form $\min_x \mathbf{E}_v [f_v\big(\mathbf{E}_w [g_w(x)]\big)]$. In order to solve this stochastic composition problem, we propose a class of stochastic compositional gradient descent (SCGD) algorithms that can be viewed as stochastic versions of quasi-gradient method. SCGD update the solutions based on noisy sample gradients of $f_v,g_{w}$ and use an auxiliary variable to track the unknown quantity $\mathbf{E}_w[g_w(x)]$. We prove that the SCGD converge almost surely to an optimal solution for convex optimization problems, as long as such a solution exists. The convergence involves the interplay of two iterations with different time scales. For nonsmooth convex problems, the SCGD achieve a convergence rate of $O(k^{-1/4})$ in the general case and $O(k^{-2/3})$ in the strongly convex case, after taking $k$ samples. For smooth convex problems, the SCGD can be accelerated to converge at a rate of $O(k^{-2/7})$ in the general case and $O(k^{-4/5})$ in the strongly convex case. For nonconvex problems, we prove that any limit point generated by SCGD is a stationary point, for which we also provide the convergence rate analysis. Indeed, the stochastic setting where one wants to optimize compositions of expected-value functions is very common in practice. The proposed SCGD methods find wide applications in learning, estimation, dynamic programming, etc.


title: Kohn-Sham Kinetic Energy Density in the Nuclear and Asymptotic Regions: Deviations from the Von Weizs\"acker Behavior and Applications to Density Functionals

abstract: We show that the Kohn-Sham positive-definite kinetic energy (KE) density significantly differs from the von Weizs\"acker (VW) one at the nuclear cusp as well as in the asymptotic region. At the nuclear cusp, the VW functional is shown to be linear and the contribution of p-type orbitals to the KE density is theoretically derived and numerically demonstrated in the limit of infinite nuclear charge, as well in the semiclassical limit of neutral large atoms. In the latter case, it reaches 12 of the KE density. In the asymptotic region we find new exact constraints for meta Generalized Gradient Approximation (meta-GGA) exchange functionals: with an exchange enhancement factor proportional to $\sqrt{\alpha}$, where $\alpha$ is the common meta-GGA ingredient, both the exchange energy density and the potential are proportional to the exact ones. In addition, this describes exactly the large-gradient limit of quasi-two dimensional systems.


title: Metric-Independent Measures for Supersymmetric Extended Object Theories on Curved Backgrounds

abstract: For Green-Schwarz superstring sigma-model on curved backgrounds, we introduce a non-metric measure $\Phi \equiv \epsilon^{i j} \epsilon^{I J} (\partial_i \varphi^I) (\partial_j \varphi^J)$ with two scalars $\varphi^I (I = 1, 2)$ used in Two Measure Theory (TMT). As in the flat-background case, the string tension $T= (2 \pi \alpha ' )^{-1}$ emerges as an integration constant for the A_i-field equation. This mechanism is further generalized to supermembrane theory, and to super p-brane theory, both on general curved backgrounds. This shows the universal applications of dynamical measure of TMT to general supersymmetric extended objects on general curved backgrounds.


title: Integrating Fuzzy and Ant Colony System for Fuzzy Vehicle Routing Problem with Time Windows

abstract: In this paper fuzzy VRPTW with an uncertain travel time is considered. Credibility theory is used to model the problem and specifies a preference index at which it is desired that the travel times to reach the customers fall into their time windows. We propose the integration of fuzzy and ant colony system based evolutionary algorithm to solve the problem while preserving the constraints. Computational results for certain benchmark problems having short and long time horizons are presented to show the effectiveness of the algorithm. Comparison between different preferences indexes have been obtained to help the user in making suitable decisions.


title: A cluster algorithm for Monte Carlo simulations of spin ice

abstract: We present an algorithm for Monte Carlo simulations of a nearest-neighbor spin ice model based on its cluster representation. To assess its performance, we estimate a relaxation time, and find that, in contrast to the Metropolis algorithm, our algorithm does not develop spin-freezing. Also, to demonstrate the efficiency, we calculate the spin and charge structure factors, and observe pinch points in a high-resolution color map. We then find that Debye screening works among defects and brings about short-range correlations, and that the deconfinement transition triggered by a fugacity of defects $z$ is dictated by a singular part of the free-energy density $f_{\rm s}\propto z$.


title: Tracking Triadic Cardinality Distributions for Burst Detection in Social Activity Streams

abstract: In everyday life, we often observe unusually frequent interactions among people before or during important events, e.g., we receive/send more greetings from/to our friends on Christmas Day, than usual. We also observe that some videos suddenly go viral through people's sharing in online social networks (OSNs). Do these seemingly different phenomena share a common structure? All these phenomena are associated with sudden surges of user activities in networks, which we call "bursts" in this work. We find that the emergence of a burst is accompanied with the formation of triangles in networks. This finding motivates us to propose a new method to detect bursts in OSNs. We first introduce a new measure, "triadic cardinality distribution", corresponding to the fractions of nodes with different numbers of triangles, i.e., triadic cardinalities, within a network. We demonstrate that this distribution changes when a burst occurs, and is naturally immunized against spamming social-bot attacks. Hence, by tracking triadic cardinality distributions, we can reliably detect bursts in OSNs. To avoid handling massive activity data generated by OSN users, we design an efficient sample-estimate solution to estimate the triadic cardinality distribution from sampled data. Extensive experiments conducted on real data demonstrate the usefulness of this triadic cardinality distribution and the effectiveness of our sample-estimate solution.


title: Parametric instability of classical Yang-Mills fields in a color magnetic background

abstract: We investigate instabilities of classical Yang-Mills fields in a time-dependent spatially homogeneous color magnetic background field in a non-expanding geometry for elucidating the earliest stage dynamics of ultra-relativistic heavy-ion collisions. The background gauge field configuration considered in this article is spatially homogeneous and temporally periodic, and is alluded by Berges-Scheffler-Schlichting-Sexty (BSSS). We discuss the whole structure of instability bands of fluctuations around the BSSS background gauge field on the basis of Floquet theory, which enables us to discuss the stability in a systematic way. We find various instability bands on the $(p_z, p_T)$-plane. These instability bands are caused by parametric resonance despite the fact that the momentum dependence of the growth rate for $|\mathbf{p}| \leq \sqrt{B}$ is similar to a Nielsen-Olesen instability. Moreover, some of instability bands are found to emerge not only in the low momentum but also in the high momentum region; typically of the order of the saturation momentum as $|\mathbf{p}| \sim \sqrt{B} \sim Q_{\rm s}$.


title: Fundamental Limits of Blind Deconvolution Part I: Ambiguity Kernel

abstract: Blind deconvolution is an ubiquitous non-linear inverse problem in applications like wireless communications and image processing. This problem is generally ill-posed, and there have been efforts to use sparse models for regularizing blind deconvolution to promote signal identifiability. Part I of this two-part paper characterizes the ambiguity space of blind deconvolution and shows unidentifiability of this inverse problem for almost every pair of unconstrained input signals. The approach involves lifting the deconvolution problem to a rank one matrix recovery problem and analyzing the rank two null space of the resultant linear operator. A measure theoretically tight (parametric and recursive) representation of the key rank two null space is stated and proved. This representation is a novel foundational result for signal and code design strategies promoting identifiability under convolutive observation models. Part II of this paper analyzes the identifiability of sparsity constrained blind deconvolution and establishes surprisingly strong negative results on scaling laws for the sparsity-ambiguity trade-off.


title: Performance Evaluation of Virtualized Hadoop Clusters

abstract: In this report we investigate the performance of Hadoop clusters, deployed with separated storage and compute layers, on top of a hypervisor managing a single physical host. We have analyzed and evaluated the different Hadoop cluster configurations by running CPU bound and I/O bound workloads. The report is structured as follows: Section 2 provides a brief description of the technologies involved in our study. An overview of the experimental platform, setup test and configurations are presented in Section 3. Our benchmark methodology is defined in Section 4. The performed experiments together with the evaluation of the results are presented in Section 5. Finally, Section 6 concludes with lessons learned.


title: Spectrally Accurate Causality Enforcement using SVD-based Fourier Continuations for High Speed Digital Interconnects

abstract: We introduce an accurate and robust technique for accessing causality of network transfer functions given in the form of bandlimited discrete frequency responses. These transfer functions are commonly used to represent the electrical response of high speed digital interconnects used on chip and in electronic package assemblies. In some cases small errors in the model development lead to non-causal behavior that does not accurately represent the electrical response and may lead to a lack of convergence in simulations that utilize these models. The approach is based on Hilbert transform relations or Kramers-Kronig dispersion relations and a construction of causal Fourier continuations using a regularized singular value decomposition (SVD) method. Given a transfer function, non-periodic in general, this procedure constructs highly accurate Fourier series approximations on the given frequency interval by allowing the function to be periodic in an extended domain. The causality dispersion relations are enforced spectrally and exactly. This eliminates the necessity of approximating the transfer function behavior at infinity and explicit computation of the Hilbert transform. We perform the error analysis of the method and take into account a possible presence of a noise or approximation errors in data. The developed error estimates can be used in verifying causality of the given data. The performance of the method is tested on several analytic and simulated examples that demonstrate an excellent accuracy and reliability of the proposed technique in agreement with the obtained error estimates. The method is capable of detecting very small localized causality violations with amplitudes close to the machine precision.


title: Matter Power Spectra in Viable $f(R)$ Gravity Models with Massive Neutrinos

abstract: We investigate the matter power spectra in the power law and exponential types of viable $f(R)$ theories along with massive neutrinos. The enhancement of the matter power spectrum is found to be a generic feature in these models. In particular, we show that in the former type, such as the Starobinsky model, the spectrum is magnified much larger than the latter one, such as the exponential model. A greater scale of the total neutrino mass, $\Sigma m_{\nu}$, is allowed in the viable $f(R)$ models than that in the $\Lambda$CDM one. We obtain the constraints on the neutrino masses by using the CosmoMC package with the modified MGCAMB. Explicitly, we get $\Sigma m_{\nu} < 0.451 \ (0.214)\ \mathrm{eV}$ at 95% C.L. in the Starobinsky (exponential) model, while the corresponding one for the $\Lambda$CDM model is $\Sigma m_{\nu} < 0.200\ \mathrm{eV}$. Furthermore, by treating the effective number of neutrino species $N_{\mathrm{eff}}$ as a free parameter along with $\Sigma m_{\nu}$, we find that $N_{\mathrm{eff}} = 3.78^{+0.64}_{-0.84} (3.47^{+0.74}_{-0.60})$ and $\Sigma m_{\nu} = 0.533^{+0.254}_{-0.411}$ ($< 0.386) \ \mathrm{eV}$ at 95% C.L. in the Starobinsky (exponential) model.


title: The Normality of Certain Varieties of Special Lattices

abstract: We begin with a short exposition of the theory of lattice varieties. This includes a description of their orbit structure and smooth locus. We construct a flat cover of the lattice variety and show that it is a complete intersection. We show that the lattice variety is locally a complete intersection and nonsingular in codimension one and hence normal. We then prove a comparison theorem showing that this theory becomes parallel to the function field case if linear algebra is replaced by $p$-linear algebra. We then compute the Lie algebra of the special linear group over the truncated Witt vectors. We conclude by applying these results to show how to describe the canonical bundle on the lattice variety and we use the description to show that lattice varieties are not isomorphic to the analogous objects in the affine Grassmannian.


title: Predictive Encoding of Contextual Relationships for Perceptual Inference, Interpolation and Prediction

abstract: We propose a new neurally-inspired model that can learn to encode the global relationship context of visual events across time and space and to use the contextual information to modulate the analysis by synthesis process in a predictive coding framework. The model learns latent contextual representations by maximizing the predictability of visual events based on local and global contextual information through both top-down and bottom-up processes. In contrast to standard predictive coding models, the prediction error in this model is used to update the contextual representation but does not alter the feedforward input for the next layer, and is thus more consistent with neurophysiological observations. We establish the computational feasibility of this model by demonstrating its ability in several aspects. We show that our model can outperform state-of-art performances of gated Boltzmann machines (GBM) in estimation of contextual information. Our model can also interpolate missing events or predict future events in image sequences while simultaneously estimating contextual information. We show it achieves state-of-art performances in terms of prediction accuracy in a variety of tasks and possesses the ability to interpolate missing frames, a function that is lacking in GBM.


title: Monte Carlo error analyses of Spearman's rank test

abstract: Spearman's rank correlation test is commonly used in astronomy to discern whether a set of two variables are correlated or not. Unlike most other quantities quoted in astronomical literature, the Spearman's rank correlation coefficient is generally quoted with no attempt to estimate the errors on its value. This is a practice that would not be accepted for those other quantities, as it is often regarded that an estimate of a quantity without an estimate of its associated uncertainties is meaningless. This manuscript describes a number of easily implemented, Monte Carlo based methods to estimate the uncertainty on the Spearman's rank correlation coefficient, or more precisely to estimate its probability distribution.