Spatially-resolved images of debris disks are necessary to determine disk morphological properties and the scattering phase function (SPF) which quantifies the brightness of scattered light as a function of phase angle. Current high-contrast imaging instruments have successfully resolved several dozens of debris disks around other stars, but few studies have investigated trends in the scattered-light, resolved population of debris disks in a uniform and consistent manner. We have combined Karhunen-Loeve Image Projection (KLIP) with radiative-transfer disk forward modeling in order to obtain the highest quality image reductions and constrain disk morphological properties of eight debris disks imaged by the Gemini Planet Imager at H-band with a consistent and uniformly-applied approach. In describing the scattering properties of our models, we assume a common SPF informed from solar system dust scattering measurements and apply it to all systems. We identify a diverse range of dust density properties among the sample, including critical radius, radial width, and vertical width. We also identify radially narrow and vertically extended disks that may have resulted from substellar companion perturbations, along with a tentative positive trend in disk eccentricity with relative disk width. We also find that using a common SPF can achieve reasonable model fits for disks that are axisymmetric and asymmetric when fitting models to each side of the disk independently, suggesting that scattering behavior from debris disks may be similar to Solar System dust.

Observations of exoplanet transits by small satellites have gained increasing attention for reducing detection biases. However, no unambiguous detection of an exoplanet has yet been demonstrated using optics with apertures smaller than 60 mm. Here, we investigated the detectability of exoplanet transits using the telescopic Optical Navigation Camera (ONC-T) onboard the Hayabusa2 spacecraft, which has an effective aperture of only 15 mm. We conducted transit observations of the hot Jupiters WASP-189 b and MASCARA-1 b, collecting data for ten and four events, respectively. The transit signal was detected with a signal-to-noise ratio (SNR) of 13 for WASP-189 b and 8 for MASCARA-1 b for each event. Stacking all events improved the SNR to 40 and 16, respectively. The transit mid-times of each event were measured with a precision of 6 minutes and were consistent with Transiting Exoplanet Survey Satellite (TESS) data to within 2 minutes. The planet-to-star radius ratio was determined with an absolute precision of 0.004 (6% relative) and agreed with TESS results to within 0.002 (3% relative). The recent ONC-T and TESS data enabled an update to the planetary ephemerides. We report a 4 sigma discrepancy between the updated orbital period of MASCARA-1 b and previously reported values. ONC-T sets a new record for the smallest-aperture instrument to detect an exoplanet transit from space, advancing the frontier of exoplanet science with miniature instrumentation. Our results suggest that optics as small as ONC-T may be capable of detecting transiting long-period Jupiters: a population that remains underrepresented in current surveys.

Mitigation of the threat from airbursting asteroids requires an understanding of the potential risk they pose for the ground. How asteroids release their kinetic energy in the atmosphere is not well understood due to the rarity of significant impacts. Ordinary chondrites, in particular L chondrites, represent a frequent type of Earth-impacting asteroids. Here, we present the first comprehensive, space-to-lab characterization of an L chondrite impact. Small asteroid 2023 CX1 was detected in space and predicted to impact over Normandy, France, on 13 February 2023. Observations from multiple independent sensors and reduction techniques revealed an unusual but potentially high-risk fragmentation behavior. The nearly spherical 650 $\pm$ 160 kg (72 $\pm$ 6 cm diameter) asteroid catastrophically fragmented around 28 km altitude, releasing 98% of its total energy in a concentrated region of the atmosphere. The resulting shockwave was spherical, not cylindrical, and released more energy closer to the ground. This type of fragmentation increases the risk of significant damage at ground level. These results warrant consideration for a planetary defense strategy for cases where a >3-4 MPa dynamic pressure is expected, including planning for evacuation of areas beneath anticipated disruption locations.

We present new VLTI/GRAVITY astrometry and updated orbit fits for the directly imaged companions YSES 1 b and HR 2562 B, substellar objects straddling the planet-brown dwarf boundary. Using high-precision astrometry, radial velocity (RV) data, and proper motions, we derive revised orbital parameters with orbitize! arXiv:1910.01756. For YSES 1 b, the inclusion of GRAVITY astrometry and a relative radial velocity measurement from arXiv:2409.16660 overcomes the traditional challenge of constraining eccentricities for distant companions, enabling the first orbit fit and yielding a constrained eccentricity of 0.44 (0.20). This represents the first full orbit fit for the system. Additionally, we calculate a median line-of-sight stellar obliquity of 12 (+11, -8) degrees, providing further insight into the system's dynamical architecture. For HR 2562 B, our analysis agrees with arXiv:2302.04893, confirming a low-eccentricity orbit (0.34 (0.20)) and an inclination of 87 (1) degrees. We find HR 2562 B's orbit to be nearly coplanar with the debris disk, with a mutual inclination of 3.7 (0.3) degrees. For both YSES 1 b and HR 2562 B the lower eccentricities favor an in situ formation scenario over extreme scattering or cloud fragmentation.

