title: Vacuum currents induced by a magnetic flux around a cosmic string with finite core

abstract: We evaluate the Hadamard function and the vacuum expectation value of the current density for a massive complex scalar field in the generalized geometry of a straight cosmic string with a finite core enclosing an arbitrary distributed magnetic flux along the string axis. For the interior geometry, a general cylindrically symmetric static metric tensor is used with finite support. In the region outside the core, both the Hadamard function and the current density are decomposed into the idealized zero-thickness cosmic string and core-induced contributions. The only nonzero component corresponds to the azimuthal current. The zero-thickness part of the latter is a periodic function of the magnetic flux inside the core, with the period equal to the quantum flux. As a consequence of the direct interaction of the quantum field with the magnetic field inside the penetrable core, the core-induced contribution, in general, is not a periodic function of the flux. In addition, the vacuum current, in general, is not a monotonic function of the distance from the string and may change the sign. For a general model of the core interior, we also evaluate the magnetic fields generated by the vacuum current. As applications of the general results, we have considered an impenetrable core modeled by Robin boundary condition, a core with the Minkowski-like interior and a core with a constant positive curvature space. Various exactly solvable distributions of the magnetic flux are discussed.


title: New Upper and Lower bounds for the Trapezoid inequality of absolutely continuous functions and Applications

abstract: In this paper, new upper and lower bounds for the Trapezoid inequality of absolutely continuous functions are obtained. Applications to some special means are provided as well.


title: The Gravitational Wave Observatory Designer: Sensitivity Limits of Spaceborne Detectors

abstract: The most promising concept for low frequency gravitational wave observatories are laser interferometric detectors in space. It is usually assumed that the noise floor for such a detector is dominated by optical shot noise in the signal readout. For this to be true, a careful balance of mission parameters is crucial to keep all other parasitic disturbances below shot noise. We developed a web application that uses over 30 input parameters and considers many important technical noise sources and noise suppression techniques. It optimizes free parameters automatically and generates a detailed report on all individual noise contributions. Thus you can easily explore the entire parameter space and design a realistic gravitational wave observatory. In this document we describe the different parameters, present all underlying calculations, and compare the final observatory's sensitivity with astrophysical sources of gravitational waves. We use as an example parameters currently assumed to be likely applied to a space mission to be launched in 2034 by the European Space Agency. The web application itself is publicly available on the Internet at http://spacegravity.org/designer.


title: An embedded active nucleus in the OH megamaser galaxy IRAS16399-0937

abstract: We present a multiwavelength study of the OH Megamaser galaxy (OHMG) IRAS16399-0937, based on new HST/ACS F814W and H$\alpha$+[NII] images and archive data from HST, 2MASS, Spitzer, Herschel and the VLA. This system has a double nucleus, whose northern (IRAS16399N) and southern (IRAS16399S) components have a projected separation of $\sim$ 6'' (3.4 kpc) and have previously been identified based on optical spectra as a Low Ionization Nuclear Emission Line Region (LINER) and starburst nucleus, respectively. The nuclei are embedded in a tidally distorted common envelope, in which star formation is mostly heavily obscured. The infrared spectrum is dominated by strong polycyclic aromatic hydrocarbon (PAH), but deep silicate and molecular absorption features are also present, and are strongest in the IRAS16399N nucleus. The 0.435 - 500$\mu$m SED was fitted with a model including stellar, ISM and AGN torus components using our new MCMC code, clumpyDREAM. The results indicate that the IRAS16399N contains an AGN (L$_{bol} \sim 10^{44}$ ergs/s) deeply embedded in a quasi-spherical distribution of optically-thick clumps with a covering fraction $\approx1$. We suggest that these clumps are the source of the OHM emission in IRAS16399-0937. The high torus covering fraction precludes AGN-photoionization as the origin of the LINER spectrum, however, the spectrum is consistent with shocks (v $\sim100-200$ km s$^{-1}$). We infer that the $\sim10^8$ M$_{\odot}$ black-hole in IRAS16399N is accreting at a small fraction ($\sim1$%) of its Eddington rate. The low accretion-rate and modest nuclear SFRs suggest that while the gas-rich major merger forming the IRAS16399-0937 system has triggered widespread star formation, the massive gas inflows expected from merger simulations have not yet fully developed.


title: Hidden Symmetries of Dynamics in Classical and Quantum Physics

abstract: This article reviews the role of hidden symmetries of dynamics in the study of physical systems, from the basic concepts of symmetries in phase space to the forefront of current research. Such symmetries emerge naturally in the description of physical systems as varied as non-relativistic, relativistic, with or without gravity, classical or quantum, and are related to the existence of conserved quantities of the dynamics and integrability. In recent years their study has grown intensively, due to the discovery of non-trivial examples that apply to different types of theories and different numbers of dimensions. Applications encompass the study of integrable systems such as spinning tops, the Calogero model, systems described by the Lax equation, the physics of higher dimensional black holes, the Dirac equation, supergravity with and without fluxes, providing a tool to probe the dynamics of non-linear systems.


title: How Well Can Graphs Represent Wireless Interference?

abstract: Efficient use of a wireless network requires that transmissions be grouped into feasible sets, where feasibility means that each transmission can be successfully decoded in spite of the interference caused by simultaneous transmissions. Feasibility is most closely modeled by a signal-to-interference-plus-noise (SINR) formula, which unfortunately is conceptually complicated, being an asymmetric, cumulative, many-to-one relationship. We re-examine how well graphs can capture wireless receptions as encoded in SINR relationships, placing them in a framework in order to understand the limits of such modelling. We seek for each wireless instance a pair of graphs that provide upper and lower bounds on the feasibility relation, while aiming to minimize the gap between the two graphs. The cost of a graph formulation is the worst gap over all instances, and the price of (graph) abstraction is the smallest cost of a graph formulation. We propose a family of conflict graphs that is parameterized by a non-decreasing sub-linear function, and show that with a judicious choice of functions, the graphs can capture feasibility with a cost of $O(\log^* \Delta)$, where $\Delta$ is the ratio between the longest and the shortest link length. This holds on the plane and more generally in doubling metrics. We use this to give greatly improved $O(\log^* \Delta)$-approximation for fundamental link scheduling problems with arbitrary power control. We explore the limits of graph representations and find that our upper bound is tight: the price of graph abstraction is $\Omega(\log^* \Delta)$. We also give strong impossibility results for general metrics, and for approximations in terms of the number of links.


title: Measurement of the Z+b-jet cross-section in pp collisions at $\sqrt{s}=7{\mathrm{\,Te\kern -0.1em V}}$ in the forward region

abstract: The associated production of a Z boson or an off-shell photon $\gamma^*$ with a bottom quark in the forward region is studied using proton-proton collisions at a centre-of-mass energy of $7{\mathrm{\,Te\kern -0.1em V}}$. The Z bosons are reconstructed in the ${\text{Z}/\gamma^*}\!\rightarrow{\mu^{+}\mu^{-}}$ final state from muons with a transverse momentum larger than $20{\mathrm{\,Ge\kern -0.1em V}}$, while two transverse momentum thresholds are considered for jets ($10{\mathrm{\,Ge\kern -0.1em V}}$ and $20{\mathrm{\,Ge\kern -0.1em V}}$). Both muons and jets are reconstructed in the pseudorapidity range $2.0 < \eta < 4.5$. The results are based on data corresponding to $1.0{\,\mbox{fb}^{-1}}$ recorded in 2011 with the LHCb detector. The measurement of the Z+b-jet cross-section is normalized to the Z+jet cross-section. The measured cross-sections are \begin{equation*} \sigma(\text{$\text{Z}/\gamma^*(\mu^{+}\mu^{-})$+b-jet}) = 295 \pm 60~(\text{stat}) \pm 51~(\text{syst}) \pm 10~(\text{lumi}) {\,\mbox{fb}} \end{equation*} for ${p_{\rm T}}$(jet)$>10{\mathrm{\,Ge\kern -0.1em V}}$, and \begin{equation*} \sigma(\text{$\text{Z}/\gamma^*(\mu^{+}\mu^{-})$+b-jet}) = 128 \pm 36~(\text{stat}) \pm 22~(\text{syst}) \pm 5~(\text{lumi}) {\,\mbox{fb}} \end{equation*} for ${p_{\rm T}}$(jet)$>20{\mathrm{\,Ge\kern -0.1em V}}$.


title: Extensive light profile fitting of galaxy-scale strong lenses

abstract: We investigate the merits of a massive forward modeling of ground-based optical imaging as a diagnostic for the strong lensing nature of Early-Type Galaxies, in the light of which blurred and faint Einstein rings can hide. We simulate several thousand mock strong lenses under ground- and space-based conditions as arising from the deflection of an exponential disk by a foreground de Vaucouleurs light profile whose lensing potential is described by a Singular Isothermal Ellipsoid. We then fit for the lensed light distribution with sl_fit after having subtracted the foreground light emission off (ideal case) and also after having fitted the deflector's light with galfit. By setting thresholds in the output parameter space, we can decide the lens/not-a-lens status of each system. We finally apply our strategy to a sample of 517 lens candidates present in the CFHTLS data to test the consistency of our selection approach. The efficiency of the fast modeling method at recovering the main lens parameters like Einstein radius, total magnification or total lensed flux, is quite comparable under CFHT and HST conditions when the deflector is perfectly subtracted off (only possible in simulations), fostering a sharp distinction between the good and the bad candidates. Conversely, for a more realistic subtraction, a substantial fraction of the lensed light is absorbed into the deflector's model, which biases the subsequent fitting of the rings and then disturbs the selection process. We quantify completeness and purity of the lens finding method in both cases. This suggests that the main limitation currently resides in the subtraction of the foreground light. Provided further enhancement of the latter, the direct forward modeling of large numbers of galaxy-galaxy strong lenses thus appears tractable and could constitute a competitive lens finder in the next generation of wide-field imaging surveys.


title: Towards a two-parameter q-deformation of AdS_3 x S^3 x M^4 superstrings

