The Dark Energy Spectroscopic Instrument (DESI) Collaboration has obtained robust measurements of baryon acoustic oscillations (BAO) in the redshift range, 0.1 < z < 4.2, based on the Lyman-$\alpha$ forest and galaxies from Data Release 2 (DR2). We combine these measurements with external cosmic microwave background (CMB) data from Planck and ACT to place our tightest constraints yet on the sum of neutrino masses. Assuming the cosmological $\Lambda$CDM model and three degenerate neutrino states, we find \sum m_\nu<0.0642 eV (95%) with a marginalized error of $\sigma(\sum m_\nu)=0.020$ eV. We also constrain the effective number of neutrino species, finding N_\rm{eff} = 3.23^{+0.35}_{-0.34} (95%), in line with the Standard Model prediction. When accounting for neutrino oscillation constraints, we find a preference for the normal mass ordering and an upper limit on the lightest neutrino mass of m_l < 0.023 eV (95%). However, we determine using frequentist and Bayesian methods that our constraints are in tension with the lower limits derived from neutrino oscillations. Correcting for the physical boundary at zero mass, we report a 95% Feldman-Cousins upper limit of \sum m_\nu<0.053 eV, breaching the lower limit from neutrino oscillations. Considering a more general Bayesian analysis with an effective cosmological neutrino mass parameter, $\sum m_{\nu,\rm{eff}}$, that allows for negative energy densities and removes unsatisfactory prior weight effects, we derive constraints that are in $3\sigma$ tension with the same oscillation limit. In the absence of unknown systematics, this finding could be interpreted as a hint of new physics not necessarily related to neutrinos. The preference of DESI and CMB data for an evolving dark energy model offers one possible solution. In the $w_0w_a$CDM model, we find \sum m_\nu<0.163 eV (95%), relaxing the neutrino tension. [Abridged]