The future space-based gravitational wave observatory LISA is expected to detect massive black hole binaries (MBHBs) with high signal-to-noise ratios (SNRs), ranging up to thousands. Such high-precision observations require accurate modeling of the detector response. However, current derivations of the response function neglect the motion of the spacecraft during light travel time, omitting velocity-dependent terms of order $\beta = v/c \sim 10^{-4}$. In this work, we derive the velocity-dependent corrections to the gravitational wave response. We analyze the contribution of the velocity-terms for MBHBs in the mass range $[10^6,10^8]\:\mathrm{M}_{\odot}$ using a modified version of the state-of-the-art response simulator lisagwresponse. We find that corrections introduce residual SNRs up to $\sim 2$ for the loudest events and fractional differences up to $0.04\%$, compared to lisagwresponse. While small, these effects are comparable to current waveform modeling uncertainties and imprint distinctive sky-localization signatures, making them potentially relevant for parameter estimation of high-mass MBHBs and simulation of mock datasets.

Scaling relations between galactic parameters represent key pieces of evidence for investigating the processes of galaxy formation and evolution. In most studies, these relations have been obtained for large portions of the galaxies (i.e., on kpc scales), but it is also important to evaluate these relations in smaller scales. In this work, we used optical data cubes of a subsample of nearby galaxies of the DIVING 3D survey. These allowed us to analyze the scaling relations involving stellar velocity dispersion, stellar population age, and stellar population metallicity at the nuclear and circumnuclear regions of galaxies. We detected correlations between the stellar velocity dispersion and the age, metallicity, and total stellar mass. These correlations are independent of galaxy inclinations, considering all morphological types, nuclear activity, and the presence or absence of galactic bars. We detected, for the first time, a correlation between the stellar velocity dispersion and stellar metallicity in the nuclear regions of galaxies. It is found to be qualitatively consistent with the well-known stellar mass-metallicity relation. We also noted that barred galaxies tend to show younger and less metal-rich stellar populations than unbarred galaxies in the central regions, which may be a consequence of the bar triggering star formation in the nuclear regions of these objects. Some active galactic nuclei (AGNs) in our sample are positioned above the observed correlation between stellar velocity dispersion and stellar population age, suggesting that their nuclear stellar populations are younger than expected. This may be a consequence of positive AGN feedback, triggering star formation. Conversely, starburst galaxies do not show nuclear stellar populations at ages over one billion years.

We present a multiphase, resolved study of the galactic wind extending from the nearby starburst galaxy NGC 4666. For this we use VLT/MUSE observations from the GECKOS program and HI data from the WALLABY survey. We identify both ionised and HI gas in a biconical structure extending to at least $z\sim$8 kpc from the galaxy disk, with increasing velocity offsets above the midplane in both phases, consistent with a multiphase wind. The measured electron density, using [SII], differs significantly from standard expectations of galactic winds. We find electron density declines from the galaxy centre to $\sim2$ kpc, then rises again, remaining high ($\sim100-300$ cm$^{-3}$) out to $\sim$5 kpc. We find that HI dominates the mass loading. The total HI mass outflow rate (above $z~>2$ kpc) is between $5-13~M_{\odot}~\rm yr^{-1}$, accounting for uncertainties from disk-blurring and group interactions. The total ionised mass outflow rate (traced by H$\alpha$) is between $0.5~M_{\odot}~\rm yr^{-1}$ and $5~M_{\odot}~\rm yr^{-1}$, depending on $n_e(z)$ assumptions. From ALMA/ACA observations, we place an upper-limit on CO flux in the outflow which correlates to $\lesssim2.9~M_{\odot}~\rm yr^{-1}$. We also show that the entire outflow is not limited to the bicone, but a secondary starburst at the edge generates a more widespread outflow, which should be included in simulations. The cool gas in NGC 4666 wind has insufficient velocity to escape the halo of a galaxy of its mass, especially because most of the mass is present in the slower atomic phase. This strong biconical wind contributes to gas cycling around the galaxy.

