title: Characterizing the finiteness of the Hausdorff distance between two algebraic curves

abstract: In this paper, we present a characterization for the Hausdorff distance between two given algebraic curves in the $n$-dimensional space (parametrically or implicitly defined) to be finite. The characterization is related with the asymptotic behavior of the two curves and it can be easily checked. More precisely, the Hausdorff distance between two curves $\cal C$ and $\overline{\cal C}$ is finite if and only if for each infinity branch of $\cal C$ there exists an infinity branch of $\overline{\cal C}$ such that the terms with positive exponent in the corresponding series are the same, and reciprocally.


title: Extension from Precoloured Sets of Edges

abstract: We consider precolouring extension problems for proper edge-colourings of graphs and multigraphs, in an attempt to prove stronger versions of Vizing's and Shannon's bounds on the chromatic index of (multi)graphs in terms of their maximum degree $\Delta$. We are especially interested in the following question: when is it possible to extend a precoloured matching to a colouring of all edges of a (multi)graph? This question turns out to be related to the notorious List Colouring Conjecture and other classic notions of choosability.


title: Kinetics of Final Degassing of Hydrogen Desorption by Metal Hydrides

abstract: The proposed model concerns the 'confluent shrinking core' scenario and reproduces the desorption kinetic after the complete decay of the stoichiometric hydride ($\beta$-phase). The exact analytical solution is obtained, the numerical values are demonstrated by the example of magnesium hydride.


title: On the Hochschild cohomology ring of the quaternion group of order eight in characteristic two

abstract: Let $k$ be an algebraically closed field of characteristic two and let $Q_8$ be the quaternion group of order $8$. We determine the Gerstenhaber Lie algebra structure and the Batalin-Vilkovisky structure on the Hochschild cohomology ring of the group algebra $kQ_8$.


title: Non-binary LDPC decoding using truncated messages in the Walsh-Hadamard domain

abstract: The Extended Min-Sum (EMS) algorithm for non-binary low-density parity-check (LDPC) defined over an alphabet of size $q$ operates on truncated messages of length $q'$ to achieve a complexity of the order $q'^2$. In contrast, Walsh-Hadamard (WH) transform based iterative decoders achieve a complexity of the order $q\log q$, which is much larger for $q'<<q$. In this paper, we demonstrate that considerable savings can be achieved by letting WH based decoders operate on truncated messages as well. We concentrate on the direct WH transform and compute the number of operations required if only $q'$ of the $q$ inputs are non-zero. Our paper does not cover the inverse WH transform and hence further research is needed to construct WH based decoders that can compete with the EMS algorithm on complexity terms.


title: Quantum Speed Limits for General Physical Processes

abstract: Quantum speed limits are relations yielding lower bounds on the evolution time of quantum systems. These results have been generalized in some ways, in particular by including evolutions to non-orthogonal states. However, there was a gap in the literature on this area, for only unitary evolutions -- closed quantum systems -- had been considered. On this Ph.D. thesis, such limitation is overcome: our main result is a bound for quantum-system evolutions in general, whether unitary or not, and correctly recovers the known bounds in the unitary case. Applications of this bound to several concrete cases of interest are herein presented. This bound is also used to extend to the non-unitary case the discussion of the role of entanglement in fast evolutions, leading to nontrivial results. For the derivation of the results, a geometric approach has been employed, which allows a clear interpretation of the bounds and a discussion of the criteria for their saturation. No previous knowledge of quantum-state geometry by the reader has been assumed.


title: Bistable switching asymptotics for the self regulating gene

abstract: A simple stochastic model of a self regulating gene that displays bistable switching is analyzed. While on, a gene transcribes mRNA at a constant rate. Transcription factors can bind to the DNA and affect the gene's transcription rate. Before an mRNA is degraded, it synthesizes protein, which in turn regulates gene activity by influencing the activity of transcription factors. Protein is slowly removed from the system through degradation. Depending on how the protein regulates gene activity, the protein concentration can exhibit noise induced bistable switching. An asymptotic approximation of the mean switching rate is derived that includes the pre exponential factor, which improves upon a previously reported logarithmically accurate approximation. With the improved accuracy, a uniformly accurate approximation of the stationary probability density, describing the gene, mRNA copy number, and protein concentration is also obtained.


title: In Quantum Computing Speedup Illusory?: The False Coin of "Counting Function Evaluations"

abstract: By using a new way to encode Boolean functions in a reversible gate, an algorithm is developed in quantum computing over Z_2, symbolized QC/2, (as opposed to QC over C) that needs only one function evaluation to solve the Grover Database Search Problem of finding a designated record among 2^m records for any m. In the usual Grover algorithm in quantum computing over C, one needs essentially Sqrt(2^m) function evaluations as opposed to the average of (2^m)/2 functions evaluations needed in the classical algorithm. The one function evaluation of the QC/2 algorithm (for any m) represents such a super speedup, even over the Grover algorithm in QC/C, that one feels something has gone awry. Indeed, our analysis of the transparent calculations of Boolean functions over Z_2 shows that the classical algorithm is just repackaged in a rather obvious way in the single function evaluation of the QC/2 algorithm--whereas the calculations are hidden and non-transparent in the Grover QC/C algorithm using C. The conclusion in both cases (which is rather obvious in the QC/2 case) is that "counting function evaluations" is a false coin to measure speedup in the comparison between quantum and classical computing.


title: Faults in Linux 2.6

abstract: In August 2011, Linux entered its third decade. Ten years before, Chou et al. published a study of faults found by applying a static analyzer to Linux versions 1.0 through 2.4.1. A major result of their work was that the drivers directory contained up to 7 times more of certain kinds of faults than other directories. This result inspired numerous efforts on improving the reliability of driver code. Today, Linux is used in a wider range of environments, provides a wider range of services, and has adopted a new development and release model. What has been the impact of these changes on code quality? To answer this question, we have transported Chou et al.'s experiments to all versions of Linux 2.6; released between 2003 and 2011. We find that Linux has more than doubled in size during this period, but the number of faults per line of code has been decreasing. Moreover, the fault rate of drivers is now below that of other directories, such as arch. These results can guide further development and research efforts for the decade to come. To allow updating these results as Linux evolves, we define our experimental protocol and make our checkers available.


title: The flat density profiles of massive, and relaxed galaxy clusters

abstract: In the present paper, we studied by means of the SIM introduced in Del Popolo (2009), the total and DM density profiles, and the correlations among different quantities, observed by Newman et al. (2012a,b), in seven massive and relaxed clusters, namely MS2137, A963, A383, A611, A2537, A2667, A2390. Similarly to Newman et al. (2012a,b), the total density profile, in the radius range 0.003 - 0.03$r_{200}$, has a mean total density profile in agreement with dissipationless simulations. The slope of the DM profiles of all clusters is flatter than -1. The slope, $\alpha$, has a maximum value (including errors) of $\alpha=-0.88$ in the case of A2390, and minimum value $\alpha=-0.14$ for A2537. The baryonic component dominates the mass distribution at radii $< 5-10$ kpc, while the outer distribution is dark matter dominated. We found an anti-correlation among the slope $\alpha$, the effective radius, $R_e$, and the BCG mass, and a correlation among the core radius $r_{core}$, and $R_e$. Moreover, the mass in 100 kpc (mainly dark matter) is correlated with the mass inside 5 kpc (mainly baryons). The behavior of the total mass density profile, the DM density profile, and the quoted correlations can be understood in a double phase scenario. In the first dissipative phase the proto-BCG forms, and in the second dissipationless phase, dynamical friction between baryonic clumps (collapsing to the center) and the DM halo flattens the inner slope of the density profile. In simple terms, the large scatter in the inner slope from cluster to cluster, and the anti-correlation among the slope, $\alpha$ and $R_e$ is due to the fact that in order to have a total mass density profile which is NFW-like, clusters having more massive BCGs at their centers must contain less DM in their center. Consequently the inner profile has a flatter slope.


title: Interacting parametrized post-Friedmann method

abstract: We apply the interacting parametrized post-Friedmann (IPPF) method to a coupled dark energy model where the interaction is proportional to dark matter density at background level. In doing so, we perform a Markov Chain Monte-Carlo analysis which combines several cosmological probes including the cosmic microwave background (WMAP9+Planck) data, baryon acoustic oscillation (BAO) measurements, JLA sample of supernovae, Hubble constant (HST), and redshift-space distortion (RSD) measurements through the ${\rm f}\sigma_{8}{\rm (z)}$ data points. The joint observational analysis of ${\rm Planck+WP+JLA+BAO+HST+RSD}$ data leads to a coupling parameter, $\xi_{c}=0.00140_{-0.00080}^{+0.00079}$ at $1\sigma$ level for vanishing momentum transfer potential; this value is reduced a when the momentum transfer potential is switched on, giving $\xi_{c}=0.00136_{-0.00073}^{+0.00080}$ at $1\sigma$ level. The CMB power spectrum shows up a correlation between the coupling parameter $\xi_{c}$ and the position of acoustic peaks or their amplitudes. The first peak's height increases when $\xi_{c}$ takes larger values and its position is shifted. We also obtain the matter power spectrum may be affected by the strength of interaction coupling over scales bigger that $10^{-2} {\rm h~ Mpc^{-1}}$, reducing its amplitude in relation to the vanilla model.


title: Neutron scattering studies on URu2Si2

abstract: This paper is aiming to review some of the neutron scattering studies performed on URu2Si2 in Grenoble. This compound has been studied for a quarter of century because of a so-called hidden order ground state visible by most of the bulk experiments but almost invisible by microscopic probes like neutrons, muons NMR or x-ray. We stress on some aspects that were not addressed previously. Firstly, the comparison of the cell parameters in the 1-2-2 systems seems to point that the magnetic properties of URu2Si2 are leading by an U4+ electronic state. Secondly, a compilation of the different studies of the tiny antiferromagnetic moment indicates that the tiny antiferromagnetic moment has a constant value which may indicate that it is not necessary extrinsic. We also present the last development on the magnetic form factor measurement in which the magnetic density rotates when entering in the hidden order state. To end, the thermal dependence of the two most intense magnetic excitation at Q0=(1,0,0) and Q1=(0.6,0,0) seems to indicate two different origins or processes for these excitations.


title: Efficient Generation of Model Bulk Heterojunction Morphologies for Organic Photovoltaic Device Modeling