Recent advances in cosmological observations have provided an unprecedented opportunity to investigate the distribution of baryons relative to the underlying matter. In this work, we show that the gas is more extended than the dark matter, and the amount of baryonic feedback at $z \lesssim 1$ disfavors low-feedback models such as that of state-of-the-art hydrodynamical simulation IllustrisTNG compared with high-feedback models such as that of the original Illustris simulation. This has important implications for bridging the gap between theory and observations and understanding galaxy formation and evolution. Furthermore, a better grasp of the baryon-dark matter link is critical to future cosmological analyses, which are currently impeded by our limited knowledge of baryonic feedback. Here, we measure the kinematic Sunyaev-Zel'dovich (kSZ) effect from the Atacama Cosmology Telescope (ACT), stacked on the luminous red galaxy (LRG) sample of the Dark Energy Spectroscopic Instrument (DESI) imaging survey. This is the first analysis to use photometric redshifts for reconstructing galaxy velocities. Due to the large number of galaxies comprising the DESI imaging survey, this is the highest signal-to-noise stacked kSZ measurement to date: we detect the signal at 13$\sigma$, finding strong evidence that the gas is more spread out than the dark matter, as well as a preference for larger feedback compared to some commonly used state-of-the-art hydrodynamical simulations. Our work opens up the possibility of recalibrating large hydrodynamical simulations using the kSZ effect. In addition, our findings point towards a way of alleviating inconsistencies between weak lensing surveys and cosmic microwave background (CMB) experiments, such as the `low $S_8$' tension, and shed light on long-standing enigmas in astrophysics, such as the `missing baryon' problem.

CARMENES is a dual-channel high-resolution spectrograph at the 3.5 m Calar Alto telescope designed to detect low-mass planets around late-type dwarfs by measuring their radial velocities (RVs). High thermal stability in both the visible (VIS) and near infrared channels is essential to achieve the precision required for these measurements. In particular, stabilising the NIR channel to the millikelvin level, which operates at cryogenic temperatures (140 K), poses significant engineering this http URL CARMENES-PLUS project was initiated to improve the instruments intrinsic RV precision. In this article, we focus on the thermal stability improvements made to the NIR channels cooling system. The NIR cooling system was originally conceived to operate with a discontinuous flow of cryogenic nitrogen gas. As part of CARMENES-PLUS, this was upgraded to a continuous flow configuration. Additional changes included the installation of an automatic vacuum system, a proportional control valve, and a pressure regulation system. These upgrades were designed to reduce thermal fluctuations and enhance long-term stability. The implemented upgrades significantly improved the intrinsic RV precision of the NIR channel. We quantified this improvement using Fabry Perot calibration spectra, obtaining an intrinsic RV precision of 0.67 ms after the interventions, an improvement of nearly 2 ms . We also assessed the stability of the nightly zero points, finding a reduced scatter of 3.9 ms post upgrade, compared to 6.1 ms before. For a sample of slowly rotating stars (vsin i below 2 kms), the median scatter decreased from 8.8 ms to 6.7 ms after the upgrades. These results demonstrate that the thermal control upgrades introduced in CARMENES PLUS have enhanced the NIR channels RV performance, bringing it closer to the VIS channels stability and reinforcing CARMENES capabilities for exoplanet detection around M dwarfs.

We investigate the stellar shape and size-mass relationship of X-ray selected Active Galactic Nuclei (AGN) host galaxies using the high-angular resolution and deep sensitivity in the near-infrared of the COSMOS-Web JWST survey field. We present the rest-frame 1-$\mu m$ size, stellar mass, Sersic index, axis-ratio, Gini-$M_{20}$ parameters of 690 moderate luminosity AGNs between redshift 0-3 and with stellar mass $\log M_s\sim 10.75$. We find that AGN host galaxies have an effective radius of 1-5 kpc, which is between star-forming (SFG) and quiescent galaxies (QGs) of the same stellar mass. AGN hosts have similar size-mass trends as SFG and QGs, being smaller at higher redshift for the same stellar mass. The slope of the size-mass relationship of AGN host galaxies is steeper than that of star-forming galaxies. Their rest-frame 1$\mu m$ stellar morphology indicates a significant spheroidal component. We observed a low merger fraction (6%) in our sample as well as substructures similar to disks, bars, and spiral arms in the residual images, which are in tension with evolutionary pathways that require major mergers. However, it may also be due to the different timescales between mergers and AGN activity.

We present James Webb Space Telescope (JWST) NIRSpec 1.7--5.5 micron observations of SN~2024ggi at +285.51 and +385.27 days post-explosion. The late-time nebular spectra are dominated by emission lines from various ionization states of H, Ca, Ar, C, Mg, Ni, Co, and Fe. We also detect strong CO emission in both the first overtone and fundamental vibrational bands. Most atomic features exhibit asymmetric line profiles, indicating an aspherical explosion. Using observed fluxes combined with non-LTE radiative-transfer simulations, we develop a data-driven method that resolves the complex molecular-emission region, constrains its 3D structure, and reproduces high-fidelity spectral profiles. We find that, CO is mostly formed prior to +285d past explosion. The subsequent evolution is dominated by the evaporation of CO with CO mass varying from M(CO) of 8.7E-3 to 1.3E-3 Mo, and with instabilities growing from almost homogeneous to highly clumped (density contrast f_c of 1.2 to 2). The minimum velocity of CO only slightly decreases between epochs (v_1 of 1200 and 1100 km/sec), with the reference temperature dropping from T_1 of 2400 and 1900K.

The ESA Euclid mission will measure the photometric redshifts of billions of galaxies in order to provide an accurate 3D view of the Universe at optical and near-infrared wavelengths. Photometric redshifts are determined by the PHZ processing function on the basis of the multi-wavelength photometry of Euclid and ground-based observations. In this paper, we describe the PHZ processing used for the Euclid Quick Data Release, the output products, and their validation. The PHZ pipeline is responsible for the following main tasks: source classification into star, galaxy, and QSO classes based on photometric colours; determination of photometric redshifts and of physical properties of galaxies. The classification is able to provide a star sample with a high level of purity, a highly complete galaxy sample, and reliable probabilities of belonging to those classes. The identification of QSOs is more problematic: photometric information seems to be insufficient to accurately separate QSOs from galaxies. The performance of the pipeline in the determination of photometric redshifts has been tested using the COSMOS2020 catalogue and a large sample of spectroscopic redshifts. The results are in line with expectations: the precision of the estimates are compatible with Euclid requirements, while, as expected, a bias correction is needed to achieve the accuracy level required for the cosmological probes. Finally, the pipeline provides reliable estimates of the physical properties of galaxies, in good agreement with findings from the COSMOS2020 catalogue, except for an unrealistically large fraction of very young galaxies with very high specific star-formation rates. The application of appropriate priors is, however, sufficient to obtain reliable physical properties for those problematic objects. We present several areas for improvement for future Euclid data releases.

This is the second paper in the HOWLS (higher-order weak lensing statistics) series exploring the usage of non-Gaussian statistics for cosmology inference within \textit{Euclid}. With respect to our first paper, we develop a full tomographic analysis based on realistic photometric redshifts which allows us to derive Fisher forecasts in the ($\sigma_8$, $w_0$) plane for a \textit{Euclid}-like data release 1 (DR1) setup. We find that the 5 higher-order statistics (HOSs) that satisfy the Gaussian likelihood assumption of the Fisher formalism (1-point probability distribution function, $\ell$1-norm, peak counts, Minkowski functionals, and Betti numbers) each outperform the shear 2-point correlation functions by a factor $2.5$ on the $w_0$ forecasts, with only marginal improvement when used in combination with 2-point estimators, suggesting that every HOS is able to retrieve both the non-Gaussian and Gaussian information of the matter density field. The similar performance of the different estimators\inlinecomment{, with a slight preference for Minkowski functionals and 1-point probability distribution function,} is explained by a homogeneous use of multi-scale and tomographic information, optimized to lower computational costs. These results hold for the $3$ mass mapping techniques of the \textit{Euclid} pipeline: aperture mass, Kaiser--Squires, and Kaiser--Squires plus, and are unaffected by the application of realistic star masks. Finally, we explore the use of HOSs with the Bernardeau--Nishimichi--Taruya (BNT) nulling scheme approach, finding promising results towards applying physical scale cuts to HOSs.

The Milky Way is known to contain a stellar bar, as are a significant fraction of disc galaxies across the universe. Our understanding of bar evolution, both theoretically and through analysis of simulations indicates that bars both grow in amplitude and slow down over time through interaction and angular momentum exchange with the galaxy's dark matter halo. Understanding the physical mechanisms underlying this coupling requires modelling of the structural deformations to the potential that are mutually induced between components. In this work we use Basis Function Expansion (BFE) in combination with multichannel Singular Spectral Analysis (mSSA) as a non-parametric analysis tool to illustrate the coupling between the bar and the dark halo in a single high-resolution isolated barred disc galaxy simulation. We demonstrate the power of mSSA to extract and quantify explicitly coupled dynamical modes, determining growth rates, pattern speeds and phase lags for different stages of evolution of the stellar bar and the dark matter response. BFE & mSSA together grant us the ability to explore the importance and physical mechanisms of bar-halo coupling, and other dynamically coupled structures across a wide range of dynamical environments.

The Large Magellanic Cloud (LMC) will induce a dynamical friction (DF) wake on infall to the Milky Way (MW). The MW's stellar halo will respond to the gravity of the LMC and the dark matter (DM) wake, forming a stellar counterpart to the DM wake. This provides a novel opportunity to constrain the properties of the DM particle. We present a suite of high-resolution, windtunnel-style simulations of the LMC's DF wake that compare the structure, kinematics, and stellar tracer response of the DM wake in cold DM (CDM), with and without self-gravity, vs. fuzzy DM (FDM) with $m_a = 10^{-23}$ eV. We conclude that the self-gravity of the DM wake cannot be ignored. Its inclusion raises the wake's density by $\sim 10\%$, and holds the wake together over larger distances ($\sim$ 50 kpc) than if self-gravity is ignored. The DM wake's mass is comparable to the LMC's infall mass, meaning the DM wake is a significant perturber to the dynamics of MW halo tracers. An FDM wake is more granular in structure and is $\sim 20\%$ dynamically colder than a CDM wake, but with comparable density. The granularity of an FDM wake increases the stars' kinematic response at the percent level compared to CDM, providing a possible avenue of distinguishing a CDM vs. FDM wake. This underscores the need for kinematic measurements of stars in the stellar halo at distances of 70-100 kpc.

We perform a frequentist analysis using the standard profile likelihood method for clustering measurements from Data Release 1 of the Dark Energy Spectroscopic Instrument (DESI). While Bayesian inferences for Effective Field Theory models of galaxy clustering can be highly sensitive to the choice of priors for extended cosmological models, frequentist inferences are not susceptible to such effects. We compare Bayesian and frequentist constraints for the parameter set $\{\sigma_8, H_0, \Omega_{\rm{m}}, w_0, w_a\}$ when fitting to the full-shape of the power spectrum multipoles, the post-reconstruction Baryon Acoustic Oscillation (BAO) measurements, as well as external datasets from the CMB and type Ia supernovae measurements. Bayesian prior effects are very significant for the $w_0w_a$CDM model; while the $1 \sigma$ frequentist confidence intervals encompass the maximum a posteriori (MAP), the Bayesian credible intervals almost always exclude the maximum likelihood estimate (MLE) and the MAP - indicating strong prior volume projection effects - unless supernovae data are included. We observe limited prior effects for the $\Lambda$CDM model, due to the reduced number of parameters. When DESI full-shape and BAO data are jointly fit, we obtain the following $1\sigma$ frequentist confidence intervals for $\Lambda$CDM ($w_0w_a$CDM): $\sigma_8 = 0.867^{+0.048}_{-0.041} , \ H_0 = 68.91^{+0.80}_{-0.79} \ \rm{km \ s^{-1}Mpc^{-1}} , \ \Omega_{\rm{m}} = 0.3038\pm0.0110$ ($\sigma_8 = 0.793^{+0.069}_{-0.048} , \ H_0 = 64.9^{+4.8}_{-2.8} \ \rm{km \ s^{-1}Mpc^{-1}} , \ \Omega_{\rm{m}} = 0.369^{+0.029}_{-0.059}$ , $w_0 = -0.24^{+0.17}_{-0.64}$ , $w_a = -2.5^{+1.9}_{}$), corresponding to 0.7$\sigma$, 0.3$\sigma$, 0.7$\sigma$ (1.9$\sigma$, 3.4$\sigma$, 5.6$\sigma$, 5.5$\sigma$, 5.6$\sigma$) shifts between the MLE relative to the Bayesian posterior mean for $\Lambda$CDM ($w_0w_a$CDM) respectively.

The Euclid mission aims to measure the positions, shapes, and redshifts of over a billion galaxies to provide unprecedented constraints on the nature of dark matter and dark energy. Achieving this goal requires a continuous reassessment of the mission's scientific performance, particularly in terms of its ability to constrain cosmological parameters, as our understanding of how to model large-scale structure observables improves. In this study, we present the first scientific forecasts using CLOE (Cosmology Likelihood for Observables in Euclid), a dedicated Euclid cosmological pipeline developed to support this endeavour. Using advanced Bayesian inference techniques applied to synthetic Euclid-like data, we sample the posterior distribution of cosmological and nuisance parameters across a variety of cosmological models and Euclid primary probes: cosmic shear, angular photometric galaxy clustering, galaxy-galaxy lensing, and spectroscopic galaxy clustering. We validate the capability of CLOE to produce reliable cosmological forecasts, showcasing Euclid's potential to achieve a figure of merit for the dark energy parameters $w_0$ and $w_a$ exceeding 400 when combining all primary probes. Furthermore, we illustrate the behaviour of the posterior probability distribution of the parameters of interest given different priors and scale cuts. Finally, we emphasise the importance of addressing computational challenges, proposing further exploration of innovative data science techniques to efficiently navigate the Euclid high-dimensional parameter space in upcoming cosmological data releases.

Many applications in transient science, gravitational wave follow-up, and galaxy population studies require all-sky galaxy catalogs with reliable distances, extents, and stellar masses. However, existing catalogs often lack completeness beyond $\sim 100$ Mpc, suffer from stellar contamination, or do not provide homogeneous stellar mass estimates and size information. Our goal is to build a high-purity, high-completeness, all-sky galaxy catalog out to 2,000 Mpc, specifically designed to support time-domain and multi-messenger astrophysics. We combined major galaxy catalogs and deep imaging surveys -- including the Legacy Surveys, Pan-STARRS, DELVE, and SDSS -- and added spectroscopic, photometric, and redshift-independent distances. We cleaned the sample using the Gaia catalog to remove stars and visually inspected all ambiguous cases below 100 Mpc through a classification platform that gathered 27,000 expert votes. Stellar masses were estimated using optical and mid-infrared profile-fit photometry, and we improved the accuracy of photometric distances by combining multiple independent estimates. The resulting catalog, REGALADE, includes nearly 80 million galaxies with distances under 2,000 Mpc. It provides stellar masses for 88% of the sample and ellipse fits for 80%. REGALADE is more than 90% complete for galaxies contributing 50% of the total $r$-band luminosity out to 360 Mpc. In science tests, it recovers 60% more known supernova hosts, doubles the number of low-luminosity transient hosts, and identifies more reliable hosts for ultraluminous and hyper-luminous X-ray sources. REGALADE is one of the most complete and reliable all-sky galaxy catalog to date for the nearby Universe, built for real-world applications in transient and multi-messenger astrophysics. The full dataset, visual classifications, and code will be released to support broad community use.

Euclid is expected to establish new state-of-the-art constraints on extensions beyond the standard LCDM cosmological model by measuring the positions and shapes of billions of galaxies. Specifically, its goal is to shed light on the nature of dark matter and dark energy. Achieving this requires developing and validating advanced statistical tools and theoretical prediction software capable of testing extensions of the LCDM model. In this work, we describe how the Euclid likelihood pipeline, Cosmology Likelihood for Observables in Euclid (CLOE), has been extended to accommodate alternative cosmological models and to refine the theoretical modelling of Euclid primary probes. In particular, we detail modifications made to CLOE to incorporate the magnification bias term into the spectroscopic two-point correlation function of galaxy clustering. Additionally, we explain the adaptations made to CLOE's implementation of Euclid primary photometric probes to account for massive neutrinos and modified gravity extensions. Finally, we present the validation of these CLOE modifications through dedicated forecasts on synthetic Euclid-like data by sampling the full posterior distribution and comparing with the results of previous literature. In conclusion, we have identified in this work several functionalities with regards to beyond-LCDM modelling that could be further improved within CLOE, and outline potential research directions to enhance pipeline efficiency and flexibility through novel inference and machine learning techniques.

The outer Galaxy is characterized by a lower metallicity than regions near the Sun, suggesting differences in the formation and survival of molecules in star-forming regions. To understand chemical evolution across the Milky Way, deriving molecular abundances in star-forming regions in the outer Galaxy is essential for refining models of sub-Solar metallicity environments. We analyzed IRAM 30 m observations at 3 and 2 mm toward 35 sources at Galactocentric distances of 9$-$24 kpc, within the "CHEMical complexity in star-forming regions of the outer Galaxy" (CHEMOUT) project. We focused on species with the highest detection rates (i.e., HCN, HCO$^+$, c-C$_3$H$_2$, H$^{13}$CO$^+$, HCO, SO) and searched for trends in column densities, abundances, and line widths with Galactocentric distance. Abundances for H$_2$CO and CH$_3$OH were updated using H$_2$ column densities from new NIKA2 dust maps. Fractional abundances relative to H$_2$ of most species (HCN, HCO$^+$, c-C$_3$H$_2$, HCO, H$_2$CO, CH$_3$OH) scale at most with the elemental carbon abundance ([C/H]) up to $\sim$24 kpc. SO shows a steeper gradient than sulfur abundance ([S/H]), while H$^{13}$CO$^+$ shows a shallower gradient than [$^{13}$C/H]. Gas turbulence, inferred from line widths, decreases with Galactocentric distance, suggesting a more quiescent environment in the outer Galaxy with respect to the inner Galaxy. In the outer Galaxy, the formation efficiency of most molecules, following the parent element availability, is comparable or higher (e.g., for H$^{13}$CO$^+$) than in the local Galaxy, whereas SO forms less efficiently. These results have significant implications for chemical models of the outermost star-forming regions and for understanding molecule formation under lower metallicity conditions.

We present reconstructed convergence maps, \textit{mass maps}, from the Dark Energy Survey (DES) third year (Y3) weak gravitational lensing data set. The mass maps are weighted projections of the density field (primarily dark matter) in the foreground of the observed galaxies. We use four reconstruction methods, each is a \textit{maximum a posteriori} estimate with a different model for the prior probability of the map: Kaiser-Squires, null B-mode prior, Gaussian prior, and a sparsity prior. All methods are implemented on the celestial sphere to accommodate the large sky coverage of the DES Y3 data. We compare the methods using realistic $\Lambda$CDM simulations with mock data that are closely matched to the DES Y3 data. We quantify the performance of the methods at the map level and then apply the reconstruction methods to the DES Y3 data, performing tests for systematic error effects. The maps are compared with optical foreground cosmic-web structures and are used to evaluate the lensing signal from cosmic-void profiles. The recovered dark matter map covers the largest sky fraction of any galaxy weak lensing map to date.

The computation of magnitudes and distances from low signal-to-noise observations is known to be problematic, in the sense that the magnitudes and distances tend to assume extreme values, or are even undefined or unphysical in the case of negative observed fluxes or parallaxes. In this work we show that magnitudes can be computed consistently at all signal-to-noise levels, and even for negative fluxes, if the prior information that the true flux or distance is non-negative is properly included. Furthermore, we derive an all-purpose estimator for distances from a prior implementing only the non-negativity of the true parallax. We apply our results to the case of combining magnitudes to colours, and magnitudes and distances to obtain absolute magnitudes. The resulting expressions are easy to compute and we show that the resulting distribution functions for magnitudes, colours, distances, and absolute magnitudes are not only consistent for all signal-to-noise levels and applicable to both, positive and negative observed fluxes and parallaxes, but also show no strong tails. While biases at very low signal-to-noise levels are unavoidable, the estimator for distances derived in this work is less biased than previously used estimators. We find that the magnitude, colour, distance, and absolute magnitude distributions for vanishing signals converge to limiting distributions, whose median values are important for assessing biases when working with data at low signal-to-noise levels.

We observed NH3 metastable inversion lines from (3, 3) to (18, 18) toward G0.66-0.13 in the Galactic center with the Shanghai Tianma 65m radio telescope and Yebes 40 m telescope. Highly-excited lines of NH3 (17, 17), (18, 18) were detected in emission for the first time in the interstellar medium, with upper energy levels up to 3100 K. Mapping observations reveal widespread hot molecular gas traced by NH3 (13, 13) toward G0.66-0.13. The rotation temperatures of hot gas traced by NH3 exceed 400 K, which amounts to five percent of the total NH3 in the Galactic Center. Hot gas (>400 K) and warm gas (100-140 K) are found in distinct clumps, with the hot gas located at the interfacing regions between different warm clouds. The theory of intermittency in turbulence reproduces the complex temperature structure in the central molecular zone, especially the hot gas observed here. The results presented here demonstrate that turbulence heating dominates the heating of the molecular gas in the Central Molecular Zone, while the turbulence is induced by the shear-motion of molecular clouds under the gravitational potential of the nuclear star clusters and the supermassive black hole. Our results suggest that shear-induced turbulence heating could be a widespread factor influencing galactic evolution.

The census of stellar streams and dwarf galaxies in the Milky Way provides direct constraints on galaxy formation models and the nature of dark matter. The DESI Milky Way survey (with a footprint of 14,000$~deg{^2}$ and a depth of r<19 mag) delivers the largest sample of distant metal-poor stars compared to previous optical fiber-fed spectroscopic surveys. This makes DESI an ideal survey to search for previously undetected streams and dwarf galaxies. We present a detailed characterization of the Cocytos stream, which was re-discovered using a clustering analysis with a catalog of giants in the DESI year 3 data, supplemented with Magellan/MagE spectroscopy. Our analysis reveals a relatively metal-rich ([Fe/H]$=-1.3$) and thick stream (width$=1.5^\circ$) at a heliocentric distance of $\approx 25$ kpc, with an internal velocity dispersion of 6.5-9 km s$^{-1}$. The stream's metallicity, radial orbit, and proximity to the Virgo stellar overdensities suggest that it is most likely a disrupted globular cluster that came in with the Gaia-Enceladus merger. We also confirm its association with the Pyxis globular cluster. Our result showcases the ability of wide-field spectroscopic surveys to kinematically discover faint disrupted dwarfs and clusters, enabling constraints on the dark matter distribution in the Milky Way.

Holographic superconductors have been studied so far in the absence of dynamical electromagnetic fields, namely in the limit in which they coincide with holographic superfluids. It is possible, however, to introduce dynamical gauge fields if a Neumann-type boundary condition is imposed on the AdS-boundary. In 3+1 dimensions, the dual theory is a 2+1 dimensional CFT whose spectrum contains a massless gauge field, signaling the emergence of a gauge symmetry. We study the impact of a dynamical gauge field in vortex configurations where it is known to significantly affect the energetics and phase transitions. We calculate the critical magnetic fields H_c1 and H_c2, obtaining that holographic superconductors are of Type II (H_c1 < H_c2). We extend the study to 4+1 dimensions where the gauge field does not appear as an emergent phenomena, but can be introduced, by a proper renormalization, as an external dynamical field. We also compare our predictions with those arising from a Ginzburg-Landau theory and identify the generic properties of Abrikosov vortices in holographic models.

WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366$-$959\,nm at $R\sim5000$, or two shorter ranges at $R\sim20\,000$. After summarising the design and implementation of WEAVE and its data systems, we present the organisation, science drivers and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy's origins by completing Gaia's phase-space information, providing metallicities to its limiting magnitude for $\sim$3 million stars and detailed abundances for $\sim1.5$ million brighter field and open-cluster stars; (ii) survey $\sim0.4$ million Galactic-plane OBA stars, young stellar objects and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey $\sim400$ neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionised gas in z&lt;0.5 cluster galaxies; (vi) survey stellar populations and kinematics in $\sim25\,000$ field galaxies at $0.3\lesssim z \lesssim 0.7$; (vii) study the cosmic evolution of accretion and star formation using &gt;1 million spectra of LOFAR-selected radio sources; (viii) trace structures using intergalactic/circumgalactic gas at z&gt;2. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator.