We produce a clean and well-characterised catalogue of objects within 100\,pc of the Sun from the \G\ Early Data Release 3. We characterise the catalogue through comparisons to the full data release, external catalogues, and simulations. We carry out a first analysis of the science that is possible with this sample to demonstrate its potential and best practices for its use. The selection of objects within 100\,pc from the full catalogue used selected training sets, machine-learning procedures, astrometric quantities, and solution quality indicators to determine a probability that the astrometric solution is reliable. The training set construction exploited the astrometric data, quality flags, and external photometry. For all candidates we calculated distance posterior probability densities using Bayesian procedures and mock catalogues to define priors. Any object with reliable astrometry and a non-zero probability of being within 100\,pc is included in the catalogue. We have produced a catalogue of \NFINAL\ objects that we estimate contains at least 92\% of stars of stellar type M9 within 100\,pc of the Sun. We estimate that 9\% of the stars in this catalogue probably lie outside 100\,pc, but when the distance probability function is used, a correct treatment of this contamination is possible. We produced luminosity functions with a high signal-to-noise ratio for the main-sequence stars, giants, and white dwarfs. We examined in detail the Hyades cluster, the white dwarf population, and wide-binary systems and produced candidate lists for all three samples. We detected local manifestations of several streams, superclusters, and halo objects, in which we identified 12 members of \G\ Enceladus. We present the first direct parallaxes of five objects in multiple systems within 10\,pc of the Sun.

We analyze the evolution of massive (log$_{10}$ [$M_\star/M_\odot$] $>10$) galaxies at $z \sim$ 4--8 selected from the JWST Cosmic Evolution Early Release Science (CEERS) survey. We infer the physical properties of all galaxies in the CEERS NIRCam imaging through spectral energy distribution (SED) fitting with dense basis to select a sample of high redshift massive galaxies. Where available we include constraints from additional CEERS observing modes, including 18 sources with MIRI photometric coverage, and 28 sources with spectroscopic confirmations from NIRSpec or NIRCam wide-field slitless spectroscopy. We sample the recovered posteriors in stellar mass from SED fitting to infer the volume densities of massive galaxies across cosmic time, taking into consideration the potential for sample contamination by active galactic nuclei (AGN). We find that the evolving abundance of massive galaxies tracks expectations based on a constant baryon conversion efficiency in dark matter halos for $z \sim$ 1--4. At higher redshifts, we observe an excess abundance of massive galaxies relative to this simple model. These higher abundances can be explained by modest changes to star formation physics and/or the efficiencies with which star formation occurs in massive dark matter halos, and are not in tension with modern cosmology.

GSN 069 is the first galactic nucleus where quasi-periodic eruptions (QPEs) have been identified. These are high-amplitude, soft X-ray bursts recurring every ~9 hr, lasting ~1 hr, and during which the X-ray count rate increases by up to two orders of magnitude with respect to an otherwise stable quiescent level. The X-ray spectral properties and the long-term evolution of GSN 069 in the first few years are consistent with a long-lived tidal disruption event (TDE). Here we derive the properties of QPEs and of the long-term X-ray evolution in GSN 069 over the past 12 yr by studying timing and spectral X-ray data from 11 XMM-Newton, one Chandra, and 34 Swift observations on timescales ranging from minutes to years. QPEs in GSN 069 are a transient phenomenon with a lifetime > 1.05 yr. The QPE intensity and recurrence time oscillate and allow for alternating strong-weak QPEs and long-short recurrence times to be defined. In observations with QPEs, the quiescent level exhibits a quasi-periodic oscillation with a period equal to the average separation between consecutive QPEs. The QPE spectral evolution is consistent with thermal emission from a very compact region that heats up quickly and subsequently cools down via X-ray emission while expanding by a factor of ~3 in radius. The long-term evolution of the quiescent level is characterised by two repeating TDEs ~9 yr apart. We detect a precursor X-ray flare prior to the second TDE that may be associated with the circularisation phase during disc formation. A similar precursor flare is tentatively detected just before the first TDE. Future X-ray observations of GSN 069 promise that the QPE origin and the relation between QPEs and repeating TDEs in this galactic nucleus will be constrained, with consequences for the other sources where QPEs have been identified. [abridged]

