Spherical collapse predicts that a single value of the turnaround density (average matter density within the scale on which a structure detaches from the Hubble flow) characterizes all cosmic structures at the same redshift. It has been recently shown by Korkidis et al. that this feature persists in complex non-spherical galaxy clusters identified in N-body simulations. Here we show that the low-redshift evolution of the turnaround density constrains the cosmological parameters, and that it can be used to derive a local constraint on $\Omega_\Lambda$ alone, independent of $\Omega_m$. The turnaround density thus provides a promising new way to exploit upcoming large cosmological datasets.

The leading order effect of non-sphericity on the maximum size $R_{\rm TA,max}$ of realistic large scale bound cosmic structures is studied in the framework of $\Lambda{\rm CDM}$. Three simple models are used, two to describe intrinsically axisymmetric static mass distributions and the Kerr-de Sitter metric for stationary ones with non-vanishing angular momentum. In all cases the fractional change $\delta R_{\rm TA,max}(\theta)/R^{(0)}_{\rm TA,max}$ of $R_{\rm TA,max}$ of a given structure, compared to a spherical one with the same mass $M$, depends on the polar angle $\theta$ and is proportional to the product of the relevant eccentricity parameter, times the square of a small quantity, which in the static examples is the ratio of two characteristic length scales, while in the spinning case it is the ratio $v_{\rm out}/c$ of the azimuthal speed of the outmost members of the structure over the speed of light. Furthermore, the anglular average $\langle \delta R_{\rm TA,max}(\theta)/R^{(0)}_{\rm TA,max}\rangle$ is zero in the two static cases, while it is negative and proportional to ${\cal O}(v^2_{\rm out}/c^2)$ in the Kerr-de Sitter example. Thus, $\delta R_{\rm TA,max}(\theta)/R^{(0)}_{\rm TA,max}$ of axisymmetric structures is very small for practically any value of the eccentricity parameter. Based on the above, we argue that the maximum turnaround size of {\it realistic} superclusters has negligible dependence on their shapes.

We use N-body simulations to examine whether a characteristic turnaround radius, as predicted from the spherical collapse model in a $\rm {\Lambda CDM}$ Universe, can be meaningfully identified for galaxy clusters, in the presence of full three-dimensional effects. We use The Dark Sky Simulations and Illustris-TNG dark-matter--only cosmological runs to calculate radial velocity profiles around collapsed structures, extending out to many times the virial radius $R_{200}$. There, the turnaround radius can be unambiguously identified as the largest non-expanding scale around a center of gravity. We find that: (a) Indeed, a single turnaround scale can meaningfully describe strongly non-spherical structures. (b) For halos of masses $M_{200}>10^{13}M_\odot$, the turnaround radius $R_{ta}$ scales with the enclosed mass $M_{ta}$ as $M_{ta}^{1/3}$, as predicted by the spherical collapse model. (c) The deviation of $R_{ta}$ in simulated halos from the spherical collapse model prediction is insensitive to halo asphericity. Rather, it is sensitive to the tidal forces due to massive neighbors when such are present. (d) Halos exhibit a characteristic average density within the turnaround scale. This characteristic density is dependent on cosmology and redshift. For the present cosmic epoch and for concordance cosmological parameters ($\Omega_m \sim 0.7$; $\Omega_\Lambda \sim 0.3$) turnaround structures exhibit an average matter density contrast with the background Universe of $\delta \sim 11$. Thus $R_{ta}$ is equivalent to $R_{11}$ -- in a way analogous to defining the "virial" radius as $R_{200}$ -- with the advantage that $R_{11}$ is shown in this work to correspond to a kinematically relevant scale in N-body simulations.

The scalar and vector cosmological perturbations at all length scales of our Universe are studied in the framework of the phantom braneworld model. The model is characterized by the parameter $\Omega_M\equiv M^3/2m^2H_0$, with $M$ and $m$ the 5- and 4-dimensional Planck scales, respectively, and $H_0$ the Hubble parameter today, while $\Omega_M\rightarrow 0$ recovers the $\Lambda\rm CDM$ model. Ignoring the backreaction due to the peculiar velocities and also the bulk cosmological constant, allows the explicit computation of the gravitational potentials, $\Phi$ and $\Psi$. They exhibit exponentially decreasing screening behaviour characterized by a screening length which is a function of the quasidensity parameter $ \Omega_M$.

