We study a new class of vector dark energy models where multi-Proca fields $A_\mu^a$ are coupled to cold dark matter by the term $f(X)\tilde{\mathcal{L}}_{m}$ where $f(X)$ is a general function of $X\equiv -\frac{1}{2}A^\mu_ a A^a_\mu$ and $\tilde{\mathcal{L}}_{m}$ is the cold dark matter Lagrangian. From here, we derive the general covariant form of the novel interaction term sourcing the field equations. This result is quite general in the sense that encompasses Abelian and non-Abelian vector fields. In particular, we investigate the effects of this type of coupling in a simple dark energy model based on three copies of canonical Maxwell fields to realize isotropic expansion. The cosmological background dynamics of the model is examined by means of a dynamical system analysis to determine the stability of the emergent cosmological solutions. As an interesting result, we find that the coupling function leads to the existence of a novel scaling solution during the dark matter dominance. Furthermore, the critical points show an early contribution of the vector field in the form of dark radiation and a stable de Sitter-type attractor at late times mimicking dark energy. The cosmological evolution of the system as well as the aforementioned features are verified by numerical computations. Observational constraints are also discussed to put the model in a more phenomenological context in the light of future observations.

The infrared source known as Orion n was detected in 1980 with observations made with the 3.8-m United Kingdom Infrared Telescope. About two decades later, sensitive observations made with the Very Large Array revealed the presence of a mJy double radio source apparently coincident in position with the infrared source n. The radio source was assumed to be the counterpart of the infrared source. However, over the years it has been concluded that the radio source shows large proper motions to the south while the infrared source n is stationary. Here we reanalyze the proper motions of the radio source adding both older and newer VLA observations than previously used. We confirm the proper motions of the radio source that at present no longer coincides positionally with the infrared source. The solution to this problem is, most probably, that the infrared source n and the radio source are not the same object: the infrared source is a stationary object in the region while the radio counterpart is moving as a result of the explosion that took place in this region some 500 years ago and that expelled large amounts of molecular gas as well as several compact sources. Considering the paper where it was first reported, we refer to this double radio source as Orion MR. In addition, we use these new observations to fully confirm the large proper motions of the sources IRc23 and Zapata 11. Together with sources BN, I, Orion MR, and x, there are at least six compact sources that recede from a point in common in Orion BN/KL. However, IRc23 is peculiar in that its ejection age appears to be only $\sim$300 years. The relatively large number of sources rules out as a possible mechanism the classic three-body scenario since then only two escaping bodies are expected: a tight binary plus the third star involved in the encounter.

Under the same spirit of the Galileon-Horndeski theories and their more modern extensions, the generalized SU(2) Proca theory was built by demanding that its action may be free of the Ostrogradski's instability. Nevertheless, the theory must also be free of other instability problems in order to ensure its viability. As a first approach to address this issue, we concentrate on a quite general variant of the theory and investigate the general conditions for the absence of ghost and gradient instabilities in the tensor sector without the need for resolving the dynamical background. The phenomenological interest of this approach as well as of the variant investigated lies on the possibility of building cosmological models driven solely by non-Abelian vector fields that may account for a successful description of both the early inflation and the late-time accelerated expansion of the universe.

Young massive clusters (YMCs) are the most compact, high-mass stellar systems still forming at the present day. The precursor clouds to such systems are, however, rare due to their large initial gas mass reservoirs and rapid dispersal timescales due to stellar feedback. Nonetheless, unlike their high-z counterparts, these precursors are resolvable down to the sites of individually forming stars, and hence represent the ideal environments in which to test the current theories of star and cluster formation. Using high angular resolution (1$^{\prime\prime}$ / 0.05pc) and sensitivity ALMA observations of two YMC progenitor clouds in the Galactic Centre, we have identified a suite of molecular line transitions -- e.g. c-C$_{3}$H$_{2} $($7-6$) -- that are believed to be optically thin, and reliably trace the gas structure in the highest density gas on star-forming core scales. We conduct a virial analysis of the identified core and proto-cluster regions, and show that half of the cores (5/10) and both proto-clusters are unstable to gravitational collapse. This is the first kinematic evidence of global gravitational collapse in YMC precursor clouds at such an early evolutionary stage. The implications are that if these clouds are to form YMCs, then they likely do so via the "conveyor-belt" mode, whereby stars continually form within dispersed dense gas cores as the cloud undergoes global gravitational collapse. The concurrent contraction of both the cluster-scale gas and embedded (proto)stars ultimately leads to the high (proto)stellar density in YMCs.

