Catalog footprint

What is connected

60works
36topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

60 published item(s)

preprint2022arXiv

Quantification of Robotic Surgeries with Vision-Based Deep Learning

Surgery is a high-stakes domain where surgeons must navigate critical anatomical structures and actively avoid potential complications while achieving the main task at hand. Such surgical activity has been shown to affect long-term patient outcomes. To better understand this relationship, whose mechanics remain unknown for the majority of surgical procedures, we hypothesize that the core elements of surgery must first be quantified in a reliable, objective, and scalable manner. We believe this is a prerequisite for the provision of surgical feedback and modulation of surgeon performance in pursuit of improved patient outcomes. To holistically quantify surgeries, we propose a unified deep learning framework, entitled Roboformer, which operates exclusively on videos recorded during surgery to independently achieve multiple tasks: surgical phase recognition (the what of surgery), gesture classification and skills assessment (the how of surgery). We validated our framework on four video-based datasets of two commonly-encountered types of steps (dissection and suturing) within minimally-invasive robotic surgeries. We demonstrated that our framework can generalize well to unseen videos, surgeons, medical centres, and surgical procedures. We also found that our framework, which naturally lends itself to explainable findings, identified relevant information when achieving a particular task. These findings are likely to instill surgeons with more confidence in our framework's behaviour, increasing the likelihood of clinical adoption, and thus paving the way for more targeted surgical feedback.

preprint2021arXiv

Control of Scanning Quantum Dot Microscopy

Scanning quantum dot microscopy is a recently developed high-resolution microscopy technique that is based on atomic force microscopy and is capable of imaging the electrostatic potential of nanostructures like molecules or single atoms. Recently, it could be shown that it not only yields qualitatively but also quantitatively cutting edge images even on an atomic level. In this paper we present how control is a key enabling element to this. The developed control approach consists of a two-degree-of-freedom control framework that comprises a feedforward and a feedback part. For the latter we design two tailored feedback controllers. The feedforward part generates a reference for the current scanned line based on the previously scanned one. We discuss in detail various aspects of the presented control approach and its implications for scanning quantum dot microscopy. We evaluate the influence of the feedforward part and compare the two proposed feedback controllers. The proposed control algorithms speed up scanning quantum dot microscopy by more than a magnitude and enable to scan large sample areas.

preprint2021arXiv

Vortical flow structures induced by red blood cells in capillaries

Knowledge about the flow field of the plasma around the red blood cells in capillary flow is important for a physical understanding of blood flow and the transport of micro- and nanoparticles and molecules in the flowing plasma. We conduct an experimental study on the flow field around red blood cells in capillary flow that are complemented by simulations of vortical flow between red blood cells. Red blood cells were injected in a 10x12 micrometer rectangular microchannel at a low hematocrit and the flow field around a single or two cells have been characterized thanks to a highspeed camera and by tracking 250 nm nanoparticles in flow behaving as tracers. While the flow field around a steady croissant shape is found to be relatively similar to that of a rigid sphere, the flow field around a slipper shape exhibits a small vortex at the rear of the red blood cell. Even more pronounced are vortex-like structures observed in the central region between two neighboring croissants. Conclusions: The rotation frequency of the vortices is to a good approximation, inversely proportional to the distance between the cells. Our experimental data are confirmed and complemented by numerical simulations.

preprint2020arXiv

A generic ensemble based deep convolutional neural network for semi-supervised medical image segmentation

Deep learning based image segmentation has achieved the state-of-the-art performance in many medical applications such as lesion quantification, organ detection, etc. However, most of the methods rely on supervised learning, which require a large set of high-quality labeled data. Data annotation is generally an extremely time-consuming process. To address this problem, we propose a generic semi-supervised learning framework for image segmentation based on a deep convolutional neural network (DCNN). An encoder-decoder based DCNN is initially trained using a few annotated training samples. This initially trained model is then copied into sub-models and improved iteratively using random subsets of unlabeled data with pseudo labels generated from models trained in the previous iteration. The number of sub-models is gradually decreased to one in the final iteration. We evaluate the proposed method on a public grand-challenge dataset for skin lesion segmentation. Our method is able to significantly improve beyond fully supervised model learning by incorporating unlabeled data.

preprint2020arXiv

Disaster Resilience and Asset Prices

This paper investigates whether security markets price the effect of social distancing on firms' operations. We document that firms that are more resilient to social distancing significantly outperformed those with lower resilience during the COVID-19 outbreak, even after controlling for the standard risk factors. Similar cross-sectional return differentials already emerged before the COVID-19 crisis: the 2014-19 cumulative return differential between more and less resilient firms is of similar size as during the outbreak, suggesting growing awareness of pandemic risk well in advance of its materialization. Finally, we use stock option prices to infer the market's return expectations after the onset of the pandemic: even at a two-year horizon, stocks of more pandemic-resilient firms are expected to yield significantly lower returns than less resilient ones, reflecting their lower exposure to disaster risk. Hence, going forward, markets appear to price exposure to a new risk factor, namely, pandemic risk.

preprint2020arXiv

Fusing Online Gaussian Process-Based Learning and Control for Scanning Quantum Dot Microscopy

Elucidating electrostatic surface potentials contributes to a deeper understanding of the nature of matter and its physicochemical properties, which is the basis for a wide field of applications. Scanning quantum dot microscopy, a recently developed technique allows to measure such potentials with atomic resolution. For an efficient deployment in scientific practice, however, it is essential to speed up the scanning process. To this end we employ a two-degree-of-freedom control paradigm, in which a Gaussian process is used as the feedforward part. We present a tailored online learning scheme of the Gaussian process, adapted to scanning quantum dot microscopy, that includes hyperparameter optimization during operation to enable fast and precise scanning of arbitrary surface structures. For the potential application in practice, the accompanying computational cost is reduced evaluating different sparse approximation approaches. The fully independent training conditional approximation, used on a reduced set of active training data, is found to be the most promising approach.

preprint2019arXiv

Flagellar number governs bacterial spreading and transport efficiency

Peritrichous bacteria synchronize and bundle their flagella to actively swim while disruption of the bundle leads to tumbling. It is still not known whether the number of flagella represents an evolutionary adaptation towards optimizing bacterial navigation. Here, we study the swimming dynamics of differentially flagellated Bacillus subtilis strains in a quasi-two-dimensional system. We find that decreasing the number of flagella N reduces the average turning angle between two successive run phases and enhances the duration and directional persistence of the run phase. As a result, having less flagella is beneficial for long-distance transport and fast spreading, while having a lot of flagella is advantageous for the processes which require localization and slow dynamics, such as biofilm formation. We develop a two-state random walk model that incorporates spontaneous switchings between the states and yields exact analytical expressions for transport properties, in remarkable agreement with experiments. The results of numerical simulations based on our two-state model suggest that the efficiency of searching and exploring the environment is optimized at intermediate values of N. The optimal choice of N, for which the search time is minimized, decreases with increasing the size of the environment in which the bacteria swim.

preprint2019arXiv

Inelastic electron tunneling spectroscopy for probing strongly correlated many-body systems by scanning tunneling microscopy

