We investigate the accuracy of weak lensing simulations by comparing the results of five independently developed lensing simulation codes run on the same input $N$-body simulation. Our comparison focuses on the lensing convergence maps produced by the codes, and in particular on the corresponding PDFs, power spectra and peak counts. We find that the convergence power spectra of the lensing codes agree to $\lesssim 2\%$ out to scales $\ell \approx 4000$. For lensing peak counts, the agreement is better than $5\%$ for peaks with signal-to-noise $\lesssim 6$. We also discuss the systematic errors due to the Born approximation, line-of-sight discretization, particle noise and smoothing. The lensing codes tested deal in markedly different ways with these effects, but they nonetheless display a satisfactory level of agreement. Our results thus suggest that systematic errors due to the operation of existing lensing codes should be small. Moreover their impact on the convergence power spectra for a lensing simulation can be predicted given its numerical details, which may then serve as a validation test.

Observations show that stellar streams originating in satellite dwarf galaxies are frequent in the Universe. While such events are predicted by theory, it is not clear how many of the streams that are generated are washed out afterwards to the point in which it is imposible to detect them. Here we study how these diffusion times are affected by the fact that typical gravitational potentials of the host galaxies can sustain chaotic orbits. We do this by comparing the behaviour of simulated stellar streams that reside in chaotic or non-chaotic regions of the phase-space. We find that chaos does reduce the time interval in which streams can be detected. By analyzing detectability criteria in configuration and velocity space, we find that the impact of these results on the observations depends on the quality of both the data and the underlying stellar halo model. For all the stellar streams, we obtain a similar upper limit to the detectable mass.

There are currently no reliable methods to measure transverse velocities of galaxies. This is an important piece of information that could allow us to probe the physics of structure formation as well as testing the underlying theory of gravity. The slingshot effect, a special case of the Integrated Sachs-Wolfe effect, is expected to create dipole signals in the temperature fluctuations of the Cosmic Microwave Background Radiation (CMB). This effect creates a hot spot behind and a cold spot in front of moving massive objects. The dipole signal created by the slingshot effect can be used to measure transverse velocities, but because the signal is expected to be weak, the effect has not been measured yet. The aim is to show that the slingshot effect can be measured by stacking the signals of galaxies falling into a collapsing cluster. We evaluate if the effect can probe modified gravity. We use data from a simulated galaxy catalogue (MDPL2) to mimic observations. We identify a massive galaxy cluster, and make maps of the slingshot effect around infalling galaxies. We add uncorrelated Gaussian noise to each map. The maps are rotated according to the direction to the cluster centre, such that the dipole signal will add up constructively when stacking. We compare each stack to a dipole stencil and we find the probability for a false positive in the absence of the slingshot signal. Each galaxy gives a signal of around $\Delta T/T \approx 10^{-9}$, while the precision of CMB experiments of today are $\Delta T/T \approx 4 \times 10^{-6}$. By stacking around 10 000 galaxies, the slingshot signal can be over the detectable threshold with experiments of today. However, future CMB experiments must be used to be certain of the strength of the observed signal.

The Quintessence model is one of the simplest and better known alternatives to Einstein's theory for gravity. The properties of the solutions have been studied in great detail in the background, linear and non-linear contexts in cosmology. Here we discuss new phenomenology that is induced by adding disformal terms to the interactions. Among other results, we show analytically and using cosmological simulations ran with the code \texttt{Isis} that the model posses a mechanism through which is it possible to obtain repulsive fifth forces, which are opposite to gravity. Although the equations are very complex, we also find that most of the new phenomenology can be explained by studying background quantities. We used our simulation data to test approximate relations that exist between the metric and scalar field perturbations as well as between the fifth force and gravity. Excellent agreement was found between exact and approximated solutions, which opens the way for running disformal gravity cosmological simulations using simply a Newtonian solver. These results could not only help us to find new ways of testing gravity, but also provide new motivations for building alternative models.