abstract: We construct a two-parameter deformation of the Metsaev-Tseytlin action for supercosets with isometry group of the form G x G. The resulting action is classically integrable and is Poisson-Lie symmetric suggesting that the symmetry of the model is q-deformed, U_q_L(G) x U_q_R(G). Focusing on the cases relevant for strings moving in AdS_3 x S^3 x T^4 and AdS_3 x S^3 x S^3 x S^1, we analyze the corresponding deformations of the AdS_3 and S^3 metrics. We also construct a two-parameter $q$-deformation of the u(1) + psu(1|1)^2 x u(1) x R^3-invariant R-matrix and closure condition, which underlie the light-cone gauge S-matrix and dispersion relation of the aforementioned string theories. With the appropriate identification of parameters, the near-BMN limit of the dispersion relation is shown to agree with that found from the deformed supercoset sigma model.


title: An Interesting Property of LPCs for Sonorant Vs Fricative Discrimination

abstract: Linear prediction (LP) technique estimates an optimum all-pole filter of a given order for a frame of speech signal. The coefficients of the all-pole filter, 1/A(z) are referred to as LP coefficients (LPCs). The gain of the inverse of the all-pole filter, A(z) at z = 1, i.e, at frequency = 0, A(1) corresponds to the sum of LPCs, which has the property of being lower (higher) than a threshold for the sonorants (fricatives). When the inverse-tan of A(1), denoted as T(1), is used a feature and tested on the sonorant and fricative frames of the entire TIMIT database, an accuracy of 99.07% is obtained. Hence, we refer to T(1) as sonorant-fricative discrimination index (SFDI). This property has also been tested for its robustness for additive white noise and on the telephone quality speech of the NTIMIT database. These results are comparable to, or in some respects, better than the state-of-the-art methods proposed for a similar task. Such a property may be used for segmenting a speech signal or for non-uniform frame-rate analysis.


title: The Peak/Dip Picture of the Cosmic Web

abstract: The initial shear field plays a central role in the formation of large-scale structures, and in shaping the geometry, morphology, and topology of the cosmic web. We discuss a recent theoretical framework for the shear tensor, termed the `peak/dip picture', which accounts for the fact that halos/voids may form from local extrema of the density field - rather than from random spatial positions; the standard Doroshkevich's formalism is generalized, to include correlations between the density Hessian and shear field at special points in space around which halos/voids may form. We then present the `peak/dip excursion-set-based' algorithm, along with its most recent applications - merging peaks theory with the standard excursion set approach.


title: Ionized gas velocity dispersion and multiple supernova explosions

abstract: Using 3D numerical simulations we study the evolution of the H$\alpha$ intensity and velocity dispersion for single and multiple supenova (SN) explosions. We find that the $I_{\rm H\alpha}-\sigma$ diagram obtained for simulated gas flows is similar in shape to that observed in dwarf galaxies. We conclude that colliding SN shells with significant difference in age are resposible for high velocity dispersion that reaches values high as $\simgt 100$kms$^{-1}$. Such a high velocity dispersion could be hardly got for a single SN remnant. Peaks of velocity dispersion on the $I_{\rm H\alpha}-\sigma$ diagram may correspond to several stand-alone or merged SN remnants with moderately different ages. The procedure of the spatial resolution degrading in the H$\alpha$ intensity and velocity dispersion maps makes the simulated $I_{\rm H\alpha}-\sigma$ diagrams close to those observed in dwarf galaxies not only in shape, but also quantitatively.


title: Adiabatic many-body state preparation and information transfer in quantum dot arrays

abstract: Quantum simulation of many-body systems are one of the most interesting tasks of quantum technology. Among them is the preparation of a many-body system in its ground state when the vanishing energy gap makes the cooling mechanisms ineffective. Adiabatic theorem, as an alternative to cooling, can be exploited for driving the many-body system to its ground state. In this paper, we study two most common disorders in quantum dot arrays, namely exchange coupling fluctuations and hyperfine interaction, in adiabatically preparation of ground state in such systems. We show that the adiabatic ground state preparation is highly robust against those disorder effects making it good analog simulator. Moreover, we also study the adiabatic classical information transfer, using singlet-triplet states, across a spin chain. In contrast to ground state preparation the transfer mechanism is highly affected by disorder and in particular, the hyperfine interaction is very destructive for the performance. This suggests that for communication tasks across such arrays adiabatic evolution is not as effective and quantum quenches could be preferable.


title: Beyond Anderson Localization in 1D: Anomalous Localization of Microwaves in Random Waveguides