We present an extension of the tunneling theory for scanning tunneling microcopy (STM) to include different types of vibrational-electronic couplings responsible for inelastic contributions to the tunnel current in the strong-coupling limit. It allows for a better understanding of more complex scanning tunneling spectra of molecules on a metallic substrate in separating elastic and inelastic contributions. The starting point is the exact solution of the spectral functions for the electronic active local orbitals in the absence of the STM tip. This includes electron-phonon coupling in the coupled system comprising the molecule and the substrate to arbitrary order including the anti-adiabatic strong coupling regime as well as the Kondo effect on a free electron spin of the molecule. The tunneling current is derived in second order of the tunneling matrix element which is expanded in powers of the relevant vibrational displacements. We use the results of an ab-initio calculation for the single-particle electronic properties as an adapted material-specific input for a numerical renormalization group approach for accurately determining the electronic properties of a NTCDA molecule on Ag(111) as a challenging sample system for our theory. Our analysis shows that the mismatch between the ab-initio many-body calculation of the tunnel current in the absence of any electron-phonon coupling to the experiment scanning tunneling spectra can be resolved by including two mechanisms: (i) a strong unconventional Holstein term on the local substrate orbital leads to reduction of the Kondo temperature and (ii) a different electron-vibrational coupling to the tunneling matrix element is responsible for inelastic steps in the $dI/dV$ curve at finite frequencies.

preprint2016arXiv

Adaptive Data Communication Interface: A User-Centric Visual Data Interpretation Framework

In this position paper, we present ideas about creating a next generation framework towards an adaptive interface for data communication and visualisation systems. Our objective is to develop a system that accepts large data sets as inputs and provides user-centric, meaningful visual information to assist owners to make sense of their data collection. The proposed framework comprises four stages: (i) the knowledge base compilation, where we search and collect existing state-ofthe-art visualisation techniques per domain and user preferences; (ii) the development of the learning and inference system, where we apply artificial intelligence techniques to learn, predict and recommend new graphic interpretations (iii) results evaluation; and (iv) reinforcement and adaptation, where valid outputs are stored in our knowledge base and the system is iteratively tuned to address new demands. These stages, as well as our overall vision, limitations and possible challenges are introduced in this article. We also discuss further extensions of this framework for other knowledge discovery tasks.

preprint2016arXiv

Applying Interval Type-2 Fuzzy Rule Based Classifiers Through a Cluster-Based Class Representation

Fuzzy Rule-Based Classification Systems (FRBCSs) have the potential to provide so-called interpretable classifiers, i.e. classifiers which can be introspective, understood, validated and augmented by human experts by relying on fuzzy-set based rules. This paper builds on prior work for interval type-2 fuzzy set based FRBCs where the fuzzy sets and rules of the classifier are generated using an initial clustering stage. By introducing Subtractive Clustering in order to identify multiple cluster prototypes, the proposed approach has the potential to deliver improved classification performance while maintaining good interpretability, i.e. without resulting in an excessive number of rules. The paper provides a detailed overview of the proposed FRBC framework, followed by a series of exploratory experiments on both linearly and non-linearly separable datasets, comparing results to existing rule-based and SVM approaches. Overall, initial results indicate that the approach enables comparable classification performance to non rule-based classifiers such as SVM, while often achieving this with a very small number of rules.

preprint2016arXiv

Cluster of red blood cells in microcapillary flow: hydrodynamic versus macromolecule induced interaction

We present experiments on RBCs that flow through microcapillaries under physiological conditions. We show that the RBC clusters form as a subtle imbrication between hydrodynamics interaction and adhesion forces because of plasma proteins. Clusters form along the capillaries and macromolecule-induced adhesion contribute to their stability. However, at high yet physiological flow velocities, shear stresses overcome part of the adhesion forces, and cluster stabilization due to hydrodynamics becomes stronger. For the case of pure hydrodynamic interaction, cell-to-cell distances have a pronounced bimodal distribution. Our 2D-numerical simulations on vesicles captures the transition between adhesive and non-adhesive clusters at different flow velocities.

preprint2016arXiv

Exploring Differences in Interpretation of Words Essential in Medical Expert-Patient Communication

In the context of cancer treatment and surgery, quality of life assessment is a crucial part of determining treatment success and viability. In order to assess it, patients completed questionnaires which employ words to capture aspects of patients well-being are the norm. As the results of these questionnaires are often used to assess patient progress and to determine future treatment options, it is important to establish that the words used are interpreted in the same way by both patients and medical professionals. In this paper, we capture and model patients perceptions and associated uncertainty about the words used to describe the level of their physical function used in the highly common (in Sarcoma Services) Toronto Extremity Salvage Score (TESS) questionnaire. The paper provides detail about the interval-valued data capture as well as the subsequent modelling of the data using fuzzy sets. Based on an initial sample of participants, we use Jaccard similarity on the resulting words models to show that there may be considerable differences in the interpretation of commonly used questionnaire terms, thus presenting a very real risk of miscommunication between patients and medical professionals as well as within the group of medical professionals.

preprint2016arXiv

Measuring Player's Behaviour Change over Time in Public Goods Game

An important issue in public goods game is whether player's behaviour changes over time, and if so, how significant it is. In this game players can be classified into different groups according to the level of their participation in the public good. This problem can be considered as a concept drift problem by asking the amount of change that happens to the clusters of players over a sequence of game rounds. In this study we present a method for measuring changes in clusters with the same items over discrete time points using external clustering validation indices and area under the curve. External clustering indices were originally used to measure the difference between suggested clusters in terms of clustering algorithms and ground truth labels for items provided by experts. Instead of different cluster label comparison, we use these indices to compare between clusters of any two consecutive time points or between the first time point and the remaining time points to measure the difference between clusters through time points. In theory, any external clustering indices can be used to measure changes for any traditional (non-temporal) clustering algorithm, due to the fact that any time point alone is not carrying any temporal information. For the public goods game, our results indicate that the players are changing over time but the change is smooth and relatively constant between any two time points.

preprint2016arXiv

Modelling Cyber-Security Experts' Decision Making Processes using Aggregation Operators

An important role carried out by cyber-security experts is the assessment of proposed computer systems, during their design stage. This task is fraught with difficulties and uncertainty, making the knowledge provided by human experts essential for successful assessment. Today, the increasing number of progressively complex systems has led to an urgent need to produce tools that support the expert-led process of system-security assessment. In this research, we use weighted averages (WAs) and ordered weighted averages (OWAs) with evolutionary algorithms (EAs) to create aggregation operators that model parts of the assessment process. We show how individual overall ratings for security components can be produced from ratings of their characteristics, and how these individual overall ratings can be aggregated to produce overall rankings of potential attacks on a system. As well as the identification of salient attacks and weak points in a prospective system, the proposed method also highlights which factors and security components contribute most to a component's difficulty and attack ranking respectively. A real world scenario is used in which experts were asked to rank a set of technical attacks, and to answer a series of questions about the security components that are the subject of the attacks. The work shows how finding good aggregation operators, and identifying important components and factors of a cyber-security problem can be automated. The resulting operators have the potential for use as decision aids for systems designers and cyber-security experts, increasing the amount of assessment that can be achieved with the limited resources available.

preprint2016arXiv

Optimising Rule-Based Classification in Temporal Data

This study optimises manually derived rule-based expert system classification of objects according to changes in their properties over time. One of the key challenges that this study tries to address is how to classify objects that exhibit changes in their behaviour over time, for example how to classify companies' share price stability over a period of time or how to classify students' preferences for subjects while they are progressing through school. A specific case the paper considers is the strategy of players in public goods games (as common in economics) across multiple consecutive games. Initial classification starts from expert definitions specifying class allocation for players based on aggregated attributes of the temporal data. Based on these initial classifications, the optimisation process tries to find an improved classifier which produces the best possible compact classes of objects (players) for every time point in the temporal data. The compactness of the classes is measured by a cost function based on internal cluster indices like the Dunn Index, distance measures like Euclidean distance or statistically derived measures like standard deviation. The paper discusses the approach in the context of incorporating changing player strategies in the aforementioned public good games, where common classification approaches so far do not consider such changes in behaviour resulting from learning or in-game experience. By using the proposed process for classifying temporal data and the actual players' contribution during the games, we aim to produce a more refined classification which in turn may inform the interpretation of public goods game data.