To date, galaxy image simulations for weak lensing surveys usually approximate the light profiles of all galaxies as a single or double Sérsic profile, neglecting the influence of galaxy substructures and morphologies deviating from such a simplified parametric characterization. While this approximation may be sufficient for previous data sets, the stringent cosmic shear calibration requirements and the high quality of the data in the upcoming Euclid survey demand a consideration of the effects that realistic galaxy substructures have on shear measurement biases. Here we present a novel deep learning-based method to create such simulated galaxies directly from HST data. We first build and validate a convolutional neural network based on the wavelet scattering transform to learn noise-free representations independent of the point-spread function of HST galaxy images that can be injected into simulations of images from Euclid's optical instrument VIS without introducing noise correlations during PSF convolution or shearing. Then, we demonstrate the generation of new galaxy images by sampling from the model randomly and conditionally. Next, we quantify the cosmic shear bias from complex galaxy shapes in Euclid-like simulations by comparing the shear measurement biases between a sample of model objects and their best-fit double-Sérsic counterparts. Using the KSB shape measurement algorithm, we find a multiplicative bias difference between these branches with realistic morphologies and parametric profiles on the order of $6.9\times 10^{-3}$ for a realistic magnitude-Sérsic index distribution. Moreover, we find clear detection bias differences between full image scenes simulated with parametric and realistic galaxies, leading to a bias difference of $4.0\times 10^{-3}$ independent of the shape measurement method. This makes it relevant for stage IV weak lensing surveys such as Euclid.

We present a new infrared survey covering the three Euclid deep fields and four other Euclid calibration fields using Spitzer's Infrared Array Camera (IRAC). We have combined these new observations with all relevant IRAC archival data of these fields in order to produce the deepest possible mosaics of these regions. In total, these observations represent nearly 11% of the total Spitzer mission time. The resulting mosaics cover a total of approximately 71.5deg$^2$ in the 3.6 and 4.5um bands, and approximately 21.8deg$^2$ in the 5.8 and 8um bands. They reach at least 24 AB magnitude (measured to sigma, in a 2.5 arcsec aperture) in the 3.6um band and up to ~ 5 mag deeper in the deepest regions. The astrometry is tied to the Gaia astrometric reference system, and the typical astrometric uncertainty for sources with 16<[3.6]<19 is <0.15 arcsec. The photometric calibration is in excellent agreement with previous WISE measurements. We have extracted source number counts from the 3.6um band mosaics and they are in excellent agreement with previous measurements. Given that the Spitzer Space Telescope has now been decommissioned these mosaics are likely to be the definitive reduction of these IRAC data. This survey therefore represents an essential first step in assembling multi-wavelength data on the Euclid deep fields which are set to become some of the premier fields for extragalactic astronomy in the 2020s.

Quasi-periodic eruptions (QPEs) are recurring soft X-ray transients emerging from the vicinity of supermassive black holes (SMBHs) in nearby, low-mass galaxy nuclei; about ten QPE hosts have been identified thus far. Here we report the \textit{NICER} discovery of QPEs in the optically-selected Tidal Disruption Event (TDE) and Extreme Coronal Line Emitter (ECLE) AT2022upj, exhibiting a large spread in recurrence times from 0.5-3.5 days, durations from 0.3-1 days, peak luminosities from $10^{42.5-43.0}$ erg s$^{-1}$, and erratic flare profiles. A wealth of evidence now links at least some QPEs to the newly-formed accretion flows emerging from TDEs; AT2022upj is the third QPE reported in an optically discovered TDE. Marginalizing over the uncertain distributions of QPE peak luminosity, recurrence time, delay after TDE peak, and lifetime, we use the burgeoning sample to make a Bayesian estimate that the fraction of optical TDEs resulting in QPEs within 5 years post-disruption is $9^{+9}_{-5}$\%. Along with AT2019qiz, AT2022upj also marks the second of the three optical TDE+X-ray QPEs showing coronal line emission, suggesting ECLEs may represent a subset of TDEs particularly efficient at forming QPEs and/or that sustained QPE X-ray emission contributes to coronal line emission in some galaxy nuclei.

