Spinning neutron stars acquire a quadrupole moment due to their own rotation. This quadratic-in-spin, self-spin effect depends on the equation of state (EOS) and affects the orbital motion and rate of inspiral of neutron star binaries. Building upon circularized post-Newtonian results, we incorporate the EOS-dependent self-spin (or monopole-quadrupole) terms in the spin-aligned effective-one-body (EOB) waveform model TEOBResumS at next-to-next-to-leading (NNLO) order, together with other (bilinear, cubic and quartic) nonlinear-in-spin effects (at leading order, LO). We point out that the structure of the Hamiltonian of TEOBResumS is such that it already incorporates, in the binary black hole case, the recently computed [Levi and Steinhoff, arXiv: 1607.04252] quartic-in-spin LO term. Using the gauge-invariant characterization of the phasing provided by the function Q(omega) = omega(2)/omega of omega = 2 pi f, where f is the gravitational wave frequency, we study the EOS dependence of the self-spin effects and show that: (i) the next-to-leading order (NLO) and NNLO monopole-quadrupole corrections yield increasingly phase-accelerating effects compared to the corresponding LO contribution; (ii) the standard TaylorF2 post-Newtonian (PN) treatment of NLO (3PN) EOS-dependent self-spin effects makes their action stronger than the corresponding EOB description; (iii) the addition to the standard 3PN TaylorF2 post-Newtonian phasing description of self-spin tail effects at LO allows one to reconcile the self-spin part of the TaylorF2 PN phasing with the corresponding TEOBResumS one up to dimensionless frequencies M omega similar or equal to 0.04-0.06. Such a tail-augmented TaylorF2 approximant then yields an analytically simplified, EOB-faithful, representation of the EOS-dependent self-spin phasing that can be useful to improve current PN-based (or phenomenological) waveform models for inspiralling neutron star binaries. Finally, by generating the inspiral dynamics using the post-adiabatic approximation, incorporated in a new implementation of TEOBResumS, one finds that the computational time needed to obtain a typical waveform (including all multipoles up to l = 8) from 10 Hz is of the order of 0.4 sec.

Using standard cylindrical-like coordinates naturally adapted to the cylindrical symmetry of the Godel spacetime, we study elliptic like geodesic motion on hyperplanes orthogonal to the symmetry axis through an eccentricity-semi-latus rectum parametrization which is familiar from the Newtonian description of a two-body system. We compute several quantities which summarize the main features of the motion, namely the coordinate time and proper time periods of the radial motion, the frequency of the azimuthal motion, the full variation of the azimuthal angle over a period, and so on. Exact as well as approximate (i.e., Taylor-expanded in the limit of small eccentricity) analytic expressions of all these quantities are obtained. Finally, we consider their application to the gyroscope precession frequency along these orbits, generalizing the existing results for the circular case.

We investigate the motion of uncharged particles scattered by a binary system consisting of extremely charged black holes in equilibrium as described by the Majumdar-Papapetrou solution. We focus on unbound orbits confined to the plane containing both black holes. We consider the two complementary situations of particles approaching the system along a direction parallel to the axis where the black holes are displaced and orthogonal to it. We numerically compute the scattering angle as a function of the particle's conserved energy parameter, which provides a gauge-invariant information of the scattering process. We also study the precession of a test gyroscope along such orbits and evaluate the accumulated precession angle after a full scattering, which is another gauge-invariant quantity.

The Detweiler-Barack-Sago redshift function for particles moving along slightly eccentric equatorial orbits around a Kerr black hole is currently known up to the second order in eccentricity, second order in spin parameter, and the 8.5 post-Newtonian order. We improve the analytical computation of such a gaugeinvariant quantity by including terms up to the fourth order in eccentricity at the same post-Newtonian approximation level. We also check that our results agrees with the corresponding post-Newtonian expectation of the same quantity, calculated by using the currently known Hamiltonian for spinning binaries.

We analytically compute the gravitational self-force correction to the gyroscope precession along slightly eccentric equatorial orbits in the Kerr spacetime, generalizing previous results for the Schwarzs-child spacetime. Our results are accurate through the 9.5 post-Newtonian order and to second order in both eccentricity and rotation parameter. We also provide a post-Newtonian check of our results based on the currently known Hamiltonian for spinning binaries.

We propose an approach to the numerical simulation of thin-film flows based on the lattice Boltzmann method. We outline the basic features of the method, show in which limits the expected thin-film equations are recovered, and perform validation tests. The numerical scheme is applied to the viscous Rayleigh-Taylor instability of a thin film and to the spreading of a sessile drop toward its equilibrium contact angle configuration. We show that the Cox-Voinov law is satisfied and that the effect of a tunable slip length on the substrate is correctly captured. We address, then, the problem of a droplet sliding on an inclined plane, finding that the Capillary number scales linearly with the Bond number, in agreement with experimental results. At last, we demonstrate the ability of the method to handle heterogenous and complex systems by showcasing the controlled dewetting of a thin film on a chemically structured substrate.