abstract: Kinetic Monte Carlo (KMC) simulations have been previously used to model and understand a wide range of behaviors in bulk heterojunction (BHJ) organic photovoltaic devices, from fundamental mechanisms to full device performance. One particularly unique and valuable aspect of this type of modeling technique is the ability to explicitly implement models for the bicontinuous nanostructured morphology present in these devices. For this purpose, an Ising-based method for creating model BHJ morphologies has become prevalent. However, this technique can be computationally expensive, and a detailed characterization of this method has not yet been published. Here, we perform a thorough characterization of this method and describe how to efficiently generate controlled model BHJ morphologies. We show how the interaction energy affects the tortuosity of the interconnected domains and the resulting charge transport behavior in KMC simulations. We also demonstrate how to dramatically reduce calculation time by several orders of magnitude without detrimentally affecting the resulting morphologies. In the end, we propose standard conditions for generating model morphologies and introduce a new open-source software tool. These developments to the Ising method provide a strong foundation for future simulation and modeling of BHJ organic photovoltaic devices that will lead to a more detailed understanding of the important link between morphological features and device performance.


title: Connectivity and Convexity Properties of the Momentum Map for Group Actions on Hilbert Manifolds

abstract: In the early $1980$s a landmark result was obtained by Atiyah and independently Guillemin and Sternberg: the image of the momentum map for a torus action on a compact symplectic manifold is a convex polyhedron. Atiyah's proof makes use of the fact that level sets of the momentum map are connected. These proofs work in the setting of finite-dimensional compact symplectic manifolds. One can ask how these results generalize. A well-known example of an infinite-dimensional symplectic manifold with a finite-dimensional torus action is the based loop group. Atiyah and Pressley proved convexity for this example, but not connectedness of level sets. A proof of connectedness of level sets for the based loop group was provided by Harada, Holm, Jeffrey and Mare in $2006$. In this thesis we study Hilbert manifolds equipped with a strong symplectic structure and a finite-dimensional group action preserving the strong symplectic structure. We prove connectedness of regular generic level sets of the momentum map. We use this to prove convexity of the image of the momentum map.


title: Angle-resolved effective potentials for disk-shaped molecules

abstract: We present an approach for calculating coarse-grained angle-resolved effective pair potentials for uniaxial molecules. For integrating out the intramolecular degrees of freedom we apply umbrella sampling and steered dynamics techniques in atomistically-resolved molecular dynamics (MD) computer simulations. Throughout this study we focus on disk-like molecules such as coronene. To develop the methods we focus on integrating out the van-der-Waals and intramolecular interactions, while electrostatic charge contributions are neglected. The resulting coarse-grained pair potential reveals a strong temperature and angle dependence. In the next step we fit the numerical data with various Gay-Berne-like potentials to be used in more efficient simulations on larger scales. The quality of the resulting coarse-grained results is evaluated by comparing their pair and many-body structure as well as some thermodynamic quantities self-consistently to the outcome of atomistic MD simulations of many-particle systems. We find that angle-resolved potentials are essential not only to accurately describe crystal structures but also for fluid systems where simple isotropic potentials start to fail already for low to moderate packing fractions. Further, in describing these states it is crucial to take into account the pronounced temperature dependence arising in selected pair configurations due to bending fluctuations.


title: Dynamics of multiple bubbles, excited by femtosecond filament in water: Role of aberrations

abstract: Using shadow photography, we observed microsecond time scale evolution of multiple cavitation bubbles, excited by tighty focused femtosecond laser pulse in water under supercritical power regime (~100 Pcr). In these extreme conditions high energy delivery into the microvolume of liquid sample leads to creation of single filament which becomes a source of cavitation region formation. When aberrations were added to the optical scheme the hot spots along the filament axis are formed. At high energies (more than 40uJ) filaments in these hot spots are fired and, as a result, complex pattern of cavitation bubbles is created. The bubbles can be isolated from each other or build exotic drop-shaped cavitation region, which evolution at the end of its life, before the final collapse, contains the jet emission. The dynamics of the cavitation pattern was investigated from pulse energy and focusing. We found that greater numerical aperture of the focusing optics leads to greater cavitation area length. The strong interaction between bubbles was also observed. This leads to a significant change of bubble evolution, which is not yet in accordance with Rayleigh model.


title: The Galactic Center cloud G2 and its gas streamer

abstract: We present new, deep near-infrared SINFONI @ VLT integral field spectroscopy of the gas cloud G2 in the Galactic Center, from late summer 2013 and spring 2014. G2 is visible in recombination line emission. The spatially resolved kinematic data track the ongoing tidal disruption. As expected for an observation near pericenter passage, roughly half of the gas in 2014 is found at the redshifted, pre-pericenter side of the orbit, while the other half is at the post-pericenter, blueshifted side. We also present an orbital solution for the gas cloud G1, which was discovered a decade ago in L'-band images when it was spatially almost coincident with Sgr A*. The orientation of the G1 orbit in the three angles is almost identical to the one of G2, but it has a lower eccentricity and smaller semi-major axis. We show that the observed astrometric positions and radial velocities of G1 are compatible with the G2 orbit, assuming that (i) G1 was originally on the G2 orbit preceding G2 by 13 years and (ii) a simple drag force acted on it during pericenter passage. Taken together with the previously described tail of G2, which we detect in recombination line emission and thermal broadband emission, we propose that G2 may be a bright knot in a much more extensive gas streamer. This matches purely gaseous models for G2, such as a stellar wind clump or the tidal debris from a partial disruption of a star.


title: Pricing for local and global WiFi markets

abstract: This paper analyzes two pricing schemes commonly used in WiFi markets: the flat-rate and the usage-based pricing. The flat-rate pricing encourages the maximum usage, while the usage-based pricing can flexibly attract more users especially those with low valuations in mobile Internet access. First, we use theoretical analysis to compare the two schemes and show that for a single provider in a market, as long as the WiFi capacity is abundant, the flat-rate pricing leads to more revenue. Second, we study how a global provider (e.g., Skype) collaborates with this monopolist in each local market to provide a global WiFi service. We formulate {the interactions between the global and local providers as a dynamic game. In Stage I, the global provider bargains with the local provider in each market to determine the global WiFi service price and revenue sharing agreement. In Stage II, local users and travelers choose local or global WiFi services. We analytically show that the global provider prefers to use the usage-based pricing to avoid a severe competition with the local provider. At the equilibrium, the global provider always shares the majority of his revenue with the local provider to incentivize the cooperation. Finally, we analytically study how the interaction changes if the local market has more than one local provider. In this case, the global provider can integrate the coverages of multiple local providers and provide a better service. Compared to the local monopoly case, local market competition enables the global provider to share less revenue with each of the local providers. However, we numerically show that the global provider's revenue could decrease, as he shares his revenue with more providers and can only charge a lower price.


title: Adiabatic theorem for bipartite quantum systems in weak coupling limit

abstract: We study the adiabatic approximation of the dynamics of a bipartite quantum system with respect to one of the components, when the coupling between its two components is perturbative. We show that the density matrix of the considered component is described by adiabatic transport formulae exhibiting operator-valued geometric and dynamical phases. The present results can be used to study the quantum control of the dynamics of qubits and of open quantum systems where the two components are the system and its environment. We treat two examples, the control of an atomic qubit interacting with another one and the control of a spin in the middle of a Heisenberg spin chain.


title: Constraints on absolute neutrino Majorana mass from current data

abstract: We present new constraints on the neutrino Majorana masses from the current data of neutrinoless double beta decay and neutrino flavour mixing. With the latest results of $0\nu\beta\beta$ progresses from various isotopes, including the recent calculations of the nuclear matrix elements, we find that the strongest constraint of the effective Majorana neutrino mass is from the $^{136}\rm{Xe}$ data of the EXO-200 and KamLAND-Zen collaborations. Further more, by combining the $0\nu\beta\beta$ experimental data with the neutrino mixing parameters from new analyses, we get the mass upper limits of neutrino mass eigenstates and flavour eigenstates and suggest several relations among these neutrino masses.


title: A theory for the zeros of Riemann Zeta and other L-functions

abstract: In these lectures we first review all of the important properties of the Riemann $\zeta$-function, necessary to understand the importance and nature of the Riemann hypothesis. In particular this first part describes the analytic continuation, the functional equation, trivial zeros, the Euler product formula, Riemann's main result relating the zeros on the critical strip to the distribution of primes, the exact counting formula for the number of zeros on the strip $N(T)$ and the GUE statistics of the zeros on the critical line. We then turn to presenting some new results obtained in the past year and describe several strategies towards proving the RH. First we describe an electro-static analogy and argue that if the electric potential along the line $\Re (z) =1$ is a regular alternating function, the RH would follow. The main new result is that the zeros on the critical line are in one-to-one correspondence with the zeros of the cosine function, and this leads to a transcendental equation for the $n$-th zero on the critical line that depends only on $n$. If there is a unique solution to this equation for every $n$, then if $N_0 (T)$ is the number of zeros on the critical line, then $N_0(T) = N(T)$, i.e. all zeros are on the critical line. These results are generalized to two infinite classes of functions, Dirichlet $L$-functions and $L$-functions based on modular forms. We present extensive numerical analysis of the solutions of these equations. We apply these methods to the Davenport-Heilbronn $L$-function, which is known to have zeros off of the line, and explain why the RH fails in this case. We also present a new approximation to the $\zeta$-function that is analogous to the Stirling approximation to the $\Gamma$-function.


title: Reducing Collective Quantum State Rotation Errors with Reversible Dephasing

abstract: We demonstrate that reversible dephasing via inhomogeneous broadening can greatly reduce collective quantum state rotation errors, and observe the suppression of rotation errors by more than 21 dB in the context of collective population measurements of the spin states of an ensemble of $2.1 \times 10^5$ laser cooled and trapped $^{87}$Rb atoms. The large reduction in rotation noise enables direct resolution of spin state populations 13(1) dB below the fundamental quantum projection noise limit. Further, the spin state measurement projects the system into an entangled state with 9.5(5) dB of directly observed spectroscopic enhancement (squeezing) relative to the standard quantum limit, whereas no enhancement would have been obtained without the suppression of rotation errors.


title: Data Assimilation by Artificial Neural Networks for an Atmospheric General Circulation Model: Conventional Observation

abstract: This paper presents an approach for employing artificial neural networks (NN) to emulate an ensemble Kalman filter (EnKF) as a method of data assimilation. The assimilation methods are tested in the Simplified Parameterizations PrimitivE-Equation Dynamics (SPEEDY) model, an atmospheric general circulation model (AGCM), using synthetic observational data simulating localization of balloon soundings. For the data assimilation scheme, the supervised NN, the multilayer perceptrons (MLP-NN), is applied. The MLP-NN are able to emulate the analysis from the local ensemble transform Kalman filter (LETKF). After the training process, the method using the MLP-NN is seen as a function of data assimilation. The NN were trained with data from first three months of 1982, 1983, and 1984. A hind-casting experiment for the 1985 data assimilation cycle using MLP-NN were performed with synthetic observations for January 1985. The numerical results demonstrate the effectiveness of the NN technique for atmospheric data assimilation. The results of the NN analyses are very close to the results from the LETKF analyses, the differences of the monthly average of absolute temperature analyses is of order 0.02. The simulations show that the major advantage of using the MLP-NN is better computational performance, since the analyses have similar quality. The CPU-time cycle assimilation with MLP-NN is 90 times faster than cycle assimilation with LETKF for the numerical experiment.


title: New tools for determining the light travel time in static, spherically symmetric spacetimes beyond the order $G^2$

abstract: This paper is mainly devoted to the determination of the travel time of a photon as a function of the positions of the emitter and the receiver in a large class of static, spherically symmetric spacetimes. Such a function - often called time transfer function - is of crucial interest for testing metric theories of gravity in the solar system. Until very recently, this function was known only up to the second order in the Newtonian gravitational constant $G$ for a 3-parameter family of static, spherically symmetric metrics generalizing the Schwarzschild metric. We present here two procedures enabling to determine - at least in principle - the time transfer function at any order of approximation when the components of the metric are expressible in power series of the Schwarzschild radius of the central body divided by the radial coordinate. These procedures exclusively work for light rays which may be described as perturbations in power series in $G$ of a Minkowskian null geodesic passing through the positions of the emitter and the receiver. It is shown that the two methodologies lead to the same expression for the time transfer function up to the third order in $G$. The second procedure presents the advantage of exclusively needing elementary integrations which may be performed with any symbolic computer program whatever the order of approximation. The vector functions characterizing the direction of light propagation at the points of emission and reception are derived up to the third order in $G$. The relevance of the third order terms in the time transfer function is briefly discussed for some solar system experiments.


title: Unextendible entangled bases with fixed Schmidt number

abstract: The unextendible product basis (UPB) is generalized to the unextendible entangled basis with any arbitrarily given Schmidt number $k$ (UEBk) for any bipartite system $\mathbb{C}^d\otimes\mathbb{C}^{d'}$ ($2\leq k<d\leq d'$), which can also be regarded as a generalization of the unextendible maximally entangled basis (UMEB). A general way of constructing such a basis with arbitrary $d$ and $d'$ is proposed. Consequently, it is shown that there are at least $k-r$ (here $r=d$ mod $k$, or $r=d'$ mod $k$) sets of UEBk when $d$ or $d'$ is not the multiple of $k$, while there are at least $2(k-1)$ sets of UEBk when both $d$ and $d'$ are the multiples of $k$.