We investigate the accuracy and range of validity of the perturbative model for the 2-point (2PCF) and 3-point (3PCF) correlation functions in real space in view of the forthcoming analysis of the Euclid mission spectroscopic sample. We take advantage of clustering measurements from four snapshots of the Flagship I N-body simulations at z = {0.9, 1.2, 1.5, 1.8}, which mimic the expected galaxy population in the ideal case of absence of observational effects such as purity and completeness. For the 3PCF we consider all available triangle configurations given a minimal separation. First, we assess the model performance by fixing the cosmological parameters and evaluating the goodness-of-fit provided by the perturbative bias expansion in the joint analysis of the two statistics, finding overall agreement with the data down to separations of 20 Mpc/h. Subsequently, we build on the state-of-the-art and extend the analysis to include the dependence on three cosmological parameters: the amplitude of scalar perturbations As, the matter density {\omega}cdm and the Hubble parameter h. To achieve this goal, we develop an emulator capable of generating fast and robust modelling predictions for the two summary statistics, allowing efficient sampling of the joint likelihood function. We therefore present the first joint full-shape analysis of the real-space 2PCF and 3PCF, testing the consistency and constraining power of the perturbative model across both probes, and assessing its performance in a combined likelihood framework. We explore possible systematic uncertainties induced by the perturbative model at small scales finding an optimal scale cut of rmin = 30 Mpc/h for the 3PCF, when imposing an additional limitation on nearly isosceles triangular configurations included in the data vector. This work is part of a Euclid Preparation series validating theoretical models for galaxy clustering.

Local Universe dwarf galaxies are both cosmological and mass assembly probes. Deep surveys have enabled the study of these objects down to the low surface brightness (LSB) regime. In this paper, we estimate Euclid's dwarf detection capabilities as well as limits of its MERge processing function (MER pipeline), responsible for producing the stacked mosaics and final catalogues. To do this, we inject mock dwarf galaxies in a real Euclid Wide Survey (EWS) field in the VIS band and compare the input catalogue to the final MER catalogue. The mock dwarf galaxies are generated with simple Sérsic models and structural parameters extracted from observed dwarf galaxy property catalogues. To characterize the detected dwarfs, we use the mean surface brightness inside the effective radius SBe (in mag arcsec-2). The final MER catalogues achieve completenesses of 91 % for SBe in [21, 24], and 54 % for SBe in [24, 28]. These numbers do not take into account possible contaminants, including confusion with background galaxies at the location of the dwarfs. After taking into account those effects, they become respectively 86 % and 38 %. The MER pipeline performs a final local background subtraction with small mesh size, leading to a flux loss for galaxies with Re > 10". By using the final MER mosaics and reinjecting this local background, we obtain an image in which we recover reliable photometric properties for objects under the arcminute scale. This background-reinjected product is thus suitable for the study of Local Universe dwarf galaxies. Euclid's data reduction pipeline serves as a test bed for other deep surveys, particularly regarding background subtraction methods, a key issue in LSB science.

We use three semi-analytic models (SAMs) of galaxy formation and evolution, run on the same 1$h^{-1}$Gpc MultiDark Planck2 cosmological simulation, to investigate the properties of [OII] emission line galaxies in the redshift range $0.6

The accretion of matter onto black holes and neutron stars often leads to the launching of outflows that can greatly affect the environments surrounding the compact object. In supermassive black holes, these outflows can even be powerful enough to dictate the evolution of the entire host galaxy, and yet, to date, we do not understand how these so-called accretion disk winds are launched - whether by radiation pressure, magnetic forces, thermal irradiation, or a combination thereof. An important means of studying disk winds produced near the central compact object is through X-ray absorption line spectroscopy, which allows us to probe outflow properties along a single line of sight, but usually provides little information about the global 3D disk wind structure that is vital for understanding the launching mechanism and total wind energy budget. Here, we study Hercules X-1, a unique, nearly edge-on X-ray binary with a warped accretion disk precessing with a period of about 35 days. This disk precession results in changing sightlines towards the neutron star, through the ionized outflow. We perform time-resolved X-ray spectroscopy over the precession phase and detect a strong decrease in the wind column density by three orders of magnitude as our sightline progressively samples the wind at greater heights above the accretion disk. The wind becomes clumpier as it rises upwards and expands away from the neutron star. Modelling the warped disk shape, we create a 2D map of wind properties. This unique measurement of the vertical structure of an accretion disk wind allows direct comparisons to 3D global simulations to reveal the outflow launching mechanism.