Single-cell RNA-seq (scRNAseq) is a powerful tool to study heterogeneity of cells. Recently, several clustering based methods have been proposed to identify distinct cell populations. These methods are based on different statistical models and usually require to perform several additional steps, such as preprocessing or dimension reduction, before applying the clustering algorithm. Individual steps are often controlled by method-specific parameters, permitting the method to be used in different modes on the same datasets, depending on the user choices. The large number of possibilities that these methods provide can intimidate non-expert users, since the available choices are not always clearly documented. In addition, to date, no large studies have invistigated the role and the impact that these choices can have in different experimental contexts. This work aims to provide new insights into the advantages and drawbacks of scRNAseq clustering methods and describe the ranges of possibilities that are offered to users. In particular, we provide an extensive evaluation of several methods with respect to different modes of usage and parameter settings by applying them to real and simulated datasets that vary in terms of dimensionality, number of cell populations or levels of noise. Remarkably, the results presented here show that great variability in the performance of the models is strongly attributed to the choice of the user-specific parameter settings. We describe several tendencies in the performance attributed to their modes of usage and different types of datasets, and identify which methods are strongly affected by data dimensionality in terms of computational time. Finally, we highlight some open challenges in scRNAseq data clustering, such as those related to the identification of the number of clusters.

Background A subset of individuals affected by imprinting disorders displays multi-locus imprinting disturbances (MLID). MLID has been associated with maternal-effect variants that alter the maintenance of methylation at germline-derived differentially methylated regions (gDMRs) in early embryogenesis. Pedigrees of individuals with MLID also include siblings with healthy phenotype. However, it is unknown if these healthy individuals have MLID themselves or if their methylation patterns differ from those associated with imprinting disorders, and in general, if MLID affects the clinical phenotype. Methods We have investigated gDMR methylation by locus-specific and whole-genome analyses in a family with multiple pregnancy losses, a child with Beckwith-Wiedemann syndrome (BWS) and a further child with no clinical diagnosis of imprinting disorder or other pathologies. Results We detected MLID with different methylation profiles in the BWS-affected and healthy siblings. Whole-exome sequencing demonstrated the presence of novel loss-of-function variants of NLRP5 in compound heterozygosity in the mother. The methylation profiles of the two siblings were compared with those of other cases with MLID and control groups by principal component analysis and unsupervised hierarchical clustering, but while their patterns were clearly separated from those of controls, we were unable to cluster those associated with specific clinical phenotypes among the MLID cases. Conclusion The identification of two novel maternal-effect variants of NLRP5 associated with poly-abortivity and MLID adds further evidence to the role of this gene in the maintenance of genomic imprinting in early embryos. Furthermore, our results demonstrate that within these pedigrees, MLID can also be present in the progeny with healthy phenotype, indicating that some sort of compensation occurs between altered imprinted loci in these individuals. The analysis of larger cohorts of patients with MLID is needed to formulate more accurate epigenotype-phenotype correlations.

We investigate nonlinear elliptic Dirichlet problems whose growth is driven by a general anisotropic N-function, which is not necessarily of power-type and need not satisfy the $\Delta_2$ nor the $\nabla_2$ -condition. Fully anisotropic, non-reflexive Orlicz-Sobolev spaces provide a natural functional framework associated with these problems. Minimal integrability assumptions are detected on the datum on the right-hand side of the equation ensuring existence and uniqueness of weak solutions. When merely integrable, or even measure, data are allowed, existence of suitably further generalized solutions--in the approximable sense--is established. Their maximal regularity in Marcinkiewicz-type spaces is exhibited as well. Uniqueness of approximable solutions is also proved in case of L^1-data.

This paper presents the results of an experiment aiming to measure the vibrational frequencies of the main structures of the medieval church of San Domenico (Matera, southern Italy) and relate them to the mechanical properties of geological stratigraphy and construction materials. Vibrational frequencies are measured by means of the ground-based radar inteferometry technique using a Ku-band radar. Time series of ground-based radar data are processed to measure displacements and vibration frequencies of the church structures. Data collected by a seismic triaxal station for the measurement of the H/V ratio are also processed to compare with radar-based frequencies measurements, providing information about the main vibration frequency ranges of the church structures and, more important, disentangle the site and structure vibration frequencies.

We develop and study stability properties of a hybrid approximation of functionals of the Bates jump model with stochastic interest rate that uses a tree method in the direction of the volatility and the interest rate and a finite-difference approach in order to handle the underlying asset price process. We also propose hybrid simulations for the model, following a binomial tree in the direction of both the volatility and the interest rate, and a space-continuous approximation for the underlying asset price process coming from a Euler-Maruyama type scheme. We test our numerical schemes by computing European and American option prices.

We present the first analytical computation of the (conservative) gravitational self-force correction to the periastron advance around a spinning black hole. Our result is accurate to the second order in the rotational parameter and through the 9.5 post-Newtonian level. It has been obtained as the circular limit of the correction to the gyroscope precession invariant along slightly eccentric equatorial orbits in the Kerr spacetime. The latter result is also new and we anticipate here the first few terms only of the corresponding post-Newtonian expansion.

In this paper we present som eremarks about the Hilbert transform on the real line and its numerical approximation, in connection with its application in signal processing.

To investigate the effects of Glatiramer Acetate (GA) on B cells by an integrated computational and experimental approach. GA is an immunomodulatory drug approved for the treatment of multiple sclerosis (MS). GA effect on B cells is yet to be fully elucidated. We compared transcriptional profiles of B cells from treatment-naive relapsing remitting MS patients, treated or not with GA for 6 hours in vitro, and of B cells before and after six months of GA administration in vivo. Microarrays were analyzed with two different computational approaches, one for functional analysis of pathways (Gene Set Enrichment Analysis) and one for the identification of new drug targets (Mode-of-action by Network Analysis). GA modulates the expression of genes involved in immune response and apoptosis. A differential expression of genes encoding ion channels, mostly regulating Ca2+ homeostasis in endoplasmic reticulum (ER) was also observed. Microfluorimetric analysis confirmed this finding, showing a specific GA effect on ER Ca2+ concentration. Our findings unveils a GA regulatory effect on the immune response by influencing B cell phenotype and function. In particular, our results highlight a new functional role for GA in modulating Ca2+ homeostasis in these cells.

Human mesenchymal/stromal stem cells (hMSC) are the most promising cell source for adult cell therapies in regenerative medicine. Many clinical trials have reported the use of autologous transplantation of hMSCs in several disorders, but with limited results. To exert their potential, hMSCs could exhibit efficient homing and migration toward lesion sites among other effects, but the underlying process is not clear enough. To further increase the knowledge, we studied the co-regulation between hypoxia-regulated genes and miRNAs. To this end, we investigated the miRNA expression profile of healthy hMSCs in low oxygen/nutrient conditions to mimic ischemia and compared with cells of patients suffering from critical limb ischemia (CLI). miRNAs are small, highly conserved, non-coding RNAs, skilled in the control of the target's expression level in a fine-tuned way. After analyzing the miRNOme in CLI-derived hMSC cells and healthy controls, and intersecting the results with the mRNA expression dataset under hypoxic conditions, we identified two miRNAs potentially relevant to the disease: miR-29b as a pathological marker of the disease and miR-638 as a therapeutic target. This study yielded a deeper understanding of stem cell biology and ischemic disorders, opening new potential treatments in the future.

Acute myeloid leukemia (AML) arises from a complex sequence of biological and finely orchestrated events that are still poorly understood. Increasingly, epigenetic studies are providing exciting findings that may be exploited in promising and personalized cutting-edge therapies. A more appropriate and broader screening of possible players in cancer could identify a master molecular mechanism in AML. Here, we build on our previously published study by evaluating a histone deacetylase (HDAC)2-mediated miRNA regulatory network in U937 leukemic cells. Following a comparative miRNA profiling analysis in genetically and enzymatically HDAC2-downregulated AML cells, we identified miR-96-5p and miR-92a-3p as potential regulators in AML etiopathology by targeting defined genes. Our findings support the potentially beneficial role of alternative physiopathological interventions.

The MagnetoEncephaloGraphy (MEG) has gained great interest in neurorehabilitation training due to its high temporal resolution. The challenge is to localize the active regions of the brain in a fast and accurate way. In this paper we use an inversion method based on random spatial sampling to solve the real-time MEG inverse problem. Several numerical tests on synthetic but realistic data show that the method takes just a few hundredths of a second on a laptop to produce an accurate map of the electric activity inside the brain. Moreover, it requires very little memory storage. For these reasons the random sampling method is particularly attractive in real-time MEG applications.

Magnetoencephalography (MEG) aims at reconstructing the unknown neuroelectric activity in the brain from non-invasive measurements of the magnetic field induced by neural sources. The solution of this ill-posed, ill-conditioned inverse problem is usually dealt with using regularization techniques that are often time-consuming, and computationally and memory storage demanding. In this paper we analyze how a slimmer procedure, random sampling, affects the estimation of the brain activity generated by both synthetic and real sources.