The Dark Energy Spectroscopic Instrument (DESI) will precisely constrain cosmic expansion and the growth of structure by collecting $\sim$40 million extra-galactic redshifts across $\sim$80\% of cosmic history and one third of the sky. The Emission Line Galaxy (ELG) sample, which will comprise about one-third of all DESI tracers, will be used to probe the Universe over the 0.6 < z < 1.6 range, which includes the $1.1

We present three separate void catalogs created using a volume-limited sample of the DESI Year 1 Bright Galaxy Survey. We use the algorithms VoidFinder and V2 to construct void catalogs out to a redshift of z=0.24. We obtain 1,461 interior voids with VoidFinder, 420 with V2 using REVOLVER pruning, and 295 with V2 using VIDE pruning. Comparing our catalog with an overlapping SDSS void catalog, we find generally consistent void properties but significant differences in the void volume overlap, which we attribute to differences in the galaxy selection and survey masks. These catalogs are suitable for studying the variation in galaxy properties with cosmic environment and for cosmological studies.

In this work we use Lagrangian perturbation theory to analyze the harmonic space galaxy clustering signal of Bright Galaxy Survey (BGS) and Luminous Red Galaxies (LRGs) targeted by the Dark Energy Spectroscopic Instrument (DESI), combined with the galaxy--galaxy lensing signal measured around these galaxies using Dark Energy Survey Year 3 source galaxies. The BGS and LRG galaxies are extremely well characterized by DESI spectroscopy and, as a result, lens galaxy redshift uncertainty and photometric systematics contribute negligibly to the error budget of our ``$2\times2$-point'' analysis. On the modeling side, this work represents the first application of the \texttt{spinosaurus} code, implementing an effective field theory model for galaxy intrinsic alignments, and we additionally introduce a new scheme (\texttt{MAIAR}) for marginalizing over the large uncertainties in the redshift evolution of the intrinsic alignment signal. Furthermore, this is the first application of a hybrid effective field theory (HEFT) model for galaxy bias based on the $\texttt{Aemulus}\, \nu$ simulations. Our main result is a measurement of the amplitude of the lensing signal, $S_8=\sigma_8 \left(\Omega_m/0.3\right)^{0.5} = 0.850^{+0.042}_{-0.050}$, consistent with values of this parameter derived from the primary CMB. This constraint is artificially improved by a factor of $51\%$ if we assume a more standard, but restrictive parameterization for the redshift evolution and sample dependence of the intrinsic alignment signal, and $63\%$ if we additionally assume the nonlinear alignment model. We show that when fixing the cosmological model to the best-fit values from Planck PR4 there is $> 5 \sigma$ evidence for a deviation of the evolution of the intrinsic alignment signal from the functional form that is usually assumed in cosmic shear and galaxy--galaxy lensing studies.

We present DESI observations of the inner halo of M31, which reveal the kinematics of a recent merger - a galactic immigration event - in exquisite detail. Of the 11,416 sources studied in 3.75 hour of on-sky exposure time, 7,438 are M31 sources with well measured radial velocities. The observations reveal intricate coherent kinematic structure in the positions and velocities of individual stars: streams, wedges, and chevrons. While hints of coherent structures have been previously detected in M31, this is the first time they have been seen with such detail and clarity in a galaxy beyond the Milky Way. We find clear kinematic evidence for shell structures in the Giant Stellar Stream, the Northeast Shelf and Western Shelf regions. The kinematics are remarkably similar to the predictions of dynamical models constructed to explain the spatial morphology of the inner halo. The results are consistent with the interpretation that much of the substructure in the inner halo of M31 is produced by a single galactic immigration event 1 - 2 Gyr ago. Significant numbers of metal-rich stars ([Fe/H]>-0.5) are present in all of the detected substructures, suggesting that the immigrating galaxy had an extended star formation history. We also investigate the ability of the shells and Giant Stellar Stream to constrain the gravitational potential of M31, and estimate the mass within a projected radius of 125 kpc to be {\rm log_{10}}\, M_{\rm NFW}(<125\,{\rm kpc})/M_\odot = 11.80_{-0.10}^{+0.12}. The results herald a new era in our ability to study stars on a galactic scale and the immigration histories of galaxies.