Two families of models compete to explain the formation of high-mass stars. The quasi-static models predict the existence of high-mass pre-stellar cores sustained by a high degree of turbulence while competitive accretion models predict that high-mass proto-stellar cores evolve from low/intermediate mass proto-stellar cores in dynamic environments. We present ALMA (1.4 mm continuum emission and $^{12}$CO emission line) and MOPRA (HCO$^{+}$, H$^{13}$CO$^{+}$ and N$_2$H$^+$ molecular line emissions) observations of a sample of 9 starless massive dense cores (MDCs) discovered in a recent Herschel/HOBYS study that have masses and sizes ($\sim$110 M$_\odot$ and $r\sim$0.1 pc, respectively) similar to the initial conditions used in the quasi-static models. The MOPRA molecular line features show that 3 of the starless MDCs are subvirialized with $\alpha_{\rm vir}\sim$0.35, and 4 MDCs show sign of collapse. Our ALMA observations, on the other hand, show very little fragmentation within the MDCs whereas the observations resolve the Jeans length ($\lambda_{\rm Jeans}\sim$0.03 pc) and are sensitive to the Jeans mass (M$_{\rm Jeans}\sim$0.65 M$_\odot$) in the 9 starless MDCs. Only two of the starless MDCs host compact continuum sources, whose fluxes correspond to $<3$ M$_\odot$ fragments. Therefore the mass reservoir of the MDCs has not yet been accreted onto compact objects, and most of the emission is filtered out by the interferometer. These observations do not support the quasi-static models for high-mass star formation since no high-mass pre-stellar core is found in NGC6334. The competitive accretion models, on the other hand, predict a level of fragmentation much higher than what we observe.

Using Galactic Plane surveys, we have selected a massive (1200 M$_\odot$), cold (14 K) 3.6-70 $\mu$m dark IRDC G331.372-00.116. This IRDC has the potential to form high-mass stars and, given the absence of current star formation signatures, it seems to represent the earliest stages of high-mass star formation. We have mapped the whole IRDC with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.1 and 1.3 mm in dust continuum and line emission. The dust continuum reveals 22 cores distributed across the IRDC. In this work, we analyze the physical properties of the most massive core, ALMA1, which has no molecular outflows detected in the CO (2-1), SiO (5-4), and H$_2$CO (3-2) lines. This core is relatively massive ($M$ = 17.6 M$_\odot$), subvirialized (virial parameter $\alpha_{vir}=M_{vir}/M=0.14$), and is barely affected by turbulence (transonic Mach number of 1.2). Using the HCO$^+$ (3-2) line, we find the first detection of infall signatures in a relatively massive, prestellar core (ALMA1) with the potential to form a high-mass star. We estimate an infall speed of 1.54 km s$^{-1}$ and a high accretion rate of 1.96 $\times$ 10$^{-3}$ M$_\odot$ yr$^{-1}$. ALMA1 is rapidly collapsing, out of virial equilibrium, more consistent with competitive accretion scenarios rather than the turbulent core accretion model. On the other hand, ALMA1 has a mass $\sim$6 times larger than the clumps Jeans mass, being in an intermediate mass regime ($M_{J}=2.7<M\lesssim$ 30 M$_\odot$), contrary to what both the competitive accretion and turbulent core accretion theories predict.