The latest observations of extensive air showers (EAS) induced by ultra-high-energy cosmic rays (UHECR) appear to indicate, prima facie, a transition to heavy primaries at the highest energies. However, this interpretation, which is based on extrapolations of the Standard Model (SM) to ultra-LHC energies, is strained from both astrophysical and particle phenomenology perspectives. We consider the alternative that after some energy threshold, the first collision of the primary in the atmosphere results in a state, the decay of which leads to a considerably increased shower particle multiplicity, so that light-primary EAS appear heavy-like. We show that a minimal implementation of such a model yields predictions for the average EAS depth and shower-to-shower fluctuations that are consistent with each other, and an excellent fit to Auger data. If such an effect indeed takes place, we predict that: (a) the center-of-momentum (CM) energy threshold for the effect is of order 50 TeV; (b) the probability with which the effect occurs is high, and it will be detected easily by next-generation accelerators; (c) the increase in multiplicity compared to the SM prediction grows with CM energy roughly as $\sim E_{\rm CM}$; (d) the cosmic-ray composition at the highest energies is light. Remarkably, if the latter is confirmed electromagnetically this would necessitate the existence of new physics by these energies.

The maximum allowable size of a spherical cosmic structure as a function of its mass is determined by the maximum turn around radius $R_{\rm TA,max}$, the distance from its centre where the attraction on a radial test particle due to the spherical mass is balanced with the repulsion due to the ambient dark energy. In this work, we extend the existing results in several directions. (a) We first show that for $w\neq -1$, the expression for $R_{\rm TA, max}$ found earlier using the cosmological perturbation theory, can be derived using a static geometry as well. (b) In the generic dark energy model with arbitrary time dependent state parameter $w(t)$, taking into account the effect of inhomogeneities upon the dark energy as well, where it is shown that the data constrain $w(t={\rm today})>-2.3$, and (c) in the quintessence and the generalized Chaplygin gas models, both of which are shown to predict structure sizes consistent with observations.

The maximum size of a cosmic structure is given by the maximum turnaround radius -- the scale where the attraction due to its mass is balanced by the repulsion due to dark energy. We derive generic formulae for the estimation of the maximum turnaround radius in any theory of gravity obeying the Einstein equivalence principle, in two situations: on a spherically symmetric spacetime and on a perturbed Friedman-Robertson-Walker spacetime. We show that the two formulae agree. As an application of our formula, we calculate the maximum turnaround radius in the case of the Brans-Dicke theory of gravity. We find that for this theory, such maximum sizes always lie above the \LCDM value, by a factor $1 + \frac{1}{3\omega}$, where $\omega\gg 1$ is the Brans-Dicke parameter, implying consistency of the theory with current data.

The maximum size of a bound cosmic structure is computed perturbatively as a function of its mass in the framework of the cubic galileon, proposed recently to model the dark energy of our Universe. Comparison of our results with observations constrains the matter-galileon coupling of the model to $0.033\lesssim \alpha \lesssim 0.17$, thus improving previous bounds based solely on solar system physics.

We demonstrate that, in the context of the $\Lambda$CDM model, it is in principle possible to measure the value of the cosmological constant by tracing, across cosmic time, the evolution of the turnaround radius of cosmic structures. The novelty of the presented method is that it is local, in the sense that it uses the effect of the cosmological constant on the relatively short scales of cosmic structures and not on the dynamics of the Universe at its largest scales. In this way, it can provide an important consistency check for the standard cosmological model and can give signs of new physics, beyond $\Lambda$CDM.

In $\Lambda$CDM cosmology, structure formation is halted shortly after dark energy dominates the mass/energy budget of the Universe. A manifestation of this effect is that in such a cosmology the turnaround radius has an upper bound. Recently, a new, local, test for the existence of dark energy in the form of a cosmological constant was proposed based on this turnaround bound. Before designing an experiment that, through high-precision determination of masses and turnaround radii, will challenge $\Lambda$CDM cosmology, we have to answer two important questions: first, when turnaround-scale structures are predicted to be close enough to their maximum size, so that a possible violation of the bound may be observable. Second, which is the best mass scale to target for possible violations of the bound. Using the Press-Schechter formalism, we find that turnaround structures have in practice already stopped forming, and consequently, the turnaround radius of structures must be very close to the maximum value today. We also find that the mass scale of $\sim 10^{13} M_\odot$ characterizes turnaround structures that start to form in a statistically important number density today. This mass scale also separates turnaround structures with different cosmological evolution: smaller structures are no longer readjusting their mass distribution inside the turnaround scale, they asymptotically approach their ultimate abundance from higher values, and they are common enough to have, at some epoch, experienced major mergers with structures of comparable mass; larger structures exhibit the opposite behavior. We call this mass scale the transitional mass scale and we argue that it is the optimal for the purpose outlined above. As a corollary result, we explain the different accretion behavior of small and larger structures observed in already conducted numerical simulations.

