We investigate the nonlinear growth of large-scale structure in the generalised Proca theory, in which a self-interacting massive vector field plays the role of driving the acceleration of the cosmic expansion. Focusing to the Proca Lagrangian at cubic order -- the cubic vector Galileon model -- we derive the simplified equations for gravity as well as the longitudinal and transverse modes of the vector field under the weak-field and quasi-static approximations, and implement them in a modified version of the ECOSMOG N-body code. Our simulations incorporate the Vainshtein screening effect, which reconciles the fifth force propagated by the longitudinal mode of the cubic vector Galileon model with local tests of gravity. The results confirm that for all scales probed by the simulation, the transverse mode has a negligible impact on structure formation in a realistic cosmological setup. It is well known that in this model the strength of the fifth force is controlled by a free model parameter, which we denote as beta_3. By running a suite of cosmological simulations for different values of beta_3, we show that this parameter also determines the effectiveness of the Vainshtein screening. The model behaves identically to the cubic scalar Galileon for beta_3 going to zero, in which the fifth force is strong in unscreened regions but is efficiently screened in high-density regions. In the opposite limit, beta_3 going to infinity, the model approaches its `quintessence' counterpart, which has a vanishing fifth force but a modified expansion history compared to LambdaCDM. This endows the model with rich phenomenology, which will be investigated in future works.

Measurements of The Hubble-Lemaitre constant from early- and local-universe observations show a significant discrepancy. In an attempt to understand the origin of this mismatch, independent techniques to measure $H_0$ are required. One such technique, strong lensing time delays, is set to become a leading contender amongst the myriad methods due to forthcoming large strong lens samples. It is therefore critical to understand the systematic effects inherent in this method. In this paper, we quantify the influence of additional structures along the line-of-sight by adopting realistic lightcones derived from the \textit{CosmoDC2} semi-analytical extra-galactic catalogue. Using multiple lens plane ray-tracing to create a set of simulated strong lensing systems, we have investigated the impact of line-of-sight structures on time-delay measurements and in turn, on the inferred value of $H_0$. We have also tested the reliability of existing procedures for correcting for line-of-sight effects. We find that if the integrated contribution of the of line-of-sight structures is close to a uniform mass sheet, the bias in $H_0$ can be adequately corrected by including a constant external convergence $\kappa_{\rm ext}$ in the lens model. However, for realistic line-of-sight structures comprising many galaxies at different redshifts, this simple correction over-estimates the bias by a factor of approximately three. We therefore conclude that lens modelling must incorporate multiple lens planes to account for line-of-sight structures for accurate and precise inference of $H_0$.

LOFT-P is a concept for a NASA Astrophysics Probe-Class (<$1B) X-ray timing mission, based on the LOFT concept originally proposed to ESAs M3 and M4 calls. LOFT-P requires very large collecting area (>6 m^2, >10x RXTE), high time resolution, good spectral resolution, broad-band spectral coverage (2-30 keV), highly flexible scheduling, and an ability to detect and respond promptly to time-critical targets of opportunity. It addresses science questions such as: What is the equation of state of ultra dense matter? What are the effects of strong gravity on matter spiraling into black holes? It would be optimized for sub-millisecond timing to study phenomena at the natural timescales of neutron star surfaces and black hole event horizons and to measure mass and spin of black holes. These measurements are synergistic to imaging and high-resolution spectroscopy instruments, addressing much smaller distance scales than are possible without very long baseline X-ray interferometry, and using complementary techniques to address the geometry and dynamics of emission regions. A sky monitor (2-50 keV) acts as a trigger for pointed observations, providing high duty cycle, high time resolution monitoring of the X-ray sky with ~20 times the sensitivity of the RXTE All-Sky Monitor, enabling multi-wavelength and multi-messenger studies. A probe-class mission concept would employ lightweight collimator technology and large-area solid-state detectors, technologies which have been recently greatly advanced during the ESA M3 study. Given the large community interested in LOFT (>800 supporters, the scientific productivity of this mission is expected to be very high, similar to or greater than RXTE (~2000 refereed publications). We describe the results of a study, recently completed by the MSFC Advanced Concepts Office, that demonstrates that LOFT-P is feasible within a NASA probe-class mission budget.

This paper discusses a new method to generate self-coherent initial conditions for young substructured stellar cluster. The expansion of a uniform system allows stellar sub-structures (clumps) to grow from fragmentation modes by adiabatic cooling. We treat the system mass elements as stars, chosen according to a Salpeter mass function, and the time-evolution is performed with a collisional N-body integrator. This procedure allows to create a fully-coherent relation between the clumps' spatial distribution and the underlying velocity field. The cooling is driven by the gravitational field, as in a cosmological Hubble-Lema\^itre flow. The fragmented configuration has a `fractal'-like geometry but with a self-grown velocity field and mass profile. We compare the characteristics of the stellar population in clumps with that obtained from hydrodynamical simulations and find a remarkable correspondence between the two in terms of the stellar content and the degree of spatial mass-segregation. In the fragmented configuration, the IMF power index is ~0.3 lower in clumps in comparison to the field stellar population, in agreement with observations in the Milky Way. We follow in time the dynamical evolution of fully fragmented and sub-virial configurations, and find a soft collapse, leading rapidly to equilibrium (timescale of 1 Myr for a ~ 10^4 Msun system). The low-concentration equilibrium implies that the dynamical evolution including massive stars is less likely to induce direct collisions and the formation of exotic objects. Low-mass stars already ejected from merging clumps are depleted in the end-result stellar clusters, which harbour a top-heavy stellar mass function.