We report on the observation and measurement of astrometry, photometry, morphology, and activity of the interstellar object 3I/ATLAS, also designated C/2025 N1 (ATLAS), with the NSF-DOE Vera C. Rubin Observatory. The third interstellar object, comet 3I/ATLAS, was first discovered on UT 2025 July 1. Serendipitously, the Rubin Observatory collected imaging in the area of the sky inhabited by the object during regular commissioning activities. We successfully recovered object detections from Rubin visits spanning UT 2025 June 21 (10 days before discovery) to UT 2025 July 7. Facilitated by Rubin's high resolution and large aperture, we report on the detection of cometary activity as early as June 21st, and observe it throughout. We measure the location and magnitude of the object on 37 Rubin images in r, i, and z bands, with typical precision of about 20 mas (100 mas, systematic) and about 10 mmag, respectively. We use these to derive improved orbit solutions, and to show there is no detectable photometric variability on hourly timescales. We derive a V-band absolute magnitude of H_V = (13.7 +/- 0.2) mag, and an equivalent effective nucleus radius of around (5.6 +/- 0.7) km. These data represent the earliest observations of this object by a large (8-meter class) telescope reported to date, and illustrate the type of measurements (and discoveries) Rubin's Legacy Survey of Space and Time (LSST) will begin to provide once operational later this year.

One of the first exoplanet hosts discovered thirty years ago, the star 55 Cnc has been constantly observed ever since. It is now known to host at least five planets with orbital periods ranging from 17 hours to 15 years. It is also one of the most extreme metal rich stars in the neighbourhood and it has a low-mass secondary star. In this article, we present data obtained at the Canada-France-Hawai'i Telescope with the SPIRou spectropolarimeter on both components of the 55 Cnc stellar system. We revisit the long-period radial-velocity signals of 55 Cnc A, with a focus on the role of the magnetic cycle, and propose the existence of a sixth planet candidate, whose period falls close to that of the magnetic cycle, or half of it. The other massive outer planet has a revised period of 13.15 years and a minimum mass of 3.8 MJup. Although some uncertainty remains on these outer planets, the characterization of the four inner planets is very robust through the combination of many different data sets, and all signals are consistent in the nIR and optical domains. In addition, the magnetic topology of the solar-type primary component of the system is observed by SPIRou at the minimum of its activity cycle, characterized by an amplitude ten times smaller than observed during its maximum in 2017. For the low-mass component 55 Cnc B, we report the discovery of two exoplanets in the system, with a period of 6.799+-0.0014 and 33.75+-0.04 days and a minimum mass of 3.5+-0.8 and 5.3+-1.4 MEarth, respectively. The secondary magnetic field is very weak and the current data set does not allow its precise characterization, setting an upper limit of 10 G. The system 55 Cnc stands out as the sixth binary system with planetary systems around both components, and the first one with non equal-mass stellar components.

We present an overview of the JWST GLIMPSE program, highlighting its survey design, primary science goals, gravitational lensing models, and first results. GLIMPSE provides ultra-deep JWST/NIRCam imaging across seven broadband filters (F090W, F115W, F200W, F277W, F356W, F444W) and two medium-band filters (F410M, F480M), with exposure times ranging from 20 to 40 hours per filter. This yields a 5$\sigma$ limiting magnitude of 30.9 AB (measured in a 0.2 arcsec diameter aperture). The field is supported by extensive ancillary data, including deep HST imaging from the Hubble Frontier Fields program, VLT/MUSE spectroscopy, and deep JWST/NIRSpec medium-resolution multi-object spectroscopy. Exploiting the strong gravitational lensing of the galaxy cluster Abell S1063, GLIMPSE probes intrinsic depths beyond 33 AB magnitudes and covers an effective source-plane area of approximately 4.4 arcmin$^2$ at $z \sim 6$. The program's central aim is to constrain the abundance of the faintest galaxies from $z \sim 6$ up to the highest redshifts, providing crucial benchmarks for galaxy formation models, which have so far been tested primarily on relatively bright systems. We present an initial sample of $\sim 540$ galaxy candidates identified at 6 < z < 16, with intrinsic UV magnitudes spanning $M_{\mathrm UV}$ = $-$20 to $-$12. This enables unprecedented constraints on the extreme faint end of the UV luminosity function at these epochs. In addition, GLIMPSE opens new windows for spatially resolved studies of star clusters in early galaxies and the detection and characterization of faint high-$z$ active galactic nuclei. This paper accompanies the first public data release, which includes reduced JWST and HST mosaics, photometric catalogs, and gravitational lensing models.

The distribution of close-in exoplanets is shaped by the interplay between atmospheric and dynamical processes. The Neptunian Desert, Ridge, and Savanna illustrate the sensitivity of these worlds to such processes, making them ideal to disentangle their roles. Determining how many Neptunes were brought close-in by early disk-driven migration (DDM; maintaining primordial spin-orbit alignment) or late high-eccentricity migration (HEM; generating large misalignments) is essential to understand how much atmosphere they lost. We propose a unified view of the Neptunian landscape to guide its exploration, speculating that the Ridge is a hot spot for evolutionary processes. Low-density Neptunes would mainly undergo DDM, getting fully eroded at shorter periods than the Ridge, while denser Neptunes would be brought to the Ridge and Desert by HEM. We embark on this exploration via ATREIDES, which relies on spectroscopy and photometry of 60 close-in Neptunes, their reduction with robust pipelines, and their interpretation through internal structure, atmospheric, and evolutionary models. We carried out a systematic RM census with VLT/ESPRESSO to measure the distribution of 3D spin-orbit angles, correlate its shape with system properties and thus relate the fraction of aligned-misaligned systems to DDM, HEM, and atmospheric erosion. Our first target, TOI-421c, lies in the Savanna with a neighboring sub-Neptune TOI-421b. We measured their 3D spin-orbit angles (Psib = 57+11-15 deg; Psic = 44.9+4.4-4.1 deg). Together with the eccentricity and possibly large mutual inclination of their orbits, this hints at a chaotic dynamical origin that could result from DDM followed by HEM. ATREIDES will provide the community with a wealth of constraints for formation and evolution models. We welcome collaborations that will contribute to pushing our understanding of the Neptunian landscape forward.

The Fermi Large Area Telescope (LAT) has revealed a mysterious extended excess of GeV gamma-ray emission around the Galactic Center, which can potentially be explained by unresolved emission from a population of pulsars, particularly millisecond pulsars (MSPs), in the Galactic bulge. We used the distributed volunteer computing system Einstein@Home to search the Fermi-LAT data for gamma-ray pulsations from sources in the inner Galaxy, to try to identify the brightest members of this putative population. We discovered four new pulsars, including one new MSP and one young pulsar whose angular separation to the Galactic Center of 0.93° is the smallest of any known gamma-ray pulsar. We demonstrate a phase-resolved difference imaging technique that allows the flux from this pulsar to be disentangled from the diffuse Galactic Center emission. No radio pulsations were detected from the four new pulsars in archival radio observations or during the MPIfR-MeerKAT Galactic Plane Survey. While the distances to these pulsars remain uncertain, we find that it is more likely that they are all foreground sources from the Galactic disk, rather than pulsars originating from the predicted bulge population. Nevertheless, our results are not incompatible with an MSP explanation for the GC excess, as only one or two members of this population would have been detectable in our searches.

We present a sample of six F200W and three F277W dropout sources identified as $1616$ candidates present mass-weighted ages around 30 Myr, and attenuations \mathrm{A(V)}<0.1 mag. Their average stellar mass is $\mathrm{M}_\bigstar\sim10^{7}\,\mathrm{M}_\odot$, implying a stellar-to-baryon mass fraction around 10% if the emissivity increases with redshift, or significantly higher otherwise. Three candidates present very blue UV spectral slopes ($\beta\sim-3$) compatible with Pop III young ($\lesssim10$ Myr) stars and/or high escape fractions of ionizing photons; the rest have $\beta\sim-2.5$ similar to $z=10-12$ samples.

The chemical composition of diffuse interstellar clouds is not fully established. They host an active chemistry despite their relatively low density and the ubiquitous presence of far-UV radiation. To further explore the chemical composition of diffuse clouds, we performed a spectral scan toward the bright radio source BL Lac in the Q band (from 32 to 50 GHz) using the Yebes 40m telescope. Yebes observations were performed interleaving Frequency Switching and Position Switching integrations toward BL Lac, using a spectral resolution of 38 kHz. The data have been reduced with the CLASS software. We achieved an unprecedented sensitivity on the continuum of 0.02 - 0.07 %, allowing the detection of very faint absorption features. We confirm previous detections of HCS+, C3H, C3H+, CH3CN and HC3N in diffuse clouds and report new detections of CCS, C4H, CH3CHO, H2CCO , HNCO and H2CS along the line of sight to BL Lac, with abundances relative to H2 from a few 10{-11} to a few 10{-10}. We compiled molecular detections toward diffuse clouds to obtain the chemical inventory of a typical diffuse interstellar cloud. The chemical inventory of diffuse interstellar clouds includes complex organic species with up to four heavy atoms. These species are efficiently formed in the diffuse interstellar gas and reach abundances similar to those measured in dense photodissociation regions, pointing to similar gas phase chemical processes.

All future lunar missions require a definition of the lunar reference system and a realization in the form of the lunar reference frame to ensure consistent products for positioning, navigation, cartography, and timing. This paper defines the origin, orientation, and scale of the Lunar Reference System (LRS), as well as provides numerical solutions for the first realization of the International Lunar Reference Frame (ILRF). ILRF is defined as the Principal Axis (PA) system, attached to the surface and co-rotating with the Moon, with its origin in the lunar center of mass (lunocenter). The ILRF realization is based on variance component estimation of the three lunar ephemeris solutions: INPOP21a, DE430, and EPM2021 for the series of the position of the lunar center of mass and rotation Euler angles -- precession, nutation, and proper rotation. The solution is valid starting with the period covered by Lunar Laser Ranging (LLR) data in 1970 and ending with extrapolated ILRF realizations in 2052 for future lunar missions. Results. The combined ILRF is characterized by the mean error of 17.6 cm for 2010-2030, where 15.3 cm comes from the origin and 8.6 cm from the orientation realization. The error in the realization of the origin is mainly caused by a poor geometry of the retroreflector network, resulting in a high correlation between the scale and the X component of the lunocenter in PA. The LLR post-fit residuals in ILRF are at the level of 2-3 cm in terms of the standard deviations of one-way ranges for best-performing LLR stations. The mean errors of the transformation between ILRF and other reference frame realizations in PA are at the level of 3 cm, whereas the mean transformation error to the DE421 Mean Earth frame equals 5 cm.

This work presents AstroPT, an autoregressive pretrained transformer developed with astronomical use-cases in mind. The AstroPT models presented here have been pretrained on 8.6 million $512 \times 512$ pixel $grz$-band galaxy postage stamp observations from the DESI Legacy Survey DR8. We train a selection of foundation models of increasing size from 1 million to 2.1 billion parameters, and find that AstroPT follows a similar saturating log-log scaling law to textual models. We also find that the models' performances on downstream tasks as measured by linear probing improves with model size up to the model parameter saturation point. We believe that collaborative community development paves the best route towards realising an open source `Large Observation Model' -- a model trained on data taken from the observational sciences at the scale seen in natural language processing. To this end, we release the source code, weights, and dataset for AstroPT under the MIT license, and invite potential collaborators to join us in collectively building and researching these models.

The multi-staged XENON program at INFN Laboratori Nazionali del Gran Sasso aims to detect dark matter with two-phase liquid xenon time projection chambers of increasing size and sensitivity. The XENONnT experiment is the latest detector in the program, planned to be an upgrade of its predecessor XENON1T. It features an active target of 5.9 tonnes of cryogenic liquid xenon (8.5 tonnes total mass in cryostat). The experiment is expected to extend the sensitivity to WIMP dark matter by more than an order of magnitude compared to XENON1T, thanks to the larger active mass and the significantly reduced background, improved by novel systems such as a radon removal plant and a neutron veto. This article describes the XENONnT experiment and its sub-systems in detail and reports on the detector performance during the first science run.

The Gaia Galactic survey mission is designed and optimized to obtain astrometry, photometry, and spectroscopy of nearly two billion stars in our Galaxy. Yet as an all-sky multi-epoch survey, Gaia also observes several million extragalactic objects down to a magnitude of G~21 mag. Due to the nature of the Gaia onboard selection algorithms, these are mostly point-source-like objects. Using data provided by the satellite, we have identified quasar and galaxy candidates via supervised machine learning methods, and estimate their redshifts using the low resolution BP/RP spectra. We further characterise the surface brightness profiles of host galaxies of quasars and of galaxies from pre-defined input lists. Here we give an overview of the processing of extragalactic objects, describe the data products in Gaia DR3, and analyse their properties. Two integrated tables contain the main results for a high completeness, but low purity (50-70%), set of 6.6 million candidate quasars and 4.8 million candidate galaxies. We provide queries that select purer sub-samples of these containing 1.9 million probable quasars and 2.9 million probable galaxies (both 95% purity). We also use high quality BP/RP spectra of 43 thousand high probability quasars over the redshift range 0.05-4.36 to construct a composite quasar spectrum spanning restframe wavelengths from 72-100 nm.