Our intuitive understanding of cosmic structure formation works best in scales small enough so that bound, relaxed gravitating systems are no longer adjusting their radius; and large enough so that space and matter follow the average expansion of the Universe. Yet one of the most robust predictions of $\Lambda$CDM cosmology concerns the scale that separates these limits: the turnaround radius, which is the non-expanding shell furthest away from the center of a bound structure. The maximum possible value of the turnaround radius within the framework of the $\Lambda$CDM model is, for a given mass $M$, equal to $(3GM/\Lambda c^2)^{1/3}$, with $G$ Newton's constant and $c$ the speed of light, independently of cosmic epoch, exact nature of dark matter, or baryonic effects. We discuss the possible use of this prediction as an observational test for $\Lambda$CDM cosmology.

For a general dark-energy equation of state, we estimate the maximum possible radius of massive structures that are not destabilized by the acceleration of the cosmological expansion. A comparison with known stable structures constrains the equation of state. The robustness of the constraint can be enhanced through the accumulation of additional astrophysical data and a better understanding of the dynamics of bound cosmic structures.

An explicitly covariant formulation is presented of a modified DGP scenario proposed recently [1], to avoid the instability of the self-accelerating branch. It is based on the introduction of a bulk scalar field with appropriate non-minimal coupling to the bulk Einstein-Hilbert term. The method is general and may be applied to other models as well.

The role of brane-bulk energy exchange and of an induced gravity term on a single braneworld of negative tension and vanishing effective cosmological constant is studied. It is shown that for the physically interesting cases of dust and radiation a unique global attractor which can realize our present universe (accelerating and 0<Omega_{m0}<1) exists for a wide range of the parameters of the model. For Omega_{m0}=0.3, independently of the other parameters, the model predicts that the equation of state for the dark energy today is w_{DE,0}=-1.4, while Omega_{m0}=0.03 leads to w_{DE,0}=-1.03. In addition, during its evolution, w_{DE} crosses the w_{DE}=-1 line to smaller values.

We present a study of the decay of an intermediate mini black hole at the first impact of a cosmic ray particle with the atmosphere, in the context of D-brane world scenarios with TeV scale gravity and large extra dimensions. We model the decay of the black hole using the semiclassical approximation and include the corrections coming from energy loss into the bulk. Extensive simulations show that mini black hole events are characterized by essentially different multiplicities and wider lateral distributions of the air showers as a function of the energy of the incoming primary, as compared to standard events. Implications for their detection and some open issues on their possible discovery are also briefly addressed.

We present a study of the multiplicities, of the lateral distributions and of the ratio of the electromagnetic to the hadronic components in the air showers, generated by the collision in the atmosphere of an incoming high energy cosmic ray and mediated by the formation of a mini black hole, predicted in TeV scale gravity models with large extra dimensions. The analysis is performed via a large scale simulation of the resulting cascades over the entire range ($10^{15}-10^{19}$eV) of ultra high initial energies, for several values of the number of large extra dimensions, for a variety of altitudes of the initial interaction and with the energy losses in the bulk taken into account. The results are compared with a representative of the standard events, namely the shower due to the collision of a primary proton with a nucleon in the atmosphere. Both the multiplicities and the lateral distribution of the showers show important differences between the two cases and, consequently, may be useful for the observational characterization of the events. The electromagnetic/hadronic ratio is strongly fluctuating and, thus, less decisive for the altitudes considered.

We present a detailed analytical and numerical study of a novel type of static, superconducting, classically stable string texture in a renormalizable topologically trivial massive U(1) gauge model with one charged and one neutral scalar. An upper bound on the mass of the charged scalar as well as on the current that the string can carry are established. A preliminary unsuccesful search for stable solutions corresponding to large superconducting loops is also reported.

We establish the existence of a static, classically stable string solution in a region of parameters of the generic two-Higgs Standard Model. In an appropriate limit of parameters, the solution reduces to the well-known soliton of the O(3) non-linear sigma model.

We derive the Feynman rules for the graviton in the presence of a flat Robertson-Walker background and give explicit expressions for the propagator in the physically interesting cases of inflation, radiation domination, and matter domination. The aforementioned background is generated by a scalar field source which should be taken to be dynamical. As an elementary application, we compute the corrections to the Newtonian gravitational force in the present matter dominated era and conclude -- as expected -- that they are negligible except for the largest scales.