title: Laboratory Characterization and Astrophysical Detection of Vibrationally Excited States of Vinyl Cyanide in Orion-KL

abstract: New laboratory data of CH$_2$CHCN (vinyl cyanide) in its ground and vibrationally excited states at the microwave to THz domain allow searching for these excited state transitions in the Orion-KL line survey. Frequency-modulated spectrometers combined into a single broadband 50-1900 GHz spectrum provided measurements of CH$_2$CHCN covering a spectral range of 18-1893 GHz, whose assignments was confirmed by Stark modulation spectra in the 18-40 GHz region and by ab-initio anharmonic force field calculations. For analyzing the emission lines of CH$_2$CHCN species detected in Orion-KL we used the excitation and radiative transfer code (MADEX) at LTE conditions. The rotational transitions of the ground state of this molecule emerge from four cloud components of hot core nature which trace the physical and chemical conditions of high mass star forming regions in the Orion-KL Nebula. The total column density of CH$_2$CHCN in the ground state is (3.0$\pm$0.9)x10$^{15}$ cm$^{-2}$. We report on the first interstellar detection of transitions in the v10=1/(v11=1,v15=1) dyad in space, and in the v11=2 and v11=3 states in Orion-KL. The lowest energy vibrationally excited states of vinyl cyanide such as v11=1 (at 328.5 K), v15=1 (at 478.6 K), v11=2 (at 657.8 K), the v10=1/(v11=1,v15=1) dyad (at 806.4/809.9 K), and v11=3 (at 987.9 K) are populated under warm and dense conditions, so they probe the hottest parts of the Orion-KL source. Column density and rotational and vibrational temperatures for CH$_2$CHCN in their ground and excited states, as well as for the isotopologues, have been constrained by means of a sample of more than 1000 lines in this survey. Moreover, we present the detection of methyl isocyanide (CH$_3$NC) for the first time in Orion-KL and a tentative detection of vinyl isocyanide (CH$_2$CHNC) and give column density ratios between the cyanide and isocyanide isomers.


title: One-Step or Two-Step Optimization and the Overfitting Phenomenon: A Case Study on Time Series Classification

abstract: For the last few decades, optimization has been developing at a fast rate. Bio-inspired optimization algorithms are metaheuristics inspired by nature. These algorithms have been applied to solve different problems in engineering, economics, and other domains. Bio-inspired algorithms have also been applied in different branches of information technology such as networking and software engineering. Time series data mining is a field of information technology that has its share of these applications too. In previous works we showed how bio-inspired algorithms such as the genetic algorithms and differential evolution can be used to find the locations of the breakpoints used in the symbolic aggregate approximation of time series representation, and in another work we showed how we can utilize the particle swarm optimization, one of the famous bio-inspired algorithms, to set weights to the different segments in the symbolic aggregate approximation representation. In this paper we present, in two different approaches, a new meta optimization process that produces optimal locations of the breakpoints in addition to optimal weights of the segments. The experiments of time series classification task that we conducted show an interesting example of how the overfitting phenomenon, a frequently encountered problem in data mining which happens when the model overfits the training set, can interfere in the optimization process and hide the superior performance of an optimization algorithm.


title: The molecular structure of the interface between water and a hydrophobic substrate is liquid-vapor like

abstract: With molecular simulation for water and a tunable hydrophobic substrate, we apply the instantaneous interface construction [A. P. Willard and D. Chandler, J. Phys. Chem. B, 114, 1954 (2010)] to examine the similarity between a water-vapor interface and a water-hydrophobic surface interface. The intrinsic interface refers to molecular structure in terms of distances from the instantaneous interface. We show that attractive interactions between a hydrophobic surface and water affect capillary wave fluctuations of the instantaneous liquid interface, but these attractive interactions have essentially no effect on the intrinsic interface. Further, the intrinsic interface of liquid water and a hydrophobic substrate differs little from that of water and its vapor.The same is not true, we show, for an interface between water and a hydrophilic substrate. In that case, strong directional substrate-water interactions disrupt the liquid-vapor-like interfacial hydrogen bonding network.


title: Geometry of flexible filament cohesion: Better contact through twist?

abstract: Cohesive interactions between filamentous molecules have broad implications for a range of biological and synthetic materials. While long-standing theoretical approaches have addressed the problem of inter-filament forces from the limit of infinitely rigid rods, the ability of flexible filaments to deform intra-filament shape in response to changes in inter-filament geometry has a profound affect on the nature of cohesive interactions. In this paper, we study two theoretical models of inter-filament cohesion in the opposite limit, in which filaments are sufficiently flexible to maintain cohesive contact along their contours, and address, in particular, the role played by helical-interfilament geometry in defining interactions. Specifically, we study models of featureless, tubular filaments interacting via 1) pair-wise Lennard-Jones (LJ) interactions between surface elements and 2) depletion-induced filament binding stabilized by electrostatic surface repulsion. Analysis of these models reveals a universal preference for cohesive filament interactions for non-zero helical skew, and further, that in the asymptotic limit of vanishing interaction range relative to filament diameter, the skew-dependence of cohesion approaches a geometrically defined limit described purely by the close-packing geometry of twisted tubular filaments. We further analyze non-universal features of the skew-dependence of cohesion at small-twist for both potentials, and argue that in the LJ model the pair-wise surface attraction generically destabilizes parallel filaments, while in the second model, pair-wise electrostatic repulsion in combination with non-pairwise additivity of depletion leads to a meta-stable parallel state.


title: Jayant Vishnu Narlikar

abstract: An essay on Jyant Narlikar in the Living Legends of Indian Science series of the journal Current Science.

categories: physics.hist-ph astro-ph.CO gr-qc hep-th


title: Differential games with asymmetric information and without Isaacs condition

abstract: We investigate a two-player zero-sum differential game with asymmetric information on the payoff and without Isaacs condition. The dynamics is an ordinary differential equation parametrised by two controls chosen by the players. Each player has a private information on the payoff of the game, while his opponent knows only the probability distribution on the information of the other player. We show that a suitable definition of random strategies allows to prove the existence of a value in mixed strategies. Moreover, the value function can be characterised in term of the unique viscosity solution in some dual sense of a Hamilton-Jacobi-Isaacs equation. Here we do not suppose the Isaacs condition which is usually assumed in differential games.


title: The Cauchy problem on large time for the Water Waves equations with large topography variations

abstract: We prove the local existence for the Water Waves equations with large bathymetric variations on a time interval of size 1/\epsilon, where $\epsilon$ measures the amplitude of the wave. We just need the presence of surface tension.


title: The Schr\"odinger-Newton equation and its foundations