abstract: Experimental evidence demonstrating that anomalous localization of waves can be induced in a controllable manner is reported. A microwave waveguide with dielectric slabs randomly placed is used to confirm the presence of anomalous localization. If the random spacing between slabs follows a distribution with a power-law tail (L\'evy-type distribution), unconventional properties in the microwave-transmission fluctuations take place revealing the presence of anomalous localization. We study both theoretically and experimentally the complete distribution of the transmission through random waveguides characterized by $\alpha=1/2$ ("L\'evy waveguides") and $\alpha=3/4$, $\alpha$ being the exponent of the power-law tail of the L\'evy-type distribution. As we show, the transmission distributions are determined by only two parameters, both of them experimentally accessible. Effects of anomalous localization on the transmission are compared with those from the standard Anderson localization.


title: Integer points on spheres and their orthogonal grids

abstract: The set of primitive vectors on large spheres in the euclidean space of dimension d>2 equidistribute when projected on the unit sphere. We consider here a refinement of this problem concerning the direction of the vector together with the shape of the lattice in its orthogonal complement. Using unipotent dynamics we obtained the desired equidistribution result in dimension d>5 and in dimension d=4,5 under a mild congruence condition on the square of the radius. The case of d=3 is considered in a separate paper.


title: Convex and concave decompositions of affine $3$-manifolds

abstract: A (flat) affine $3$-manifold is a $3$-manifold with an atlas of charts to an affine space $\mathbb{R}^3$ with transition maps in the affine transformation group $\mathrm{Aff}(\mathbb{R}^3)$. We will show that a connected closed affine $3$-manifold is either an affine Hopf $3$-manifold or decomposes canonically to concave affine submanifolds with incompressible boundary, toral $\pi$-submanifolds and $2$-convex affine manifolds, each of which is an irreducible $3$-manifold. It follows that if there is no toral $\pi$-submanifold, then $M$ is prime. Finally, we prove that if a closed affine manifold is covered by a domain in $\mathbb{R}^{n}$, then $M$ is irreducible or is an affine Hopf manifold.


title: Anatomy and efficiency of urban multimodal mobility

abstract: The growth of transportation networks and their increasing interconnections, although positive, has the downside effect of an increasing complexity which make them difficult to use, to assess, and limits their efficiency. On average in the UK, 23% of travel time is lost in connections for trips with more than one mode, and the lack of synchronization decreases very slowly with population size. This lack of synchronization between modes induces differences between the theoretical quickest trip and the `time-respecting' path, which takes into account waiting times at interconnection nodes. We analyse here the statistics of these paths on the multilayer, temporal network of the entire, multimodal british public transportation system. We propose a statistical decomposition -- the `anatomy' -- of trips in urban areas, in terms of riding, waiting and walking times, and which shows how the temporal structure of trips varies with distance and allows us to compare different cities. Weaknesses in systems can be either insufficient transportation speed or service frequency, but the key parameter controlling their global efficiency is the total number of stop events per hour for all modes. This analysis suggests the need for better optimization strategies, adapted to short, long unimodal or multimodal trips.


title: Mapping cone formula in Heegaard Floer homology and Dehn surgery on knots in $S^3$

abstract: We write down an explicit formula for the $+$ version of the Heegaard Floer homology (as an absolutely graded vector space over an arbitrary field) of the results of Dehn surgery on a knot $K$ in $S^3$ in terms of homological data derived from $CFK^{\infty}(K)$. This allows us to prove some results about Dehn surgery on knots in $S^3$. In particular, we show that for a fixed manifold there are only finitely many alternating knots that can produce it by surgery. This is an improvement on a recent result by Lackenby and Purcell. We also derive a lower bound on the genus of knots depending on the manifold they give by surgery. Some new restrictions on Seifert fibred surgery are also presented.


title: The cosmic web in CosmoGrid void regions

abstract: We study the formation and evolution of the cosmic web, using the high-resolution CosmoGrid \$\Lambda\$CDM simulation. In particular, we investigate the evolution of the large-scale structure around void halo groups, and compare this to observations of the VGS-31 galaxy group, which consists of three interacting galaxies inside a large void. The structure around such haloes shows a great deal of tenuous structure, with most of such systems being embedded in intra-void filaments and walls. We use the Nexus+ algorithm to detect walls and filaments in CosmoGrid, and find them to be present and detectable at every scale. The void regions embed tenuous walls, which in turn embed tenuous filaments. We hypothesize that the void galaxy group of VGS-31 formed in such an environment.


title: The All Sky Automated Survey. The Catalog of Bright Variable Stars in the I-band, South of Declination +28

abstract: This paper presents the results of our extensive search for the bright variable stars in approximately 30000 square degrees of the south sky in the I-band data collected by 9 deg x 9 deg camera of the All Sky Automated Survey between 2002 and 2009. Lists of over 27000 variable stars brighter than 9 mag at maximum light, with amplitudes ranging from 0.02 mag to 7 mag and variability time-scales from hours to years, as well as corresponding light curves are provided. Automated classification algorithm based on stellar properties (period, Fourier coefficients, 2MASS J, H, K, colors, ASAS V-band data) was used to roughly classify objects. Despite low spatial resolution of the ASAS data (~15 arcs) we cross-identified all objects with other available data sources. Coordinates of the most probable 2MASS counterparts are provided. 27705 stars brighter than I=9 mag were found to be variable, of which 7842 objects were detected to be variable for the first time. Brief statistics and discussion of the presented data is provided. All the photometric data is available over the Internet at http://www.astrouw.edu.pl/~gp/asas/AsasBrightI.html


title: Notes on the Infinity-Laplace Equation

abstract: These notes are written up after my lectures at the University of Pittsburgh in March 2014 and at Tsinghua University in May 2014. My objective is the $\infty$-Laplace Equation, a marvellous kin to the ordinary Laplace Equation. The $\infty$-Laplace Equation has delightful counterparts to the Dirichlet integral, the Mean Value Theorem, the Brownian Motion, Harnack's Inequality and so on. It has applications to image processing and to mass transfer problems and provides optimal Lipschitz extensions of boundary values. My treaty of this "fully non-linear" degenerate equation is far from complete and generalizations are deliberately avoided.


title: Streaming, Memory Limited Algorithms for Community Detection

abstract: In this paper, we consider sparse networks consisting of a finite number of non-overlapping communities, i.e. disjoint clusters, so that there is higher density within clusters than across clusters. Both the intra- and inter-cluster edge densities vanish when the size of the graph grows large, making the cluster reconstruction problem nosier and hence difficult to solve. We are interested in scenarios where the network size is very large, so that the adjacency matrix of the graph is hard to manipulate and store. The data stream model in which columns of the adjacency matrix are revealed sequentially constitutes a natural framework in this setting. For this model, we develop two novel clustering algorithms that extract the clusters asymptotically accurately. The first algorithm is {\it offline}, as it needs to store and keep the assignments of nodes to clusters, and requires a memory that scales linearly with the network size. The second algorithm is {\it online}, as it may classify a node when the corresponding column is revealed and then discard this information. This algorithm requires a memory growing sub-linearly with the network size. To construct these efficient streaming memory-limited clustering algorithms, we first address the problem of clustering with partial information, where only a small proportion of the columns of the adjacency matrix is observed and develop, for this setting, a new spectral algorithm which is of independent interest.


title: Storms in Mobile Networks

abstract: Mobile networks are vulnerable to signalling attacks and storms that are caused by traffic patterns that overload the control plane, and differ from distributed denial of service (DDoS) attacks in the Internet since they directly attack the control plane, and also reserve wireless bandwidth without actually using it. Such attacks can result from malware and mobile botnets, as well as from poorly designed applications, and can cause service outages in 3G and 4G networks which have been experienced by mobile operators. Since the radio resource control (RRC) protocol in 3G and 4G networks is particularly susceptible to such attacks, we analyze their effect with a mathematical model that helps to predict the congestion that is caused by an attack. A detailed simulation model of a mobile network is used to better understand the temporal dynamics of user behavior and signalling in the network and to show how RRC based signalling attacks and storms cause significant problems in the control plane and the user plane of the network. Our analysis also serves to identify how storms can be detected, and to propose how system parameters can be chosen to mitigate their effect.


title: The artificial retina for track reconstruction at the LHC crossing rate

abstract: We present the results of an R&D study for a specialized processor capable of precisely reconstructing events with hundreds of charged-particle tracks in pixel and silicon strip detectors at $40\,\rm MHz$, thus suitable for processing LHC events at the full crossing frequency. For this purpose we design and test a massively parallel pattern-recognition algorithm, inspired to the current understanding of the mechanisms adopted by the primary visual cortex of mammals in the early stages of visual-information processing. The detailed geometry and charged-particle's activity of a large tracking detector are simulated and used to assess the performance of the artificial retina algorithm. We find that high-quality tracking in large detectors is possible with sub-microsecond latencies when the algorithm is implemented in modern, high-speed, high-bandwidth FPGA devices.


title: Multiplexity and multireciprocity in directed multiplexes

abstract: Real-world multi-layer networks feature nontrivial dependencies among links of different layers. Here we argue that, if links are directed, dependencies are twofold. Besides the ordinary tendency of links of different layers to align as the result of `multiplexity', there is also a tendency to anti-align as the result of what we call `multireciprocity', i.e. the fact that links in one layer can be reciprocated by \emph{opposite} links in a different layer. Multireciprocity generalizes the scalar definition of single-layer reciprocity to that of a square matrix involving all pairs of layers. We introduce multiplexity and multireciprocity matrices for both binary and weighted multiplexes and validate their statistical significance against maximum-entropy null models that filter out the effects of node heterogeneity. We then perform a detailed empirical analysis of the World Trade Multiplex (WTM), representing the import-export relationships between world countries in different commodities. We show that the WTM exhibits strong multiplexity and multireciprocity, an effect which is however largely encoded into the degree or strength sequences of individual layers. The residual effects are still significant and allow to classify pairs of commodities according to their tendency to be traded together in the same direction and/or in opposite ones. We also find that the multireciprocity of the WTM is significantly lower than the usual reciprocity measured on the aggregate network. Moreover, layers with low (high) internal reciprocity are embedded within sets of layers with comparably low (high) mutual multireciprocity. This suggests that, in the WTM, reciprocity is inherent to groups of related commodities rather than to individual commodities. We discuss the implications for international trade research focusing on product taxonomies, the product space, and fitness/complexity metrics.


title: Central exclusive production in the ALICE experiment at the LHC

abstract: The ALICE experiment at the Large Hadron Collider (LHC) at CERN consists of a central barrel, a muon spectrometer and additional detectors for trigger and event classification purposes. The low transverse momentum threshold of the central barrel gives ALICE a unique opportunity to study the low mass sector of central exclusive production at the LHC.


title: Kullback-Leibler divergence for interacting multiple model estimation with random matrices

abstract: This paper studies the problem of interacting multiple model (IMM) estimation for jump Markov linear systems with unknown measurement noise covariance. The system state and the unknown covariance are jointly estimated in the framework of Bayesian estimation, where the unknown covariance is modeled as a random matrix according to an inverse-Wishart distribution. For the IMM estimation with random matrices, one difficulty encountered is the combination of a set of weighted inverse-Wishart distributions. Instead of using the moment matching approach, this difficulty is overcome by minimizing the weighted Kullback-Leibler divergence for inverse-Wishart distributions. It is shown that a closed form solution can be derived for the optimization problem and the resulting solution coincides with an inverse-Wishart distribution. Simulation results show that the proposed filter performs better than the previous work using the moment matching approach.


title: Controlling false discoveries in high-dimensional situations: Boosting with stability selection

abstract: Modern biotechnologies often result in high-dimensional data sets with much more variables than observations (n $\ll$ p). These data sets pose new challenges to statistical analysis: Variable selection becomes one of the most important tasks in this setting. We assess the recently proposed flexible framework for variable selection called stability selection. By the use of resampling procedures, stability selection adds a finite sample error control to high-dimensional variable selection procedures such as Lasso or boosting. We consider the combination of boosting and stability selection and present results from a detailed simulation study that provides insights into the usefulness of this combination. Limitations are discussed and guidance on the specification and tuning of stability selection is given. The interpretation of the used error bounds is elaborated and insights for practical data analysis are given. The results will be used to detect differentially expressed phenotype measurements in patients with autism spectrum disorders. All methods are implemented in the freely available R package stabs.


title: Extremal Radii, Diameter, and Minimum Width in Generalized Minkowski Spaces

abstract: We discuss the notions of circumradius, inradius, diameter, and minimum width in generalized Minkowski spaces (that is, with respect to gauges), i.e., we measure the "size" of a given convex set in a finite-dimensional real vector space with respect to another convex set. This is done via formulating some kind of containment problem incorporating homothetic bodies of the latter set or strips bounded by parallel supporting hyperplanes thereof. The paper can be seen as a theoretical starting point for studying metrical problems of sets in generalized Minkowski spaces.


title: Cosmic-ray driven winds

abstract: The theory of Galactic Winds, driven by the cosmic-ray pressure gradient, is reviewed both on the magnetohydrodynamic and on the kinetic level. In this picture the magnetic field of the Galaxy above the dense gas disk is assumed to have a flux tube geometry, the flux tubes rising locally perpendicular out of the disk to become radially directed at large distances, with the cosmic-ray sources located deep within the Galactic disk. At least above the gas disk, the magnetic fluctuations which resonantly scatter the cosmic rays are selfconsistently excited as Alf{`e}n waves by the escaping cosmic rays. The fluctuation amplitudes remain finite through nonlinear wave dissipation. The spatially increasing speed of the resulting outflow results in a diffusion-convection boundary whose position depends on particle momentum. It replaces the escape boundary of static diffusion models. New effects like overall Galactic mass and angular momentum loss as well as gas heating beyond the disk appear. Also particle re-acceleration in the distant wind halo suggests itself. The resulting magnetohydrodynamic flow properties and the cosmic-ray transport properties are compared with observations. On the whole they show remarkable agreement. General limitations and generalisations of the basic model arise due to the expected simultaneous infall of matter from the environment of the Galaxy. On an intergalactic scale the combined winds from the Local Group galaxies should form a ``Local Group Bubble``. Its properties remain to be studied in detail.


title: The High Time Resolution Universe Survey - XI. Discovery of five recycled pulsars and the optical detectability of survey white dwarf companions

abstract: We present the discovery of a further five recycled pulsar systems in the mid-Galactic latitude portion of the High Time Resolution Universe (HTRU) Survey. The pulsars have rotational periods ranging from 2 ms to 66 ms, and four are in binary systems with orbital periods between 10.8 hours and 9.0 days. Three of these binary systems are particularly interesting; PSR J1227-6208 has a pulse period of 34.5 ms and the highest mass function of all pulsars with near-circular orbits. The circular orbit suggests that the companion is not another neutron star, so future timing experiments may reveal one of the heaviest white dwarfs ever found ($>$ 1.3 M$_\odot$). Timing observations of PSR J1431$-$4715 indicate that it is eclipsed by its companion which has a mass indicating it belongs to the redback class of eclipsing millisecond pulsars. PSR J1653-2054 has a companion with a minimum mass of only $0.08$ M$_\odot$, placing it among the class of pulsars with low-mass companions. Unlike the majority of such systems, however, no evidence of eclipses is seen at 1.4 GHz.


title: The 3D MHD code GOEMHD3 for large-Reynolds-number astrophysical plasmas

abstract: The numerical simulation of turbulence and flows in almost ideal, large-Reynolds-number astrophysical plasmas motivates the implementation of almost conservative MHD computer codes. They should efficiently calculate, use highly parallelized schemes scaling well with large numbers of CPU cores, allows to obtain a high grid resolution over large simulation domains and which can easily be adapted to new computer architectures as well as to new initial and boundary conditions, allow modular extensions. The new massively parallel simulation code GOEMHD3 enables efficient and fast simulations of almost ideal, large-Reynolds-number astrophysical plasma flows, well resolved and on huge grids covering large domains. Its abilities are validated by major tests of ideal and weakly dissipative plasma phenomena. The high resolution ($2048^3$ grid points) simulation of a large part of the solar corona above an observed active region proved the excellent parallel scalability of the code using more than 30.000 processor cores.


title: On the sum of squared logarithms inequality and related inequalities

abstract: We consider the sum of squared logarithms inequality and investigate possible connections with the theory of majorization. We also discuss alternative sufficient conditions on two sets of vectors $a,b\in\mathbb{R}_+^n$ so that $\sum_{i=1}^n(\log a_i)^2\ \leq\ \sum_{i=1}^n(\log b_i)^2\,.\notag $ Generalizations of some inequalities from information theory are obtained, including a generalized information inequality and a generalized log sum inequality, which states for $a,b\in\mathbb{R}_+^n$ and $k_1,...,k_n\in [0,\infty)$: $ \sum_{i=1}^na_i\,\log\prod_{s=1}^m(\frac{a_i}{b_i} + k_s)\ \geq\ \log\prod_{s=1}^m(1+k_s)\,.\notag $


title: Fragmentation of fractal random structures

abstract: We analyze the fragmentation behavior of random clusters on the lattice under a process where bonds between neighboring sites are successively broken. Modeling such structures by configurations of a generalized Potts or random-cluster model allows us to discuss a wide range of systems with fractal properties including trees as well as dense clusters. We present exact results for the densities of fragmenting edges and the distribution of fragment sizes for critical clusters in two dimensions. Dynamical fragmentation with a size cutoff leads to broad distributions of fragment sizes. The resulting power laws are shown to encode characteristic fingerprints of the fragmented objects.


title: Monitoring Count Time Series in R: Aberration Detection in Public Health Surveillance

abstract: Public health surveillance aims at lessening disease burden, e.g., in case of infectious diseases by timely recognizing emerging outbreaks. Seen from a statistical perspective, this implies the use of appropriate methods for monitoring time series of aggregated case reports. This paper presents the tools for such automatic aberration detection offered by the R package surveillance. We introduce the functionality for the visualization, modelling and monitoring of surveillance time series. With respect to modelling we focus on univariate time series modelling based on generalized linear models (GLMs), multivariate GLMs, generalized additive models and generalized additive models for location, shape and scale. This ranges from illustrating implementational improvements and extensions of the well-known Farrington algorithm, e.g, by spline-modelling or by treating it in a Bayesian context. Furthermore, we look at categorical time series and address overdispersion using beta-binomial or Dirichlet-Multinomial modelling. With respect to monitoring we consider detectors based on either a Shewhart-like single timepoint comparison between the observed count and the predictive distribution or by likelihood-ratio based cumulative sum methods. Finally, we illustrate how surveillance can support aberration detection in practice by integrating it into the monitoring workflow of a public health institution. Altogether, the present article shows how well surveillance can support automatic aberration detection in a public health surveillance context.


title: Multiscale Universal Interface: A Concurrent Framework for Coupling Heterogeneous Solvers

abstract: Concurrently coupled numerical simulations using heterogeneous solvers are powerful tools for modeling multiscale phenomena. However, major modifications to existing codes are often required to enable such simulations, posing significant difficulties in practice. In this paper we present a C++ library, i.e. the Multiscale Universal Interface (MUI), which is capable of facilitating the coupling effort for a wide range of multiscale simulations. The library adopts a header-only form with minimal external dependency and hence can be easily dropped into existing codes. A data sampler concept is introduced, combined with a hybrid dynamic/static typing mechanism, to create an easily customizable framework for solver-independent data interpretation. The library integrates MPI MPMD support and an asynchronous communication protocol to handle inter-solver information exchange irrespective of the solvers' own MPI awareness. Template metaprogramming is heavily employed to simultaneously improve runtime performance and code flexibility. We validated the library by solving three different multiscale problems, which also serve to demonstrate the flexibility of the framework in handling heterogeneous models and solvers. In the first example, a Couette flow was simulated using two concurrently coupled Smoothed Particle Hydrodynamics (SPH) simulations of different spatial resolutions. In the second example, we coupled the deterministic SPH method with the stochastic Dissipative Particle Dynamics (DPD) method to study the effect of surface grafting on the hydrodynamics properties on the surface. In the third example, we consider conjugate heat transfer between a solid domain and a fluid domain by coupling the particle-based energy-conserving DPD (eDPD) method with the Finite Element Method (FEM).

categories: physics.comp-ph cs.CE cs.DC physics.flu-dyn


title: Far-infrared line spectra of Seyfert galaxies from the Herschel-PACS Spectrometer

abstract: We present spectroscopic observations of FIR fine-structure lines of 26 Seyfert galaxies obtained with the Herschel-PACS spectrometer. These observations are complemented by spectroscopy with Spitzer-IRS and Herschel-SPIRE. The ratios of the OIII, NII, SIII and NeV lines have been used to determine electron densities in the ionised gas regions. The CI lines, observed with SPIRE, have been used to measure the densities in the neutral gas, while the OI lines provide a measure of the gas temperature, at densities below 10000 cm-3. Using the OI145/63um and SIII33/18um line ratios we find an anti-correlation of the temperature with the gas density. Using various fine-structure line ratios, we find that density stratification is common in these active galaxies. On average, the electron densities increase with the ionisation potential of the ions producing the NII, SIII and NeV emission. The infrared emission lines arise partly in the Narrow Line Region (NLR) photoionised by the AGN central engine, partly in HII regions photo ionised by hot stars and partly in neutral gas in photo-dissociated regions (PDRs). We attempt to separate the contributions to the line emission produced in these different regions by comparing our emission line ratios to empirical and theoretical values. In particular, we tried to separate the contribution of AGN and star formation by using a combination of Spitzer and Herschel lines, and we found that, besides the well known mid-IR line ratios, the mixed mid-IR/far-IR line ratio of OIII88um/OIV26um can reliably discriminate the two emission regimes, while the far-IR line ratio of CII157um/OI63um is only able to mildly separate the two regimes. By comparing the observed CII157um/NII205um ratio with photoionisation models, we also found that most of the CII emission in the galaxies we examined is due to PDRs.


title: Well-posedness for dislocation based gradient visco-plasticity with isotropic hardening

abstract: In this work we establish the well-posedness for infinitesimal dislocation based gradient viscoplasticity with isotropic hardening for general gradient monotone plastic flows. We assume an additive split of the displacement gradient into non-symmetric elastic distortion and non-symmetric plastic distortion. The thermodynamic potential is augmented with a term taking the dislocation density tensor $\operatorname{Curl} p$ into account. The constitutive equations in the models we study are assumed to be of self-controlling type. Based on the generalized version of Korn's inequality for incompatible tensor fields (the non-symmetric plastic distortion) due to Neff/Pauly/Witsch the existence of solutions of quasi-static initial-boundary value problems under consideration is shown using a time-discretization technique and a monotone operator method.


title: CW posets after the Poincare Conjecture

abstract: Anders Bjorner characterized which finite graded partially ordered sets arise as the posets of closure relations on cells of a finite, regular CW complex. His characterization of these "CW posets" required each open interval $(\hat{0},u)$ to have order complex homeomorphic to a sphere of dimension $rk(u)-2$. Work of Danaraj and Klee showed that sufficient conditions were for the poset to be thin and shellable. The proof of the Poincare Conjecture enables the requirement of shellability to be replaced by the homotopy Cohen-Macaulay property. This expands the range of tools that may be used to prove a poset is a CW poset.


title: Edge Detection based on Kernel Density Estimation

abstract: Edges of an image are considered a crucial type of information. These can be extracted by applying edge detectors with different methodology. Edge detection is a vital step in computer vision tasks, because it is an essential issue for pattern recognition and visual interpretation. In this paper, we propose a new method for edge detection in images, based on the estimation by kernel of the probability density function. In our algorithm, pixels in the image with minimum value of density function are labeled as edges. The boundary between two homogeneous regions is defined in two domains: the spatial/lattice domain and the range/color domain. Extensive experimental evaluations proved that our edge detection method is significantly a competitive algorithm.


title: The Milnor fiber of the singularity f(x,y) + zg(x,y) = 0

abstract: We give a description of the Milnor fiber and the monodromy of a singularity of the form f+zg = 0 where f and g define plane curves and have no common components. The description depends only on the topological type of the two plane curve germs defined by f and g. In particular, this gives a description of the boundary of the Milnor fiber. As a corollary, we give a simple formula for the monodromy zeta function and the Euler characteristic of the fiber.


title: Modifications of Tutte-Grothendieck invariants and Tutte polynomials

abstract: We transform Tutte-Grothedieck invariants thus also Tutte polynomials on matroids so that the contraction-deletion rule for loops (isthmuses) coincides with the general case.


title: Probing the Solar Wind Acceleration Region with the Sun--Grazing Comet C/2002 S2

abstract: Comet C/2002 S2, a member of the Kreutz family of Sungrazing comets, was discovered in white light images of the SOHO/LASCO coronagraph on 2002 September 18 and observed in \hi\, \lya\, emission by the SOHO/UVCS instrument at four different heights as it approached the Sun. The \hi\, \lya\, line profiles detected by UVCS are analyzed to determine the spectral parameters: line intensity, width and Doppler shift with respect to the coronal background. Two dimensional comet images of these parameters are reconstructed at the different heights. A novel aspect of the observations of this sungrazing comet data is that, whereas the emission from the most of the tail is blue--shifted, that along one edge of the tail is red--shifted. We attribute these shifts to a combination of solar wind speed and interaction with the magnetic field. In order to use the comet to probe the density, temperature and speed of the corona and solar wind through which it passes, as well as to determine the outgassing rate of the comet, we develop a Monte Carlo simulation of the \hi\, \lya\, emission of a comet moving through a coronal plasma. From the outgassing rate, we estimate a nucleus diameter of about 9 meters. This rate steadily increases as the comet approaches the Sun while the optical brightness decreases by more than a factor of ten and suddenly recovers. This indicates that the optical brightness is determined by the lifetimes of the grains, sodium atoms and molecules produced by the comet.


title: The a-function for gauge theories

abstract: The a-function is a proposed quantity defined for quantum field theories which has a monotonic behaviour along renormalisation group flows, being related to the beta-functions via a gradient flow equation involving a positive definite metric. We construct the a-function at four loop order for a general gauge theory with fermions and scalars, using only one and two loop beta-functions; we are then able to provide a stringent consistency check on the general three-loop gauge beta-function. In the case of an N=1 supersymmetric gauge theory, we present a general condition on the chiral field anomalous dimension which guarantees an exact all-orders expression for the a-function; and we verify this up to fifth order (corresponding to the three-loop anomalous dimension).


title: Analytic Dependence is an Unnecessary Requirement in Renormalization of Locally Covariant QFT

abstract: Finite renormalization freedom in locally covariant quantum field theories on curved spacetime is known to be tightly constrained, under certain standard hypotheses, to the same terms as in flat spacetime up to finitely many curvature dependent terms. These hypotheses include, in particular, locality, covariance, scaling, microlocal regularity and continuous and analytic dependence on the metric and coupling parameters. The analytic dependence hypothesis is somewhat unnatural, because it requires that locally covariant observables (which are simultaneously defined on all spacetimes) depend continuously on an arbitrary metric, with the dependence strengthened to analytic on analytic metrics. Moreover the fact that analytic metrics are globally rigid makes the implementation of this requirement at the level of local $*$-algebras of observables rather technically cumbersome. We show that the conditions of locality, covariance, scaling and a naturally strengthened microlocal spectral condition, are actually sufficient to constrain the allowed finite renormalizations equally strongly, thus eliminating both the continuity and the somewhat unnatural analyticity hypotheses. The key step in the proof uses the Peetre-Slov\'ak theorem on the characterization of (in general non-linear) differential operators by their locality and regularity properties.


title: Convex polygons in geometric triangulations

abstract: We show that the maximum number of convex polygons in a triangulation of $n$ points in the plane is $O(1.5029^n)$. This improves an earlier bound of $O(1.6181^n)$ established by van Kreveld, L\"offler, and Pach (2012) and almost matches the current best lower bound of $\Omega(1.5028^n)$ due to the same authors. Given a planar straight-line graph $G$ with $n$ vertices, we show how to compute efficiently the number of convex polygons in $G$.


abstract: A function of positive type can be defined as a positive functional on a convolution algebra of a locally compact group. In the case where the group is abelian, by Bochner's theorem a function of positive type is, up to normalization, the Fourier transform of a probability measure. Therefore, considering the group of translations on phase space, a suitably normalized phase-space function of positive type can be regarded as a realization of a classical state. Thus, it may be called a function of classical positive type. Replacing the ordinary convolution on phase space with the twisted convolution, one obtains a noncommutative algebra of functions whose positive functionals we may call functions of quantum positive type. In fact, by a quantum version of Bochner's theorem, a continuous function of quantum positive type is, up to normalization, the (symplectic) Fourier transform of a Wigner quasi-probability distribution; hence, it can be regarded as a phase-space realization of a quantum state. Playing with functions of positive type, classical and quantum, one is led in a natural way to consider a class of semigroups of operators, the classical-quantum semigroups. The physical meaning of these mathematical objects is unveiled via quantization, so obtaining a class of quantum dynamical semigroups that, borrowing terminology from quantum information science, may be called classical-noise semigroups.


title: Biclosed sets in real hyperplane arrangements

abstract: The set of chambers of a real hyperplane arrangement may be ordered by separation from some fixed chamber. When this poset is a lattice, Bjorner, Edelman, and Ziegler proved that the chambers are in natural bijection with the biconvex sets of the arrangement. Two families of examples of arrangements with a lattice of chambers are simplicial and supersolvable arrangements. For these arrangements, we prove that the chambers correspond to biclosed sets, a weakening of the biconvex property.


title: On the radio properties of the intermediate-mass black hole candidate ESO 243-49 HLX-1

abstract: We present follow-up radio observations of ESO 243-49 HLX-1 from 2012 using the Australia Telescope Compact Array (ATCA) and the Karl G. Jansky Very Large Array (VLA). We report the detection of radio emission at the location of HLX-1 during its hard X-ray state using the ATCA. Assuming that the `Fundamental Plane' of accreting black holes is applicable, we provide an independent estimate of the black hole mass of $M_{\rm{BH}}\leq2.8^{+7.5}_{-2.1} \times 10^{6}$ M$_{\odot}$ at 90% confidence. However, we argue that the detected radio emission is likely to be Doppler-boosted and our mass estimate is an upper limit. We discuss other possible origins of the radio emission such as being due to a radio nebula, star formation, or later interaction of the flares with the large-scale environment. None of these were found adequate. The VLA observations were carried out during the X-ray outburst. However, no new radio flare was detected, possibly due to a sparse time sampling. The deepest, combined VLA data suggests a variable radio source and we briefly discuss the properties of the previously detected flares and compare them with microquasars and active galactic nuclei.


title: An Open Distributed Architecture for Flexible Hybrid Assembly Systems: A Model Driven Engineering Approach

abstract: Assembly systems constitute one of the most important fields in today industry. In this paper we propose an open distributed architecture for the engineering of evolvable flexible hybrid assembly systems. The proposed architecture is based on the model driven development paradigm. Models are used to represent structure and behavior and a domain specific engineering tool is defined to facilitate the assembly system engineer in the engineering process of the assembly system. Specific meta models are defined to capture domain knowledge to guide the engineer in the construction of the models required to construct the assembly system. This work is a specialization of our previous work that defined a SOA based framework for embedded industrial automation systems. It adapts and extends, in the assembly systems domain, the 3+1 SysML-view model architecture defined for the engineering of mechatronics Manufacturing systems. The proposed architecture can be used to develop a framework for evolvable flexible and reconfigurable assembly systems that would exploit the benefits the Cyber Physical paradigm utilizing web technologies, the IoT, the Cloud computing and Big Data.


title: Adaptive Error Covariances Estimation Methods for Ensemble Kalman Filters

abstract: This paper presents a computationally fast algorithm for estimating, both, the system and observation noise covariances of nonlinear dynamics, that can be used in an ensemble Kalman filtering framework. The new method is a modification of Belanger's recursive method, to avoid an expensive computational cost in inverting error covariance matrices of product of innovation processes of different lags when the number of observations becomes large. When we use only product of innovation processes up to one-lag, the computational cost is indeed comparable to a recently proposed method by Berry-Sauer's. However, our method is more flexible since it allows for using information from product of innovation processes of more than one-lag. Extensive numerical comparisons between the proposed method and both the original Belanger's and Berry-Sauer's schemes are shown in various examples, ranging from low-dimensional linear and nonlinear systems of SDE's and 40-dimensional stochastically forced Lorenz-96 model. Our numerical results suggest that the proposed scheme is as accurate as the original Belanger's scheme on low-dimensional problems and has a wider range of more accurate estimates compared to Berry-Sauer's method on L-96 example.


title: Analogies between nuclear physics and Dark Matter

abstract: A fermionic description of dark matter using analogies with nuclear physics is developed. At tree level, scalar and vector processes are considered and the two-body potential are explicitly calculated using the Breit approximation. We show that the total cross sections in both cases exhibit Sommerfeld enhancement.


title: Entanglement swapping between discrete and continuous variables

abstract: We experimentally realize "hybrid" entanglement swapping between discrete-variable (DV) and continuous-variable (CV) optical systems. DV two-mode entanglement as obtainable from a single photon split at a beam splitter is robustly transferred by means of efficient CV entanglement and operations, using sources of squeezed light and homodyne detections. The DV entanglement after the swapping is verified without post-selection by the logarithmic negativity of up to 0.28$\pm$0.01. Furthermore, our analysis shows that the optimally transferred state can be post-selected into a highly entangled state that violates a Clauser-Horne-Shimony-Holt inequality by more than four standard deviations, and thus it may serve as resources for quantum teleportation and quantum cryptography.


title: Consistency of equations of motion in conformal frames

abstract: Four dimensional scalar-tensor theory is considered within two conformal frames, the Jordan frame (JF) and the Einstein frame (EF). The actions for the theory are equivalent and equations of motion can be obtained from each action. It is found that the JF equations of motion, expressed in terms of EF variables, translate directly into and agree with the EF equations of motion obtained from the EF action, provided that certain simple consistency conditions are satisfied, which is always the case. The implication is that a solution set obtained in one conformal frame can be reliably translated into a solution set for the other frame, and therefore the two frames are, at least, mathematically equivalent.


title: A Generalization of an Integral Arising in the Theory of Distance Correlation

abstract: We generalize an integral which arises in several areas in probability and statistics and which is at the core of the field of distance correlation, a concept developed by Sz\'ekely, Rizzo and Bakirov (2007) to measure dependence between random variables. Let $m$ be a positive integer and let ${\cos_m}(u)$, $u \in \mathbb{R}$, be the truncated Maclaurin expansion of ${\cos}(u)$, where the expansion is truncated at the $m$th summand. For $t, x \in \mathbb{R}^d$, let $\langle t,x\rangle$ and $\|x\|$ denote the standard Euclidean inner product and norm, respectively. We establish the integral formula: For $\alpha \in \mathbb{C}$ and $x \in \mathbb{R}^d$, $\int_{{\mathbb{R}}^d} [\cos_m(\langle t,x\rangle) - \cos(\langle t,x\rangle)] \,{\rm d}t/{\|t\|^{d+\alpha}} = C(d,\alpha) \, \|x\|^{\alpha}$, with absolute convergence if and only if $2(m-1) < \Re(\alpha) < 2m$. Moreover, the constant $C(d,\alpha)$ does not depend on $m$.


title: Lattice Three-Dimensional Skyrmions Revisited

abstract: In the continuum a skyrmion is a topological nontrivial map between Riemannian manifolds, and a stationary point of a particular energy functional. This paper describes lattice analogues of the aforementioned skyrmions, namely a natural way of using the topological properties of the three-dimensional continuum Skyrme model to achieve topological stability on the lattice. In particular, using fixed point iterations, numerically exact lattice skyrmions are constructed; and their stability under small perturbations is explored by means of linear stability analysis. While stable branches of such solutions are identified, it is also shown that they possess a particularly delicate bifurcation structure, especially so in the vicinity of the continuum limit. The corresponding bifurcation diagram is elucidated and a prescription for selecting the branch asymptoting to the well-known continuum limit is given. Finally, the robustness of the spectrally stable solutions is corroborated by virtue of direct numerical simulations .


title: Computation of Gaussian orthant probabilities in high dimension

abstract: We study the computation of Gaussian orthant probabilities, i.e. the probability that a Gaussian falls inside a quadrant. The Geweke-Hajivassiliou-Keane (GHK) algorithm [Genz, 1992; Geweke, 1991; Hajivassiliou et al., 1996; Keane, 1993], is currently used for integrals of dimension greater than 10. In this paper we show that for Markovian covariances GHK can be interpreted as the estimator of the normalizing constant of a state space model using sequential importance sampling (SIS). We show for an AR(1) the variance of the GHK, properly normalized, diverges exponentially fast with the dimension. As an improvement we propose using a particle filter (PF). We then generalize this idea to arbitrary covariance matrices using Sequential Monte Carlo (SMC) with properly tailored MCMC moves. We show empirically that this can lead to drastic improvements on currently used algorithms. We also extend the framework to orthants of mixture of Gaussians (Student, Cauchy etc.), and to the simulation of truncated Gaussians.


title: Lax operator for Macdonald symmetric functions

abstract: Using the Lax operator formalism, we construct a family of pairwise commuting operators such that the Macdonald symmetric functions of infinitely many variables and of two parameters $q,t$ are their eigenfunctions. We express our operators in terms of the Hall-Littlewood symmetric functions of the same variables and of the parameter $t$ corresponding to the partitions with one part only. Our expression is based on the notion of Baker-Akhiezer function.


title: Rapid Skill Capture in a First-Person Shooter

abstract: Various aspects of computer game design, including adaptive elements of game levels, characteristics of 'bot' behavior, and player matching in multiplayer games, would ideally be sensitive to a player's skill level. Yet, while difficulty and player learning have been explored in the context of games, there has been little work analyzing skill per se, and how it pertains to a player's input. To this end, we present a data set of 476 game logs from over 40 players of a first-person shooter game (Red Eclipse) as a basis of a case study. We then analyze different metrics of skill and show that some of these can be predicted using only a few seconds of keyboard and mouse input. We argue that the techniques used here are useful for adapting games to match players' skill levels rapidly, perhaps more rapidly than solutions based on performance averaging such as TrueSkill.


title: Radiative properties of multi-carrier bound excitons in GaAs

abstract: Excitons in semiconductors can have multiple lifetimes due to spin dependent oscillator strengths and interference between different recombination pathways. In addition, strain and symmetry effects can further modify lifetimes via the removal of degeneracies. We present a convenient formalism for predicting the optical properties of ${k=0}$ excitons with an arbitrary number of charge carriers in different symmetry environments. Using this formalism, we predict three distinct lifetimes for the neutral acceptor bound exciton in GaAs, and confirm this prediction through polarization dependent and time-resolved photoluminescence experiments. We find the acceptor bound-exciton lifetimes to be ${T_o (1,3,3/4)}$ where ${T_o = (0.61 \pm 0.12) \text{ns}}$. Furthermore, we provide an estimate of the intra-level and inter-level exciton spin-relaxation rates.


title: Radial and orbital excitation energies of charmonium

abstract: The charmonium system has several excited states below the energy threshold for decay into $D$ and $\bar{D}$ mesons, which can in principle be studied accurately in lattice QCD. Studies that include many states in the spectrum have typically only been done at one value of the lattice spacing and with relatively heavy light quarks in the sea. Here we give preliminary results for radial and orbital excitation energies for charmonium from a calculation on 2+1+1 MILC configurations at multiple lattice spacings and including physical values for $u/d$ quark masses. We use the HISQ formulation for $c$ to obtain small discretisation errors and smeared operators to improve excited state overlap.


title: Election vs. Selection: Two Ways of Finding the Largest Node in a Graph

abstract: Finding the node with the largest label in a network, modeled as an undirected connected graph, is one of the fundamental problems in distributed computing. This is the way in which $\textit{leader election}$ is usually solved. We consider two distinct tasks in which the largest-labeled node is found deterministically. In $\textit{selection}$, this node must output 1 and all other nodes must output 0. In $\textit{election}$, the other nodes must additionally learn the largest label. Our aim is to compare the difficulty of these two tasks executed under stringent running time constraints. The measure of difficulty is the amount of information that nodes of the network must initially possess in order to solve the given task in an imposed amount of time. Following the standard framework of $\textit{algorithms with advice}$, this information (a single binary string) is provided to all nodes at the start by an oracle knowing the entire graph. The length of this string is called the $\textit{size of advice}$. Consider the class of $n$-node graphs with any diameter $diam \leq D$. If time is larger than $diam$, then both tasks can be solved without advice. For the task of $\textit{election}$, we show that if time is smaller than $diam$, then the optimal size of advice is $\Theta(\log n)$, and if time is exactly $diam$, then the optimal size of advice is $\Theta(\log D)$. For the task of $\textit{selection}$, the situation changes dramatically, even within the class of rings. Indeed, for the class of rings, we show that, if time is $O(diam^{\epsilon})$, for any $\epsilon <1$, then the optimal size of advice is $\Theta(\log D)$, and, if time is $\Theta(diam)$ (and at most $diam$) then this optimal size is $\Theta(\log \log D)$.


title: Twisted spectral triple for the Standard Model and spontaneous breaking of the Grand Symmetry

abstract: Grand symmetry models in noncommutative geometry have been introduced to explain how to generate minimally (i.e. without adding new fermions) an extra scalar field beyond the standard model, which both stabilizes the electroweak vacuum and makes the computation of the mass of the Higgs compatible with its experimental value. In this paper, we use Connes-Moscovici twisted spectral triples to cure a technical problem of the grand symmetry, that is the appearance together with the extra scalar field of unbounded vectorial terms. The twist makes these terms bounded and - thanks to a twisted version of the first-order condition that we introduce here - also permits to understand the breaking to the standard model as a dynamical process induced by the spectral action. This is a spontaneous breaking from a pre-geometric Pati-Salam model to the almost-commutative geometry of the standard model, with two Higgs-like fields: scalar and vector.


title: On the geometry behind a recurrent relation

abstract: We consider a certain linear recursive relation with integer parameters and study some of its algebraic and geometric properties, with the purpose of estimating the number of chains of valences in the Farey series.


title: Water and ice structure in the range 220 - 365K from radiation total scattering experiments

abstract: The past decade or so has witnessed a large number of articles about water structure. The most incisive experiments involve radiation with a wavelength compatible with the observed inter-molecular separations found in water, of order $\sim 3$\AA, in other words mostly $<$1eV neutrons and $>$10keV x-rays. Because x-rays are scattered by electrons while neutrons are scattered by nuclei, the two probes give complementary information about the three site-site radial distribution functions for water, namely O-O, O-H and H-H. Here a version of Monte Carlo simulation called Empirical Potential Structure Refinement is used to devise an empirical intermolecular potential which attempts to drive the simulated radial distribution functions as close as possible to the data. New x-ray and neutron scattering data on water in the temperature range 280 - 365K are presented for the first time, alongside a new analysis of some much older neutron data on ice 1h at 220K. This temperature analysis, above and below the water freezing point of water, reveals some non-intuitive water properties in the liquid and solid states.


title: Fast cooling for a system of stochastic oscillators

abstract: We study feedback control of coupled nonlinear stochastic oscillators in a force field. We first consider the problem of asymptotically driving the system to a desired {\em steady state} corresponding to reduced thermal noise. Among the feedback controls achieving the desired asymptotic transfer, we find that the most efficient one {from an energy point of view} is characterized by {\em time-reversibility}. We also extend the theory of Schr\"{o}dinger bridges to this model, thereby steering the system in {\em finite time} and with minimum effort to a target steady-state distribution. The system can then be maintained in this state through the optimal steady-state feedback control. The solution, in the finite-horizon case, involves a space-time harmonic function $\varphi$, and $-\log\varphi$ plays the role of an artificial, time-varying potential in which the desired evolution occurs. This framework appears extremely general and flexible and can be viewed as a considerable generalization of existing active control strategies such as macromolecular cooling. In the case of a quadratic potential, the results assume a form particularly attractive from the algorithmic viewpoint as the optimal control can be computed via deterministic matricial differential equations. An example involving inertial particles illustrates both transient and steady state optimal feedback control.


title: Hybrid optical-electrical detection of donor electron spins with bound excitons in silicon

abstract: Electrical detection of spins is an essential tool in understanding the dynamics of spins in semiconductor devices, providing valuable insights for applications ranging from optoelectronics and spintronics to quantum information processing. For electron spins bound to shallow donors in silicon, bulk electrically-detected magnetic resonance has relied on coupling to spin readout partners such as paramagnetic defects or conduction electrons which fundamentally limits spin coherence times. Here we demonstrate electrical detection of phosphorus donor electron spin resonance by transport through a silicon device, using optically-driven donor-bound exciton transitions. We use this method to measure electron spin Rabi oscillations, and, by avoiding use of an ancillary spin for readout, we are able to obtain long intrinsic electron spin coherence times, limited only by the donor concentration. We go on to experimentally address critical issues for adopting this scheme for single spin measurement in silicon nanodevices, including the effects of strain, electric fields, and non-resonant excitation. This lays the foundations for realising a versatile readout method for single spin readout with relaxed magnetic field and temperature requirements compared with spin-dependent tunneling.


title: Making sense of the local Galactic escape speed estimates in direct dark matter searches

abstract: Direct detection (DD) of dark matter (DM) candidates in the $\lesssim$10 GeV mass range is very sensitive to the tail of their velocity distribution. The important quantity is the maximum WIMP speed in the observer's rest frame, i.e. in average the sum of the local Galactic escape speed $v_{\rm esc}$ and of the circular velocity of the Sun $v_c$. While the latter has been receiving continuous attention, the former is more difficult to constrain. The RAVE Collaboration has just released a new estimate of $v_{\rm esc}$ (Piffl {\em et al.}, 2014 --- P14) that supersedes the previous one (Smith {\em et al.}, 2007), which is of interest in the perspective of reducing the astrophysical uncertainties in DD. Nevertheless, these new estimates cannot be used blindly as they rely on assumptions in the dark halo modeling which induce tight correlations between the escape speed and other local astrophysical parameters. We make a self-consistent study of the implications of the RAVE results on DD assuming isotropic DM velocity distributions, both Maxwellian and ergodic. Taking as references the experimental sensitivities currently achieved by LUX, CRESST-II, and SuperCDMS, we show that: (i) the exclusion curves associated with the best-fit points of P14 may be more constraining by up to $\sim 40$% with respect to standard limits, because the underlying astrophysical correlations induce a larger local DM density; (ii) the corresponding relative uncertainties inferred in the low WIMP mass region may be moderate, down to 10-15% below 10 GeV. We finally discuss the level of consistency of these results with other independent astrophysical constraints. This analysis is complementary to others based on rotation curves.

categories: astro-ph.CO astro-ph.GA astro-ph.HE hep-ph


title: Why is the electrocaloric effect so small in ferroelectrics?

abstract: Ferroelectrics are attractive candidate materials for environmentally friendly solid state refrigeration free of greenhouse gases. Their thermal response upon variations of external electric fields is largest in the vicinity of their phase transitions, which may occur near room temperature. The magnitude of the effect, however, is too small for useful cooling applications even when they are driven close to dielectric breakdown. Insight from microscopic theory is therefore needed to characterize materials and provide guiding principles to search for new ones with enhanced electrocaloric performance. Here, we derive from well-known microscopic models of ferroelectricity meaningful figures of merit which provide insight into the relation between the strength of the effect and the characteristic interactions of ferroelectrics such as dipole forces. We find that the long range nature of these interactions results in a small effect. A strategy is proposed to make it larger by shortening the correlation lengths of fluctuations of polarization.


title: Variational Monte Carlo study of chiral spin liquid in the extended Heisenberg model on the Kagome lattice

abstract: We investigate the extended Heisenberg model on the Kagome lattice by using Gutzwiller projected fermionic states and the variational Monte Carlo technique. In particular, when both second- and third-neighbor super-exchanges are considered, we find that a gapped spin liquid described by non-trivial magnetic fluxes and long-range chiral-chiral correlations is energetically favored compared to the gapless U(1) Dirac state. Furthermore, the topological Chern number, obtained by integrating the Berry curvature, and the degeneracy of the ground state, by constructing linearly independent states, lead us to identify this flux state as the chiral spin liquid with $C=1/2$ fractionalized Chern number.


title: Discrete-Time Models Resulting From Dynamic Continuous-Time Perturbations In Phase-Amplitude Modulation-Demodulation Schemes

abstract: We consider discrete-time (DT) systems S in which a DT input is first tranformed to a continuous-time (CT) format by phase-amplitude modulation, then modified by a non-linear CT dynamical transformation F, and finally converted back to DT output using an ideal de-modulation scheme. Assuming that F belongs to a special class of CT Volterra series models with fixed degree and memory depth, we provide a complete characterization of S as a series connection of a DT Volterra series model of fixed degree and memory depth, and an LTI system with special properties. The result suggests a new, non-obvious, analytically motivated structure of digital compensation of analog nonlinear distortions (for example, those caused by power amplifiers) in digital communication systems. Results from a MATLAB simulation are used to demonstrate effectiveness of the new compensation scheme, as compared to the standard Volterra series approach.


title: Multiple Comparisons using Composite Likelihood in Clustered Data

abstract: We study the problem of multiple hypothesis testing for multidimensional data when inter-correlations are present. The problem of multiple comparisons is common in many applications. When the data is multivariate and correlated, existing multiple comparisons procedures based on maximum likelihood estimation could be prohibitively computationally intensive. We propose to construct multiple comparisons procedures based on composite likelihood statistics. We focus on data arising in three ubiquitous cases: multivariate Gaussian, probit, and quadratic exponential models. To help practitioners assess the quality of our proposed methods, we assess their empirical performance via Monte Carlo simulations. It is shown that composite likelihood based procedures maintain good control of the familywise type I error rate in the presence of intra-cluster correlation, whereas ignoring the correlation leads to erratic performance. Using data arising from a diabetic nephropathy study, we show how our composite likelihood approach makes an otherwise intractable analysis possible.


title: Anomalous CDW ground state in Cu$_2$Se: a wave-like fluctuation of $\it{dc}$ I-V curve near 50 K

abstract: A charge density wave (CDW) ground state is observed in polycrystalline Cu$_2$Se below 125 K, which corresponds to an energy gap of 40.9 meV and an electron-phonon coupling constant of 0.6. Due to the polycrystalline structure, the Peierls transition process has been expanded to a wide temperature range from 90 to 160 K. The Hall carrier concentration shows a continuous decrease from 2.1$\times$10$^{20}$ to 1.6$\times$10$^{20}$ cm$^{-3}$ in the temperature range from 160 K to 80 K, while almost unchanged above 160 K and below 90 K. After entering the CDW ground state, a wave-like fluctuation was observed in the I-V curve near 50 K, which exhibits as a periodic negative differential resistivity in an applied electric field due to the current. We also investigated the doping effect of Zn, Ni, and Te on the CDW ground state. Both Zn and Ni doped Cu$_2$Se show a CDW character with increased energy gap and electron-phonon coupling constant, but no notable Peierls transition was observed in Te doped Cu$_2$Se. Similar wave-like I-V curve was also seen in Cu$_{1.98}$Zn$_{0.02}$Se near 40 K. The regular fluctuation in $\it{dc}$ I-V curve was not magnetic field sensitive, but temperature and sample size sensitive.


title: Four loop anomalous dimension of the third and fourth moments of the non-singlet twist-2 operator in QCD

abstract: We present the result of a calculation for the third and fourth moments of the non-singlet four-loop anomalous dimension of Wilson twist-2 operators in QCD with full color and flavour structures. We discuss also a general expressions for some contributions to the full four-loop anomalous dimension obtained with the help of the method, based on LLL-algorithm, which was proposed earlier by us for the reconstruction of a general form of the anomalous dimension from a fixed values.


title: Strong C+ emission in galaxies at z~1-2: Evidence for cold flow accretion powered star formation in the early Universe

abstract: We have recently detected the [CII] 157.7 micron line in eight star forming galaxies at redshifts 1 to 2 using the redshift(z) Early Universe Spectrometer (ZEUS). Our sample targets star formation dominant sources detected in PAH emission. This represents a significant addition to [CII] observations during the epoch of peak star formation. We have augmented this survey with observations of the [OI] 63 micron line and far infrared photometry from the PACS and SPIRE Herschel instruments as well as Spitzer IRS spectra from the literature showing PAH features. Our sources exhibit above average gas heating efficiency, many with both [OI]/FIR and [CII]/FIR ~1% or more. The relatively strong [CII] emission is consistent with our sources being dominated by star formation powered PDRs, extending to kpc scales. We suggest that the star formation mode in these systems follows a Schmidt-Kennicutt law similar to local systems, but at a much higher rate due to molecular gas surface densities 10 to 100 times that of local star forming systems. The source of the high molecular gas surface densities may be the infall of neutral gas from the cosmic web. In addition to the high [CII]/FIR values, we also find high [CII]/PAH ratios and, in at least one source, a cool dust temperature. This source, SWIRE 4-5, bears a resemblance in these diagnostics to shocked regions of Stephan's Quintet, suggesting that another mode of [CII] excitation in addition to normal photoelectric heating may be contributing to the observed [CII] line.


title: Parabolic theory as a high-dimensional limit of elliptic theory

abstract: The aim of this article is to show how certain parabolic theorems follow from their elliptic counterparts. This technique is demonstrated through new proofs of five important theorems in parabolic unique continuation and the regularity theory of parabolic equations and geometric flows. Specifically, we give new proofs of an $L^2$ Carleman estimate for the heat operator, and the monotonicity formulas for the frequency function associated to the heat operator, the two-phase free boundary problem, the flow of harmonic maps, and the mean curvature flow. The proofs rely only on the underlying elliptic theorems and limiting procedures belonging essentially to probability theory. In particular, each parabolic theorem is proved by taking a high-dimensional limit of the related elliptic result.


title: A game with divisors and absolute differences of exponents

abstract: In this work, we discuss a number game that develops in a manner similar to that on which Gilbreath's conjecture on iterated absolute differences between consecutive primes is formulated. In our case the action occurs at the exponent level and there, the evolution is reminiscent of that in a final Ducci game. We present features of the whole field of the game created by the successive generations, prove an analog of Gilbreath's conjecture, and raise some open questions.


title: Discrimination of dark matter models in future experiments

abstract: Phenomenological aspects of simple dark matter models are studied. We discuss ways to discriminate the dark matter models in future experiments. We find that the measurements of the branching fraction of the Higgs boson into two photons and the electric dipole moment of the electron as well as the direct detection experiments are quite useful in discriminating particle models of dark matter. We also discuss the prospects of finding new particles in dark sector at the LHC/ILC.


title: Hard probes and the event generator EPOS

abstract: After a short presentation of the event generator EPOS, we discuss the production of heavy quarks and prompt photons which has been recently implemented. Whereas we have satisfying results for the charm, work on photons is still in progress.


title: Entanglement-enhanced time-continuous quantum control in optomechanics

abstract: The cavity-optomechanical radiation pressure interaction provides the means to create entanglement between a mechanical oscillator and an electromagnetic field interacting with it. Here we show how we can utilize this entanglement within the framework of time-continuous quantum control, in order to engineer the quantum state of the mechanical system. Specifically, we analyze how to prepare a low-entropy mechanical state by (measurement-based) feedback cooling operated in the blue detuned regime, the creation of bipartite mechanical entanglement via time-continuous entanglement swapping, and preparation of a squeezed mechanical state by time-continuous teleportation. The protocols presented here are feasible in optomechanical systems exhibiting a cooperativity larger than 1.


title: A Proof for Poisson Bracket in Non-commutative Algebra of Quantum Mechanics

abstract: The widely accepted approach to the foundation of quantum mechanics is that the Poisson bracket, governing the non-commutative algebra of operators, is taken as a postulate with no underlying physics. In this manuscript, it is shown that this postulation is in fact unnecessary and may be replaced by a few deeper concepts, which ultimately lead to the derivation of Poisson bracket. One would only need to use Fourier transform pairs and Kramers-Kronig identities in the complex domain. We present a definition of Hermitian time-operator and discuss some of its basic properties.


title: On Match Lengths, Zero Entropy and Large Deviations - with Application to Sliding Window Lempel-Ziv Algorithm

abstract: The Sliding Window Lempel-Ziv (SWLZ) algorithm that makes use of recurrence times and match lengths has been studied from various perspectives in information theory literature. In this paper, we undertake a finer study of these quantities under two different scenarios, i) \emph{zero entropy} sources that are characterized by strong long-term memory, and ii) the processes with weak memory as described through various mixing conditions. For zero entropy sources, a general statement on match length is obtained. It is used in the proof of almost sure optimality of Fixed Shift Variant of Lempel-Ziv (FSLZ) and SWLZ algorithms given in literature. Through an example of stationary and ergodic processes generated by an irrational rotation we establish that for a window of size $n_w$, a compression ratio given by $O(\frac{\log n_w}{{n_w}^a})$ where $a$ depends on $n_w$ and approaches 1 as $n_w \rightarrow \infty$, is obtained under the application of FSLZ and SWLZ algorithms. Also, we give a general expression for the compression ratio for a class of stationary and ergodic processes with zero entropy. Next, we extend the study of Ornstein and Weiss on the asymptotic behavior of the \emph{normalized} version of recurrence times and establish the \emph{large deviation property} (LDP) for a class of mixing processes. Also, an estimator of entropy based on recurrence times is proposed for which large deviation principle is proved for sources satisfying similar mixing conditions.