We study the [OII] profiles of emission line galaxies (ELGs) from the Early Data Release of the Dark Energy Spectroscopic Instrument (DESI). To this end, we decompose and classify the shape of [OII] profiles with the first two eigenspectra derived from Principal Component Analysis. Our results show that DESI ELGs have diverse line profiles which can be categorized into three main types: (1) narrow lines with a median width of ~50 km/s, (2) broad lines with a median width of ~80 km/s, and (3) two-redshift systems with a median velocity separation of ~150 km/s, i.e., double-peak galaxies. To investigate the connections between the line profiles and galaxy properties, we utilize the information from the COSMOS dataset and compare the properties of ELGs, including star-formation rate (SFR) and galaxy morphology, with the average properties of reference star-forming galaxies with similar stellar mass, sizes, and redshifts. Our findings show that on average, DESI ELGs have higher SFR and more asymmetrical/disturbed morphology than the reference galaxies. Moreover, we uncover a relationship between the line profiles, the excess SFR and the excess asymmetry parameter, showing that DESI ELGs with broader [OII] line profiles have more disturbed morphology and higher SFR than the reference star-forming galaxies. Finally, we discuss possible physical mechanisms giving rise to the observed relationship and the implications of our findings on the galaxy clustering measurements, including the halo occupation distribution modeling of DESI ELGs and the observed excess velocity dispersion of the satellite ELGs.

The acceleration process of massive particles as well as the production of very high energy (VHE) photons and neutrinos remains a fundamental challenge in astrophysics. We investigate the parsec-scale jet structure and magnetic field of the blazar PKS 1424+240, that has been selected on the basis of strong VHE gamma-ray emission and identified with one of the highest peaks in the IceCube 9-year neutrino sky. We analyze 15 GHz VLBA observations of this BL Lac object, stacking 42 polarization-sensitive images collected in 2009-2025 to enhance the signal and reveal persistent parsec-scale structure. Our observations uncover a rare scenario. The object is viewed inside the jet cone, very close to the axis of its relativistic jet, with a viewing angle of <0.6 deg. This effectively maximizes Doppler boosting to values ~30 and enhances both electromagnetic and neutrino emission in the observer's direction. Based on polarimetric observations, we unambiguously detect a net toroidal component of the jet's magnetic field, indicating a current carrying jet flowing almost directly towards our line of sight. Blazars with very small jet viewing angles offer a solution to the longstanding mismatch between Doppler factors inferred from low VLBI apparent jet speed and those derived from VHE observations -- the so-called "Doppler factor crisis". We show that relativistic beaming plays the critical role in the gamma-ray and neutrino emission of blazars, with direct implications for models of their multi-messenger emission.

We measure the clustering of Lyman Alpha Emitting galaxies (LAEs) selected from the One-hundred-square-degree DECam Imaging in Narrowbands (ODIN) survey, with spectroscopic follow-up from Dark Energy Spectroscopic Instrument (DESI). We use DESI spectroscopy to optimize our selection and to constrain the interloper fraction and redshift distribution of our narrow-band selected sources. We select samples of 4000 LAEs at z=2.45 and 3.1 in 9 sq.deg. centered on the COSMOS field with median LyA fluxes of 10^{-16}erg/s/cm2. Covariances and cosmological inferences are obtained from a series of mock catalogs built upon high-resolution N-body simulations that match the footprint, number density, redshift distribution and observed clustering of the sample. We find that both samples have a correlation length of r_0=(3.0\pm 0.2)Mpc/h. Within our fiducial cosmology these correspond to 3D number densities of 10^{-3} h^3/Mpc^3 and, from our mock catalogs, biases of 1.7 and 2.0 at z=2.45 and 3.1, respectively. We discuss the implications of these measurements for the use of LAEs as large-scale structure tracers for high-redshift cosmology.

We reproduce the galaxy clustering catalogue from the SDSS-III Baryon Oscillation Spectroscopic Survey Final Data Release (BOSS DR11&DR12) with high fidelity on all relevant scales in order to allow a robust analysis of baryon acoustic oscillations and redshift space distortions. We have generated (6,000) 12,288 MultiDark PATCHY BOSS (DR11) DR12 light-cones corresponding to an effective volume of $\sim192,000\,[h^{-1}\,{\rm Gpc}]^3$ (the largest ever simulated volume), including cosmic evolution in the redshift range from 0.15 to 0.75. The mocks have been calibrated using a reference galaxy catalogue based on the halo abundance matching modelling of the BOSS DR11&DR12 galaxy clustering data and on the data themselves. The production follows three steps. First, we apply the PATCHY code to generate a dark matter field and an object distribution including nonlinear stochastic galaxy bias. Secondly, we run the halo/stellar distribution reconstruction HADRON code to assign masses to the various objects. This step uses the mass distribution as a function of local density and non-local indicators (i.e., tidal field tensor eigenvalues and relative halo exclusion separation for massive objects) from the reference simulation applied to the corresponding PATCHY dark matter and galaxy distribution. Finally, we apply the SUGAR code to build the light cones. The resulting MultiDark PATCHY mock light cones reproduce the number density, selection function, survey geometry, and in general within 1 $\sigma$, for arbitrary stellar mass bins, the power spectrum up to $k=0.3\,h\,{\rm Mpc}^{-1}$, the two-point correlation functions down to a few Mpc scales, and the three-point statistics of the BOSS DR11&DR12 galaxy samples.

