We investigate analytically and numerically the orbits of spinning particles around black holes in the post Newtonian limit and in the presence of cosmic expansion. We show that orbits that are circular in the absence of spin, get deformed when the orbiting particle has spin. We show that the origin of this deformation is twofold: a. the background expansion rate which induces an attractive (repulsive) interaction due to the cosmic background fluid when the expansion is decelerating (accelerating) and b. a spin-orbit interaction which can be attractive or repulsive depending on the relative orientation between spin and orbital angular momentum and on the expansion rate.

A recent analysis by one of the authors\cite{Perivolaropoulos:2016ucs} has indicated the presence of a $2\sigma$ signal of spatially oscillating new force residuals in the torsion balance data of the Washington experiment. We extend that study and analyse the data of the Stanford Optically Levitated Microsphere Experiment (SOLME) \cite{Rider:2016xaq} (kindly provided by the authors of \cite{Rider:2016xaq}) searching for sub-mm spatially oscillating new force signals. We find a statistically significant oscillating signal for a force residual of the form $F(z)=\alpha \; cos(\frac{2\pi}{\lambda}\; z +c)$ where $z$ is the distance between the macroscopic interacting masses (levitated microsphere and cantilever). The best fit parameter values are $\alpha=(1.1 \pm 0.4)\times 10^{-17}N$, $\lambda=(35.2\pm 0.6)\mu m$. Monte Carlo simulation of the SOLME data under the assumption of zero force residuals has indicated that the statistical significance of this signal is at about $2\sigma$ level. Private communication with members of the SOLME group has indicated that this previously unnoticed signal is most probably due to a systematic effect (diffraction of non-Gaussian tails of the laser beam). Thus it can only be useful as an upper bound of the amplitude of new spatially oscillating forces on sub-mm scales. Assuming gravitational origin emerging from a fundamental modification of the Newtonian potential of the form $V_{eff}(r)=-G\frac{M}{r}(1+\alpha_O \cos(\frac{2\pi}{\lambda}\; r+\theta))\equiv V_{N}(r)+V_{osc}(r)$, we evaluate the source integral and obtain a bound $\alpha_O < 10^7$ for $\lambda \simeq 35 \mu m$. Thus, we get no useful constraints on the modified gravity parameter $\alpha_O$. However, the constraints on the general phenomenological parameter $\alpha$ can be useful in constraining other fifth force models related to dark energy (chameleon oscillating potentials etc).

We review a possible non-minimal coupling (dilatonic) of a scalar field (axion like particle) to electromagnetism, through experimental and observational constraints. Such a coupling is motivated from recent quasar spectrum observations that indicate a possible spatial and/or temporal variation of the fine-structure constant. We consider a dilatonic coupling of the form $B_F (\phi) = 1 + g\phi$. The strongest bound on the coupling parameter $g$ is derived from weak equivalence principle tests, which impose $g < 1.6 \times 10^{-17}GeV^{-1}$. This constraint is strong enough to rule out this class of models as a cause for an observable cosmological variation of the fine structure constant unless a chameleon mechanism is implemented. Also, we argue that a similar coupling occurs in chameleon cosmology, another candidate dark mater particle and we estimate the cosmological consequences by both effects. It should be clarified that this class of models is not necessarily ruled out in the presence of a chameleon mechanism which can freeze the dynamics of the scalar field in high density laboratory regions.

We investigate the geodesics of a Schwarzschild spacetime embedded in an isotropic expanding cosmological background (McVittie metric). We focus on bound particle geodesics in a background including matter and phantom dark energy with constant dark energy equation of state parameter $w<-1$ involving a future Big Rip singularity at a time $t_{\ast}$. Such geodesics have been previously studied in the Newtonian approximation and found to lead to dissociation of bound systems at a time $t_{rip}<t_{\ast}$ which for fixed background $w$, depends on a single dimensionless parameter $\bar{\omega}_0$ related to the angular momentum and depending on the mass and the size of the bound system. We extend this analysis to large massive bound systems where the Newtonian approximation is not appropriate and we compare the derived dissociation time with the corresponding time in the context of the Newtonian approximation. By identifying the time when the effective potential minimum disappears due to the repulsive force of dark energy we find that the dissociation time of bound systems occurs earlier than the prediction of the Newtonian approximation. However, the effect is negligible for all existing cosmological bound systems and it would become important only in hypothetical bound extremely massive ($10^{20} M_\odot$) and large ($100Mpc$) bound systems. We verify this result by explicit solution of the geodesic equations. This result is due to an interplay between the repulsive phantom dark energy effects and the existence of the well known innermost stable orbits of Schwarzschild spacetimes.

We consider Lagrangians in 3+1 dimensions admitting topological defects where there is an additional coupling between the defect scalar field $\Phi$ and the gauge field kinetic term (eg $B(\vert \Phi \vert^2) F_{\mu \nu}F^{\mu \nu}$). Such a {\it dilatonic} coupling in the context of a static defect, induces a spatially dependent effective gauge charge and effective mass for the scalar field which leads to modified properties of the defect core. In particular, the scale of the core gets modified while the stability properties of the corresponding embedded defects are also affected. These modifications are illustrated for gauged (Nielsen-Olesen) vortices and for gauged ('t Hooft-Polyakov) monopoles. The corresponding dilatonic global defects are also studied in the presence of an external gauge field.

We review, compare and extend recent studies searching for evidence for a preferred cosmological axis. We start from the Union2 SnIa dataset and use the hemisphere comparison method to search for a preferred axis in the data. We find that the hemisphere of maximum accelerating expansion rate is in the direction $(l,b)=({309^\circ}^{+23^\circ}_{-3^\circ}, {18^\circ}^{+11^\circ}_{-10^\circ})$ ($\omm=0.19$) while the hemisphere of minimum acceleration is in the opposite direction $(l,b)=({129^\circ}^{+23^\circ}_{-3^\circ},{-18^\circ}^{+10^\circ}_{-11^\circ})$ ($\omm=0.30$). The level of anisotropy is described by the normalized difference of the best fit values of $\omm$ between the two hemispheres in the context of \lcdm fits. We find a maximum anisotropy level in the Union2 data of $\frac{\Delta \ommax}{\bomm}=0.43\pm 0.06$. Such a level does not necessarily correspond to statistically significant anisotropy because it is reproduced by about $30%$ of simulated isotropic data mimicking the best fit Union2 dataset. However, when combined with the axes directions of other cosmological observations (bulk velocity flow axis, three axes of CMB low multipole moments and quasar optical polarization alignment axis), the statistical evidence for a cosmological anisotropy increases dramatically. We estimate the probability that the above independent six axes directions would be so close in the sky to be less than $1%$. Thus either the relative coincidence of these six axes is a very large statistical fluctuation or there is an underlying physical or systematic reason that leads to their correlation.