We have extended the time baseline for observations of the proper motions of radio sources in the Orion BN/KL region from 14.7 to 22.5 years. We present improved determinations for the sources BN and I. In addition, we address the proper motions of the double radio source n, that have been questioned in the literature. We confirm that all three sources are moving away at transverse velocities of tens of km s$^{-1}$ from a region in-between them, where they were located about 500 years ago. Source n exhibits a new component that we interpret as due to a one-sided ejection of free-free emitting plasma that took place after 2006.36. We used the highly accurate relative proper motions between sources BN and I to determine that their closest separation took place in the year 1475$\pm$6, when they were within $\sim$100 AU or less from each other in the plane of the sky.

Observations using the 7 mm receiver system on the Australia Telescope Compact Array have revealed large reservoirs of molecular gas in two high-redshift radio galaxies: HATLAS J090426.9+015448 (z = 2.37) and HATLAS J140930.4+003803 (z = 2.04). Optically the targets are very faint, and spectroscopy classifies them as narrow-line radio galaxies. In addition to harbouring an active galactic nucleus the targets share many characteristics of sub-mm galaxies. Far-infrared data from Herschel-ATLAS suggest high levels of dust (>10^9 M_solar) and a correspondingly large amount of obscured star formation (~1000 M_solar / yr). The molecular gas is traced via the J = 1-0 transition of 12CO, its luminosity implying total H_2 masses of (1.7 +/- 0.3) x 10^11 and (9.5 +/- 2.4) x 10^10 (alpha_CO/0.8) M_solar in HATLAS J090426.9+015448 and HATLAS J140930.4+003803 respectively. Both galaxies exhibit molecular line emission over a broad (~1000 km/s) velocity range, and feature double-peaked profiles. We interpret this as evidence of either a large rotating disk or an on-going merger. Gas depletion timescales are ~100 Myr. The 1.4 GHz radio luminosities of our targets place them close to the break in the luminosity function. As such they represent `typical' z > 2 radio sources, responsible for the bulk of the energy emitted at radio wavelengths from accretion-powered sources at high redshift, and yet they rank amongst the most massive systems in terms of molecular gas and dust content. We also detect 115 GHz rest-frame continuum emission, indicating a very steep high-radio-frequency spectrum, possibly classifying the targets as compact steep spectrum objects.

We will expose in this paper our advances towards a proof of the equivalence between FRW background expansion, during some period of time that contains primordial inflation, and the statistical isotropy of the primordial curvature perturbation $\zeta$ at the end of this period of time. Our motivation rests on the growing interest in the existence of a preferred direction in the Universe hinted by the continuous presence of anomalies in the CMB data.

We combine previously published interferometric and single-dish data of relatively nearby massive dense cores that are actively forming stars to test whether their `fragmentation level' is controlled by turbulent or thermal support. We find no clear correlation between the fragmentation level and velocity dispersion, nor between the observed number of fragments and the number of fragments expected when the gravitationally unstable mass is calculated including various prescriptions for `turbulent support'. On the other hand, the best correlation is found for the case of pure thermal Jeans fragmentation, for which we infer a core formation efficiency around 13 per cent, consistent with previous works. We conclude that the dominant factor determining the fragmentation level of star-forming massive dense cores at 0.1 pc scale seems to be thermal Jeans fragmentation.

We present an analysis of the dust continuum emission at 870 um in order to investigate the mass distribution of clumps within infrared dark clouds (IRDCs). We map six IRDCs with the Large APEX BOlometer CAmera (LABOCA) at APEX, reaching an rms noise level of 28-44 mJy/beam. The dust continuum emission coming from these IRDCs was decomposed by using two automated algorithms, Gaussclumps and Clumpfind. Moreover, we carried out single-pointing observations of the N_2H^+ (3-2) line toward selected positions to obtain kinematic information. The mapped IRDCs are located in the range of kinematic distances of 2.7-3.2 kpc. We identify 510 and 352 sources with Gaussclumps and Clumpfind, respectively, and estimate masses and other physical properties assuming a uniform dust temperature. The mass ranges are 6-2692 Msun (Gaussclumps) and 7-4254 Msun (Clumpfind) and the ranges in effective radius are around 0.10-0.74 pc (Gaussclumps) and 0.16-0.99 pc (Clumpfind). The mass distribution, independent of the decomposition method used, is fitted by a power law, dN/dM propto M^alpha, with an index of -1.60 +/- 0.06, consistent with the CO mass distribution and other high-mass star-forming regions.