preprint2016arXiv

Precision measurement of the local bias of dark matter halos

We present accurate measurements of the linear, quadratic, and cubic local bias of dark matter halos, using curved "separate universe" N-body simulations which effectively incorporate an infinite-wavelength overdensity. This can be seen as an exact implementation of the peak-background split argument. We compare the results with the linear and quadratic bias measured from the halo-matter power spectrum and bispectrum, and find good agreement. On the other hand, the standard peak-background split applied to the Sheth & Tormen (1999) and Tinker et al. (2008) halo mass functions matches the measured linear bias parameter only at the level of 10%. The prediction from the excursion set-peaks approach performs much better, which can be attributed to the stochastic moving barrier employed in the excursion set-peaks prediction. We also provide convenient fitting formulas for the nonlinear bias parameters $b_2(b_1)$ and $b_3(b_1)$, which work well over a range of redshifts.

preprint2015arXiv

On the rheology of red blood cell suspensions with different amounts of dextran: separating the effect of aggregation and increase in viscosity of the suspending phase

We investigate the shear thinning of red blood cell - dextran suspensions. Microscopic images show that at low polymer concentration, aggregation increases with increasing concentration until it reaches a maximum and then decreases again to non-aggregation. This bell shape dependency is also deduced from the rheological measurements, if the data are correctly normalized by the viscosity of the suspending phase since a significant amount of polymers adsorb to the cell surfaces. We find that the position of the maximum of this shear rate dependent bell shape increases with increasing viscosity of the suspending phase, which indicates a that the dynamic process of aggregation and disaggregation is coupled via hydrodynamic interactions. This hydrodynamic coupling can be suppressed by characterizing a suspension of 80% hematrocrit which yields good agreement with the results from the microscopical images.

preprint2015arXiv

Position-dependent correlation function from the SDSS-III Baryon Oscillation Spectroscopic Survey Data Release 10 CMASS Sample

We report on the first measurement of the three-point function with the position-dependent correlation function from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 10 CMASS sample. This new observable measures the correlation between two-point functions of galaxy pairs within different subvolumes, $\hatξ({\rm r},{\rm r}_L)$, where ${\rm r}_L$ is the location of a subvolume, and the corresponding mean overdensities, $\barδ({\rm r}_L)$. This correlation, which we call the "integrated three-point function", $iζ(r)=\langle\hatξ({\rm r},{\rm r}_L)\barδ({\rm r}_L)\rangle$, measures a three-point function of two short- and one long-wavelength modes, and is generated by nonlinear gravitational evolution and possibly also by the physics of inflation. The $iζ(r)$ measured from the BOSS data lies within the scatter of those from the mock galaxy catalogs in redshift space, yielding a ten-percent-level determination of the amplitude of $iζ(r)$. The tree-level perturbation theory in redshift space predicts how this amplitude depends on the linear and quadratic nonlinear galaxy bias parameters ($b_1$ and $b_2$), as well as on the amplitude and linear growth rate of matter fluctuations ($σ_8$ and $f$). Combining $iζ(r)$ with the constraints on $b_1σ_8$ and $fσ_8$ from the global two-point correlation function and that on $σ_8$ from the weak lensing signal of BOSS galaxies, we measure $b_2=0.41\pm0.41$ (68% C.L.) assuming standard perturbation theory at the tree level and the local bias model.

preprint2015arXiv

Scanning Quantum Dot Microscopy

Interactions between atomic and molecular objects are to a large extent defined by the nanoscale electrostatic potentials which these objects produce. We introduce a scanning probe technique that enables three-dimensional imaging of local electrostatic potential fields with sub-nanometer resolution. Registering single electron charging events of a molecular quantum dot attached to the tip of a (qPlus tuning fork) atomic force microscope operated at 5 K, we quantitatively measure the quadrupole field of a single molecule and the dipole field of a single metal adatom, both adsorbed on a clean metal surface. Because of its high sensitivity, the technique can record electrostatic potentials at large distances from their sources, which above all will help to image complex samples with increased surface roughness.

preprint2015arXiv

Single particles accelerate final stages of capillary break up

Droplet formation of suspensions is present in many industrial and technological processes such as coating and food engineering. Whilst the finite time singularity of the minimum neck diameter in capillary break-up of simple liquids can be described by well known self-similarity solutions, the pinching of non-Brownian suspension depends in a complex way on the particle dynamics in the thinning thread. Here we focus on the very dilute regime where the filament contains only isolated beads to identify the physical mechanisms leading to the pronounced acceleration of the filament thinning observed. This accelerated regime is characterized by an asymmetric shape of the filament with an enhanced curvature that depends on the size and the spatial distribution of the particles within the capillary thread.

preprint2015arXiv

The angle-averaged squeezed limit of nonlinear matter N-point functions

We show that in a certain, angle-averaged squeezed limit, the $N$-point function of matter is related to the response of the matter power spectrum to a long-wavelength density perturbation, $P^{-1}d^nP(k|δ_L)/dδ_L^n|_{δ_L=0}$, with $n=N-2$. By performing N-body simulations with a homogeneous overdensity superimposed on a flat Friedmann-Robertson-Lemaître-Walker (FRLW) universe using the \emph{separate universe} approach, we obtain measurements of the nonlinear matter power spectrum response up to $n=3$, which is equivalent to measuring the fully nonlinear matter $3-$ to $5-$point function in this squeezed limit. The sub-percent to few percent accuracy of those measurements is unprecedented. We then test the hypothesis that nonlinear $N$-point functions at a given time are a function of the linear power spectrum at that time, which is predicted by standard perturbation theory (SPT) and its variants that are based on the ideal pressureless fluid equations. Specifically, we compare the responses computed from the separate universe simulations and simulations with a rescaled initial (linear) power spectrum amplitude. We find discrepancies of 10\% at $k\simeq 0.2 - 0.5 \,h\,{\rm Mpc}^{-1}$ for $5-$ to $3-$point functions at $z=0$. The discrepancy occurs at higher wavenumbers at $z=2$. Thus, SPT and its variants, carried out to arbitrarily high order, are guaranteed to fail to describe matter $N$-point functions ($N>2$) around that scale.

preprint2015arXiv

The power spectrum and bispectrum of SDSS DR11 BOSS galaxies I: bias and gravity

We analyse the anisotropic clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS Data Release 11 sample, which consists of $690 827$ galaxies in the redshift range $0.43 < z < 0.70$ and has a sky coverage of $8 498$ deg$^2$ corresponding to an effective volume of $\sim6\,\rm{Gpc}^3$. We fit the Fourier space statistics, the power spectrum and bispectrum monopoles to measure the linear and quadratic bias parameters, $b_1$ and $b_2$, for a non-linear non-local bias model, the growth of structure parameter $f$ and the amplitude of dark matter density fluctuations parametrised by $σ_8$. We obtain $b_1(z_{\rm eff})^{1.40}σ_8(z_{\rm eff})=1.672\pm 0.060$ and $b_2^{0.30}(z_{\rm eff})σ_8(z_{\rm eff})=0.579\pm0.082$ at the effective redshift of the survey, $z_{\rm eff}=0.57$. The main cosmological result is the constraint on the combination $f^{0.43}(z_{\rm eff})σ_8(z_{\rm eff})=0.582\pm0.084$, which is complementary to $fσ_8$ constraints obtained from 2-point redshift space distortion analyses. A less conservative analysis yields $f^{0.43}(z_{\rm eff})σ_8(z_{\rm eff})=0.584\pm0.051$. We ensure that our result is robust by performing detailed systematic tests using a large suite of survey galaxy mock catalogs and N-body simulations. The constraints on $f^{0.43}σ_8$ are useful for setting additional constrains on neutrino mass, gravity, curvature as well as the number of neutrino species from galaxy surveys analyses (as presented in a companion paper).