Upcoming and existing large-scale surveys of galaxies require accurate theoretical predictions of the dark matter clustering statistics for thousands of mock galaxy catalogs. We demonstrate that this goal can be achieve with our new Parallel Particle-Mesh (PM) Nbody code (PPM-GLAM) at a very low computational cost. We run about 15,000 simulations with ~2 billion particles that provide ~1% accuracy of the dark matter power spectra P(k) for wave-numbers up to k~ 1h/Mpc. Using this large data-set we study the power spectrum covariance matrix, the stepping stone for producing mock catalogs. In contrast to many previous analytical and numerical results, we find that the covariance matrix normalised to the power spectrum C(k,k')/P(k)P(k') has a complex structure of non-diagonal components. It has an upturn at small k, followed by a minimum at k=0.1-0.2h/Mpc. It also has a maximum at k=0.5-0.6h/Mpc. The normalised covariance matrix strongly evolves with redshift: C(k,k')~delta(t)^alpha P(k)P(k'), where delta is the linear growth factor and alpha ~ 1-1.25, which indicates that the covariance matrix depends on cosmological parameters. We also show that waves longer than 1Gpc have very little impact on the power spectrum and covariance matrix. This significantly reduces the computational costs and complexity of theoretical predictions: relatively small volume ~ (1Gpc)^3 simulations capture the necessary properties of dark matter clustering statistics. All the power spectra obtained from many thousands of our simulations are publicly available.

The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) gamma-ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE; E > 100 MeV) gamma-rays, poses intriguing questions on the location of the gamma-ray emitting region in this flat spectrum radio quasar. We present multifrequency data of 4C +21.35 collected from centimeter to VHE during 2010 to investigate the properties of this source and discuss a possible emission model. The first hint of detection at VHE was observed by MAGIC on 2010 May 3, soon after a gamma-ray flare detected by Fermi-LAT that peaked on April 29. The same emission mechanism may therefore be responsible for both the HE and VHE emission during the 2010 flaring episodes. Two optical peaks were detected on 2010 April 20 and June 30, close in time but not simultaneous with the two gamma-ray peaks, while no clear connection was observed between the X-ray an gamma-ray emission. An increasing flux density was observed in radio and mm bands from the beginning of 2009, in accordance with the increasing gamma-ray activity observed by Fermi-LAT, and peaking on 2011 January 27 in the mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C +21.35 for the two periods of the VHE detection and a quiescent state, using a one-zone model with the emission coming from a very compact region outside the broad line region. The three SEDs can be fit with a combination of synchrotron self-Compton and external Compton emission of seed photons from a dust torus, changing only the electron distribution parameters between the epochs. The fit of the optical/UV part of the spectrum for 2010 April 29 seems to favor an inner disk radius of <6 gravitational radii, as one would expect from a prograde-rotating Kerr black hole.

We present RadioAstron Space VLBI imaging observations of the BL Lac object S5 0716+71 made on January 3-4 2015 at a frequency of 22 GHz (wavelength $\lambda=1.3$ cm). The observations were made in the framework of the AGN Polarization Key Science Program. The source was detected on projected space-ground baselines up to 70 833 km (5.6 Earth diameters) for both, parallel hand and cross-hand interferometric visibilities. We have used these detections to obtain a full-polarimetric image of the blazar at an unprecedented angular resolution of 24 $\mu$as, the highest for this source to date. This enabled us to estimate the size of the radio core to be &lt;12\times5~\muas and to reveal a complex structure and a significant curvature of the blazar jet in the inner 100 $\mu$as, which is an indication that the jet viewing angle lies inside the opening angle of the jet conical outflow. Fairly highly (15%) linearly polarized emission is detected in a jet region of 19 $\mu$as in size, located 58 $\mu$as downstream from the core. The highest brightness temperature in the source frame is estimated to be &gt;2.2\times10^{13} K for the blazar core. This implies that the inverse Compton limit must be violated in the rest frame of the source, even for the largest Doppler factor $\delta\thicksim25$ reported for 0716+714.