We study the most general contributions due to scalar field perturbations, vector field perturbations, and anisotropic expansion to the generation of statistical anisotropy in the primordial curvature perturbation \zeta. Such a study is done using the \delta N formalism where only linear terms are considered. Here, we consider two specific cases that lead to determine the power spectrum P_\zeta(k) of the primordial curvature perturbation. In the first one, we consider the possibility that the n-point correlators of the field perturbations in real space are invariant under rotations in space (statistical isotropy); as a result, we obtain as many levels of statistical anisotropy as vector fields present and, therefore, several preferred directions. The second possibility arises when we consider anisotropic expansion, which leads us to obtain I+a additional contributions to the generation of statistical anisotropy of \zeta compared with the former case, being I and a the number of scalar and vector fields involved respectively.

We present an astrometry study of the radio source VLA 1 at the core of the HH 111 outflow using new data (2007) as well as archival observations (1992-1996). All data were taken at 3.6 cm with the Very Large Array in its most extended (A) configuration. The source VLA 1 has undergone a dramatic morphological change, showing a one-sided knot ejection in the 2007 epoch. We also report on the detection of a 3.6 cm compact continuum source (VLA 3) located at (-10.6",98.7") from VLA 1. No significant absolute proper motions were found for VLA 1 and VLA 3 and the upper limits are consistent with those found for (embedded) radio sources in the Orion Nebula. We favor the interpretation that in the continuum at 3.6 cm we are observing two nearly perpendicular jets. HH 111 presents a new case of one-sided jet ejection in a young stellar object. The Galactic (or extragalactic) nature of VLA 3 remains unclear.

Recent studies suggest that infrared dark clouds (IRDCs) have the potential of harboring the earliest stages of massive star formation and indeed evidence for this is found toward distinct regions within them. We present a study with the Plateau de Bure Interferometer of a core in the archetypal filamentary IRDC G11.11-0.12 at few arcsecond resolution to determine its physical and chemical structure. The data consist of continuum and line observations covering the C34S 2-1 line and the methanol 2_k-1_k v_t=0 lines at 3mm and the methanol 5_k-4_k v_t =0 lines at 1mm. Our observations show extended emission in the continuum at 1 and 3 mm. The methanol 2_k-1_k v_t=0 emission presents three maxima extending over 1 pc scale (when merged with single-dish short-spacing observations); one of the maxima is spatially coincident with the continuum emission. The fitting results show enhanced methanol fractional abundance (~3x10^-8) at the central peak with respect to the other two peaks, where it decreases by about an order of magnitude (~4-6x10^-9). Evidence of extended 4.5 microns emission, "wings" in the CH3OH 2_k-1_k spectra, and CH3OH abundance enhancement point to the presence of an outflow in the East-West direction. In addition, we find a gradient of ~4 km/s in the same direction, which we interpret as being produced by an outflow(s)-cloud interaction.

We present Very Large Array observations of the 7_0-6_1 A^+ methanol maser transition at 44 GHz towards NGC 6334F, G8.67-0.36, and M17. These arcsecond resolution observations complete a previous, larger VLA survey of this maser transition in high-mass star-forming regions reported by Kurtz et al. We confirm the presence of 44 GHz methanol maser emission in all three sources, detecting eight distinct maser components in NGC 6334F, twelve components in G8.67-0.36 and one in M17.

We present two Very Large Array observations of the pre-main-sequence star PV Cephei, taken with a separation of 10.5 years. These data show that the proper motions of this star are $\mu_\alpha \cos \delta = +10.9 \pm 3.0$ mas yr$^{-1}$; $\mu_\delta = +0.2 \pm 1.8$ mas yr$^{-1}$, very similar to those -previously known- of HD 200775, the B2Ve star that dominates the illumination of the nearby reflection nebula NGC 7023. This result suggests that PV Cephei is not a rapidly moving run-away star as suggested by previous studies. The large velocity of PV Cephei had been inferred from the systematic eastward displacement of the bisectors of successive pairs of Herbig Haro knots along its flow. These systematic shifts might instead result from an intrinsic dissymmetry in the ejection mechanisms, or from an asymmetric distribution of the circumstellar material.