We report on NuSTAR observations of transient pulsations in the neutron star X-ray binary SMC X-1. The transition from non-pulsing to pulsing states was not accompanied by a large change in flux. Instead, both pulsing and non-pulsing states were observed in a single observation during the low-flux super-orbital state. During the high-state, we measure a pulse period of $P = 0.70117(9)\,\mathrm{s}$at$T_{ref} = 56145\,\mathrm{MJD}$. Spectral analysis during non-pulsing and pulsing states reveals that the observations can be consistently modeled by an absorbed power law with a phenomenological cutoff resembling a Fermi-Dirac distribution, or by a partially obscured cutoff power law. The shapes of the underlying continua show little variability between epochs, while the covering fraction and column density vary between super-orbital states. The strength of pulsations also varies, leading us to infer that the absence and reemergence of pulsations are related to changing obscuration, such as by a warped accretion disk. SMC X-1 is accreting near or above its Eddington limit, reaching an unabsorbed X-ray luminosity of $L_{\rm X}({\rm 2-10~keV}) \approx 5 \times 10^{38}\, {\rm erg}\, {\rm s}^{-1}$. This suggests that SMC X-1 may be a useful local analog to ultraluminous X-ray pulsars (ULXPs), which likewise exhibit strong variability in their pulsed fractions, as well as flux variability on similar timescales. In particular, the gradual pulse turn-on which has been observed in M82 X-2 is similar to the behavior we observe in SMC X-1. Thus we propose that pulse fraction variability of ULXPs may also be due to variable obscuration.

Ultraluminous X-ray sources (ULXs) are X-ray binary systems containing an accreting neutron star (NS) or black hole emitting at luminosities above the Eddington limit of a $10M_{\odot}$ black hole. Approximately 1900 (either confirmed or candidate) ULXs have been identified to date. Three systems have been confirmed to exhibit coherent signals consistent with NS spin frequencies and quasi-periodic oscillations (QPOs) in the mHz range. Several interpretations for generating such QPOs have been proposed, including general relativistic frame-dragging effects. In this work, we test if an alternative model in which magnetically-driven precession of the inner accretion flow can self-consistently reproduce the observed NS spin and QPO frequencies for reasonable values for accretion rates and NS magnetic field strengths. For a range of parameters, we recover family of solutions with accretion rates $\approx10^{-7}$--$10^{-5}$\,M$_{\odot}$\,yr$^{-1}$ and surface magnetic fields $\gtrsim10^{12}$\,G, in agreement with previous estimates. If validated, this interpretation could reconcile several observed properties of pulsating ULXs, including QPO frequencies and the observed high luminosities of these systems, in a self-consistent framework without requiring general relativistic effects and/or strong beaming due to specific viewing angles. Although the predictive power of the model is currently limited by parameter degeneracies and uncertainties, searching for and discovering more pulsating ULX systems will allow to further test or refute the proposed model.

Planetary systems orbiting close binaries are valuable testing grounds for planet formation and migration models. More detections with good mass measurements are needed. We present a new planet discovered during the BEBOP survey for circumbinary exoplanets using radial velocities. We use data taken with the SOPHIE spectrograph at the Observatoire de Haute-Provence, and perform a spectroscopic analysis to obtain high precision radial velocities. This planet is the first radial velocity detection of a previously unknown circumbinary system. The planet has a mass of $0.56$ $M_{Jup}$ and orbits its host binary in 550 days with an eccentricity of 0.25. Compared to most of the previously known circumbinary planets, BEBOP-3b has a long period (relative to the binary) and a high eccentricity. There also is a candidate outer planet with a $\sim1400$ day orbital period. We test the stability of potential further candidate signals inside the orbit of BEBOP-3b, and demonstrate that there are stable orbital solutions for planets near the instability region which is where the Kepler circumbinary planets are located. We also use our data to obtain independent dynamical masses for the two stellar components of the eclipsing binary using High Resolution Cross-Correlation Spectroscopy (HRCCS), and compare those results to a more traditional approach, finding them compatible with one another.