preprint2015arXiv

The power spectrum and bispectrum of SDSS DR11 BOSS galaxies II: cosmological interpretation

We examine the cosmological implications of the measurements of the linear growth rate of cosmological structure obtained in a companion paper from the power spectrum and bispectrum monopoles of the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey Data, Release 11, CMASS galaxies. This measurement was of $f^{0.43}σ_8$, where $σ_8$ is the amplitude of dark matter density fluctuations, and $f$ is the linear growth rate, at the effective redshift of the survey, $z_{\rm eff}=0.57$. In conjunction with Cosmic Microwave Background (CMB) data, interesting constraints can be placed on models with non-standard neutrino properties and models where gravity deviates from general relativity on cosmological scales. In particular, the sum of the masses of the three species of the neutrinos is constrained to $m_ν<0.49\,{\rm eV}$ (at 95\% confidence level) when the $f^{0.43}σ_8$ measurement is combined with state-of-the-art CMB measurements. Allowing the effective number of neutrinos to vary as a free parameter does not significantly change these results. When we combine the measurement of $f^{0.43}σ_8$ with the complementary measurement of $fσ_8$ from the monopole and quadrupole of the two-point correlation function we are able to obtain an independent measurements of $f$ and $σ_8$. We obtain $f=0.63\pm0.16$ and $σ_8=0.710\pm0.086$ (68\% confidence level). This is the first time when these parameters have been able to be measured independently using the redshift-space power spectrum and bispectrum measurements from galaxy clustering data only.

preprint2014arXiv

A Fuzzy Directional Distance Measure

The measure of distance between two fuzzy sets is a fundamental tool within fuzzy set theory, however, distance measures currently within the literature use a crisp value to represent the distance between fuzzy sets. A real valued distance measure is developed into a fuzzy distance measure which better reflects the uncertainty inherent in fuzzy sets and a fuzzy directional distance measure is presented, which accounts for the direction of change between fuzzy sets. A multiplicative version is explored as a full maximal assignment is computationally intractable so an intermediate solution is offered.

preprint2014arXiv

Analysing Fuzzy Sets Through Combining Measures of Similarity and Distance

Reasoning with fuzzy sets can be achieved through measures such as similarity and distance. However, these measures can often give misleading results when considered independently, for example giving the same value for two different pairs of fuzzy sets. This is particularly a problem where many fuzzy sets are generated from real data, and while two different measures may be used to automatically compare such fuzzy sets, it is difficult to interpret two different results. This is especially true where a large number of fuzzy sets are being compared as part of a reasoning system. This paper introduces a method for combining the results of multiple measures into a single measure for the purpose of analysing and comparing fuzzy sets. The combined measure alleviates ambiguous results and aids in the automatic comparison of fuzzy sets. The properties of the combined measure are given, and demonstrations are presented with discussions on the advantages over using a single measure.

preprint2014arXiv

Comparison of Distance Metrics for Hierarchical Data in Medical Databases

Distance metrics are broadly used in different research areas and applications, such as bio-informatics, data mining and many other fields. However, there are some metrics, like pq-gram and Edit Distance used specifically for data with a hierarchical structure. Other metrics used for non-hierarchical data are the geometric and Hamming metrics. We have applied these metrics to The Health Improvement Network (THIN) database which has some hierarchical data. The THIN data has to be converted into a tree-like structure for the first group of metrics. For the second group of metrics, the data are converted into a frequency table or matrix, then for all metrics, all distances are found and normalised. Based on this particular data set, our research question: which of these metrics is useful for THIN data? This paper compares the metrics, particularly the pq-gram metric on finding the similarities of patients' data. It also investigates the similar patients who have the same close distances as well as the metrics suitability for clustering the whole patient population. Our results show that the two groups of metrics perform differently as they represent different structures of the data. Nevertheless, all the metrics could represent some similar data of patients as well as discriminate sufficiently well in clustering the patient population using $k$-means clustering algorithm.

preprint2014arXiv

Dark matter and halo bispectrum in redshift space: theory and applications

We present a phenomenological modification of the standard perturbation theory prediction for the bispectrum in redshift space that allows us to extend the model to mildly non-linear scales over a wide range of redshifts, $z\leq1.5$. We find that we can describe the bispectrum of dark matter particles with $\sim5%$ accuracy for $k_i\lesssim0.10\,h/{\rm Mpc}$ at $z=0$, for $k_i\lesssim0.15\,h/{\rm Mpc}$ at $z=0.5$, for $k_i\lesssim0.17\,h/{\rm Mpc}$ at $z=1.0$ and for $k_i\lesssim0.20\,h/{\rm Mpc}$ at $z=1.5$. We also test that the fitting formula is able to describe with similar accuracy the bispectrum of cosmologies with different $Ω_m$, in the range $0.2\lesssim Ω_m \lesssim 0.4$, and consequently with different values of the logarithmic grow rate $f$ at $z=0$, $0.4\lesssim f(z=0) \lesssim 0.6$. We apply this new formula to recover the bias parameters, $f$ and $σ_8$, by combining the redshift space power spectrum monopole and quadrupole with the bispectrum monopole for both dark matter particles and haloes. We find that the combination of these three statistics can break the degeneracy between $b_1$, $f$ and $σ_8$. For dark matter particles the new model can be used to recover $f$ and $σ_8$ with $\sim1%$ accuracy. For dark matter haloes we find that $f$ and $σ_8$ present larger systematic shifts, $\sim10%$. The systematic offsets arise because of limitations in the modelling of the interplay between bias and redshift space distortions, and represent a limitation as the statistical errors of forthcoming surveys reach this level. Conveniently, we find that these residual systematics are mitigated for combinations of parameters. The improvement on the modelling of the bispectrum presented in this paper will be useful for extracting information from current and future galaxy surveys. [abridged]

preprint2014arXiv

Position-dependent power spectrum of the large-scale structure: a novel method to measure the squeezed-limit bispectrum

The influence of large-scale density fluctuations on structure formation on small scales is described by the three-point correlation function (bispectrum) in the so-called "squeezed configurations," in which one wavenumber, say $k_3$, is much smaller than the other two, i.e., $k_3\ll k_1\approx k_2$. This bispectrum is generated by non-linear gravitational evolution and possibly also by inflationary physics. In this paper, we use this fact to show that the bispectrum in the squeezed configurations can be measured without employing three-point function estimators. Specifically, we use the "position-dependent power spectrum," i.e., the power spectrum measured in smaller subvolumes of the survey (or simulation box), and correlate it with the mean overdensity of the corresponding subvolume. This correlation directly measures an integral of the bispectrum dominated by the squeezed configurations. Measuring this correlation is only slightly more complex than measuring the power spectrum itself, and sidesteps the considerable complexity of the full bispectrum estimation. We use cosmological $N$-body simulations of collisionless particles with Gaussian initial conditions to show that the measured correlation between the position-dependent power spectrum and the long-wavelength overdensity agrees with the theoretical expectation. The position-dependent power spectrum thus provides a new, efficient, and promising way to measure the squeezed-limit bispectrum from large-scale structure observations such as galaxy redshift surveys.

preprint2014arXiv

Separate Universe Simulations

