Discontinuous collocation methods and gravitational self-force applications

Numerical simulations of extereme mass ratio inspirals, the mostimportant sources for the LISA detector, face several computational challenges. We present a new approach to evolving partial differential equations occurring in black hole perturbation theory and calculations of the self-force acting on point particles orbiting supermassive black holes. Such equations are distributionally sourced, and standard numerical methods, such as finite-difference or spectral methods, face difficulties associated with approximating discontinuous functions. However, in the self-force problem we typically have access to full a-priori information about the local structure of the discontinuity at the particle. Using this information, we show that high-order accuracy can be recovered by adding to the Lagrange interpolation formula a linear combination of certain jump amplitudes. We construct discontinuous spatial and temporal discretizations by operating on the corrected Lagrange formula. In a method-of-lines framework, this provides a simple and efficient method of solving time-dependent partial differential equations, without loss of accuracy near moving singularities or discontinuities. This method is well-suited for the problem of time-domain reconstruction of the metric perturbation via the Teukolsky or Regge-Wheeler-Zerilli formalisms. Parallel implementations on modern CPU and GPU architectures are discussed.Comment: 29 pages, 5 figure

Inferring the neutron star equation of state from binary inspiral waveforms

The properties of neutron star matter above nuclear density are not precisely known. Gravitational waves emitted from binary neutron stars during their late stages of inspiral and merger contain imprints of the neutron-star equation of state. Measuring departures from the point-particle limit of the late inspiral waveform allows one to measure properties of the equation of state via gravitational wave observations. This and a companion talk by J. S. Read reports a comparison of numerical waveforms from simulations of inspiraling neutron-star binaries, computed for equations of state with varying stiffness. We calculate the signal strength of the difference between waveforms for various commissioned and proposed interferometric gravitational wave detectors and show that observations at frequencies around 1 kHz will be able to measure a compactness parameter and constrain the possible neutron-star equations of state.Comment: Talk given at the 12th Marcel Grossman Meeting, Paris, France, 12-18 Jul 200

Efficient integration of the variational equations of multi-dimensional Hamiltonian systems: Application to the Fermi-Pasta-Ulam lattice

We study the problem of efficient integration of variational equations in multi-dimensional Hamiltonian systems. For this purpose, we consider a Runge-Kutta-type integrator, a Taylor series expansion method and the so-called `Tangent Map' (TM) technique based on symplectic integration schemes, and apply them to the Fermi-Pasta-Ulam $\beta$ (FPU-$\beta$) lattice of $N$ nonlinearly coupled oscillators, with $N$ ranging from 4 to 20. The fast and accurate reproduction of well-known behaviors of the Generalized Alignment Index (GALI) chaos detection technique is used as an indicator for the efficiency of the tested integration schemes. Implementing the TM technique--which shows the best performance among the tested algorithms--and exploiting the advantages of the GALI method, we successfully trace the location of low-dimensional tori.Comment: 14 pages, 6 figure

Noise reduction in protein-protein interaction graphs by the implementation of a novel weighting scheme

<p>Abstract</p> <p>Background</p> <p>Recent technological advances applied to biology such as yeast-two-hybrid, phage display and mass spectrometry have enabled us to create a detailed map of protein interaction networks. These interaction networks represent a rich, yet noisy, source of data that could be used to extract meaningful information, such as protein complexes. Several interaction network weighting schemes have been proposed so far in the literature in order to eliminate the noise inherent in interactome data. In this paper, we propose a novel weighting scheme and apply it to the <it>S. cerevisiae </it>interactome. Complex prediction rates are improved by up to 39%, depending on the clustering algorithm applied.</p> <p>Results</p> <p>We adopt a two step procedure. During the first step, by applying both novel and well established protein-protein interaction (PPI) weighting methods, weights are introduced to the original interactome graph based on the confidence level that a given interaction is a true-positive one. The second step applies clustering using established algorithms in the field of graph theory, as well as two variations of Spectral clustering. The clustered interactome networks are also cross-validated against the confirmed protein complexes present in the MIPS database.</p> <p>Conclusions</p> <p>The results of our experimental work demonstrate that interactome graph weighting methods clearly improve the clustering results of several clustering algorithms. Moreover, our proposed weighting scheme outperforms other approaches of PPI graph weighting.</p

Geological modelling for investigating CO2 emissions in Florina Basin, Greece

Published version also available at http://dx.doi.org/10.1515/geo-2015-0039This paper presents an investigation of naturally occurring CO2 emissions from the Florina natural analogue site in Greece. The main objective was to interpret previously collected depth sounding data, convert them into surfaces, and use them as input to develop, for the rst time, 3D geological models of the Florina basin. By also locating the extent of the aquifer, the location of the CO2 source, the location of other natural CO2 accumulations, and the points where CO2 reaches the surface, we were able to assess the potential for CO2 leakage. Geological models provided an estimate of the lithological composition of the Florina Basin and allowed us to determine possible directions of groundwater ow and pathways of CO2 ow throughout the basin. Important modelling parameters included the spatial positions of boundaries, faults, and major stratigraphic units (which were subdivided into layers of cells). We used various functions in Petrel software to rst construct a structural model describing the main rock boundaries. We then de ned a 3D mesh honouring the structural model, and nally we populated each cell in the mesh with geologic properties, such as rock type and relative permeability. According to the models, the thickest deposits are located around Mesochorion village where we estimate that around 1000 m of sediments were deposited above the basement. Initiation of CO2 ow at Florina Basin could have taken place between 6.5 Ma and 1.8 Ma ago. The NESW oriented faults, which acted as uid ow pathways, are still functioning today, allowing for localised leakage at the surface. CO2 leakage may be spatially variable and episodic in rate. The episodicity can be linked to the timing of Almopia volcanic activity in the area

The clinical evidence and the role of imiquimod in the extramammary Paget disease

Extramammary Paget disease is a rare disease. Different treatment options are used for its management, among them treatment with imiquimod. Our aim is to clarify the efficacy and safety of the use of topical imiquimod cream in the treatment of extramammary Paget disease.In this article the current literature in Pubmed and Scopus on this topic is reviewed.According to the literature the topical use of imiquimod can be used as first-line treatment in different dosages and durations with excellent results, reaching a 87.5% cure rate with rare complications. The available evidence is currently limited to case reports and a single case series, but is definitely in favor of imiquimod use as an alternative treatment for extramammary Paget disease.</p

Neutron star equation of state via gravitational wave observations

Gravitational wave observations can potentially measure properties of neutron star equations of state by measuring departures from the point-particle limit of the gravitational waveform produced in the late inspiral of a neutron star binary. Numerical simulations of inspiraling neutron star binaries computed for equations of state with varying stiffness are compared. As the stars approach their final plunge and merger, the gravitational wave phase accumulates more rapidly if the neutron stars are more compact. This suggests that gravitational wave observations at frequencies around 1 kHz will be able to measure a compactness parameter and place stringent bounds on possible neutron star equations of state. Advanced laser interferometric gravitational wave observatories will be able to tune their frequency band to optimize sensitivity in the required frequency range to make sensitive measures of the late-inspiral phase of the coalescence.Comment: Talk given at the 13th Conference on Recent Developments in Gravity (NEB XIII), Thessaloniki, Greece, 4-6 Jun 200