The inherently unstable nature of domain walls makes their detection in laboratory experiments extremely challenging. We propose a method to stabilise domain walls in a particular modified gravity model inside a cavity. We suggest two ways in which the walls could be detected once stabilized: studying the trajectories of Ultra Cold Neutrons (UCN's) either via the deflection angle of a neutron beam induced by the attraction towards the wall or through the time difference of these particles passing through the wall. We give realistic estimates for these effects and expect that they should be detectable experimentally.

Most statistical inference from cosmic large-scale structure relies on two-point statistics, i.e.\ on the galaxy-galaxy correlation function (2PCF) or the power spectrum. These statistics capture the full information encoded in the Fourier amplitudes of the galaxy density field but do not describe the Fourier phases of the field. Here, we quantify the information contained in the line correlation function (LCF), a three-point Fourier phase correlation function. Using cosmological simulations, we estimate the Fisher information (at redshift $z=0$) of the 2PCF, LCF and their combination, regarding the cosmological parameters of the standard $\Lambda$CDM model, as well as a Warm Dark Matter (WDM) model and the $f(R)$ and Symmetron modified gravity models. The galaxy bias is accounted for at the level of a linear bias. The relative information of the 2PCF and the LCF depends on the survey volume, sampling density (shot noise) and the bias uncertainty. For a volume of $1h^{-3}\rm Gpc^3$, sampled with points of mean density $\bar{n} = 2\times10^{-3} h^{3}\ \rm Mpc^{-3}$ and a bias uncertainty of 13\%, the LCF improves the parameter constraints by about 20\% in the $\Lambda$CDM cosmology and potentially even more in alternative models. Finally, since a linear bias only affects the Fourier amplitudes (2PCF), but not the phases (LCF), the combination of the 2PCF and the LCF can be used to break the degeneracy between the linear bias and $\sigma_8$, present in 2-point statistics.

The dynamics of globular clusters has been studied in great detail in the context of general relativity as well as with modifications of gravity that strongly depart from the standard paradigm such as MOND. However, at present there are no studies that aim to test the impact that less extreme modifications of gravity (e.g. models constructed as alternatives to dark energy) have on the behaviour of globular clusters. This Letter presents fits to the velocity dispersion profile of the cluster NGC 2419 under the symmetron modified gravity model. The data shows an increase in the velocity dispersion towards the centre of the cluster which could be difficult to explain within general relativity. By finding the best fitting solution associated with the symmetron model, we show that this tension does not exist in modified gravity. However, the best fitting parameters give a model that is inconsistent with the dynamics of the Solar System. Exploration of different screening mechanisms should give us the chance to understand if it is possible to maintain the appealing properties of the symmetron model when it comes to globular clusters and at the same time recover the Solar System dynamics properly.

The advent of the new generation of wide field galaxy surveys makes large N-body cosmological simulations a necessary evil. While the cosmological simulation codes have evolved a lot since the first calculations in the 80s, the computational requirements for generating data that is relevant for large surveys remain challenging. We propose an alternative approach that can speed up these simulations. The framework is based on the idea of reducing the size of the integration region following the lightcone of an observer at redshift zero and thus simulating only the parts of the Universe that can be observed. A possible implementation of this framework is presented, as well as tests of its accuracy and performance. These simple tests, based on conservative assumptions, show that the new framework gives a factor of three speed up with respect to the usual approach.

When it comes to searches for extensions to general relativity, large efforts are being dedicated to accurate predictions for the power spectrum of density perturbations. While this observable is known to be sensitive to the gravitational theory, its efficiency as a diagnostic for gravity is significantly reduced when Solar System constraints are strictly adhered to. We show that this problem can be overcome by studying weigthed density fields. We propose a transformation of the density field for which the impact of modified gravity on the power spectrum can be increased by more than a factor of three. The signal is not only amplified, but the modified gravity features are shifted to larger scales which are less affected by baryonic physics. Furthermore, the overall signal-to-noise increases, which in principle makes identifying signatures of modified gravity with future galaxy surveys more feasible. While our analysis is focused on modified gravity, the technique can be applied to other problems in cosmology, such as the detection of neutrinos, the effects of baryons or baryon acoustic oscillations.

Extending general relativity by adding extra degrees of freedom is a popular approach for explaining the accelerated expansion of the Universe and to build high energy completions of the theory of gravity. The presence of such new degrees of freedom is, however, tightly constrained from several observations and experiments that aim to test general relativity in a wide range of scales. The viability of a given modified theory of gravity, therefore, strongly depends on the existence of a screening mechanism that suppresses the extra degrees of freedom. We perform simulations, and find that waves propagating in the new degrees of freedom can significantly impact the efficiency of some screening mechanisms, thereby threatening the viability of these modified gravity theories. Specifically, we show that the waves produced in the symmetron model can increase the amplitude of the fifth force and the parametrized post-Newtonian parameters by several orders of magnitude.

We describe and demonstrate the potential of a new and very efficient method for simulating certain classes of modified gravity theories, such as the widely studied $f(R)$ gravity models. High resolution simulations for such models are currently very slow due to the highly nonlinear partial differential equation that needs to be solved exactly to predict the modified gravitational force. This nonlinearity is partly inherent, but is also exacerbated by the specific numerical algorithm used, which employs a variable redefinition to prevent numerical instabilities. The standard Newton-Gauss-Seidel iterative method used to tackle this problem has a poor convergence rate. Our new method not only avoids this, but also allows the discretised equation to be written in a form that is analytically solvable. We show that this new method greatly improves the performance and efficiency of $f(R)$ simulations. For example, a test simulation with $512^3$ particles in a box of size $512 \, \mathrm{Mpc}/h$ is now 5 times faster than before, while a Millennium-resolution simulation for $f(R)$ gravity is estimated to be more than 20 times faster than with the old method. Our new implementation will be particularly useful for running very high resolution, large-sized simulations which, to date, are only possible for the standard model, and also makes it feasible to run large numbers of lower resolution simulations for covariance analyses. We hope that the method will bring us to a new era for precision cosmological tests of gravity.

We study the imprints that theories of gravity beyond GR can leave on the lensing signal around line of sight directions that are predominantly halo-underdense (called troughs) and halo-overdense. To carry out our investigations, we consider the normal branch of DGP gravity, as well as a phenomenological variant thereof that directly modifies the lensing potential. The predictions of these models are obtained with N-body simulation and ray-tracing methods using the ECOSMOG and Ray-Ramses codes. We analyse the stacked lensing convergence profiles around the underdense and overdense lines of sight, which exhibit, respectively, a suppression and a boost w.r.t. the mean in the field of view. The modifications to gravity in these models strengthen the signal w.r.t. $\Lambda{\rm CDM}$ in a scale-independent way. We find that the size of this effect is the same for both underdense and overdense lines of sight, which implies that the density field along the overdense directions on the sky is not sufficiently evolved to trigger the suppression effects of the screening mechanism. These results are robust to variations in the minimum halo mass and redshift ranges used to identify the lines of sight, as well as to different line of sight aperture sizes and criteria for their underdensity and overdensity thresholds.

Alternative theories of gravity typically invoke an environment-dependent screening mechanism to allow phenomenologically interesting deviations from general relativity (GR) to manifest on larger scales, while reducing to GR on small scales. The observation of the transition from screened to unscreened behavior would be compelling evidence for beyond-GR physics. We show that pairwise peculiar velocity statistics, in particular the relative radial velocity dispersion, $\sigma_\parallel$ , can be used to observe this transition when they are binned by some measure of halo environment. We established this by measuring the radial velocity dispersion between pairs of halos in N-body simulations for three $f(R)$ gravity and four symmetron models. We developed an estimator involving only line-of-sight velocities to show that this quantity is observable, and binned the results in halo mass, ambient density, and the isolatedness of halos. Ambient density is found to be the most relevant measure of environment; it is distinct from isolatedness, and correlates well with theoretical expectations for the symmetron model. By binning $\sigma_\parallel$ in ambient density, we find a strong environment-dependent signature for the symmetron models, with the velocities showing a clear transition from GR to non-GR behavior. No such transition is observed for $f(R)$, as the relevant scales are deep in the unscreened regime. Observations of the relative radial velocity dispersion in forthcoming peculiar velocity surveys, if binned appropriately by environment, therefore offer a valuable way of detecting the screening signature of modified gravity.

We use a spherical model and an extended excursion set formalism with drifting diffusive barriers to predict the abundance of cosmic voids in the context of general relativity as well as f(R) and symmetron models of modified gravity. We detect spherical voids from a suite of N-body simulations of these gravity theories and compare the measured void abundance to theory predictions. We find that our model correctly describes the abundance of both dark matter and galaxy voids, providing a better fit than previous proposals in the literature based on static barriers. We use the simulation abundance results to fit for the abundance model free parameters as a function of modified gravity parameters, and show that counts of dark matter voids can provide interesting constraints on modified gravity. For galaxy voids, more closely related to optical observations, we find that constraining modified gravity from void abundance alone may be significantly more challenging. In the context of current and upcoming galaxy surveys, the combination of void and halo statistics including their abundances, profiles and correlations should be effective in distinguishing modified gravity models that display different screening mechanisms.

We present a ray tracing code to compute integrated cosmological observables on the fly in AMR N-body simulations. Unlike conventional ray tracing techniques, our code takes full advantage of the time and spatial resolution attained by the N-body simulation by computing the integrals along the line of sight on a cell-by-cell basis through the AMR simulation grid. Moroever, since it runs on the fly in the N-body run, our code can produce maps of the desired observables without storing large (or any) amounts of data for post-processing. We implemented our routines in the RAMSES N-body code and tested the implementation using an example of weak lensing simulation. We analyse basic statistics of lensing convergence maps and find good agreement with semi-analytical methods. The ray tracing methodology presented here can be used in several cosmological analysis such as Sunyaev-Zel'dovich and integrated Sachs-Wolfe effect studies as well as modified gravity. Our code can also be used in cross-checks of the more conventional methods, which can be important in tests of theory systematics in preparation for upcoming large scale structure surveys.

We present measurements of the spatial clustering statistics in redshift space of various scalar field modified gravity simulations. We utilise the two-point and the three-point correlation functions to quantify the spatial distribution of dark matter halos within these simulations and thus discern between the models. We compare $\Lambda$CDM simulations to various modified gravity scenarios and find consistency with previous work in terms of 2-point statistics in real and redshift-space. However using higher order statistics such as the three-point correlation function in redshift space we find significant deviations from $\Lambda$CDM hinting that higher order statistics may prove to be a useful tool in the hunt for deviations from General Relativity.

Despite its continued observational successes, there is a persistent (and growing) interest in extending cosmology beyond the standard model, $\Lambda$CDM. This is motivated by a range of apparently serious theoretical issues, involving such questions as the cosmological constant problem, the particle nature of dark matter, the validity of general relativity on large scales, the existence of anomalies in the CMB and on small scales, and the predictivity and testability of the inflationary paradigm. In this paper, we summarize the current status of $\Lambda$CDM as a physical theory, and review investigations into possible alternatives along a number of different lines, with a particular focus on highlighting the most promising directions. While the fundamental problems are proving reluctant to yield, the study of alternative cosmologies has led to considerable progress, with much more to come if hopes about forthcoming high-precision observations and new theoretical ideas are fulfilled.

We investigate the propagation of scalar waves induced by matter sources in the context of scalar-tensor theories of gravity which include screening mechanisms for the scalar degree of freedom. The usual approach when studying these theories in the non-linear regime of cosmological perturbations is based on the assumption that scalar waves travel at the speed of light. Within General Relativity such approximation is good and leads to no loss of accuracy in the estimation of observables. We find, however, that mass terms and non-linearities in the equations of motion lead to propagation and dispersion velocities significantly different from the speed of light. As the group velocity is the one associated to the propagation of signals, a reduction of its value has direct impact on the behavior and dynamics of nonlinear structures within modified gravity theories with screening. For instance, the internal dynamics of galaxies and satellites submerged in large dark matter halos could be affected by the fact that the group velocity is smaller than the speed of light. It is therefore important, within such framework, to take into account the fact that different part of a galaxy will see changes in the environment at different times. A full non-static analysis may be necessary under those conditions.

We investigate statistical properties of the distribution of matter at redshift zero in disformal gravity by using N-body simulations. The disformal model studied here consists of a conformally coupled symmetron field with an additional exponential disformal term. We conduct cosmological simulations to discover the impact of the new disformal terms in the matter power spectrum, halo mass function, and radial profile of the scalar field. We calculated the disformal geodesic equation and the equation of motion for the scalar field. We then implemented these equations into the N-body code ISIS, which is a modified gravity version of the code RAMSES. The presence of a conformal symmetron field increases both the power spectrum and mass function compared to standard gravity on small scales. Our main finding is that the newly added disformal terms tend to counteract these effects and can make the evolution slightly closer to standard gravity. We finally show that the disformal terms give rise to oscillations of the scalar field in the centre of the dark matter haloes.

Self-consistent ${\it N}$-body simulations of modified gravity models are a key ingredient to obtain rigorous constraints on deviations from General Relativity using large-scale structure observations. This paper provides the first detailed comparison of the results of different ${\it N}$-body codes for the $f(R)$, DGP, and Symmetron models, starting from the same initial conditions. We find that the fractional deviation of the matter power spectrum from $\Lambda$CDM agrees to better than $1\%$ up to $k \sim 5-10~h/{\rm Mpc}$ between the different codes. These codes are thus able to meet the stringent accuracy requirements of upcoming observational surveys. All codes are also in good agreement in their results for the velocity divergence power spectrum, halo abundances and halo profiles. We also test the quasi-static limit, which is employed in most modified gravity ${\it N}$-body codes, for the Symmetron model for which the most significant non-static effects among the models considered are expected. We conclude that this limit is a very good approximation for all of the observables considered here.

In this paper we present the first results from implementing two scalar-tensor modified gravity theories, the symmetron and the Hu-Sawicki $f(R)$-gravity model, into a hydrodynamic N-body code with dark matter particles and a baryonic ideal gas. The study is a continuation of previous work where the symmetron and $f(R)$ have been successfully implemented in the RAMSES code, but for dark matter only. By running simulations, we show that the deviation from $\Lambda$CDM in these models for the gas density profiles are significantly lower than the dark matter equivalents. When it comes to the matter power-spectrum we find that hydrodynamic simulations agree very well with dark matter only simulations as long as we consider scales larger than $k\sim 0.5$ h/Mpc. In general the effects of modified gravity on the baryonic gas is found to not always mirror the effects it has on the dark matter. The largest signature is found when considering temperature profiles. We find that the gas temperatures in the modified gravity model studied here show deviations, when compared to $\Lambda$CDM, that can be a factor of a few larger than the deviations found in density profiles and power spectra.

Screened modified gravity predicts potentially large signatures in the peculiar velocity field that makes it an interesting probe to test gravity on cosmological scales. We investigate the signatures induced by the Symmetron and a Chameleon $f(R)$ model in the peculiar velocity field using $N$-body simulations. By studying fifth force and halo velocity profiles we identify three general categories of effects found in screened modified gravity models: a fully screened regime where we recover $\Lambda$CDM to high precision, an unscreened regime where the fifth force is in full operation, and, a partially screened regime where screening occurs in the inner part of a halo, but the fifth force is active at larger radii. These three regimes can be pointed out very clearly by analyzing the deviation in the maximum cluster velocity. Observationally, the partially screened regime is of particular interest since an uniform increase of the gravitational force - as present in the unscreened regime - is degenerate with the (dynamical) halo mass estimate, and, thus, hard to detect.

We study properties of domain walls in the symmetron model, in which the scalar gravitational degree of freedom decouples from matter in regions of high density, and exhibits a spontaneously broken $Z_2$ symmetry at low densities. The non-minimal coupling of the scalar to matter leads to a host of interesting properties of the domain walls that are not present in minimally coupled theories. We estimate the cosmological energy fraction in domain walls and find that this leads to an upper bound on the redshift of the symmetry breaking. We also show that a spherical symmetron wall can remain stable if it is "pinned" on matter halos and derive a criterion for the stability. In addition, we present results of numerical simulations of representative interactions between domain walls and matter over-densities.

In the context of scalar tensor theories for gravity, there is a universally adopted hypothesis when running N-body simulations that time derivatives in the equation of motion for the scalar field are negligible. In this work we propose to test this assumption for one specific scalar-tensor model with a gravity screening mechanism: the symmetron. To this end, we implemented the necessary modifications to include the non-static terms in the N-body code Ramses. We present test cases and results from cosmological simulations. Our main finding when comparing static vs. non-static simulations is that the global power spectrum is only slightly modified when taking into account the inclusion of non-static terms. On the contrary, we find that the calculation of the local power spectrum gives different measurements. Such results imply one must be careful when assuming the quasi-static approximation when investigating the environmental effects of modified gravity and screening mechanisms in structure formation of halos and voids distributions.

Scalar fields are crucial components in high energy physics and extensions of General Relativity. The fact they are not observed in the solar system may be due to a mechanism which screens their presence in high dense regions. We show how observations of the ellipticity of galaxy clusters can discriminate between models with and without scalar fields and even between different screening mechanisms. Using nowadays X-ray observations we put novel constraints on the different models.

Several extensions of General Relativity and high energy physics include scalar fields as extra degrees of freedom. In the search for predictions in the non-linear regime of cosmological evolution, the community makes use of numerical simulations in which the quasi-static limit is assumed when solving the equation of motion of the scalar field. In this Letter, we propose a method to solve the full equations of motion for scalar degrees of freedom coupled to matter. We run cosmological simulations which track the full time and space evolution of the scalar field, and find striking differences with respect to the commonly used quasi-static approximation. This novel procedure reveals new physical properties of the scalar field and uncovers concealed astrophysical phenomena which were hidden in the old approach.

Several extensions of the standard cosmological model include scalar fields as new degrees of freedom in the underlying gravitational theory. A particular class of these scalar field theories include screening mechanisms intended to hide the scalar field below observational limits in the solar system, but not on galactic scales, where data still gives freedom to find possible signatures of their presence. In order to make predictions to compare with observations coming from galactic and clusters scales (i.e. in the non-linear regime of cosmological evolution), cosmological N-body simulations are needed, for which codes that can solve for the scalar field must be developed. We present a new implementation of scalar-tensor theories of gravity which include screening mechanisms. The code is based in the already existing code RAMSES, to which we have added a non-linear multigrid solver that can treat a large class of scalar tensor theories of modified gravity. We present details of the implementation and the tests that we made to it. As application of the new code, we have studied the influence that two particular modified gravity theories, the symmetron and $f(R)$ gravity, have on the shape of cluster sized dark matter halos and found consistent results with previous estimations made with a static analysis.

Aims. We investigate the gravitational redshift in clusters of galaxies in the symmetron and Hu-Sawicky $f(R)$ models. The characteristic feature of both models is the screening mechanism that hides the fifth force in dense environments recovering general relativity. Methods. We use $N$-body simulations that were run using the code \texttt{Isis}, which includes scalar fields, to analyse the deviation of observables in modified gravity models with respect to $\Lambda$CDM. Results. We find that due to the presence of the screening, the deviation is highly dependent on the halo mass. For instance, the $f(R)$ parameters $|f_{R0}|=10^{-5}$, $n=1$ cause an enhancement of the gravitational signal by up to 50% for halos with masses between $10^{13}\,M_\odot h^{-1}$ and $10^{14}\,M_\odot h^{-1}$. The characteristic mass range where the fifth force is most active varies with the model parameters. The usual assumption is that the presence of a fifth force leads to a deeper potential well and thus, a stronger gravitational redshift. However, we find that in cases in which only the central regions of the halos are screened, there could also be a weaker gravitational redshift.

The recent years have seen combined measurements of X-ray and (weak) lensing contours for colliding galaxy clusters such as, for instance, the famous "Bullet" cluster. These observations have revealed offsets in the peaks of the baryonic and (dominant) gravitational matter component of order ~(100-200) kpc. Such discrepancies are difficult to explain using modified theories for gravity other than dark matter. Or are they not? Here we use the concept of "phantom dark matter" that is based upon a Newtonian interpretation of the MONDian gravitational potential. We show that this idea is in fact capable of producing substantial offsets in idealistic density configurations, involving a uniform external field. However, when analysed in a MONDian cosmological framework we deduce that the size (and probablity) of the effect is too small to explain the observed offsets found in the most recent observations, at least in the simplest incarnation of phantom dark matter as discussed here. The lensing centers in merging galaxy clusters are likely very close to the centers of true mass even in a MONDian cosmology. This gives the support to the idea that neutrino-like non-collisional matter might be responsible for the observed offsets of lensing and X-ray peaks.

We perform a stability test of triaxial models in MOdified Newtonian Dynamics (MOND) using N-body simulations. The triaxial models considered here have densities that vary with $r^{-1}$ in the center and $r^{-4}$ at large radii. The total mass of the model varies from $10^8\Msun$ to $10^{10}\Msun$, representing the mass scale of dwarfs to medium-mass elliptical galaxies, respectively, from deep MOND to quasi-Newtonian gravity. We build triaxial galaxy models using the Schwarzschild technique, and evolve the systems for 200 Keplerian dynamical times (at the typical length scale of 1.0 kpc). We find that the systems are virial overheating, and in quasi-equilibrium with the relaxation taking approximately 5 Keplerian dynamical times (1.0 kpc). For all systems, the change of the inertial (kinetic) energy is less than 10% (20%) after relaxation. However, the central profile of the model is flattened during the relaxation and the (overall) axis ratios change by roughly 10% within 200 Keplerian dynamical times (at 1.0kpc) in our simulations. We further find that the systems are stable once they reach the equilibrium state.

The statistics of high speed satellite galaxies, as reported in the recent literature, can be a powerful diagnosis of the depth of the potential well of the host halo, and hence discriminate between competing gravitational theories. Naively one expects that high speed satellites are more common in Modified Newtonian Dynamics (MOND) than in cold dark matter (CDM) since an isolated MONDian system has an infinite potential well, while CDM halos have finite potential wells. In this \textit{Letter} we report on an initial test of this hypothesis in the context of the first generation of cosmological simulations utilising a rigorous MONDian Poisson solver. We find that such high speed encounters are approximately a factor of four more common in MOND than in the concordance $\Lambda$CDM model of cosmic structure formation.

The periodicity of 5.5 years for some observational events occurring in Eta Carinae manifests itself across a large wavelength range and has been associated with its binary nature. These events are supposed to occur when the binary components are close to periastron. To detect the previous periastron passage of Eta Car in 2003, we started an intensive, ground-based, optical, photometric observing campaign. We continued observing the object to monitor its photometric behavior and variability across the entire orbital cycle. Our observation program consisted of daily differential photometry from CCD images, which were acquired using a 0.8 m telescope and a standard BVRI filter set at La Plata Observatory. The photometry includes the central object and the surrounding Homunculus nebula. We present up-to-date results of our observing program, including homogeneous photometric data collected between 2003 and 2008. Our observations demonstrated that Eta Car has continued increasing in brightness at a constant rate since 1998. In 2006, it reached its brightest magnitude (V ~ 4.7) since about 1860s. The object then suddenly reverted its brightening trend, fading to V = 5.0 at the beginning of 2007, and has maintained a quite steady state since then. We continue the photometric monitoring of Eta Car in anticipation of the next "periastron passage", predicted to occur at the beginning of 2009.

We present a novel solver for an analogue to Poisson's equation in the framework of modified Newtonian dynamics (MOND). This equation is highly non-linear and hence standard codes based upon tree structures and/or FFT's in general are not applicable; one needs to defer to multi-grid relaxation techniques. After a detailed description of the necessary modifications to the cosmological N-body code AMIGA (formerly known as MLAPM) we utilize the new code to revisit the issue of cosmic structure formation under MOND. We find that the proper (numerical) integration of a MONDian Poisson's equation has some noticable effects on the final results when compared against simulations of the same kind but based upon rather ad-hoc assumptions about the properties of the MONDian force field. Namely, we find that the large-scale structure evolution is faster in our revised MOND model leading to an even stronger clustering of galaxies, especially when compared to the standard LCDM paradigm.

After the Observing Campaign of Eta Carinae carried out during 2003, we have continued the CCD observations using the 0.8 m Reflector telescope at La Plata Observatory, in the BVRI bands. Here, we present the results obtained since November 2003 until August 2004. In addition, we present differential photometry of other objects included in the Eta Car's frames, belonging to the open cluster Trumpler 16.

In 2003.5 Eta Carinae was expected to undergo an X-ray eclipse (Damineli et al., 2000). In the framework of an international campaign to obtain multi-wavelength observations of this event, we have obtained optical CCD images of Eta Carinae. Here, we present the B, V, R, I, data of Eta Car obtained before and during the X-ray eclipse.