The large-scale statistics of observables such as the galaxy density are chiefly determined by their dependence on the local coarse-grained matter density. This dependence can be measured directly and efficiently in N-body simulations by using the fact that a uniform density perturbation with respect to some fiducial background cosmology is equivalent to modifying the background and including curvature, i.e., by simulating a "separate universe". We derive this mapping to fully non-linear order, and provide a step-by-step description of how to perform and analyse the separate universe simulations. This technique can be applied to a wide range of observables. As an example, we calculate the response of the non-linear matter power spectrum to long-wavelength density perturbations, which corresponds to the angle-averaged squeezed limit of the matter bispectrum and higher $n$-point functions. Using only a modest simulation volume, we obtain results with percent-level precision over a wide range of scales.

preprint2014arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Signs of neutrino mass in current cosmological datasets

We investigate the cosmological implications of the latest growth of structure measurement from the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS Data Release 11 with particular focus on the sum of the neutrino masses, $\sum m_ν$. We examine the robustness of the cosmological constraints from the Baryon Acoustic Oscillation (BAO) scale, the Alcock-Paczynski effect and redshift-space distortions ($D_V/r_s$, $F_{\rm AP}$, $fσ_8$) of \citet{Beutler:2013b}, when introducing a neutrino mass in the power spectrum template. We then discuss how the neutrino mass relaxes discrepancies between the Cosmic Microwave Background (CMB) and other low-redshift measurements within $Λ$CDM. Combining our cosmological constraints with WMAP9 yields $\sum m_ν = 0.36\pm0.14\,$eV ($68\%$ c.l.), which represents a $2.6σ$ preference for non-zero neutrino mass. The significance can be increased to $3.3σ$ when including weak lensing results and other BAO constraints, yielding $\sum m_ν = 0.35\pm0.10\,$eV ($68\%$ c.l.). However, combining CMASS with Planck data reduces the preference for neutrino mass to $\sim 2σ$. When removing the CMB lensing effect in the Planck temperature power spectrum (by marginalising over $A_{\rm L}$), we see shifts of $\sim 1σ$ in $σ_8$ and $Ω_m$, which have a significant effect on the neutrino mass constraints. In case of CMASS plus Planck without the $A_{\rm L}$-lensing signal, we find a preference for a neutrino mass of $\sum m_ν = 0.34\pm 0.14\,$eV ($68\%$ c.l.), in excellent agreement with the WMAP9+CMASS value. The constraint can be tightened to $3.4σ$ yielding $\sum m_ν = 0.36\pm 0.10\,$eV ($68\%$ c.l.) when weak lensing data and other BAO constraints are included.

preprint2014arXiv

Vesicle Dynamics in a Confined Poiseuille Flow: From Steady-State to Chaos

Red blood cells (RBCs) are the major component of blood and the flow of blood is dictated by that of RBCs. We employ vesicles, which consist of closed bilayer membranes enclosing a fluid, as a model system to study the behavior of RBCs under a confined Poiseuille flow. We extensively explore two main parameters: i) the degree of confinement of vesicles within the channel, and ii) the flow strength. Rich and complex dynamics for vesicles are revealed ranging from steady-state shapes (in the form of parachute and slipper) to chaotic dynamics of shape. Chaos occurs through a cascade of multiple periodic oscillations of the vesicle shape. We summarize our results in a phase diagram in the parameter plane (degree of confinement, flow strength). This finding highlights the level of complexity of a flowing vesicle in the small Reynolds number where the flow is laminar in the absence of vesicles and can be rendered turbulent due to elasticity of vesicles.

preprint2013arXiv

Creating Personalised Energy Plans. From Groups to Individuals using Fuzzy C Means Clustering

Changes in the UK electricity market mean that domestic users will be required to modify their usage behaviour in order that supplies can be maintained. Clustering allows usage profiles collected at the household level to be clustered into groups and assigned a stereotypical profile which can be used to target marketing campaigns. Fuzzy C Means clustering extends this by allowing each household to be a member of many groups and hence provides the opportunity to make personalised offers to the household dependent on their degree of membership of each group. In addition, feedback can be provided on how user's changing behaviour is moving them towards more "green" or cost effective stereotypical usage.

preprint2013arXiv

Effective viscosity of non-gravitactic Chlamydomonas Reinhardtii microswimmer suspensions

Active microswimmers are known to affect the macroscopic viscosity of suspensions in a more complex manner than passive particles. For puller-like microswimmers an increase in the viscosity has been observed. It has been suggested that the persistence of the orientation of the microswimmers hinders the rotation that is normally caused by the vorticity. It was previously shown that some sorts of algaes are bottom-heavy swimmers, i.e. their centre of mass is not located in the centre of the body. In this way, the algae affects the vorticity of the flow when it is perpendicular oriented to the axis of gravity. This orientation of gravity to vorticity is given in a rheometer that is equipped with a cone-plate geometry. Here we present measurements of the viscosity both in a cone-plate and a Taylor-Couette cell. The two set-ups yielded the same increase in viscosity although the axis of gravitation in the Taylor-Couette cell is parallel to the direction of vorticity. In a complementary experiment we tested the orientation of the direction of swimming through microscopic observation of single \textit{Chlamydomonas reinhardtii} and could not identify a preferred orientation, i. e. our specific strain of \textit{Chlamydomonas reinhardtii} are not bottom-heavy swimmers. We thus conclude that bottom heaviness is not a prerequisite for the increase of viscosity and that the effect of gravity on the rheology of our strain of \textit{Chlamydomonas reinhardtii} is negligible. This finding reopens the question of whether origin of persistence in the orientation of cells is actually responsible for the increased viscosity of the suspension.

preprint2013arXiv

Extending Similarity Measures of Interval Type-2 Fuzzy Sets to General Type-2 Fuzzy Sets

Similarity measures provide one of the core tools that enable reasoning about fuzzy sets. While many types of similarity measures exist for type-1 and interval type-2 fuzzy sets, there are very few similarity measures that enable the comparison of general type-2 fuzzy sets. In this paper, we introduce a general method for extending existing interval type-2 similarity measures to similarity measures for general type-2 fuzzy sets. Specifically, we show how similarity measures for interval type-2 fuzzy sets can be employed in conjunction with the zSlices based general type-2 representation for fuzzy sets to provide measures of similarity which preserve all the common properties (i.e. reflexivity, symmetry, transitivity and overlapping) of the original interval type-2 similarity measure. We demonstrate examples of such extended fuzzy measures and provide comparisons between (different types of) interval and general type-2 fuzzy measures.

preprint2013arXiv

Measuring the Directional Distance Between Fuzzy Sets

The measure of distance between two fuzzy sets is a fundamental tool within fuzzy set theory. However, current distance measures within the literature do not account for the direction of change between fuzzy sets; a useful concept in a variety of applications, such as Computing With Words. In this paper, we highlight this utility and introduce a distance measure which takes the direction between sets into account. We provide details of its application for normal and non-normal, as well as convex and non-convex fuzzy sets. We demonstrate the new distance measure using real data from the MovieLens dataset and establish the benefits of measuring the direction between fuzzy sets.

preprint2013arXiv

Serpentine channels: micro -- rheometers for fluid relaxation times

We propose a novel device capable of measuring the relaxation time of viscoelastic fluids as small as 1\,ms. In contrast to most rheometers, which by their very nature are concerned with producing viscometric or nearly-viscometric flows, here we make use of an elastic instability which occurs in the flow of viscoelastic fluids with curved streamlines. To calibrate the rheometer we combine simple scaling arguments with relaxation times obtained from first normal-stress difference data measured in a classical shear rheometer. As an additional check we also compare these relaxation times to those obtained from Zimm theory and good agreement is observed. Once calibrated, we show how the serpentine rheometer can be used to access smaller polymer concentrations and lower solvent viscosities where classical measurements become difficult or impossible to use due to inertial and/or resolution limitations. In the absence of calibration the serpentine channel can still be a very useful comparative or index device.