title: Synchronization by noise

abstract: We provide sufficient conditions for synchronization by noise, i.e. under these conditions we prove that weak random attractors for random dynamical systems consist of single random points. In the case of SDE with additive noise, these conditions are also essentially necessary. In addition, we provide sufficient conditions for the existence of a minimal weak point random attractor consisting of a single random point. As a result, synchronization by noise is proven for a large class of SDE with additive noise. In particular, we prove that the random attractor for an SDE with drift given by a (multidimensional) double-well potential and additive noise consists of a single random point. All examples treated in [Tearne, PTRF; 2008] are also included.


title: Consistent mass matrix of ten nodes tetrahedral element based on analytical integration

abstract: Currently, components of consistent mass matrix are computed using various numerical integration schemes, each one alters in number of integration (Gauss) points, requires different amount of computations and possess different level of accuracy. We discuss the closed-form mass matrix based on analytical integration. Curved-sided and straight-sided elements are considered. For a straight-sided element we derive an exact analytical easy to implement consistent mass matrix. For a curved-sided element an exact analytical mass matrix is derived, however it is rather lengthy, hence approximations are proposed. Three systematic approximations to the metric (jacobian determinant) are suggested; constant metric (CM), linearly varying metric (LM) and quadratic metric (QM). CM requires evaluation of the metric at the centroid, LM requires metric evaluations at the four corner nodes and QM uses metric values at all the ten nodes. Analytical integration together with approximated metric models yields closed-form semi-analytical mass matrices. The accuracy of the schemes is studied numerically using randomly generated coarse mesh. Our findings reveal significant superiority in accuracy and computations over equivalent schemes. An important implication of this study is that based on the results, it is superior to use our CM, LM and QM semi-analytical mass matrices over mass matrices based on numerical integration schemes which involve four, five and fifteen point Gauss quadrature. For a straight-sided element, CM, LM and QM admit an exact consistent mass matrix.