The Circinus galaxy is one of the nearest obscured AGN, making it an ideal target for detailed study. Combining archival Chandra and XMM-Newton data with new NuSTAR observations, we model the 2-79 keV spectrum to constrain the primary AGN continuum and to derive physical parameters for the obscuring material. Chandra's high angular resolution allows a separation of nuclear and off-nuclear galactic emission. In the off-nuclear diffuse emission we find signatures of strong cold reflection, including high equivalent-width neutral Fe lines. This Compton-scattered off-nuclear emission amounts to 18% of the nuclear flux in the Fe line region, but becomes comparable to the nuclear emission above 30 keV. The new analysis no longer supports a prominent transmitted AGN component in the observed band. We find that the nuclear spectrum is consistent with Compton-scattering by an optically-thick torus, where the intrinsic spectrum is a powerlaw of photon index $\Gamma = 2.2-2.4$, the torus has an equatorial column density of $N_{\rm H} = (6-10)\times10^{24}$cm$^{-2}$ and the intrinsic AGN $2-10$ keV luminosity is $(2.3-5.1)\times 10^{42}$ erg/s. These values place Circinus along the same relations as unobscured AGN in accretion rate-vs-$\Gamma$ and $L_X$-vs-$L_{IR}$ phase space. NuSTAR's high sensitivity and low background allow us to study the short time-scale variability of Circinus at X-ray energies above 10 keV for the first time. The lack of detected variability favors a Compton-thick absorber, in line with the the spectral fitting results.

We report on the Swift discovery of a second high-amplitude (factor 100) outburst of the Seyfert 1.9 galaxy IC 3599, and discuss implications for outburst scenarios. Swift detected this active galactic nucleus (AGN) again in February 2010 in X-rays at a level of (1.50\plm0.11)$\times 10^{36}$ W (0.2-2.0 keV), which is nearly as luminous as the first outburst detected with ROSAT in 1990. Optical data from the Catalina sky survey show that the optical emission was already bright two years before the Swift X-ray high-state. Our new Swift observations performed between 2013 and 2015 show that IC 3599 is currently again in a very low X-ray flux state. This repeat optical and X-ray outburst, and the long optical duration, suggest that IC 3599 is likely not a tidal disruption event (TDE). Instead, variants of AGN-related variability are explored. The data are consistent with an accretion disk instability around a black hole of mass on the order 10$^6$--10$^7$ M$_{\odot}$; a value estimated using several different methods.

The recent gravitational wave measurements have demonstrated the existence of stellar mass black hole binaries. It is essential for our understanding of massive star evolution to identify the contribution of binary evolution to the formation of double black holes. A promising way to progress is investigating the progenitors of double black hole systems and comparing predictions with local massive star samples such as the population in 30 Doradus in the Large Magellanic Cloud (LMC). Methods. To this purpose, we analyse a large grid of detailed binary evolution models at LMC metallicity with initial primary masses between 10 and 40 Msun, and identify which model systems potentially evolve into a binary consisting of a black hole and a massive main sequence star. We then derive the observable properties of such systems, as well as peculiarities of the OB star component. We find that about 3% of the LMC late O and early B stars in binaries are expected to possess a black hole companion, when assuming stars with a final helium core mass above 6.6 M to form black holes. While the vast majority of them may be X-ray quiet, our models suggest that these may be identified in spectroscopic binaries, either by large amplitude radial velocity variations ( > 50 km s ) and simultaneous nitrogen surface enrichment, or through a moderate radial velocity ( > 10 km/s ) and simultaneously rapid rotation of the OB star. The predicted mass ratios are such that main sequence companions could be excluded in most cases. A comparison to the observed OB+WR binaries in the LMC, Be/X-ray binaries, and known massive BH binaries supports our conclusion. We expect spectroscopic observations to be able to test key assumptions in our models, with important implications for massive star evolution in general, and for the formation of double-black hole mergers in particular.