preprint2013arXiv

The Coyote Universe III: Simulation Suite and Precision Emulator for the Nonlinear Matter Power Spectrum

Many of the most exciting questions in astrophysics and cosmology, including the majority of observational probes of dark energy, rely on an understanding of the nonlinear regime of structure formation. In order to fully exploit the information available from this regime and to extract cosmological constraints, accurate theoretical predictions are needed. Currently such predictions can only be obtained from costly, precision numerical simulations. This paper is the third in a series aimed at constructing an accurate calibration of the nonlinear mass power spectrum on Mpc scales for a wide range of currently viable cosmological models, including dark energy. The first two papers addressed the numerical challenges, and the scheme by which an interpolator was built from a carefully chosen set of cosmological models. In this paper we introduce the "Coyote Univers"' simulation suite which comprises nearly 1,000 N-body simulations at different force and mass resolutions, spanning 38 wCDM cosmologies. This large simulation suite enables us to construct a prediction scheme, or emulator, for the nonlinear matter power spectrum accurate at the percent level out to k~1 h/Mpc. We describe the construction of the emulator, explain the tests performed to ensure its accuracy, and discuss how the central ideas may be extended to a wider range of cosmological models and applications. A power spectrum emulator code is released publicly as part of this paper.

preprint2013arXiv

Wide band Fresnel super-resolution applied to capillary break up of viscoelastic fluids

We report a technique based on Fresnel diffraction with white illumination that permits the resolution of capillary surface patterns of less than 100 nanometers. We investigate Rayleigh Plateaux like instability on a viscoelastic capillary bridge and show that we can overcome the resolution limit of optical microscopy. The viscoelastic filaments are approximately 20 microns thick at the end of the thinning process when the instability sets in. The wavy distortions grow exponentially in time and the pattern is resolved by an image treatment that is based on an approximation of the measured rising flank of the first Fresnel peak.

preprint2012arXiv

An improved fitting formula for the dark matter bispectrum

In this paper we present an improved fitting formula for the dark matter bispectrum motivated by the previous phenomenological approach of Scoccimarro & Couchman (2001). We use a set of LCDM simulations to calibrate the fitting parameters in the k-range of 0.03 h/Mpc<k<0.4 h/Mpc and in the redshift range of 0<z<1.5. This new proposed fit describes well the BAO-features although it was not designed to. The deviation between the simulations output and our analytic prediction is typically less than 5% and in the worst case is never above 10%. We envision that this new analytic fitting formula will be very useful in providing reliable predictions for the non-linear dark matter bispectrum for LCDM models.

preprint2012arXiv

Bayesian Analysis of Inflation III: Slow Roll Reconstruction Using Model Selection

We implement Slow Roll Reconstruction -- an optimal solution to the inverse problem for inflationary cosmology -- within ModeCode, a publicly available solver for the inflationary dynamics. We obtain up-to-date constraints on the reconstructed inflationary potential, derived from the WMAP 7-year dataset and South Pole Telescope observations, combined with large scale structure data derived from SDSS Data Release 7. Using ModeCode in conjunction with the MultiNest sampler, we compute Bayesian evidence for the reconstructed potential at each order in the truncated slow roll hierarchy. We find that the data are well-described by the first two slow roll parameters, εand η, and that there is no need to include a nontrivial ξparameter.

preprint2012arXiv

Effects of the neutrino mass splitting on the non-linear matter power spectrum

We have performed cosmological N-body simulations which include the effect of the masses of the individual neutrino species. The simulations were aimed at studying the effect of different neutrino hierarchies on the matter power spectrum. Compared to the linear theory predictions, we find that non-linearities enhance the effect of hierarchy on the matter power spectrum at mildly non-linear scales. The difference between the different hierarchies is about 0.5% for a sum of neutrino masses of 0.1eV. Albeit this is a small effect, it is potentially measurable from upcoming surveys. In combination with neutrinoless double-beta decay experiments, this opens up the possibility of using the sky to determine if neutrinos are Majorana or Dirac fermions.

preprint2012arXiv

Elasto-inertial turbulence

Turbulence is ubiquitous in nature yet even for the case of ordinary Newtonian fluids like water our understanding of this phenomenon is limited. Many liquids of practical importance however are more complicated (e.g. blood, polymer melts or paints), they exhibit elastic as well as viscous characteristics and the relation between stress and strain is nonlinear. We here demonstrate for a model system of such complex fluids that at high shear rates turbulence is not simply modified as previously believed but it is suppressed and replaced by a new type of disordered motion, elasto-inertial turbulence (EIT). EIT is found to occur at much lower Reynolds numbers than Newtonian turbulence and the dynamical properties differ significantly. In particular the drag is strongly reduced and the observed friction scaling resolves a longstanding puzzle in non-Newtonian fluid mechanics regarding the nature of the so-called maximum drag reduction asymptote. Theoretical considerations imply that EIT will arise in complex fluids if the extensional viscosity is sufficiently large.

preprint2012arXiv

Mirror-Descent Methods in Mixed-Integer Convex Optimization

In this paper, we address the problem of minimizing a convex function f over a convex set, with the extra constraint that some variables must be integer. This problem, even when f is a piecewise linear function, is NP-hard. We study an algorithmic approach to this problem, postponing its hardness to the realization of an oracle. If this oracle can be realized in polynomial time, then the problem can be solved in polynomial time as well. For problems with two integer variables, we show that the oracle can be implemented efficiently, that is, in O(ln(B)) approximate minimizations of f over the continuous variables, where B is a known bound on the absolute value of the integer variables.Our algorithm can be adapted to find the second best point of a purely integer convex optimization problem in two dimensions, and more generally its k-th best point. This observation allows us to formulate a finite-time algorithm for mixed-integer convex optimization.

preprint2012arXiv

Perturbation theory approach for the power spectrum: from dark matter in real space to massive haloes in redshift space

We investigate the accuracy of Eulerian perturbation theory for describing the matter and galaxy power spectra in real and redshift space in light of future observational probes for precision cosmology. Comparing the analytical results with a large suite of N-body simulations (160 independent boxes of 13.8 (Gpc/h)^3 volume each, which are publicly available), we find that re-summing terms in the standard perturbative approach predicts the real-space matter power spectrum with an accuracy of < 2% for k < 0.20 h/Mpc at redshifts z < 1.5. This is obtained following the widespread technique of writing the resummed propagator in terms of 1-loop contributions. We show that the accuracy of this scheme increases by considering higher-order terms in the resummed propagator. By combining resummed perturbation theories with several models for the mappings from real to redshift space discussed in the literature, the multipoles of the dark-matter power spectrum can be described with sub-percent deviations from N-body results for k < 0.15h/Mpc at z < 1. As a consequence, the logarithmic growth rate, f, can be recovered with sub-percent accuracy on these scales. Extending the models to massive dark-matter haloes in redshift space, our results describe the monopole term from N-body data within 2% accuracy for scales k < 0.15 h/Mpc at z < 0.5; here f can be recovered within < 5% when the halo bias is known. We conclude that these techniques are suitable to extract cosmological information from future galaxy surveys.

preprint2012arXiv

Quantification of depletion-induced adhesion of Red Blood Cells

Red blood cells (RBC) are known to form aggregates in the forms of rouleaux due to the presence of plasma proteins under physiological conditions. Rouleaux formation can be also induced in vitro by the addition of macromolecules to the RBC solution. Current data on the adhesion strength between red blood cells in their natural discocyte shapes mostly rely on indirect measurements like flow chamber experiments, but on the single cell level data is lacking. Here we present measurements on the dextran induced aggregation of red blood cells by use of atomic force microscopy based single cell force spectroscopy (SCFS). The effects of dextran concentration and molecular weight on the interaction energy of adhering RBCs was determined. The results are in good agreement with a model based on the depletion effect and former experimental studies.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the Data Release 9 Spectroscopic Galaxy Sample

We present measurements of galaxy clustering from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). These use the Data Release 9 (DR9) CMASS sample, which contains 264,283 massive galaxies covering 3275 square degrees with an effective redshift z=0.57 and redshift range 0.43 < z < 0.7. Assuming a concordance Lambda-CDM cosmological model, this sample covers an effective volume of 2.2 Gpc^3, and represents the largest sample of the Universe ever surveyed at this density, n = 3 x 10^-4 h^-3 Mpc^3. We measure the angle-averaged galaxy correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of 5σin both the correlation function and power spectrum. Combining with the SDSS-II Luminous Red Galaxy Sample, the detection significance increases to 6.7σ. Fitting for the position of the acoustic features measures the distance to z=0.57 relative to the sound horizon DV /rs = 13.67 +/- 0.22 at z=0.57. Assuming a fiducial sound horizon of 153.19 Mpc, which matches cosmic microwave background constraints, this corresponds to a distance DV(z=0.57) = 2094 +/- 34 Mpc. At 1.7 per cent, this is the most precise distance constraint ever obtained from a galaxy survey. We place this result alongside previous BAO measurements in a cosmological distance ladder and find excellent agreement with the current supernova measurements. We use these distance measurements to constrain various cosmological models, finding continuing support for a flat Universe with a cosmological constant.

preprint2012arXiv

Thinking Outside the Box: Effects of Modes Larger than the Survey on Matter Power Spectrum Covariance