abstract: The necessity of quantising the gravitational field is still subject to an open debate. In this paper we compare the approach of quantum gravity, with that of a fundamentally semi-classical theory of gravity, in the weak-field non-relativistic limit. We show that, while in the former case the Schr\"odinger equation stays linear, in the latter case one ends up with the so-called Schr\"odinger-Newton equation, which involves a nonlinear, non-local gravitational contribution. We further discuss that the Schr\"odinger-Newton equation does not describe the collapse of the wave-function, although it was initially proposed for exactly this purpose. Together with the standard collapse postulate, fundamentally semi-classical gravity gives rise to superluminal signalling. A consistent fundamentally semi-classical theory of gravity can therefore only be achieved together with a suitable prescription of the wave-function collapse. We further discuss, how collapse models avoid such superluminal signalling and compare the nonlinearities appearing in these models with those in the Schr\"odinger-Newton equation.


title: Iwasawa theory for symmetric powers of CM modular forms at nonordinary primes, II

abstract: Continuing the study of the Iwasawa theory of symmetric powers of CM modular forms at supersingular primes begun by the first author and Antonio Lei, we prove a Main Conjecture equating the "admissible" $p$-adic $L$-functions to characteristic ideals of "finite-slope" Selmer modules constructed by the second author. As a key ingredient, we improve Rubin's result on the Main Conjecture of Iwasawa theory for imaginary quadratic fields to an equality at inert primes.


title: Stein estimation of the intensity of a spatial homogeneous Poisson point process

abstract: In this paper, we revisit the original ideas of Stein and propose an estimator of the intensity parameter of a homogeneous Poisson point process defined in $\R^d$ and observed in a bounded window. The procedure is based on a new general integration by parts formula for Poisson point processes. We show that our Stein estimator outperforms the maximum likelihood estimator in terms of mean squared error. In particular, we show that in many practical situations we have a gain larger than 30\%.


title: Online Asynchronous Distributed Regression

abstract: Distributed computing offers a high degree of flexibility to accommodate modern learning constraints and the ever increasing size of datasets involved in massive data issues. Drawing inspiration from the theory of distributed computation models developed in the context of gradient-type optimization algorithms, we present a consensus-based asynchronous distributed approach for nonparametric online regression and analyze some of its asymptotic properties. Substantial numerical evidence involving up to 28 parallel processors is provided on synthetic datasets to assess the excellent performance of our method, both in terms of computation time and prediction accuracy.


title: Therapeutic target discovery using Boolean network attractors: avoiding pathological phenotypes

abstract: Target identification, one of the steps of drug discovery, aims at identifying biomolecules whose function should be therapeutically altered in order to cure the considered pathology. This work proposes an algorithm for in silico target identification using Boolean network attractors. It assumes that attractors of dynamical systems, such as Boolean networks, correspond to phenotypes produced by the modeled biological system. Under this assumption, and given a Boolean network modeling a pathophysiology, the algorithm identifies target combinations able to remove attractors associated with pathological phenotypes. It is tested on a Boolean model of the mammalian cell cycle bearing a constitutive inactivation of the retinoblastoma protein, as seen in cancers, and its applications are illustrated on a Boolean model of Fanconi anemia. The results show that the algorithm returns target combinations able to remove attractors associated with pathological phenotypes and then succeeds in performing the proposed in silico target identification. However, as with any in silico evidence, there is a bridge to cross between theory and practice, thus requiring it to be used in combination with wet lab experiments. Nevertheless, it is expected that the algorithm is of interest for target identification, notably by exploiting the inexpensiveness and predictive power of computational approaches to optimize the efficiency of costly wet lab experiments.


title: Molecular dynamics simulation of growth of Cu nanoclusters from Cu-ions in a plasma

abstract: A recently developed method of nanoclusters growth in a pulsed plasma is studied by means of molecular dynamics. A model that allows one to consider high-energy charged particles in classical molecular dynamics is suggested, and applied for studies of single impact events in nanoclusters growth. In particular, we provide a comparative analysis of the well-studied inert gas aggregation method and the growth from ions in a plasma. The importance to consider of the angular distribution of incoming ions in the simulations of the nanocluster growth is underlined. A detailed study of the energy transfer from the incoming ions to a nanocluster, as well as the diffusion of incoming ions on the cluster surface is carried out. Our results are important for understanding and control of the nanocluster growth process.


title: Common price and volatility jumps in noisy high-frequency data

abstract: We introduce a statistical test for simultaneous jumps in the price of a financial asset and its volatility process. The proposed test is based on high-frequency data and is robust to market microstructure frictions. For the test, local estimators of volatility jumps at price jump arrival times are designed using a nonparametric spectral estimator of the spot volatility process. A simulation study and an empirical example with NASDAQ order book data demonstrate the practicability of the proposed methods and highlight the important role played by price volatility co-jumps.


title: Recovery-Based Error Estimators for Diffusion Problems: Explicit Formulas

abstract: We introduced and analyzed robust recovery-based a posteriori error estimators for various lower order finite element approximations to interface problems in [9, 10], where the recoveries of the flux and/or gradient are implicit (i.e., requiring solutions of global problems with mass matrices). In this paper, we develop fully explicit recovery-based error estimators for lower order conforming, mixed, and non- conforming finite element approximations to diffusion problems with full coefficient tensor. When the diffusion coefficient is piecewise constant scalar and its distribution is local quasi-monotone, it is shown theoretically that the estimators developed in this paper are robust with respect to the size of jumps. Numerical experiments are also performed to support the theoretical results.


title: PaPy: Parallel and Distributed Data-processing Pipelines in Python

abstract: PaPy, which stands for parallel pipelines in Python, is a highly flexible framework that enables the construction of robust, scalable workflows for either generating or processing voluminous datasets. A workflow is created from user-written Python functions (nodes) connected by 'pipes' (edges) into a directed acyclic graph. These functions are arbitrarily definable, and can make use of any Python modules or external binaries. Given a user-defined topology and collection of input data, functions are composed into nested higher-order maps, which are transparently and robustly evaluated in parallel on a single computer or on remote hosts. Local and remote computational resources can be flexibly pooled and assigned to functional nodes, thereby allowing facile load-balancing and pipeline optimization to maximize computational throughput. Input items are processed by nodes in parallel, and traverse the graph in batches of adjustable size -- a trade-off between lazy-evaluation, parallelism, and memory consumption. The processing of a single item can be parallelized in a scatter/gather scheme. The simplicity and flexibility of distributed workflows using PaPy bridges the gap between desktop -> grid, enabling this new computing paradigm to be leveraged in the processing of large scientific datasets.


title: Intrinsic and extrinsic decay of edge magnetoplasmons in graphene

abstract: We investigate intrinsic and extrinsic decay of edge magnetoplasmons (EMPs) in graphene quantum Hall (QH) systems by high-frequency electronic measurements. From EMP resonances in disk shaped graphene, we show that the dispersion relation of EMPs is nonlinear due to interactions, giving rise to intrinsic decay of EMP wavepacket. We also identify extrinsic dissipation mechanisms due to interaction with localized states in bulk graphene from the decay time of EMP wavepackets. We indicate that, owing to the unique linear and gapless band structure, EMP dissipation in graphene can be lower than that in GaAs systems.


title: Hierarchy of neural organization in the embryonic spinal cord: Granger-causality graph analysis of in vivo calcium imaging data

abstract: The recent development of genetically encoded calcium indicators enables monitoring in vivo the activity of neuronal populations. Most analysis of these calcium transients relies on linear regression analysis based on the sensory stimulus applied or the behavior observed. To estimate the basic properties of the functional neural circuitry, we propose a network-based approach based on calcium imaging recorded at single cell resolution. Differently from previous analysis based on cross-correlation, we used Granger-causality estimates to infer activity propagation between the activities of different neurons. The resulting functional networks were then modeled as directed graphs and characterized in terms of connectivity and node centralities. We applied our approach to calcium transients recorded at low frequency (4 Hz) in ventral neurons of the zebrafish spinal cord at the embryonic stage when spontaneous coiling of the tail occurs. Our analysis on population calcium imaging data revealed a strong ipsilateral connectivity and a characteristic hierarchical organization of the network hubs that supported established propagation of activity from rostral to caudal spinal cord. Our method could be used for detecting functional defects in neuronal circuitry during development and pathological conditions.


title: Rippling transition from electron-induced condensation of curvature field in graphene

abstract: A quantum field theory approach is applied to investigate the dynamics of flexural phonons in a metallic membrane like graphene, looking for the effects deriving from the strong interaction between the electronic excitations and elastic deformations. Relying on a self-consistent screening approximation to the phonon self-energy, we show that the theory has a critical point characterized by the vanishing of the effective bending rigidity of the membrane at low momentum. We also check that the instability in the sector of flexural phonons takes place without the development of an in-plane static distortion, which is avoided due to the significant reduction of the electron-phonon couplings for in-plane phonons at large momenta. Furthermore, we analyze the scaling properties of the many-body theory to identify the order parameter that opens up at the point of the transition. We find that the vanishing of the effective bending rigidity and the onset of a nonzero expectation value of the mean curvature are concurrent manifestations of the critical behavior. The results presented here imply that, even in the absence of tension, the theory has a critical point at which the flat geometry becomes an unstable configuration of the metallic membrane, with a condensation of the mean curvature field that may well reproduce the smooth distribution of ripples observed in free-standing graphene.


title: Classical bulges, supermassive blackholes and AGN feedback: Extension to low-mass galaxies

abstract: The empirical model of Lu et al. 2014a for the relation between star formation rate and halo mass growth is adopted to predict the classical bulge mass ($M_{\rm cb}$) - total stellar mass ($M_\star$) relation for central galaxies. The assumption that the supermassive black hole (SMBH) mass ($M_{\rm BH}$) is directly proportional to the classical bulge mass, with the proportionality given by that for massive galaxies, predicts a $M_{\rm BH}$ - $M_\star$ relation that matches well the observed relation for different types of galaxies. In particular, the model reproduces the strong transition at $M_\star=10^{10.5}$ - $10^{11}M_{\odot}$, below which $M_{\rm BH}$ drops rapidly with decreasing $M_\star$. Our model predicts a new sequence at $M_\star <10^{10.5}M_{\odot}$, where $M_{\rm BH} \propto M_\star$ but the amplitude is a factor of $\sim 50$ lower than the amplitude of the sequence at $M_\star>10^{11}M_{\odot}$. If all SMBH grow through similar quasar modes with a feedback efficiency of a few percent, the energy produced in low-mass galaxies at redshift $z\gtrsim 2$ can heat the circum-galactic medium up to a specific entropy level that is required to prevent excessive star formation in low-mass dark matter halos.