The diversity of structures in the Universe (from the smallest galaxies to the largest superclusters) has formed under the pull of gravity from the tiny primordial perturbations that we see imprinted in the cosmic microwave background. A quantitative description of this process would require description of motion of zillions of dark matter particles. This impossible task is usually circumvented by coarse-graining the problem: one either considers a Newtonian dynamics of "particles" with macroscopically large masses or approximates the dark matter distribution with a continuous density field. There is no closed system of equations for the evolution of the matter density field alone and instead it should still be discretized at each timestep. In this work we describe a method of solving the full 6-dimensional Vlasov-Poisson equation via a system of auxiliary Schroedinger-like equations. The complexity of the problem gets shifted into the choice of the number and shape of the initial wavefunctions that should only be specified at the beginning of the computation (we stress that these wavefunctions have nothing to do with quantum nature of the actual dark matter particles). We discuss different prescriptions to generate the initial wave functions from the initial conditions and demonstrate the validity of the technique on two simple test cases. This new simulation algorithm can in principle be used on an arbitrary distribution function, enabling the simulation of warm and hot dark matter structure formation scenarios.

In the scope of the star formation process, it is unclear how the environment shapes the initial mass function (IMF). While observations of open clusters propose a universal picture for the IMF from the substellar domain up to a few solar masses, the young association eta Chamaeleontis presents an apparent lack of low mass objects (m<0.1 Msun). Another unusual feature of this cluster is the absence of wide binaries with a separation > 50 AU. We aim to test whether dynamical evolution alone can reproduce the peculiar properties of the association assuming a universal IMF. We use a pure N-body code to simulate the dynamical evolution of the cluster for 10 Myr, and compare the results with observations. A wide range of values for the initial parameters are tested in order to identify the initial state that would most likely lead to observations. In this context we also investigate the influence of the initial binary population on the dynamics and the possibility of having a discontinuous single IMF near the transition to the brown dwarf regime. We consider as an extreme case an IMF with no low mass systems (m<0.1 Msun). The initial configurations cover a wide range of initial density, from 10^2 to 10^8 stars/pc^3, in virialized, hot and cold dynamical state. We do not find any initial state that would evolve from a universal single IMF to fit the observations. Only when starting with a truncated IMF without any very low mass systems and no wide binaries, can we reproduce the cluster core properties with a success rate of 10% at best. Pure dynamical evolution alone cannot explain the observed properties of eta Cha from universal initial conditions. The lack of brown dwarfs and very low mass stars, and the peculiar binary properties (low binary fraction and lack of wide binaries), are probably the result of the star formation process in this association. (abridged)

From 1994 to 1997 we monitored the optical spectrum of RXJ0019.8+2156. This supersoft X-ray source is one of only two accreting white dwarfs in the Galaxy thought to be burning hydrogen on their surface as a consequence of a high rate of mass transfer from a binary companion. Accurate orbital ephemerides are derived from radial velocity measure- ments of HeII emission lines, which exhibit a stable velocity semi- amplitude of 71.2 +- 3.6 km/s. We report the discovery of transient, low-velocity, bipolar jets. These jets are respresented by redshifted and blueshifted pairs of emission lines with an outflow velocity of v cos(i) = 815 km/s, where i is the binary inclination angle. When present, the jet lines also exhibit an orbital modulation of 71 km/s, which indicates that the jets are oriented nearly perpendicular to the orbital plane. On most occasions, the Hline profiles are further altered by P Cygni absorption. We show that the jets and the P Cygni features have very different temporal characteristics and binary phase dependence. The jet material and the absorbing wind must therefore occupy very different geometries. Finally, the measured mass function is combined with binary evolution models to suggest a limit, i < 40 deg. A particular model invoked to explain a high rate of mass transfer requires 16 < i < 25. However, at such small inclination it is difficult to explain the large amplitude of the orbital light curve. By contrast, a simple model fit to the jet outflow lines indicates an orbital inclination angle of 35 < i < 60.

Optical photometry is presented for the quadruple gravitational lens PG1115+080. A preliminary reduction of data taken from November 1995 to June 1996 gives component ``C'' leading component ``B'' by 23.7+/-3.4 days and components ``A1'' and ``A2'' by 9.4 days. A range of models has been fit to the image positions, none of which gives an adequate fit. The best fitting and most physically plausible of these, taking the lensing galaxy and the associated group of galaxies to be singular isothermal spheres, gives a Hubble constant of 42 km/s/Mpc for Omega=1, with an observational uncertainty of 14%, as computed from the B-C time delay measurement. Taking the lensing galaxy to have an approximately E5 isothermal mass distribution yields H0=64 km/sec/Mpc while taking the galaxy to be a point mass gives H0=84 km/sec/Mpc. The former gives a particularly bad fit to the position of the lensing galaxy, while the latter is inconsistent with measurements of nearby galaxy rotation curves. Constraints on these and other possible models are expected to improve with planned HST observations.