Gaia-CRF3 is the celestial reference frame for positions and proper motions in the third release of data from the Gaia mission, Gaia DR3 (and for the early third release, Gaia EDR3, which contains identical astrometric results). The reference frame is defined by the positions and proper motions at epoch 2016.0 for a specific set of extragalactic sources in the (E)DR3 catalogue. We describe the construction of Gaia-CRF3, and its properties in terms of the distributions in magnitude, colour, and astrometric quality. Compact extragalactic sources in Gaia DR3 were identified by positional cross-matching with 17 external catalogues of quasars (QSO) and active galactic nuclei (AGN), followed by astrometric filtering designed to remove stellar contaminants. Selecting a clean sample was favoured over including a higher number of extragalactic sources. For the final sample, the random and systematic errors in the proper motions are analysed, as well as the radio-optical offsets in position for sources in the third realisation of the International Celestial Reference Frame (ICRF3). The Gaia-CRF3 comprises about 1.6 million QSO-like sources, of which 1.2 million have five-parameter astrometric solutions in Gaia DR3 and 0.4 million have six-parameter solutions. The sources span the magnitude range G = 13 to 21 with a peak density at 20.6 mag, at which the typical positional uncertainty is about 1 mas. The proper motions show systematic errors on the level of 12 ${\mu}$as yr${}^{-1}$ on angular scales greater than 15 deg. For the 3142 optical counterparts of ICRF3 sources in the S/X frequency bands, the median offset from the radio positions is about 0.5 mas, but exceeds 4 mas in either coordinate for 127 sources. We outline the future of the Gaia-CRF in the next Gaia data releases.

Radio observations allow us to identify a wide range of active galactic nuclei (AGN), which play a significant role in the evolution of galaxies. Amongst AGN at low radio-luminosities is the 'radio-quiet' quasar (RQQ) population, but how they contribute to the total radio emission is under debate, with previous studies arguing that it is predominantly through star formation. In this talk, SVW summarised the results of recent papers on RQQs, including the use of far-infrared data to disentangle the radio emission from the AGN and that from star formation. This provides evidence that black-hole accretion, instead, dominates the radio emission in RQQs. In addition, we find that this accretion-related emission is correlated with the optical luminosity of the quasar, whilst a weaker luminosity-dependence is evident for the radio emission connected with star formation. What remains unclear is the process by which this accretion-related emission is produced. Understanding this for RQQs will then allow us to investigate how this type of AGN influences its surroundings. Such studies have important implications for modelling AGN feedback, and for determining the accretion and star-formation histories of the Universe.

After the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered a significant brightening of the inner region of NGC 2617, we began a ~70 day photometric and spectroscopic monitoring campaign from the X-ray through near-infrared (NIR) wavelengths. We report that NGC 2617 went through a dramatic outburst, during which its X-ray flux increased by over an order of magnitude followed by an increase of its optical/ultraviolet (UV) continuum flux by almost an order of magnitude. NGC 2617, classified as a Seyfert 1.8 galaxy in 2003, is now a Seyfert 1 due to the appearance of broad optical emission lines and a continuum blue bump. Such "changing look Active Galactic Nuclei (AGN)" are rare and provide us with important insights about AGN physics. Based on the Hbeta line width and the radius-luminosity relation, we estimate the mass of central black hole to be (4 +/- 1) x 10^7 M_sun. When we cross-correlate the light curves, we find that the disk emission lags the X-rays, with the lag becoming longer as we move from the UV (2-3 days) to the NIR (6-9 days). Also, the NIR is more heavily temporally smoothed than the UV. This can largely be explained by a simple model of a thermally emitting thin disk around a black hole of the estimated mass that is illuminated by the observed, variable X-ray fluxes.