We report on atomic gas (HI) and molecular gas (as traced by CO(2-1)) redshifted absorption features toward the nuclear regions of the closest powerful radio galaxy, Centaurus A (NGC 5128). Our HI observations using the Very Long Baseline Array allow us to discern with unprecedented sub-parsec resolution HI absorption profiles toward different positions along the 21 cm continuum jet emission in the inner 0."3 (or 5.4 pc). In addition, our CO(2-1) data obtained with the Submillimeter Array probe the bulk of the absorbing molecular gas with little contamination by emission, not possible with previous CO single-dish observations. We shed light with these data on the physical properties of the gas in the line of sight, emphasizing the still open debate about the nature of the gas that produces the broad absorption line (~55 km/s). First, the broad H I line is more prominent toward the central and brightest 21 cm continuum component than toward a region along the jet at a distance ~ 20 mas (or 0.4 pc) further from it. This suggests that the broad absorption line arises from gas located close to the nucleus, rather than from diffuse and more distant gas. Second, the different velocity components detected in the CO(2-1) absorption spectrum match well other molecular lines, such as those of HCO+(1-0), except the broad absorption line that is detected in HCO+(1-0) (and most likely related to that of the H I). Dissociation of molecular hydrogen due to the AGN seems to be efficient at distances <= 10 pc, which might contribute to the depth of the broad H I and molecular lines.

The Daniel K. Inouye Solar Telescope (DKIST) will revolutionize our ability to measure, understand and model the basic physical processes that control the structure and dynamics of the Sun and its atmosphere. The first-light DKIST images, released publicly on 29 January 2020, only hint at the extraordinary capabilities which will accompany full commissioning of the five facility instruments. With this Critical Science Plan (CSP) we attempt to anticipate some of what those capabilities will enable, providing a snapshot of some of the scientific pursuits that the Daniel K. Inouye Solar Telescope hopes to engage as start-of-operations nears. The work builds on the combined contributions of the DKIST Science Working Group (SWG) and CSP Community members, who generously shared their experiences, plans, knowledge and dreams. Discussion is primarily focused on those issues to which DKIST will uniquely contribute.

We report on high spatial resolution observations, using the Australia Telescope Compact Array (ATCA), of ground-state OH masers. These observations were carried out toward 196 pointing centres previously identified in the Southern Parkes Large-Area Survey in Hydroxyl (SPLASH) pilot region, between Galactic longitudes of $334^{\circ}$ and $344^{\circ}$ and Galactic latitudes of $-2^{\circ}$ and $+2^{\circ}$. Supplementing our data with data from the MAGMO (Mapping the Galactic Magnetic field through OH masers) survey, we find maser emission towards 175 of the 196 target fields. We conclude that about half of the 21 non-detections were due to intrinsic variability. Due to the superior sensitivity of the follow-up ATCA observations, and the ability to resolve nearby sources into separate sites, we have identified 215 OH maser sites towards the 175 fields with detections. Among these 215 OH maser sites, 111 are new detections. After comparing the positions of these 215 maser sites to the literature, we identify 122 (57 per cent) sites associated with evolved stars (one of which is a planetary nebula), 64 (30 per cent) with star formation, two sites with supernova remnants and 27 (13 per cent) of unknown origin. The infrared colors of evolved star sites with symmetric maser profiles tend to be redder than those of evolved star sites with asymmetric maser profiles, which may indicate that symmetric sources are generally at an earlier evolutionary stage.

Since the early 1970s, inversion techniques have become the most useful tool for inferring the magnetic, dynamic, and thermodynamic properties of the solar atmosphere. The intrinsic model dependence makes it necessary to formulate specific means that include the physics in a properly quantitative way. The core of this physics lies in the radiative transfer equation (RTE), where the properties of the atmosphere are assumed to be known while the unknowns are the four Stokes profiles. The solution of the (differential) RTE is known as the direct or forward problem. From an observational point of view, the problem is rather the opposite: the data are made up of the observed Stokes profiles and the unknowns are the solar physical quantities. Inverting the RTE is therefore mandatory. Indeed, the formal solution of this equation can be considered an integral equation. The solution of such an integral equation is called the inverse problem. Inversion techniques are automated codes aimed at solving the inverse problem. The foundations of inversion techniques are critically revisited with an emphasis on making explicit the many assumptions underlying each of them. An incremental complexity procedure is advised for the implementation in practice. Coarse details of the profiles or coarsely sam- pled profiles should be reproduced first with simple model atmospheres (with, for example, a few physical quantities that are constant with optical depth). If the Stokes profiles are well sampled and differences between synthetic and observed ones are greater than the noise, then the inversion should proceed by using more complex models (that is, models where physical quantities vary with depth or, eventually, with more than one component). Significant im- provements are expected as well from the use of new inversion techniques that take the spatial degradation by the instruments into account.