title: Correlated starting points for the functional renormalization group

abstract: We present a general frame to extend functional renormalization group (fRG) based computational schemes by using an exactly solvable interacting reference problem as starting point for the RG flow. The systematic expansion around this solution accounts for a non-perturbative inclusion of correlations. Introducing auxiliary fermionic fields by means of a Hubbard-Stratonovich transformation, we derive the flow equations for the auxiliary fields and determine the relation to the conventional weak-coupling truncation of the hierarchy of flow equations. As a specific example we consider the dynamical mean-field theory (DMFT) solution as reference system, and discuss the relation to the recently introduced DMF$^2$RG and the dual-fermion formalism.


title: Analysis and calibration techniques for superconducting resonators

abstract: A method is proposed and experimentally explored for in-situ calibration of complex transmission data for superconducting microwave resonators. This cryogenic calibration method accounts for the instrumental transmission response between the vector network analyzer reference plane and the device calibration plane. Once calibrated, the observed resonator response is analyzed in detail by two approaches. The first, a phenomenological model based on physically realizable rational functions, enables the extraction of multiple resonance frequencies and widths for coupled resonators without explicit specification of the circuit network. In the second, an ABCD-matrix representation for the distributed transmission line circuit is used to model the observed response from the characteristic impedance and propagation constant. When used in conjunction with electromagnetic simulations, the kinetic inductance fraction can be determined with this method with an accuracy of 2%. Datasets for superconducting microstrip and coplanar-waveguide resonator devices were investigated and a recovery within 1% of the observed complex transmission amplitude was achieved with both analysis approaches. The experimental configuration used in microwave characterization of the devices and self-consistent constraints for the electromagnetic constitutive relations for parameter extraction are also presented.


title: Tuning of the electronic and optical properties of single layer black phosphorus by strain

abstract: Using first principles calculations we showed that the electronic and optical properties of single layer black phosphorus (BP) depend strongly on the applied strain. Due to the strong anisotropic atomic structure of BP, its electronic conductivity and optical response are sensitive to the magnitude and the orientation of the applied strain. We found that the inclusion of many body effects is essential for the correct description of the electronic properties of monolayer BP; for example while the electronic gap of strainless BP is found to be 0.90 eV by using semilocal functionals, it becomes 2.31 eV when many-body effects are taken into account within the G0W0 scheme. Applied tensile strain was shown to significantly enhances electron transport along zigzag direction of BP. Furthermore, biaxial strain is able to tune the optical band gap of monolayer BP from 0.38 eV (at -8% strain) to 2.07 eV (at 5.5%). The exciton binding energy is also sensitive to the magnitude of the applied strain. It is found to be 0.40 eV for compressive biaxial strain of -8%, and it becomes 0.83 eV for tensile strain of 4%. Our calculations demonstrate that the optical response of BP can be significantly tuned using strain engineering which appears as a promising way to design novel photovoltaic devices that capture a broad range of solar spectrum.


title: Parametric resonance in spherical immersed elastic shells

abstract: We perform a stability analysis for a fluid-structure interaction problem in which a spherical elastic shell or membrane is immersed in a 3D viscous, incompressible fluid. The shell is an idealised structure having zero thickness, and has the same fluid lying both inside and outside. The problem is formulated mathematically using the immersed boundary framework in which Dirac delta functions are employed to capture the two-way interaction between fluid and immersed structure. The elastic structure is driven parametrically via a time-periodic modulation of the elastic membrane stiffness. We perform a Floquet stability analysis, considering the case of both a viscous and inviscid fluid, and demonstrate that the forced fluid-membrane system gives rise to parametric resonances in which the solution becomes unbounded even in the presence of viscosity. The analytical results are validated using numerical simulations with a 3D immersed boundary code for a range of wavenumbers and physical parameter values. Finally, potential applications to biological systems are discussed, with a particular focus on the human heart and investigating whether or not fluid-structure interaction-mediated instabilities could play a role in cardiac fluid dynamics.


title: Noncommutative minisuperspace, gravity-driven acceleration and kinetic inflation

abstract: In this paper, we introduce a noncommutative version of the Brans-Dicke (BD) theory and obtain the Hamiltonian equations of motion for a spatially flat Friedmann--Lema\^{\i}tre--Robertson--Walker universe filled with a perfect fluid. We focus on the case where the scalar potential as well as the ordinary matter sector are absent. Then, we investigate gravity-driven acceleration and kinetic inflation in this noncommutative BD cosmology. In contrast to the commutative case, in which the scale factor and BD scalar field are in a power-law form, in the noncommutative case the power-law scalar factor is multiplied by a dynamical exponential warp factor. This warp factor depends on the noncommutative parameter as well as the momentum conjugate associated to the BD scalar field. We show that the BD scalar field and the scale factor effectively depend on the noncommutative parameter. For very small values of this parameter, we obtain an appropriate inflationary solution, which can overcome problems within BD standard cosmology in a more efficient manner. Furthermore, a graceful exit from an early acceleration epoch towards a decelerating radiation epoch is provided. For late times, due to the presence of the noncommutative parameter, we obtain a zero acceleration epoch, which can be interpreted as the coarse-grained explanation.


title: Weyl-Pedersen calculus for some semidirect products of nilpotent Lie groups

abstract: For certain nilpotent real Lie groups constructed as semidirect products, algebras of invariant differential operators on some coadjoint orbits are used in the study of boundedness properties of the Weyl-Pedersen calculus of their corresponding unitary irreducible representations. Our main result is applicable to all unitary irreducible representations of arbitrary 3-step nilpotent Lie groups.


title: Modelling cross-border systemic risk in the European banking sector: a copula approach

abstract: We propose a new methodology based on the Marshall-Olkin (MO) copula to model cross-border systemic risk. The proposed framework estimates the impact of the systematic and idiosyncratic components on systemic risk. Initially, we propose a maximum-likelihood method to estimate the parameter of the MO copula. In order to use the data on non-distressed banks for these estimates, we consider times to bank failures as censored samples. Hence, we propose an estimation procedure for the MO copula on censored data. The empirical evidence from European banks shows that the proposed censored model avoid possible underestimation of the contagion risk.


title: On the unit of a monoidal model category

abstract: In this paper we show how to modify cofibrations in a monoidal model category so that the tensor unit becomes cofibrant while keeping the same weak equivalences. We obtain aplications to enriched categories and coloured operads in stable homotopy theory.


title: Geometric Network Comparison

abstract: Network analysis has a crucial need for tools to compare networks and assess the significance of differences between networks. We propose a principled statistical approach to network comparison that approximates networks as probability distributions on negatively curved manifolds. We outline the theory, as well as implement the approach on simulated networks.