title: Reissner-Nordstrom Black Holes in the Inverse Electrodynamics Model

abstract: We study electric and magnetic monopoles in static, spherically symmetric and constant curvature geometries in the context of the inverse electrodynamics model. We prove that this U(1) invariant Lagrangian density is able to support the standard metric of a Reissner-Nordstrom Black Hole, but with more complex thermodynamical properties than in the standard case. By employing the Euclidean Action approach we perform a complete analysis of its phase space depending on the sign and singularities of the heat capacity and the Helmholtz free energy.


title: Improved Estimates for the Parameters of the Heavy Quark Expansion

abstract: We give improved estimates for the non-perturbative parameters appearing in the heavy quark expansion for inclusive decays. While the parameters appearing in low orders of this expansion can be extracted from data, the number of parameters in higher orders proliferates strongly, making a determination of these parameters from data impossible. Thus, one has to rely on theoretical estimates which may be obtained from an insertion of intermediate states. In this paper we refine this method and attempt to estimate the uncertainties of this approach.


title: LOFAR tied-array imaging of Type III solar radio bursts

abstract: The Sun is an active source of radio emission which is often associated with energetic phenomena such as solar flares and coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), the Sun has not been imaged extensively because of the instrumental limitations of previous radio telescopes. Here, the combined high spatial, spectral and temporal resolution of the Low Frequency Array (LOFAR) was used to study solar Type III radio bursts at 30-90 MHz and their association with CMEs. The Sun was imaged with 126 simultaneous tied-array beams within 5 solar radii of the solar centre. This method offers benefits over standard interferometric imaging since each beam produces high temporal (83 ms) and spectral resolution (12.5 kHz) dynamic spectra at an array of spatial locations centred on the Sun. LOFAR's standard interferometric output is currently limited to one image per second. Over a period of 30 minutes, multiple Type III radio bursts were observed, a number of which were found to be located at high altitudes (4 solar radii from the solar center at 30 MHz) and to have non-radial trajectories. These bursts occurred at altitudes in excess of values predicted by 1D radial electron density models. The non-radial high altitude Type III bursts were found to be associated with the expanding flank of a CME. The CME may have compressed neighbouring streamer plasma producing larger electron densities at high altitudes, while the non-radial burst trajectories can be explained by the deflection of radial magnetic fields as the CME expanded in the low corona.


title: Bipartite separability of one parameter families of states using conditional quantum relative Tsallis entropy

abstract: In any bipartition of a quantum state, it is proved that the negative values of the conditional version of sandwiched Tsallis relative entropy necessarily implies quantum entanglement. For any N, the separability ranges in the $1:N-1$ partition of symmetric one parameter families of noisy $N$-qubit W, GHZ, states are determined using the conditional quantum relative Tsallis entropy approach. The 1:N-1 separability range matches exactly with the range obtained through positive partial transpose criterion, for all N. The advantages of using non-commuting version of $q$-conditional relative Tsallis entropy is brought out through this and other one-parameter families of states.


title: Tunable Light Filtering in a Bragg Mirror/Heavily-Doped Semiconducting Nanocrystal Composite

abstract: Tunable light filters are critical components for many optical applications in which light in-coupling, out-coupling or rejection is crucial, such as lasing, sensing, photovoltaics and information and communication technology. For this purpose, Bragg mirrors, band-pass filters with high reflectivity represent good candidates. However, their optical characteristics are determined at the stage of fabrication. Heavily doped semiconductor nanocrystals (NCs) on the other hand deliver a high degree of optical tunability through the active modulation of their carrier density ultimately influencing their plasmonic absorption properties. Here, we propose the design of a tunable light filter composed of a Bragg mirror and a layer of plasmonic semiconductor NCs. We demonstrate that the filtering properties of the coupled device can be tuned to cover a wide range of frequencies from the visible to the near infrared (vis-NIR) spectral region when employing varying carrier densities. As tunable component we implemented a dispersion of copper selenide (Cu2-xSe) NCs and a film of indium tin oxide (ITO) NCs, which are known to show optical tunablility by chemical or electrochemical treatments. We utilized the Mie theory to describe the carrier dependent plasmonic properties of the Cu2-xSe NC dispersion and the effective medium theory to describe the optical characteristics of the ITO film. The transmission properties of the Bragg mirror have been modelled with the transfer matrix method. We foresee ease of experimental realization of the coupled device, where filtering modulation is achieved upon chemical and electrochemical post-fabrication treatment of the heavily doped semiconductor NC component, eventually resulting in tunable transmission properties of the coupled device.


title: Delivery of focused short pulses through a multimode fiber

abstract: Light propagation through multimode fibers suffers from spatial distortions that lead to a scrambled intensity profile. In previous work, the correction of such distortions using various wavefront control methods has been demonstrated in the continuous wave case. However, in the ultra-fast pulse regime, modal dispersion temporally broadens a pulse after propagation. Here, we present a method that compensates for spatial distortions and mitigates temporal broadening due to modal dispersion by a selective phase conjugation process in which only modes of similar group velocities are excited. The selectively excited modes are forced to follow certain paths through the multimode fiber and interfere constructively at the distal tip to form a focused spot with minimal temporal broadening. We demonstrate the delivery of focused 500 fs pulses through a 30 cm long step-index multimode fiber. The achieved pulse duration corresponds to approximately 1/30th of the duration obtained if modal dispersion was not controlled. Moreover, we measured a detailed two-dimensional map of the pulse duration at the output of the fiber and confirmed that the focused spot produces a two-photon absorption effect. This work opens new possibilities for ultra-thin multiphoton imaging through multimode fibers.


title: Electronic and magnetic properties of the graphene/Eu/Ni(111) hybrid system

abstract: The electronic and magnetic properties of the graphene/Eu/Ni(111) intercalation-like system are studied in the framework of the GGA+U approach with dispersive interactions taken into account. Intercalation of monoatomic Eu layer underneath graphene on Ni(111) leads to the drastic changes of the electronic structure of graphene compared to free-standing graphene as well as graphene/Ni(111). The strong influence of the spin-polarised Eu 4f states, crossing the graphene-derived $\pi$ states, on magnetic properties of graphene and on spin-filtering properties of the graphene/Eu/Ni(111) trilayer is discussed.


title: A preliminary threshold model of parasitism in the Cockle\emph{Cerastoderma edule} using delayed exchange of stability

abstract: Thresholds occur in the dynamics of many biological communities. Here we model a persistence type threshold which has been shown experimentally to exist in hyperparasitised flukes in the cockle, a shellfish. Our model consists of a periodically driven slow-fast host-parasite system of equations for a slow flukes population and a fast Unikaryon hyperparasite population. The model exhibits two branches of the critical curve crossing in a transcritical bifurcation scenario. We discuss two thresholds due to immediate and delayed exchange of stability effects; and we derive algebraic relationships for parameters of the periodic solution in the limit of the infinite ratio of the time scales. Flukes parasitise cockles and in turn are hyperparasitised by the microsporidian Unikaryon legeri; the life cycle of flukes includes several life stages and a number of different hosts. That is, the flukes-hyperparasite system in a cockle is part of a larger estuarine ecosystem of species involving parasites, shellfish and birds which prey on shellfish. A population dynamics model which accounts for such multi-species interactions and includes the fluke-hyperparasite model in a cockle as a subsystem is presented. We provide evidence that the threshold effect we observed in the flukes-hyperparasite subsystem remains apparent in the multi-species system. Assuming that flukes damage cockles, and taking into account that the hyperparasite is detrimental to flukes, it is natural to suggest that the hyperparasitism may support the abundance of cockles and, thereby, the persistence of the ecosystem, including shellfish and birds. We confirm the possibility of this scenario in our model by removing the hyperparasite and demonstrating that this may result in a substantial drop in cockle numbers. The result indicates a possible significant role for the microparasite in this estuarine ecosystem.


title: Panchromatic Hubble Andromeda Treasury IX: A Photometric Survey of Planetary Nebulae in M31

abstract: We search Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) and Wide Field Camera 3 (WFC3) broadband imaging data from the Panchromatic Hubble Andromeda Treasury (PHAT) survey to identify detections of cataloged planetary nebulae (PNe). Of the 711 PNe currently in the literature within the PHAT footprint, we find 467 detected in the broadband. For these 467 we are able to refine their astrometric accuracy from ~0."3 to 0."05. Using the resolution of HST, we are able to show that 152 objects currently in the catalogs are definitively not PNe, and we show that 32 objects thought to be extended in ground-based images are actually point-like and therefore good PN candidates. We also find one PN candidate that is marginally resolved. If this is a PN, it is up to 0.7 pc in diameter. With our new photometric data, we develop a method of measuring the level of excitation in individual PNe by comparing broadband and narrowband imaging and describe the effects of excitation on a PN's photometric signature. Using the photometric properties of the known PNe in the PHAT catalogs, we search for more PN, but do not find any new candidates, suggesting that ground-based emission-line surveys are complete in the PHAT footprint to F475W $\simeq$ 24.


title: Exponential Convergence Towards Stationary States for the 1D Porous Medium Equation with Fractional Pressure

abstract: We analyse the asymptotic behaviour of solutions to the one dimensional fractional version of the porous medium equation introduced by Caffarelli and V\'azquez, where the pressure is obtained as a Riesz potential associated to the density. We take advantage of the displacement convexity of the Riesz potential in one dimension to show a functional inequality involving the entropy, entropy dissipation, and the Euclidean transport distance. An argument by approximation shows that this functional inequality is enough to deduce the exponential convergence of solutions in self-similar variables to the unique steady states.


title: EXCEDE Technology Development III: First Vacuum Tests

abstract: This paper is the third in the series on the technology development for the EXCEDE (EXoplanetary Circumstellar Environments and Disk Explorer) mission concept, which in 2011 was selected by NASA's Explorer program for technology development (Category III). EXCEDE is a 0.7m space telescope concept designed to achieve raw contrasts of 1e6 at an inner working angle of 1.2 l/D and 1e7 at 2 l/D and beyond. This will allow it to directly detect and spatially resolve low surface brightness circumstellar debris disks as well as image giant planets as close as in the habitable zones of their host stars. In addition to doing fundamental science on debris disks, EXCEDE will also serve as a technological and scientific precursor for any future exo-Earth imaging mission. EXCEDE uses a Starlight Suppression System (SSS) based on the PIAA coronagraph, enabling aggressive performance. We report on our continuing progress of developing the SSS for EXCEDE, and in particular (a) the reconfiguration of our system into a more flight-like layout, with an upstream deformable mirror and an inverse PIAA system, as well as a LOWFS, and (b) testing this system in a vacuum chamber, including IWA, contrast, and stability performance. The results achieved so far are 2.9e-7 contrast between 1.2-2.0 l/D and 9.7e-8 contrast between 2.0-6.0 l/D in monochromatic light; as well as 1.4e-6 between 2.0-6.0 l/D in a 10% band, all with a PIAA coronagraph operating at an inner working angle of 1.2 l/D. This constitutes better contrast than EXCEDE requirements (in those regions) in monochromatic light, and progress towards requirements in broadband light. Even though this technology development is primarily targeted towards EXCEDE, it is also germane to any exoplanet direct imaging space-based telescopes because of the many challenges common to different coronagraph architectures and mission requirements.


title: Parameterized Verification of Graph Transformation Systems with Whole Neighbourhood Operations

abstract: We introduce a new class of graph transformation systems in which rewrite rules can be guarded by universally quantified conditions on the neighbourhood of nodes. These conditions are defined via special graph patterns which may be transformed by the rule as well. For the new class for graph rewrite rules, we provide a symbolic procedure working on minimal representations of upward closed sets of configurations. We prove correctness and effectiveness of the procedure by a categorical presentation of rewrite rules as well as the involved order, and using results for well-structured transition systems. We apply the resulting procedure to the analysis of the Distributed Dining Philosophers protocol on an arbitrary network structure.


title: PresenceSense: Zero-training Algorithm for Individual Presence Detection based on Power Monitoring

abstract: Non-intrusive presence detection of individuals in commercial buildings is much easier to implement than intrusive methods such as passive infrared, acoustic sensors, and camera. Individual power consumption, while providing useful feedback and motivation for energy saving, can be used as a valuable source for presence detection. We conduct pilot experiments in an office setting to collect individual presence data by ultrasonic sensors, acceleration sensors, and WiFi access points, in addition to the individual power monitoring data. PresenceSense (PS), a semi-supervised learning algorithm based on power measurement that trains itself with only unlabeled data, is proposed, analyzed and evaluated in the study. Without any labeling efforts, which are usually tedious and time consuming, PresenceSense outperforms popular models whose parameters are optimized over a large training set. The results are interpreted and potential applications of PresenceSense on other data sources are discussed. The significance of this study attaches to space security, occupancy behavior modeling, and energy saving of plug loads.


title: The CP-conserving 2HDM after the 8 TeV run

abstract: We confront the most common CP-conserving 2HDM with the LHC data analysed so far while taking into account all previously available experimental data. A special allowed corner of the parameter space is analysed - the so-called wrong-sign scenario where the Higgs coupling to down-type quarks changes sign relative to the Standard Model while the coupling to the massive vector bosons does not.


title: Influence of Specific Surface Area of Powder on Hydrogen Desorption Kinetics for Metal Hydrides

abstract: The observable results for desorption kinetics by powder of metal hydride on the example of mangesium hydride are reproduced with the model formulated in terms of specific surface of powder. A volumetric measurement of hydrogen desorption process is evaluated on an example of wet ball milled magnesium hydride, and can be applied generally for any metal hydride. The exact solution of the model reproduces the shape of experimental curves for desorption process providing a satisfying agreement with experimental data.


title: Constructive Geometry and the Parallel Postulate

abstract: Euclidean geometry consists of straightedge-and-compass constructions and reasoning about the results of those constructions. We show that Euclidean geometry can be developed using only intuitionistic logic. We consider three versions of Euclid's parallel postulate: Euclid's own formulation in his Postulate 5; Playfair's 1795 version, and a new version we call the strong parallel postulate. These differ in that Euclid's version and the new version both assert the existence of a point where two lines meet, while Playfair's version makes no existence assertion. Classically, the models of Euclidean (straightedge-and-compass) geometry are planes over Euclidean fields. We prove a similar theorem for constructive Euclidean geometry, by showing how to define addition and multiplication without a case distinction about the sign of the arguments. With intuitionistic logic, there are two possible definitions of Euclidean fields, which turn out to correspond to the different versions of the parallel axiom. In this paper, we completely settle the questions about implications between the three versions of the parallel postulate: the strong parallel postulate easily implies Euclid 5, and in fact Euclid 5 also implies the strong parallel postulate, although the proof is lengthy, depending on the verification that Euclid 5 suffices to define multiplication geometrically. We show that Playfair does not imply Euclid 5, and we also give some other independence results. Our independence proofs are given without discussing the exact choice of the other axioms of geometry; all we need is that one can interpret the geometric axioms in Euclidean field theory. The proofs use Kripke models of Euclidean field theories based on carefully constructed rings of real-valued functions.


title: A Constructive Version of Tarski's Geometry

abstract: Euclid's reasoning is essentially constructive. Tarski's elegant and concise first-order theory of Euclidean geometry, on the other hand, is essentially non-constructive, even if we restrict attention (as we do here) to the theory with line-circle and circle-circle continuity in place of first-order Dedekind completeness. Here we exhibit three constructive versions of Tarski's theory. One, like Tarski's theory, has existential axioms and no function symbols. We then consider a version in which function symbols are used instead of existential quantifiers. The third version has a function symbol for the intersection point of two non-parallel, non-coincident lines, instead of only for intersection points produced by Pasch's axiom and the parallel axiom; this choice of function symbols connects directly to ruler-and-compass constructions. All three versions have this in common: the axioms have been modified so that the points they assert to exist are unique and depend continuously on parameters. This modification of Tarski's axioms, with classical logic, has the same theorems as Tarski's theory, but we obtain results connecting it with ruler-and-compass constructions as well. In particular, points constructively proved to exist can be constructed with ruler and compass, uniformly in parameters; the same is true with non-constructive proofs if several constructions are allowed for different cases.


title: Star formation relations and CO SLEDs across the J-ladder and redshift

abstract: We present FIR-CO luminosity relations ($\log L_{\rm FIR} = \alpha \log L'_{\rm CO} + \beta$) for the full CO rotational ladder from J=1-0 to J=13-12 for 62 local (z < 0.1) (Ultra) Luminous Infrared Galaxies (LIRGs) using data from Herschel SPIRE-FTS and ground-based telescopes. We extend our sample to high redshifts (z > 1) by including 35 (sub)-millimeter selected dusty star forming galaxies from the literature with robust CO observations. The addition of luminous starbursts at high redshifts enlarge the range of the FIR-CO luminosity relations towards the high-IR-luminosity end while also significantly increasing the small amount of mid-/high-J CO line data available prior to Herschel. This new data-set (both in terms of IR luminosity and J-ladder) reveals linear FIR-CO luminosity relations ($\alpha \sim 1$) for J=1-0 up to J=5-4, with a nearly constant normalisation ($\beta \sim 2$). This is expected from the (also) linear FIR-(molecular line) relations found for the dense gas tracer lines (HCN and CS), as long as the dense gas mass fraction does not vary strongly within our (merger/starburst)-dominated sample. However from J=6-5 and up to J=13-12 we find an increasingly sub-linear slope and higher normalization constant with increasing J. We argue that these are caused by a warm (~100K) and dense ($>10^4{\rm cm^{-3}}$) gas component whose thermal state is unlikely to be maintained by star formation powered far-UV radiation fields (and thus is no longer directly tied to the star formation rate). We suggest that mechanical heating (e.g., supernova driven turbulence and shocks), and not cosmic rays, is the more likely source of energy for this component. The global CO spectral line energy distributions (SLEDs), which remain highly excited from J=6-5 up to J=13-12, are found to be a generic feature of the (U)LIRGs in our sample, and further support the presence of this gas component.


title: Corroborating the equivalence between the Duffin-Kemmer-Petiau and the Klein-Gordon and Proca equations

abstract: It is shown that the Hamiltonian version of the Duffin-Kemmer-Petiau theory with electromagnetic coupling brings about a source term at the current. It is also shown that such a source term disappears from the scenario if one uses the correct physical form for the Duffin-Kemmer-Petiau field, regardless the choice for representing the Duffin-Kemmer-Petiau matrices. This result is used to fix the ambiguity in the electromagnetic coupling in the Duffin-Kemmer-Petiau theory. Moreover, some widespread misconceptions about the Hermiticity in the Duffin-Kemmer-Petiau theory are discussed.


title: New classes of systematic effects in gas spin co-magnetometers

abstract: Atomic co-magnetometers are widely used in precision measurements searching for spin interactions beyond the Standard Model. We describe a new $^3$He-$^{129}$Xe co-magnetometer probed by Rb atoms and use it to identify two general classes of systematic effects in gas co-magnetometers, one associated with diffusion in second-order magnetic field gradients and another due to temperature gradients. We also develop a general and practical approach for calculating spin relaxation and frequency shifts due to arbitrary magnetic field gradients and confirm it experimentally.


title: Lie groups as 3-dimensional almost contact B-metric manifolds

abstract: The object of investigations are almost contact B-metric structures on 3-dimensional Lie groups considered as smooth manifolds. There are established the existence and some geometric characteristics of these manifolds in all basic classes. An example is given as a support of obtained results.


title: The Chromospheric Solar Millimeter-wave Cavity; a Common Property in the Semi-empirical Models

abstract: The semi-empirical models of the solar chromosphere are useful in the study of the solar radio emission at millimeter - infrared wavelengths. However, current models do not reproduce the observations of the quiet sun. In this work we present a theoretical study of the radiative transfer equation for four semi- empirical models at these wavelengths. We found that the Chromospheric Solar Milimeter-wave Cavity (CSMC), a region where the atmosphere becomes locally optically thin at millimeter wavelengths, is present in the semi-empirical models under study. We conclude that the CSMC is a general property of the solar chromosphere where the semi-empirical models shows temperature minimum.


title: Gravitational wave recoils in non-axisymmetric Robinson-Trautman spacetimes

abstract: We examine the gravitational wave recoil waves and the associated net kick velocities in non-axisymmetric Robinson-Trautman spacetimes. We use characteristic initial data for the dynamics corresponding to non-head-on collisions of black holes. We make a parameter study of the kick distributions, corresponding to an extended range of the incidence angle $\rho_0$ in the initial data. For the range of $\rho_0$ examined ($3^{\circ} \leq \rho_0 \leq 110^{\circ}$) the kick distributions as a function of the symmetric mass parameter $\eta$ satisfy a law obtained from an empirical modification of the Fitchett law, with a parameter $C$ that accounts for the non-zero net gravitational momentum wave fluxes for the equal mass case. The law fits accurately the kick distributions for the range of $\rho_0$ examined, with a rms normalized error of the order of $5 \%$. For the equal mass case the nonzero net gravitational wave momentum flux increases as $\rho_0$ increases, up to $\rho_0 \simeq 55^{\circ}$ beyond which it decreases. The maximum net kick velocity is about $190 {\rm km/s}$ for for the boost parameter considered. For $\rho_0 \geq 50^{\circ}$ the distribution is a monotonous function of $\eta$. The angular patterns of the gravitational waves emitted are examined. Our analysis includes the two polarization modes present in wave zone curvature.


title: Uniqueness for Some Higher-Order Geometric Flows

abstract: We show that solutions to certain higher-order intrinsic geometric flows on a compact manifold, including some flows generated by the ambient obstruction tensor, are unique. With the goal of providing a complete self-contained proof, details surrounding map covariant derivatives and a careful application of the DeTurck trick are provided.


title: Tunnel Magnetoresistance scan of a pristine three-dimensional topological insulator

abstract: Though the Fermi surface of surface states of a 3D topological insulator (TI) has zero magnetization, an arbitrary segment of the full Fermi surface has a unique magnetic moment consistent with the type of spin-momentum locking in hand. We propose a three-terminal set up, which directly couples to the magnetization of a chosen segment of a Fermi surface hence leading to a finite tunnel magnetoresistance (TMR) response of the nonmagnetic TI surface states, when coupled to spin polarized STM probe. This multiterminal TMR not only provides a unique signature of spin-momentum locking for a pristine TI but also provides a direct measure of momentum resolved out of plane polarization of hexagonally warped Fermi surfaces relevant for $Bi_2Te_3$, which could be as comprehensive as spin-resolved ARPES. Implication of this unconventional TMR is also discussed in the broader context of 2D spin-orbit (SO) materials.


title: Improved limit on a temporal variation of $m_p/m_e$ from comparisons of Yb$^+$ and Cs atomic clocks

abstract: Accurate measurements of different transition frequencies between atomic levels of the electronic and hyperfine structure over time are used to investigate temporal variations of the fine structure constant $\alpha$ and the proton-to-electron mass ratio $\mu$. We measure the frequency of the $^2S_{1/2}\rightarrow {^2F_{7/2}}$ electric octupole (E3) transition in $^{171}$Yb$^+$ against two caesium fountain clocks as $f(E3) = 642\,121\,496\,772\,645.36(25)$~Hz with an improved fractional uncertainty of $3.9\times 10^{-16}$. This transition frequency shows a strong sensitivity to changes of $\alpha$. Together with a number of previous and recent measurements of the $^2S_{1/2}\rightarrow {^2D_{3/2}}$ electric quadrupole transition in $^{171}$Yb$^+$ and with data from other elements, a least-squares analysis yields $(1/\alpha)(d\alpha/dt)=-0.20(20)\times 10^{-16}/\mathrm{yr}$ and $(1/\mu)(d\mu/dt)=-0.5(1.6)\times 10^{-16}/\mathrm{yr}$, confirming a previous limit on $d\alpha/dt$ and providing the most stringent limit on $d \mu/dt$ from laboratory experiments.


title: SoundLoc: Acoustic Method for Indoor Localization without Infrastructure

abstract: Identifying locations of occupants is beneficial to energy management in buildings. A key observation in indoor environment is that distinct functional areas are typically controlled by separate HVAC and lighting systems and room level localization is sufficient to provide a powerful tool for energy usage reduction by occupancy-based actuation of the building facilities. Based upon this observation, this paper focuses on identifying the room where a person or a mobile device is physically present. Existing room localization methods, however, require special infrastructure to annotate rooms. SoundLoc is a room-level localization system that exploits the intrinsic acoustic properties of individual rooms and obviates the needs for infrastructures. As we show in the study, rooms' acoustic properties can be characterized by Room Impulse Response (RIR). Nevertheless, obtaining precise RIRs is a time-consuming and expensive process. The main contributions of our work are the following. First, a cost-effective RIR measurement system is implemented and the Noise Adaptive Extraction of Reverberation (NAER) algorithm is developed to estimate room acoustic parameters in noisy conditions. Second, a comprehensive physical and statistical analysis of features extracted from RIRs is performed. Also, SoundLoc is evaluated using the dataset consisting of ten (10) different rooms. The overall accuracy of 97.8% achieved demonstrates the potential to be integrated into automatic mapping of building space.


title: The limit of binomial means of a sequence

abstract: For a sequence $\{a_n\}_{n\geq 0}$ of real numbers and for a parameter $0<p<1$, we define the sequence of its arithmetic means $\{a^*_n\}_{n\geq 0}$ and the sequence of its $p$-binomial means $\{a^p_n\}_{n\geq 0}$ as \begin{align*} a^*_n=\frac{1}{n+1}\sum_{i=0}^n a_i & & \textrm{and} && a^p_n=\sum_{i=0}^n\binom{n}{i}p^i(1-p)^{n-i} a_i. \end{align*} We compare the convergence of sequences $\{a_n\}_{n\geq 0}$, $\{a_n^*\}_{n\geq 0}$ and $\{a_n^p\}_{n\geq 0}$ for various $0<p<1$, i.e. we analyze when the convergence of one sequence implies the convergence of the other. While the sequence $\{a^*_n\}_{n\geq 0}$, known also as the sequence of Ces\`{a}ro means of a sequence, is well studied in the literature, the results about $\{a^p_n\}_{n\geq 0}$ are hard to find. Our main result shows that, if $\{a_n\}_{n\geq 0}$ is a sequence of non-negative real numbers such that $\{a^p_n\}_{n\geq 0}$ converges to $a\in\mathbb{R}\cup\{\infty\}$ for some $0<p<1$, then $\{a^*_n\}_{n\geq 0}$ also converges to $a$. We give an application of this result on finite Markov chains.


title: Hole Spin Pumping and Re-pumping in a p-type {\delta}-doped InAs Quantum Dot

abstract: We have grown high quality p-type {\delta}-doped InAs quantum dots and have demonstrated coherent spin pumping and repumping of a hole spin in a positively charged quantum dot by means of a single-laser driving scheme under a high magnetic field in the Voigt configuration. Modeling of our system shows excellent qualitative agreement with the experimental findings and further explores the performance of the single-laser scheme for spin pumping and re-pumping.


title: A CUSUM type change detection test based on martingale differences

abstract: The paper is about detecting changes in the parameters of certain parameterized stochastic models. We apply CUSUM (Cumulated Sums) type test statistics that are based on martingale difference sequences.


title: Control Function Assisted IPW Estimation with a Secondary Outcome in Case-Control Studies

abstract: Case-control studies are designed towards studying associations between risk factors and a single, primary outcome. Information about additional, secondary outcomes is also collected, but association studies targeting such secondary outcomes should account for the case-control sampling scheme, or otherwise results may be biased. Often, one uses inverse probability weighted (IPW) estimators to estimate population effects in such studies. However, these estimators are inefficient relative to estimators that make additional assumptions about the data generating mechanism. We propose a class of estimators for the effect of risk factors on a secondary outcome in case-control studies, when the mean is modeled using either the identity or the log link. The proposed estimator combines IPW with a mean zero control function that depends explicitly on a model for the primary disease outcome. The efficient estimator in our class of estimators reduces to standard IPW when the model for the primary disease outcome is unrestricted, and is more efficient than standard IPW when the model is either parametric or semiparametric.


title: Particle Filtering and Smoothing Using Windowed Rejection Sampling

abstract: "Particle methods" are sequential Monte Carlo algorithms, typically involving importance sampling, that are used to estimate and sample from joint and marginal densities from a collection of a, presumably increasing, number of random variables. In particular, a particle filter aims to estimate the current state $X_{n}$ of a stochastic system that is not directly observable by estimating a posterior distribution $\pi(x_{n}|y_{1},y_{2}, \ldots, y_{n})$ where the $\{Y_{n}\}$ are observations related to the $\{X_{n}\}$ through some measurement model $\pi(y_{n}|x_{n})$. A particle smoother aims to estimate a marginal distribution $\pi(x_{i}|y_{1},y_{2}, \ldots, y_{n})$ for $1 \leq i < n$. Particle methods are used extensively for hidden Markov models where $\{X_{n}\}$ is a Markov chain as well as for more general state space models. Existing particle filtering algorithms are extremely fast and easy to implement. Although they suffer from issues of degeneracy and "sample impoverishment", steps can be taken to minimize these problems and overall they are excellent tools for inference. However, if one wishes to sample from a posterior distribution of interest, a particle filter is only able to produce dependent draws. Particle smoothing algorithms are complicated and far less robust, often requiring cumbersome post-processing, "forward-backward" recursions, and multiple passes through subroutines. In this paper we introduce an alternative algorithm for both filtering and smoothing that is based on rejection sampling "in windows" . We compare both speed and accuracy of the traditional particle filter and this "windowed rejection sampler" (WRS) for several examples and show that good estimates for smoothing distributions are obtained at no extra cost.


title: The Effect of Doppler Broadening on the $6.3 \ PeV$ $W^-$ Resonance in $\bar{\nu}_e e^-$ Collisions

abstract: We calculate the Doppler broadening of the $W^-$ resonance produced in $\bar{\nu}_e e^-$ collisions of cosmic anti-neutrinos with $E_{\nu}\approx 6.3 \ PeV$ with electrons in atoms up to Iron. Revisiting this issue is prompted by recent observations of PeV neutrinos by Ice-Cube. Despite its poor energy resolution, the $20\%$ Doppler broadening of the resonance due to electronic motions can produce observable effects via non-linear neutrino absorption near the resonance. The attendant suppression of the peak cross section allows $\bar{\nu}_e$ to travel correspondingly longer distances. While this effect is unlikely to be directly detected in the near future, it may facilitate terrestrial tomography at depths of $\sim 10 \ km$, complementing deeper explorations using the more frequent nuclear interactions at lower energies.


title: In Defense of MinHash Over SimHash

abstract: MinHash and SimHash are the two widely adopted Locality Sensitive Hashing (LSH) algorithms for large-scale data processing applications. Deciding which LSH to use for a particular problem at hand is an important question, which has no clear answer in the existing literature. In this study, we provide a theoretical answer (validated by experiments) that MinHash virtually always outperforms SimHash when the data are binary, as common in practice such as search. The collision probability of MinHash is a function of resemblance similarity ($\mathcal{R}$), while the collision probability of SimHash is a function of cosine similarity ($\mathcal{S}$). To provide a common basis for comparison, we evaluate retrieval results in terms of $\mathcal{S}$ for both MinHash and SimHash. This evaluation is valid as we can prove that MinHash is a valid LSH with respect to $\mathcal{S}$, by using a general inequality $\mathcal{S}^2\leq \mathcal{R}\leq \frac{\mathcal{S}}{2-\mathcal{S}}$. Our worst case analysis can show that MinHash significantly outperforms SimHash in high similarity region. Interestingly, our intensive experiments reveal that MinHash is also substantially better than SimHash even in datasets where most of the data points are not too similar to each other. This is partly because, in practical data, often $\mathcal{R}\geq \frac{\mathcal{S}}{z-\mathcal{S}}$ holds where $z$ is only slightly larger than 2 (e.g., $z\leq 2.1$). Our restricted worst case analysis by assuming $\frac{\mathcal{S}}{z-\mathcal{S}}\leq \mathcal{R}\leq \frac{\mathcal{S}}{2-\mathcal{S}}$ shows that MinHash indeed significantly outperforms SimHash even in low similarity region. We believe the results in this paper will provide valuable guidelines for search in practice, especially when the data are sparse.


title: Variable Accretion Processes in the Young Binary-Star System UY Aur

abstract: We present new K-band spectroscopy of the UY Aur binary star system. Our data are the first to show H$_{2}$ emission in the spectrum of UY Aur A and the first to spectrally resolve the Br{\gamma} line in the spectrum of UY Aur B. We see an increase in the strength of the Br{\gamma} line in UY Aur A and a decrease in Br{\gamma} and H$_{2}$ line luminosity for UY Aur B compared to previous studies. Converting Br{\gamma} line luminosity to accretion rate, we infer that the accretion rate onto UY Aur A has increased by $2 \times 10^{-9}$ M$_{\odot}$ yr$^{-1}$ per year since a rate of zero was observed in 1994. The Br{\gamma} line strength for UY Aur B has decreased by a factor of 0.54 since 1994, but the K-band flux has increased by 0.9 mags since 1998. The veiling of UY Aur B has also increased significantly. These data evince a much more luminous disk around UY Aur B. If the lower Br{\gamma} luminosity observed in the spectrum of UY Aur B indicates an intrinsically smaller accretion rate onto the star, then UY Aur A now accretes at a higher rate than UY Aur B. However, extinction at small radii or mass pile-up in the circumstellar disk could explain decreased Br{\gamma} emission around UY Aur B even when the disk luminosity implies an increased accretion rate. In addition to our scientific results for the UY Aur system, we discuss a dedicated pipeline we have developed for the reduction of echelle-mode data from the ARIES spectrograph.


title: On Gaussian multiplicative chaos

abstract: We propose a new definition of the Gaussian multiplicative chaos (GMC) and an approach based on the relation of subcritical GMC to randomized shifts of a Gaussian measure. Using this relation we prove general uniqueness and convergence results for subcritical GMC that hold for Gaussian fields with arbitrary covariance kernels.


title: Irreversibility and Entanglement Spectrum Statistics in Quantum Circuits

abstract: We show that in a quantum system evolving unitarily under a stochastic quantum circuit the notions of irreversibility, universality of computation, and entanglement are closely related. As the state evolves from an initial product state, it gets asymptotically maximally entangled. We define irreversibility as the failure of searching for a disentangling circuit using a Metropolis-like algorithm. We show that irreversibility corresponds to Wigner-Dyson statistics in the level spacing of the entanglement eigenvalues, and that this is obtained from a quantum circuit made from a set of universal gates for quantum computation. If, on the other hand, the system is evolved with a non-universal set of gates, the statistics of the entanglement level spacing deviates from Wigner-Dyson and the disentangling algorithm succeeds. These results open a new way to characterize irreversibility in quantum systems.


title: Kernel Nonnegative Matrix Factorization Without the Curse of the Pre-image - Application to Unmixing Hyperspectral Images

abstract: The nonnegative matrix factorization (NMF) is widely used in signal and image processing, including bio-informatics, blind source separation and hyperspectral image analysis in remote sensing. A great challenge arises when dealing with a nonlinear formulation of the NMF. Within the framework of kernel machines, the models suggested in the literature do not allow the representation of the factorization matrices, which is a fallout of the curse of the pre-image. In this paper, we propose a novel kernel-based model for the NMF that does not suffer from the pre-image problem, by investigating the estimation of the factorization matrices directly in the input space. For different kernel functions, we describe two schemes for iterative algorithms: an additive update rule based on a gradient descent scheme and a multiplicative update rule in the same spirit as in the Lee and Seung algorithm. Within the proposed framework, we develop several extensions to incorporate constraints, including sparseness, smoothness, and spatial regularization with a total-variation-like penalty. The effectiveness of the proposed method is demonstrated with the problem of unmixing hyperspectral images, using well-known real images and results with state-of-the-art techniques.


title: Global simulations of axisymmetric radiative black hole accretion disks in general relativity with a sub-grid magnetic dynamo

abstract: We present a sub-grid model that emulates the magnetic dynamo operating in magnetized accretion disks. We have implemented this model in the general relativisic radiation magnetohydrodynamic (GRRMHD) code \koral, using results from local shearing sheet simulations of the magnetorotational instability to fix the parameters of the dynamo. With the inclusion of this dynamo, we are able to run 2D axisymmetric GRRMHD simulations of accretion disks for arbitrarily long times. The simulated disks exhibit sustained turbulence, with the poloidal and toroidal magnetic field components driven towards a state similar to that seen in 3D studies. Using this dynamo code, we present a set of long-duration global simulations of super-Eddington, optically-thick disks around non-spinning and spinning black holes. Super-Eddington disks around non-rotating black holes exhibit a surprisingly large efficiency, $\eta\approx0.04$, independent of the accretion rate, where we measure efficiency in terms of the total energy output, both radiation and mechanical, flowing out to infinity. Super-Eddington disks around spinning black holes are even more efficient, and appear to extract black hole rotational energy through a process similar to the Blandford-Znajek mechanism. All the simulated models are characterized by highly super-Eddington radiative fluxes collimated along the rotation axis. We also present a set of simulations that were designed to have Eddington or slightly sub-Eddington accretion rates ($\dot{M} \lesssim 2\dot M_{\rm Edd}$). None of these models reached a steady state. Instead, the disks collapsed as a result of runaway cooling, presumably because of a thermal instability.


title: Subspace Restricted Boltzmann Machine

abstract: The subspace Restricted Boltzmann Machine (subspaceRBM) is a third-order Boltzmann machine where multiplicative interactions are between one visible and two hidden units. There are two kinds of hidden units, namely, gate units and subspace units. The subspace units reflect variations of a pattern in data and the gate unit is responsible for activating the subspace units. Additionally, the gate unit can be seen as a pooling feature. We evaluate the behavior of subspaceRBM through experiments with MNIST digit recognition task, measuring reconstruction error and classification error.


title: Conditioning and covariance on caterpillars

abstract: Let $X_1, \dots, X_n$ be joint $\{ \pm 1\}$-valued random variables. It is known that conditioning on a random subset of $O(1/\epsilon^2)$ of them reduces their average pairwise covariance to below $\epsilon$ (in expectation). We conjecture that $O(1/\epsilon^2)$ can be improved to $O(1/\epsilon)$. The motivation for the problem and our conjectured improvement comes from the theory of global correlation rounding for convex relaxation hierarchies. We suggest attempting the conjecture in the case that $X_1, \dots, X_n$ are the leaves of an information flow tree. We prove the conjecture in the case that the information flow tree is a caterpillar graph (similar to a two-state hidden Markov model).


title: $\alpha$-Molecules

abstract: Within the area of applied harmonic analysis, various multiscale systems such as wavelets, ridgelets, curvelets, and shearlets have been introduced and successfully applied. The key property of each of those systems are their (optimal) approximation properties in terms of the decay of the $L^2$-error of the best $N$-term approximation for a certain class of functions. In this paper, we introduce the general framework of $\alpha$-molecules, which encompasses most multiscale systems from applied harmonic analysis, in particular, wavelets, ridgelets, curvelets, and shearlets as well as extensions of such with $\alpha$ being a parameter measuring the degree of anisotropy, as a means to allow a unified treatment of approximation results within this area. Based on an $\alpha$-scaled index distance, we first prove that two systems of $\alpha$-molecules are almost orthogonal. This leads to a general methodology to transfer approximation results within this framework, provided that certain consistency and time-frequency localization conditions of the involved systems of $\alpha$-molecules are satisfied. We finally utilize these results to enable the derivation of optimal sparse approximation results \msch{for} a specific class of cartoon-like functions by sufficient conditions on the 'control' parameters of a system of $\alpha$-molecules.


title: Corecursive Algebras, Corecursive Monads and Bloom Monads

abstract: An algebra is called corecursive if from every coalgebra a unique coalgebra-to-algebra homomorphism exists into it. We prove that free corecursive algebras are obtained as coproducts of the terminal coalgebra (considered as an algebra) and free algebras. The monad of free corecursive algebras is proved to be the free corecursive monad, where the concept of corecursive monad is a generalization of Elgot's iterative monads, analogous to corecursive algebras generalizing completely iterative algebras. We also characterize the Eilenberg-Moore algebras for the free corecursive monad and call them Bloom algebras.


title: Schur indices in GAP: wedderga 4.6+

abstract: We describe a algorithms and their implementations that calculate local and global Schur indices of ordinary irreducible characters of finite groups, cyclotomic algebras over abelian number fields, and rational quaternion algebras. These functions are available with the latest release of the GAP package wedderga, versions 4.6 and higher.