Type-2 Low-ionization Narrow Emission-line Regions (LINERs) have been optically classified with the Palomar data as not presenting a broad component in the Balmer emission lines associated to the Broad Line Region (BLR) of the Active Galactic Nuclei (AGN). We unveil the presence of different kinematic components of emission lines in the nuclear region of a sample of local (z$\leq$0.022) type-2 LINERs. We focus on the analysis of their true nature by means of the detection (or not) of a broad component originated in the BLR. Additionally, we search for the possible presence of non-rotational motions such as outflows. We have applied a decomposition of the nuclear emission lines using a spectroscopic analysis of the optical spectra of 9 type-2 LINERs of intermediate-resolution data retrieved from the Hubble Space Telescope (HST) archive. The study is completed with archival spectra from the Double Spectrograph from the Palomar Observatory. The emission line fitting reveals the presence of a broad component associated to the BLR in 6 out of the 9 galaxies for the space-based data, and for 2 out of the 8 from the ground-based spectra. The velocity dispersion of two galaxies (NGC 4486 and NGC 4594) measured in HST/STIS data suggest the presence of outflows. The results indicate that the spatial resolution plays a major role in the detection of the BLR, as it appears diluted in the ground-based data (even after removing stellar contribution). This is also true for the emission line diagnostics, as the contaminant light contributes to lower emission line ratios towards the star-forming area of standard BPTs. We propose to reclassify NGC4594 as a type-1 LINER, since a BLR component is seen in both space- and ground-based spectra. We find ambiguous results for the BLR component of NGC 4486. The modest outflow detection in our sample may indicate that they are not as frequent as for type-1 LINERs.

Magneto-static models may overcome some of the issues facing force-free magnetic field extrapolations. So far they have seen limited use and have faced problems when applied to quiet-Sun data. Here we present a first application to an active region. We use solar vector magnetic field measurements gathered by the IMaX polarimeter during the flight of the \sunrise{} balloon-borne solar observatory in June 2013 as boundary condition for a magneto-static model of the higher solar atmosphere above an active region. The IMaX data are embedded in active region vector magnetograms observed with SDO/HMI. This work continues our magneto-static extrapolation approach, which has been applied earlier ({\it Paper I}) to a quiet Sun region observed with \sunrise{} I. In an active region the signal-to-noise-ratio in the measured Stokes parameters is considerably higher than in the quiet Sun and consequently the IMaX measurements of the horizontal photospheric magnetic field allow us to specify the free parameters of the model in a special class of linear magneto-static equilibria. The high spatial resolution of IMaX (110-130 km, pixel size 40 km) enables us to model the non-force-free layer between the photosphere and the mid chromosphere vertically by about 50 grid points. In our approach we can incorporate some aspects of the mixed beta layer of photosphere and chromosphere, e.g., taking a finite Lorentz force into account, which was not possible with lower resolution photospheric measurements in the past. The linear model does not, however, permit to model intrinsic nonlinear structures like strongly localized electric currents.

We present two supernovae (SNe) discovered with the Hubble Space Telescope (HST) in the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS), an HST multi-cycle treasury program. We classify both objects as Type Ia SNe and find redshifts of z = 1.80+-0.02 and 2.26 +0.02 -0.10, the latter of which is the highest redshift Type Ia SN yet seen. Using light curve fitting we determine luminosity distances and find that both objects are consistent with a standard Lambda-CDM cosmological model. These SNe were observed using the HST Wide Field Camera 3 infrared detector (WFC3-IR), with imaging in both wide- and medium-band filters. We demonstrate that the classification and redshift estimates are significantly improved by the inclusion of single-epoch medium-band observations. This medium-band imaging approximates a very low resolution spectrum (lambda/delta lambda ~ 100) which can isolate broad spectral absorption features that differentiate Type Ia SNe from their most common core collapse cousins. This medium-band method is also insensitive to dust extinction and (unlike grism spectroscopy) it is not affected by contamination from the SN host galaxy or other nearby sources. As such, it can provide a more efficient - though less precise - alternative to IR spectroscopy for high-z SNe.