Considering the matter power spectrum covariance matrix, it has recently been found that there is a potentially dominant effect on mildly non-linear scales due to power in modes of size equal to and larger than the survey volume. This {\it beat coupling} effect has been derived analytically in perturbation theory and while it has been tested with simulations, some questions remain unanswered. Moreover, there is an additional effect of these large modes, which has so far not been included in analytic studies, namely the effect on the estimated {\it average} density which enters the power spectrum estimate. In this article, we work out analytic, perturbation theory based expressions including both the beat coupling and this {\it local average effect} and we show that while, when isolated, beat coupling indeed causes large excess covariance in agreement with the literature, in a realistic scenario this is compensated almost entirely by the local average effect, leaving only $\sim 10 %$ of the excess. We test our analytic expressions by comparison to a suite of large N-body simulations. For the variances, we find excellent agreement with the analytic expressions for $k < 0.2 h$Mpc$^{-1}$ at $z=0.5$, while the correlation coefficients agree to beyond $k=0.4 h$Mpc$^{-1}$. As expected, the range of agreement increases towards higher redshift and decreases slightly towards $z=0$. We finish by including the large-mode effects in a full covariance matrix description for arbitrary survey geometry and confirming its validity using simulations. This may be useful as a stepping stone towards building an actual galaxy (or other tracer's) power spectrum covariance matrix. [abridged]

preprint2011arXiv

A probabilistic comparison of the strength of split, triangle, and quadrilateral cuts (extended version)

We consider mixed integer linear sets defined by two equations involving two integer variables and any number of non-negative continuous variables. The non-trivial valid inequalities of such sets can be classified into split, type 1, type 2, type 3, and quadrilateral inequalities. We use a strength measure of Goemans to analyze the benefit from adding a non-split inequality on top of the split closure. Applying a probabilistic model, we show that the importance of a type 2 inequality decreases with decreasing lattice width, on average. Our results suggest that this is also true for type 3 and quadrilateral inequalities.

preprint2011arXiv

Evidence Feed Forward Hidden Markov Model: A New Type of Hidden Markov Model

The ability to predict the intentions of people based solely on their visual actions is a skill only performed by humans and animals. The intelligence of current computer algorithms has not reached this level of complexity, but there are several research efforts that are working towards it. With the number of classification algorithms available, it is hard to determine which algorithm works best for a particular situation. In classification of visual human intent data, Hidden Markov Models (HMM), and their variants, are leading candidates. The inability of HMMs to provide a probability in the observation to observation linkages is a big downfall in this classification technique. If a person is visually identifying an action of another person, they monitor patterns in the observations. By estimating the next observation, people have the ability to summarize the actions, and thus determine, with pretty good accuracy, the intention of the person performing the action. These visual cues and linkages are important in creating intelligent algorithms for determining human actions based on visual observations. The Evidence Feed Forward Hidden Markov Model is a newly developed algorithm which provides observation to observation linkages. The following research addresses the theory behind Evidence Feed Forward HMMs, provides mathematical proofs of their learning of these parameters to optimize the likelihood of observations with a Evidence Feed Forwards HMM, which is important in all computational intelligence algorithm, and gives comparative examples with standard HMMs in classification of both visual action data and measurement data; thus providing a strong base for Evidence Feed Forward HMMs in classification of many types of problems.

preprint2011arXiv

Mass Function Predictions Beyond LCDM

The mass distribution of halos, as specified by the halo mass function, is a key input for several cosmological probes. The sizes of $N$-body simulations are now such that, for the most part, results need no longer be statistics-limited, but are still subject to various systematic uncertainties. We investigate and discuss some of the reasons for these differences. Quantifying error sources and compensating for them as appropriate, we carry out a high-statistics study of dark matter halos from 67 $N$-body simulations to investigate the mass function and its evolution for a reference $Λ$CDM cosmology and for a set of $w$CDM cosmologies. For the reference $Λ$CDM cosmology (close to WMAP5), we quantify the breaking of universality in the form of the mass function as a function of redshift, finding an evolution of as much as 10% away from the universal form between redshifts $z=0$ and $z=2$. For cosmologies very close to this reference we provide a fitting formula to our results for the (evolving) $Λ$CDM mass function over a mass range of $6\cdot 10^{11}-3\cdot 10^{15}$ M$_{\odot}$ to an estimated accuracy of about 2%. The set of $w$CDM cosmologies is taken from the Coyote Universe simulation suite. The mass functions from this suite (which includes a $Λ$CDM cosmology and others with $w\simeq-1$) are described by the fitting formula for the reference $Λ$CDM case at an accuracy level of 10%, but with clear systematic deviations. We argue that, as a consequence, fitting formulae based on a universal form for the mass function may have limited utility in high precision cosmological applications.

preprint2011arXiv

Maximal lattice-free polyhedra: finiteness and an explicit description in dimension three

A convex set with nonempty interior is maximal lattice-free if it is inclusion-maximal with respect to the property of not containing integer points in its interior. Maximal lattice-free convex sets are known to be polyhedra. The precision of a rational polyhedron $P$ in $\mathbb{R}^d$ is the smallest integer $s$ such that $sP$ is an integral polyhedron. In this paper we show that, up to affine mappings preserving $\mathbb{Z}^d$, the number of maximal lattice-free rational polyhedra of a given precision $s$ is finite. Furthermore, we present the complete list of all maximal lattice-free integral polyhedra in dimension three. Our results are motivated by recent research on cutting plane theory in mixed-integer linear optimization.

preprint2011arXiv

N-body simulations with generic non-Gaussian initial conditions II: Halo bias

We present N-body simulations for generic non-Gaussian initial conditions with the aim of exploring and modelling the scale-dependent halo bias. This effect is evident at very large scales requiring large simulation boxes. In addition, the previously available prescription to implement generic non-Gaussian initial conditions has been improved to keep under control higher-order terms which were spoiling the power spectrum on large scales. We pay particular attention to the differences between physical, inflation-motivated primordial bispectra and their factorizable templates, and to the operational definition of the non-Gaussian halo bias (which has both a scale-dependent and an approximately scale-independent contributions). We find that analytic predictions for both the non-Gaussian halo mass function and halo bias work well once a calibration factor (which was introduced before) is calibrated on simulations. The halo bias remains therefore an extremely promising tool to probe primordial non-Gaussianity and thus to give insights into the physical mechanism that generated the primordial perturbations. The simulation outputs and tables of the analytic predictions will be made publicly available via the non-Gaussian comparison project web site http://icc.ub.edu/~liciaverde/NGSCP.html

preprint2011arXiv

Stimulation of human red blood cells leads to Ca2+-mediated intercellular adhesion

Red blood cells (RBCs) are a major component of blood clots, which form physiologically as a response to injury or pathologically in thrombosis. The active participation of RBCs in thrombus solidification has been previously proposed but not yet experimentally proven. Holographic optical tweezers and single-cell force spectroscopy were used to study potential cell-cell adhesion between RBCs. Irreversible intercellular adhesion of RBCs could be induced by stimulation with lysophosphatidic acid (LPA), a compound known to be released by activated platelets. We identified Ca2+ as an essential player in the signaling cascade by directly inducing Ca2+ influx using A23187. Elevation of the internal Ca2+ concentration leads to an intercellular adhesion of RBCs similar to that induced by LPA stimulation. Using single-cell force spectroscopy, the adhesion of the RBCs was identified to be approximately 100 pN, a value large enough to be of significance inside a blood clot or in pathological situations like the vasco-occlusive crisis in sickle cell disease patients.

preprint2010arXiv

Inequalities for the lattice width of lattice-free convex sets in the plane

A closed, convex set $K$ in $\mathbb{R}^2$ with non-empty interior is called lattice-free if the interior of $K$ is disjoint with $\mathbb{Z}^2$. In this paper we study the relation between the area and the lattice width of a planar lattice-free convex set in the general and centrally symmetric case. A correspondence between lattice width on the one hand and covering minima on the other, allows us to reformulate our results in terms of covering minima introduced by Kannan and Lovász. We obtain a sharp upper bound for the area for any given value of the lattice width. The lattice-free convex sets satisfying the upper bound are characterized. Lower bounds are studied as well. Parts of our results are applied in a paper by the authors and Weismantel for cutting plane generation in mixed integer linear optimization, which was the original inducement for this paper. We further rectify a result of Kannan and Lovász with a new proof.

preprint2010arXiv

N-body simulations with generic non-Gaussian initial conditions I: Power Spectrum and halo mass function

We address the issue of setting up generic non-Gaussian initial conditions for N-body simulations. We consider inflationary-motivated primordial non-Gaussianity where the perturbations in the Bardeen potential are given by a dominant Gaussian part plus a non-Gaussian part specified by its bispectrum. The approach we explore here is suitable for any bispectrum, i.e. it does not have to be of the so-called separable or factorizable form. The procedure of generating a non-Gaussian field with a given bispectrum (and a given power spectrum for the Gaussian component) is not univocal, and care must be taken so that higher-order corrections do not leave a too large signature on the power spectrum. This is so far a limiting factor of our approach. We then run N-body simulations for the most popular inflationary-motivated non-Gaussian shapes. The halo mass function and the non-linear power spectrum agree with theoretical analytical approximations proposed in the literature, even if they were so far developed and tested only for a particular shape (the local one). We plan to make the simulations outputs available to the community via the non-Gaussian simulations comparison project web site http://icc.ub.edu/~liciaverde/NGSCP.html.

preprint2010arXiv

The Coyote Universe I: Precision Determination of the Nonlinear Matter Power Spectrum

Near-future cosmological observations targeted at investigations of dark energy pose stringent requirements on the accuracy of theoretical predictions for the clustering of matter. Currently, N-body simulations comprise the only viable approach to this problem. In this paper we demonstrate that N-body simulations can indeed be sufficiently controlled to fulfill these requirements for the needs of ongoing and near-future weak lensing surveys. By performing a large suite of cosmological simulation comparison and convergence tests we show that results for the nonlinear matter power spectrum can be obtained at 1% accuracy out to k~1 h/Mpc. The key components of these high accuracy simulations are: precise initial conditions, very large simulation volumes, sufficient mass resolution, and accurate time stepping. This paper is the first in a series of three, with the final aim to provide a high-accuracy prediction scheme for the nonlinear matter power spectrum.

preprint2010arXiv

The Coyote Universe II: Cosmological Models and Precision Emulation of the Nonlinear Matter Power Spectrum

The power spectrum of density fluctuations is a foundational source of cosmological information. Precision cosmological probes targeted primarily at investigations of dark energy require accurate theoretical determinations of the power spectrum in the nonlinear regime. To exploit the observational power of future cosmological surveys, accuracy demands on the theory are at the one percent level or better. Numerical simulations are currently the only way to produce sufficiently error-controlled predictions for the power spectrum. The very high computational cost of (precision) N-body simulations is a major obstacle to obtaining predictions in the nonlinear regime, while scanning over cosmological parameters. Near-future observations, however, are likely to provide a meaningful constraint only on constant dark energy equation of state 'wCDM' cosmologies. In this paper we demonstrate that a limited set of only 37 cosmological models -- the "Coyote Universe" suite -- can be used to predict the nonlinear matter power spectrum at the required accuracy over a prior parameter range set by cosmic microwave background observations. This paper is the second in a series of three, with the final aim to provide a high-accuracy prediction scheme for the nonlinear matter power spectrum for wCDM cosmologies.

preprint2009arXiv

On the starting redshift for cosmological simulations: Focusing on halo properties

We systematically study the effects of varying the starting redshift z_i for cosmological simulations in the highly non-linear regime. Our primary focus lies with the (individual) properties of dark matter halos -- namely the mass, spin, triaxiality, and concentration -- where we find that even substantial variations in z_i leave only a small imprint, at least for the probed mass range M \in [10^{10}, 10^{13}] Msun/h and when investigated at redshift z=0. We further compare simulations started by using the standard Zel'dovich approximation to runs based upon initial conditions produced with second order Lagrangian perturbation theory. Here we observe the same phenomenon, i.e. that differences in the studied (internal) properties of dark matter haloes are practically undetectable. These findings are (for the probed mass range) in agreement with other work in the literature. We therefore conclude that the commonly used technique for setting up cosmological simulations leads to stable results at redshift z=0 for the mass, the spin parameter, the triaxiality, and the concentration of dark matter haloes.

preprint2008arXiv

Dynamics of polymers in elongational flow investigated by inelastic neutron scattering

The nanoscale fluctuation dynamics of semi dilute high molecular weight polymer solutions of Polyethylenoxide (PEO) in D2O under non-equilibrium flow conditions were studied by the neutron spin-echo technique. The sample cell was in contraction flow geometry and provided a pressure driven flow with a high elongational component that stretched the polymers most efficiently. The experiments suggest that the mobility on the scale of a few monomers, comparable to the Kuhn segment length, becomes highly anisotropic and is enhanced perpendicular to the flow direction. Diffraction data show a weak structural correlation along the chains on a length scale of about 17 Angstroems, which might be related to the Kuhn length in this system.