Type 2 quasars are luminous Active Galactic Nuclei whose central engines are seen through large amounts of gas and dust. We present Spitzer spectra of twelve type 2 quasars selected on the basis of their optical emission line properties. Within this sample, we find a surprising diversity of spectra, from those that are featureless to those showing strong PAH emission, deep silicate absorption at 10 micron, hydrocarbon absorption, high-ionization emission lines and H_2 rotational emission lines. About half of the objects in the sample are likely Compton-thick, including the two with the deepest Si absorption. The median star-formation luminosity of the objects in our sample measured from the strength of the PAH features is 5x10^11 L_sun, much higher than for field galaxies or for any other AGN sample, but similar to other samples of type 2 quasars. This suggests an evolutionary link between obscured quasars and peak star formation activity in the host galaxy. Despite the high level of star formation, the bolometric output is dominated by the quasar in all cases. For a given strength of 10 micron Si absorption, ULIRGs are significantly colder than are type 2 quasars (their F_nu[14.5 micron]/F_nu[27.5 micron] ratio is 0.5 dex lower), perhaps reflecting different obscuration geometries in these sources. We find that the appearance of the 10 micron feature (i.e., whether it shows in emission or in absorption) is well-correlated with the optical classification in type 1 and type 2 quasars, contrary to some models of clumpy obscuration. Furthermore, this correlation is significantly stronger in quasars (L_bol>10^45 erg/s) than it is in Seyfert galaxies (L_bol<<10^45 erg/s).

We present absolute astrometry of four radio sources in the Becklin-Neugebauer/Kleinman-Low (BN/KL) region, derived from archival data (taken in 1991, 1995, and 2000) as well as from new observations (taken in 2006). All data consist of 3.6 cm continuum emission and were taken with the Very Large Array in its highest angular resolution A configuration. We confirm the large proper motions of the BN object, the radio source I (GMR I) and the radio counterpart of the infrared source n (Orion-n), with values from 15 to 26 km/s. The three sources are receding from a point between them from where they seem to have been ejected about 500 years ago, probably via the disintegration of a multiple stellar system. We present simulations of very compact stellar groups that provide a plausible dynamical scenario for the observations. The radio source Orion-n appeared as a double in the first three epochs, but as single in 2006. We discuss this morphological change. The fourth source in the region, GMR D, shows no statistically significant proper motions. We also present new, accurate relative astrometry between BN and radio source I that restrict possible dynamical scenarios for the region. During the 2006 observations, the radio source GMR A, located about 1' to the NW of the BN/KL region, exhibited an increase in its flux density of a factor of ~3.5 over a timescale of one hour. This rapid variability at cm wavelengths is similar to that previously found during a flare at millimeter wavelengths that took place in 2003.

We present absolute astrometry of 35 radio sources in the Orion Trapezium and Becklin-Neugebauer/Kleinman-Low regions, obtained from Very Large Array archival observations collected over a period of 15 years. By averaging the results for all the sources, we estimate the mean absolute proper motion of Orion to be --in Galactic coordinates-- $\mu_\ell \cos b$ = +2.1 $\pm$ 0.2 mas yr$^{-1}$; $\mu_b$ = $-$0.1 $\pm$ 0.2 mas yr$^{-1}$. These values agree remarkably well with those expected from the differential rotation of the Milky Way. Subtraction of this mean motion from the individual measurements allows us to register all proper motions to the rest frame of the Orion nebula, and identify radio sources with large residual velocities. In the KL region, we find three sources in this situation: the BN object, the radio source I, and the radio counterpart of the infrared source n. All three objects appear to be moving away from a common point where they must all have been located about 500 years ago. This suggests that all three sources were originally part of a multiple massive stellar system that recently disintegrated as a result of a close dynamical