Papers for Monday, Jul 15 2024

We present two uniformly observed spectroscopic samples of Ly-alpha emitters (LAEs) (127 at z = 5.7 and 82 at z = 6.6), which we use to investigate the evolution of the LAE population at these redshifts. The observations cover a large field (44 sq. deg) in the North Ecliptic Pole (HEROES), as well as several smaller fields. We have a small number of exotic LAEs in the samples: double-peaked Ly-alpha profiles; very extended red wings; and one impressive lensed LAE cross. We also find three broad-line AGNs. We compare the Ly-alpha line width measurements at the two redshifts, finding that the lower-luminosity LAEs show a strong evolution of decreasing line width with increasing redshift, while the high-luminosity LAEs do not, with a transition luminosity of log L(Ly-alpha) = 43.25 erg s-1 . Thus, at z = 6.6, the high-luminosity LAEs may be producing large ionized bubbles themselves, or they may be residing in overdense galaxy sites that are producing such bubbles. In order to avoid losses in the red wing, the radius of the ionized bubble must be larger than 1 pMpc. The double-peaked LAEs also require transmission on the blue side. For the four at z = 6.6, we use models to estimate the proximity radii, Ra , where the ionizing flux of the galaxy is sufficient to make the surroundings have a low enough neutral fraction to pass the blue light. Since the required Ra are large, multiple ionizing sources in the vicinity may be needed.

The chameleon screening mechanism has been constrained many a time using dynamic and kinematic galaxy cluster observables. Current constraints are, however, insensitive to different mass components within galaxy clusters and have been mainly focused on a single mass density profile, the Navarro-Frenk-While mass density model. In this work, we extend the study of the Chameleon screening mechanism in galaxy clusters by considering a series of mass density models, namely: generalized-Navarro-Frenk-While, b-Navarro-Frenk-While, Burket, Isothermal and Einasto. The coupling strength ($\beta$) and asymptotic value of the chameleon field ($\phi_\infty$) are constrained by using kinematics analyses of simulated galaxy clusters, generated both assuming General Relativity and a strong chameleon scenario. By implementing a Bayesian analysis we comprehensively show that the biases introduced due to an incorrect assumption of the mass model are minimal. Similarly, we also demonstrate that a spurious detection of evidence for modifications to gravity is highly unlikely when utilizing the kinematics of galaxy clusters.

Large area observations of extragalactic deep fields with the James Webb Space Telescope (JWST) have provided a wealth of candidate low-mass L- and T-class brown dwarfs. The existence of these sources, which are at derived distances of hundreds of parsecs to several kiloparsecs from the Sun, has strong implications for the low-mass end of the stellar initial mass function, and the link between stars and planets at low metallicities. In this letter, we present a JWST/NIRSpec PRISM spectrum of brown dwarf JADES-GS-BD-9, confirming its photometric selection from observations taken as part of the JWST Advanced Deep Extragalactic Survey (JADES) program. Fits to this spectrum indicate that the brown dwarf has an effective temperature of 800-900K (T5 - T6) at a distance of $1.8 - 2.3$kpc from the Sun, with evidence of the source being at low metallicity ([M/H] $\leq -0.5$). Finally, because of the cadence of JADES NIRCam observations of this source, we additionally uncover a proper motion between the 2022 and 2023 centroids, and we measure a proper motion of $20 \pm 4$ mas yr$^{-1}$ (a transverse velocity of 214 km s$^{-1}$ at 2.25 kpc). At this predicted metallicity, distance, and transverse velocity, it is likely that this source belongs either to the edge of the Milky Way thick disk or the galactic halo. This spectral confirmation demonstrates the efficacy of photometric selection of these important sources across deep extragalactic JWST imaging.

Nulling interferometry is a powerful observing technique to reach exoplanets and circumstellar dust at separations too small for direct imaging with single-dish telescopes and too large for indirect methods. With near-future instrumentation, it bears the potential to detect young, hot planets near the snow lines of their host stars. A future space mission could detect and characterize a large number of rocky, habitable-zone planets around nearby stars at thermal-infrared wavelengths. The null self-calibration is a method aiming at modelling the statistical distribution of the nulled signal. It has proven to be more sensitive and accurate than average-based data reduction methods in nulling interferometry. This statistical approach opens the possibility of designing a GPU-based Python package to reduce the data from any of these instruments, by simply providing the data and a simulator of the instrument. GRIP is a toolbox to reduce nulling and interferometric data based on the statistical self-calibration method. In this article, we present the main features of GRIP as well as applications on real data.

Magnetars in quiescent states continue to emit hard X-rays with a power far exceeding the loss of rotational energy. It has recently been noted that this hard X-ray continuum may bear a direct signature of quantum electrodynamic (QED) effects in magnetic fields stronger than the Schwinger field ($B_{\rm Q} = 4.4\times 10^{13}$ G). When the current flowing into the magnetosphere is driven by narrow structures in the solid crust, the $e^\pm$ pair plasma supporting the current relaxes to a collisional and trans-relativistic state. The decay of a pair into two photons produces a broad, bremsstrahlung-like spectrum of hard X-rays, similar to that observed and extending up to $0.5-1$ MeV. The conversion of two gamma rays to a pair is further enhanced by a factor $\sim B/B_{\rm Q}$. Monte Carlo calculations of pair creation in a dipole magnetic field are presented. Non-local particle injection is found to be strong enough to suppress the high voltage that otherwise would accompany polar magnetic twist; the hard X-rays are mostly emitted away from the magnetic poles. Some of the pairs annihilate in an optically thin surface layer. The prototypical anomalous X-ray pulsar 1E 2259$+$586, which shows a hard X-ray continuum but relatively weak torque noise, slow spindown, and no radio emission, is a Rosetta Stone for understanding the magnetar circuit, consistent with the picture advanced here. For a $15-60$ keV luminosity as low as $10^{34}$ erg s$^{-1}$, the polar flux of sub-relativistic pairs produces an optical depth $3-30$ to electron cyclotron scattering in the $1-10$ keV band, reducing the net X-ray polarization.

Fifth forces are ubiquitous in modified gravity theories, and must be screened to evade stringent local tests. This can introduce unusual behaviour in galaxy phenomenology by affecting galaxies' components differently. Here we use the SDSS-IV-MaNGA dataset to search for a systematic excess of gas circular velocity over stellar circular velocity, expected in thin-shell-screened theories in the partially screened regime. Accounting for asymmetric drift and calibrating our model on screened subsamples, we find no significant evidence for a screened fifth force. We bound the fifth-force strength to $\Delta G/G_\text{N} < 0.1$ for all astrophysical ranges, strengthening to $\sim$0.01 at Compton wavelength of 3 Mpc for the Hu-Sawicki model, for instance. This implies a stringent constraint on scalar-tensor theories, for example $f_{\mathcal{R}0} \lesssim 10^{-8}$ in Hu-Sawicki $f(\mathcal{R})$ gravity.

Discovering new actively-accreting protoplanets is crucial to answering open questions about planet formation. However, identifying such planets at orbital distances where they are expected to be abundant is extremely challenging, both due to the technical requirements and large distances to star-forming regions. Here we use the kernel phase interferometry (KPI) technique to search for companions around the $\sim$6 and $\sim$8 Myr old Herbig Ae stars MWC 758 and MWC 480. KPI is a data analysis technique which is sensitive to moderate asymmetries, arising from eg. a circumstellar disk or companions with contrasts of up to 6-8 mags, at separations down to and even below the classical Rayleigh diffraction limit ($\sim 1.2\lambda / D$). Using the high spectral resolution K-band mode of the SCExAO/CHARIS integral field spectrograph, we search for both excess Br-$\gamma$ line emission and continuum emission from companions around MWC 480 and MWC 758. We are able to set limits on the presence of rapidly accreting protoplanets and brown dwarfs between 4 and 16 au, well interior to those of previous studies. In Br-$\gamma$, we set limits on excess line emission equivalent to accretion rates ranging from $10^{-5} M_{j}^{2}.yr^{-1}$ to $10^{-6}M_{j}^{2}.yr^{-1}$. Our achievable contrasts demonstrate that KPI using SCExAO/CHARIS is a promising technique to search for giant accreting protoplanets at smaller separations compared to conventional imaging.

The Solar eruptioN Integral Field Spectrograph (SNIFS) is a solar-gazing spectrograph scheduled to fly in the summer of 2025 on a NASA sounding rocket. Its goal is to view the solar chromosphere and transition region at a high cadence (1s) both spatially (0.5") and spectrally (33 mÅ) viewing wavelengths around Lyman Alpha (1216 Å), Si iii (1206 Å) and O v (1218 Å) to observe spicules, nanoflares, and possibly a solar flare. This time cadence will provide yet-unobserved detail about fast-changing features of the Sun. The instrument is comprised of a Gregorian-style reflecting telescope combined with a spectrograph via a specialized mirrorlet array that focuses the light from each spatial location in the image so that it may be spectrally dispersed without overlap from neighboring locations. This paper discusses the driving science, detailed instrument and subsystem design, and pre-integration testing of the SNIFS instrument.

We investigate the influence of bulk viscosity on late-time cosmic acceleration within an extended $f(Q, L_m)$ gravity framework, where the non-metricity $Q$ is non-minimally coupled with the matter Lagrangian $L_m$. Analyzing the function $f(Q, L_m) = \alpha Q + \beta L_m$, we derive exact solutions under non-relativistic matter domination. Using observational datasets ($H(z)$, Pantheon supernovae, and their combination), we constrain the model parameters $H_0$, $\alpha$, $\beta$, and $\zeta$. The deceleration parameter transitions from positive to negative values around redshifts $z_t \approx 0.80$ to $0.99 $, indicating current accelerated expansion. Moreover, the effective equation of state parameter, $\omega_{eff}$, resembles quintessence dark energy ($-1 < \omega_{eff} < -\frac{1}{3}$), with corresponding values from respective datasets. Finally, we use the $Om(z)$ diagnostic, which confirms that our model demonstrates quintessence-like behavior. Our findings underscore the significant role of bulk viscosity in understanding accelerated expansion in the universe within alternative gravity theories.

Using an established classification technique, we leverage standard observations and analyses to predict the progenitors of gamma-ray bursts (GRBs). This technique, utilizing support vector machine (SVM) statistics, provides a more nuanced prediction than the previous two-component Gaussian mixture in duration of the prompt gamma-ray emission. Based on further covariance testing from \textit{Fermi}-GBM, \textit{Swift}-BAT, and \textit{Swift}-XRT data, we find that our classification based only on prompt emission properties gives perspective on the recent evidence that mergers and collapsars exist in both long and short GRB populations.

We extend our previous work (Proukakis {\em et al.}, Phys.~Rev.~D~108,~083513 (2023)) on the dynamics of bosonic, non-relativistic and self-interacting dark matter that simultaneously contains both a ``fuzzy'' low-momentum component and one with higher momenta that may be well approximated as a collection of distinct particles and described by a corresponding phase-space distribution. Starting from the non-relativistic Schwinger-Keldysh action and working beyond leading-order in the Keldysh basis fields, encoding stochastic fluctuations of the slow modes and all fluctuations of the fast modes, we obtain stochastic self-consistently coupled Gross-Pitaevskii, collisional Boltzmann kinetic and Poisson equations. Our final set of equations, which feature various collisional (dissipative and scattering) contributions and two corresponding independent stochastic force terms, are consistent with generalized fluctuation-dissipation type relations in the limit of thermal equilibrium between the particles.

Front-end polarization modulation enables improved polarization measurement stability by modulating the targeted signal above the low-frequency $1/f$ drifts associated with atmospheric and instrumental instabilities and diminishes the impact of instrumental polarization. In this work, we present the design and characterization of a new 60-cm diameter Reflective Half-Wave Plate (RHWP) polarization modulator for the 90 GHz band telescope of the Cosmology Large Angular Scale Surveyor (CLASS) project. The RHWP consists of an array of parallel wires (diameter $50~\mathrm{\mu m}$, $175~\mathrm{\mu m}$ pitch) positioned $0.88~\mathrm{mm}$ from an aluminum mirror. In lab tests, it was confirmed that the wire resonance frequency ($f_\mathrm{res}$) profile is consistent with the target, $139~\mathrm{Hz}<f_\mathrm{res}<154~\mathrm{Hz}$ in the optically active region (diameter smaller than $150~\mathrm{mm}$), preventing the wire vibration during operation and reducing the RHWP deformation under the wire tension. The mirror tilt relative to the rotating axis was controlled to be $<15''$, corresponding to an increase in beam width due to beam smearing of $<0.6''$, negligible compared to the beam's full-width half-maximum of $36'$. The median and 16/84th percentile of the wire--mirror separation residual was $0.048^{+0.013}_{-0.014}~\mathrm{mm}$ in the optically active region, achieving a modulation efficiency $\epsilon=96.2_{+0.5}^{-0.4}\%$ with an estimated bandpass of 34 GHz. The angular velocity of the RHWP was maintained to an accuracy of within $0.005\%$ at the nominal rotation frequency ($2.5~\mathrm{Hz}$). The RHWP has been successfully integrated into the CLASS 90 GHz telescope and started taking data in June 2024, replacing the previous modulator that has been in operation since June 2018.

Understanding the surface temperature and interior structure of cold-to-temperate sub-Neptunes is critical for assessing their habitability, yet direct observations are challenging. In this study, we investigate the impact of water condensation on the atmospheric compositions of sub-Neptunes, focusing on the implications for JWST spectroscopic observations. By modeling the atmospheric photochemistry of two canonical sub-Neptunes, K2-18 b and LHS 1140 b, both with and without water condensation and with and without thick atmospheres, we demonstrate that water condensation can significantly affect the predicted atmospheric compositions. This effect is driven by oxygen depletion from the condensation of water vapor and primarily manifests as an increase in the C/O ratio within the photochemically active regions of the atmosphere. This change in composition particularly affects planets with thin H$_2$-dominated atmospheres, leading to a transition in dominant nitrogen and carbon carriers from N$_2$ and oxygen-rich species like CO/CO$_2$ towards heavier hydrocarbons and nitriles. While our models do not fully account for the loss mechanisms of these higher-order species, such molecules can go on to form more refractory molecules or hazes. Planets with thin H2-rich atmospheres undergoing significant water condensation are thus likely to exhibit very hazy atmospheres. The relatively flat JWST spectra observed for LHS 1140 b could be consistent with such a scenario, suggesting a shallow surface with extensive water condensation or a high atmospheric C/O ratio. Conversely, the JWST observations of K2-18 b are better aligned with a volatile-rich mini-Neptune with a thick atmosphere.

The distribution and evolution of the magnetic field at the solar poles through a solar cycle is an important parameter in understanding the solar dynamo. The accurate observations of the polar magnetic flux is very challenging from the ecliptic view, mainly due to (a) geometric foreshortening which limits the spatial resolution, and (b) the oblique view of predominantly vertical magnetic flux elements, which presents rather small line-of-sight component of the magnetic field towards the ecliptic. Due to these effects the polar magnetic flux is poorly measured. Depending upon the measurement technique, longitudinal versus full vector field measurement, where the latter is extremely sensitive to the SNR achieved and azimuth disamiguation problem, the polar magnetic flux measurements could be underestimated or overestimated. To estimate the extent of systematic errors in magetic flux measurements at the solar poles due to aforementioned projection effects we use MHD simulations of quiet sun network as a reference solar atmosphere. Using the numerical model of the solar atmosphere we simulate the observations from the ecliptic as well as from out-of-ecliptic vantage points, such as from a solar polar orbit at various heliographic latitudes. Using these simulated observations we make an assessment of the systematic errors in our measurements of the magnetic flux due to projection effects and the extent of under- or over estimation. We suggest that such effects could contribute to reported missing open magnetic flux in the heliosphere and that the multi-viewpoint observations from out-of-the-ecliptic plane together with innovative Compact Doppler Magnetographs provide the best bet for the future measurements.

In Sun and solar-type stars, there is a critical dynamo number for the operation of a large-scale dynamo, below which the dynamo ceases to operate. This region is known as the subcritical region. Previous studies showed the possibility of operating the solar-like large-scale (global) dynamo in the subcritical region without a small-scale dynamo. As in the solar convection zone, both large- and small-scale dynamos are expected to operate at the same time and location, we check the robustness of the previously identified subcritical dynamo branch in a numerical model in which both large- and small-scale dynamos are excited. For this, we use the {\sc Pencil Code} and set up an $\alpha\Omega$ dynamo model with uniform shear and helically forced turbulence. We have performed a few sets of simulations at different relative helicity to explore the generation of large-scale oscillatory fields in the presence of small-scale dynamo. We find that in some parameter regimes, the dynamo shows hysteresis behavior, i.e., two dynamo solutions are possible depending on the initial parameters used. A decaying solution when the dynamo was started with a weak field and a strong oscillatory solution if the dynamo was initialized with a strong field. Thus, the existence of the sub-critical branch of the large-scale dynamo in the presence of small-scale dynamo is established. However, the regime of hysteresis is quite narrow with respect to the case without the small-scale dynamo. Our work supports the possible existence of large-scale dynamo in the sub-critical regime of slowly rotating stars.

We present MeerKAT Fornax Survey HI observations of NGC 1427A, a blue irregular galaxy with a stellar mass of 2e+9 Msun located near the centre of the Fornax galaxy cluster. Thanks to the excellent resolution (1 to 6 kpc spatially, 1.4 km/s in velocity) and HI column density sensitivity (4e+19/cm^2 to 1e+18/cm^2 depending on resolution), our data deliver new insights on the long-debated interaction of this galaxy with the cluster environment. We confirm the presence of a broad, one-sided, starless HI tail stretching from the outer regions of the stellar body and pointing away from the cluster centre. We find the tail to have 50% more HI (4e+8 Msun) and to be 3 times longer (70 kpc) than in previous observations. In fact, we detect scattered HI clouds out to 300 kpc from the galaxy in the direction of the tail -- possibly the most ancient remnant of the passage of NGC 1427A through the intracluster medium of Fornax. Both the velocity gradient along the HI tail and the peculiar kinematics of HI in the outer region of the stellar body are consistent with the effect of ram pressure given the line-of-sight motion of the galaxy within the cluster. However, several properties cannot be explained solely by ram pressure and suggest an ongoing tidal interaction. This includes: the close match between dense HI and stars within the disturbed stellar body; the abundant kinematically-anomalous HI; and the inversion of the HI velocity gradient near the base of the HI tail. We rule out an interaction with the cluster tidal field, and conclude that NGC 1427A is the result of a high-speed galaxy encounter or of a merger started at least 300 Myr ago, where ram pressure shapes the distribution and kinematics of the HI in the perturbed outer stellar body and in the tidal tails.

MICADO is the first-light camera of the ESO ELT, allowing NIR imaging and long-slit spectroscopy assisted by adaptive optics. MICADO is now entering its construction phase, and the software for data reduction is reaching an adequate maturity level. The PSF Reconstruction (PSF-R) of MICADO is a software tool for the blind derivation of the PSF, only using adaptive optics telemetry data. An update of the status of the PSF-R service is provided here. The PSF-R prototype has been tested on ERIS@VLT data in order to check the reconstruction of on- and off-axis PSFs. The on-axis PSF-R is accurate at a few percent level on Strehl, FWHM, Encircled Energy, and half light radius, while for the off-axis case the match is within 10-15 percent at a distance of half isoplanatic angle. The first version of the workflow for the PSF-R pipeline has been developed and verified using the latest release of the ESO data processing system. A set of simulations has been implemented on the morphological analysis of distant galaxies, showing that the accuracy of the PSF-R matches the goals needed to study their morphology. In summary, the PSF-R team is on the right track towards the ELT first light.

We report results from numerical simulations assessing astrometry measurements with the Multiconjugate Adaptive Optics Relay for ELT Observations (MORFEO) instrument on the Extremely Large Telescope (ELT). Using the Advanced Exposure Time Calculator (AETC), we evaluate MORFEO astrometric accuracy in moderately crowded fields. Our simulations account for spatially variable Point Spread Function (PSF), geometric distortion, and rotation-dependent variations. We computed focal plane coordinates using observed stellar distribution and computed population synthesis with the SPISEA tool, generating stellar magnitude distributions for MICADO filters at selected metallicities and stellar ages. Our analysis shows that MORFEO can achieve high-precision astrometry in the galaxy neighborhood (within $\mu < 24$ mag) by minimizing PSF enlargement and optimizing calibration strategies. These results inform future observational campaigns and contribute to the development of astrometric science cases for the ELT.

The FlyEye design makes its debut in the ESA's NEOSTEL developed by OHB-Italia. This pioneering FlyEye telescope integrates a monolithic 1-meter class primary mirror feeding 16 CCD cameras for discovering Near-Earth Object (NEO) and any class of transient phenomena. OHB-Italia is the prime contractor, receiving extended support from the Italian National Institute for Astrophysics (INAF) in the ESA's NEOSTED program's integration and testing. The FlyEye distinctive design splits the Field of View into 16 channels, creating a unique multi-telescope system with a panoramic 44 square degree Field of View and a seeing-size pixel-scale, enabling NEOs detection down to apparent magnitudes 21.5 insisting on a 1m diameter spherical mirror. The scientific products of a similar FlyEye telescope can complement facilities such as Vera Rubin (former LSST) and ZTF. The FlyEye has the ability to survey two-thirds of the visible sky about three times per night can revolutionize time-domain astronomy, enabling comprehensive studies of transient phenomena, placing FlyEye in a new era of exploration of the dynamic universe. Efforts to develop automated calibration and testing procedures are keys to realizing this transformative potential.

We present a silicon drift detector (SDD) system for the spectroscopy focusing array (SFA) of the enhanced X-ray timing and polarimetry (eXTP) mission. The SFA focuses on fast timing (time resolution below 10 {\mu}s) and good spectroscopy capabilities (energy resolution better than 180 eV @ 6 keV). The sensor, consisting of 19 hexagonally shaped pixels with a total sensitive area of ${5.05}\, cm^{2}$, is connected to three high time resolution spectroscopy (HTRS) ASICs, allowing a fast readout of the detector signals. The detector works in a Charge- Sensitive Amplifier configuration. We assembled a prototype detector module and present here its mechanical design, describe the used sensor, and report about its performance.

We report results from 8 hours of JWST/MIRI LRS spectroscopic monitoring directly followed by 7 hours of JWST/NIRSpec prism spectroscopic monitoring of the benchmark binary brown dwarf WISE 1049AB, the closest, brightest brown dwarfs known. We find water, methane, and CO absorption features in both components, including the 3.3 $\mu$m methane absorption feature and a tentative detection of small grain ($<$ 1$\mu$m) silicate absorption at $>$8.5 $\mu$m in WISE 1049A. Both components vary significantly ($>$1$\%$), with WISE 1049B displaying larger variations than WISE 1049A. Using K-means clustering, we find three main transition points in wavelength for both components of the binary: 1) change in behavior at $\sim$2.3 $\mu$m coincident with a CO absorption bandhead, 2) change in behavior at 4.2 $\mu$m, close to the CO fundamental band at $\lambda >$ 4.4 $\mu$m, and 3) change in behavior at 8.3-8.5 $\mu$m, potentially corresponding to silicate absorption. We interpret the lightcurves observed with both NIRSpec and MIRI as likely stemming from 1) a deep pressure level driving the double-peaked variability seen in WISE 1049B at wavelengths $<$2.3 $\mu$m and $>$8.5 $\mu$m, 2) an intermediate pressure level shaping the lightcurve morphology between 2.3 and 4.2 $\mu$m, and 3) a higher-altitude pressure level producing single-peaked and plateaued lightcurve behavior between 4.2 and 8.5 $\mu$m.

Context: Helioseismic inversions rely largely on sensitivity kernels, 3-D spatial functions describing how do the changes in the solar interior translate into the changes of helioseismic observables. These sensitivity kernels come in most cases from the forward modelling utilising state-of-the-art solar models. Aims: We aim to test the sensitivity kernels by comparing their volume integrals with measured values from helioseismic travel times. Methods: By manipulating the tracking rate, we mimic the additional zonal velocity present in the Dopplergram datacubes. These datacubes are then processed by a standard travel-time measurements pipeline. We investigate the dependence of the east--west travel time averaged over a box around the disc centre on the implanted tracking velocity. The slope of this dependence is directly proportional to the total volume integral of the sensitivity kernel corresponding to the used travel-time geometry. Results: We find a very good to acceptable agreement between measurements and models for travel times with ridge filtering applied. The sought dependence indeed resembles a linear function and the slope of it agrees with the expected volume integral from the forward-modelled sensitivity kernel. The agreement is less optimistic for the phase-speed filtered datacubes. A disagreement is particularly large for smallest phase speeds (filters td1--td4), for the larger phase speeds, our result indicate that the measured kernel integrals are systematically larger than expected from the forward modelling. We admit that for large phase-speeds and higher radial modes our testing procedure does not have to be appropriate.

The highest scientific return, for adaptive optics (AO) observations, is achieved with a reliable reconstruction of the PSF. This is especially true for MICADO@ELT. In this presentation, we will focus on extending the MICADO PSF reconstruction (PSF-R) method to the off-axis case. Specifically, a novel approach based on temporal-based tomography of AO telemetry data has been recently implemented. Results from the PSF-R of both simulated and real data show that, at half isoplanatic angle distances, a precision of about 10-15% is achievable in both Strehl ratio and full-width at half maximum, paving the way to extend the MICADO PSF-R tool also to the multi-conjugated AO case.

The discovery of pulsations in (at least) six ultraluminous X-ray sources (ULXs) has shown that neutron stars can accrete at (highly) super-Eddington rates, challenging the standard accretion theories. M51 ULX-7, with a spin signal of $P\simeq2.8$ s, is the pulsating ULX (PULX) with the shortest known orbital period ($P_\mathrm{orb}\simeq2$ d) and has been observed multiple times by XMM-Newton, Chandra, and NuSTAR. We report on the timing and spectral analyses of three XMM-Newton observations of M51 ULX-7 performed between the end of 2021 and the beginning of 2022, together with a timing re-analysis of XMM-Newton, Chandra, and NuSTAR archival observations. We investigated the spin signal by applying accelerated search techniques and studied the power spectrum through the fast Fourier transform, looking for (a)periodic variability in the source flux. We analysed the energy spectra of the 2021-2022 observations and compared them to the older ones. We report the discovery of a recurrent, significant ($>$3$\sigma$) broad complex at mHz frequencies in the power spectra of M51 ULX-7. We did not detect the spin signal, setting a 3$\sigma$ upper limit on the pulsed fraction of $\lesssim10\%$ for the single observation. The complex is significantly detected also in five Chandra observations performed in 2012. M51 ULX-7 represents the second PULX for which we have a significant detection of mHz-QPOs at super-Eddington luminosities. These findings suggest that one should avoid using the observed QPO frequency to infer the mass of the accretor in a ULX. The absence of spin pulsations when the broad complex is detected suggests that the mechanism responsible for the aperiodic modulation also dampens the spin signal's pulsed fraction. If true, this represents an additional obstacle in the detection of new PULXs, suggesting an even larger occurrence of PULXs among ULXs.

Modified theories of gravity predict deviations from General Relativity (GR) in the propagation of gravitational waves (GW) across cosmological distances. A key prediction is that the GW luminosity distance will vary with redshift, differing from the electromagnetic (EM) luminosity distance due to varying effective Planck mass. We introduce a model-independent, data-driven approach to explore these deviations using multi-messenger observations of dark standard sirens (Binary Black Holes, BBH). By combining GW luminosity distance measurements from dark sirens with Baryon Acoustic Oscillation (BAO) measurements, BBH redshifts inferred from cross-correlation with spectroscopic or photometric galaxy surveys, and sound horizon measurements from the Cosmic Microwave Background (CMB), we can make a data-driven test of GR (jointly with the Hubble constant) as a function of redshift. Using the multi-messenger technique with the spectroscopic DESI galaxy survey, we achieve precise measurements of deviations in the effective Planck mass variation with redshift. For the Cosmic Explorer and Einstein Telescope (CEET), the best precision is approximately 3.6\%, and for LIGO-Virgo-KAGRA (LVK), it is 7.4\% at a redshift of $\rm{z = 0.425}$. Additionally, we can measure the Hubble constant with a precision of about 1.1\% from CEET and 7\% from LVK over five years of observation with a 75\% duty cycle. We also explore the potential of cross-correlation with photometric galaxy surveys from the Rubin Observatory, extending measurements up to a redshift of $\rm{z \sim 2.5}$. This approach can reveal potential deviations from models affecting GW propagation using numerous dark standard sirens in synergy with DESI and the Rubin Observatory.

The fascination with the blazar OJ 287 stems not only from its status as a prominent candidate for a close supermassive black hole (SMBH) binary, but also because of the thermal bremsstrahlung origin proposed for its giant optical outbursts. These outbursts arrive as pairs, quasi-periodically every $\sim 12$ years, based on the unique 130-year-long, well-sampled optical light curve available for this blazar. For its three well-known, large quasi-periodic optical outbursts (QPOOs), observed in 1983, 2007 and 2015, optical photo-polarimetric monitoring has been reported in the literature. For these initially radio-undetected QPOOs, widely acclaimed as `bremsstrahlung flares', we have scrutinised the available measurements of optical polarisation and spectral index during the rising phase. Several inconsistencies of these data with the optical bremsstrahlung interpretation are noted, which point towards a synchrotron-dominated alternative interpretation for all these prominent QPOOs, just as for the optical emission observed between the outbursts. Possible reasons for the radio non-detection of the QPOOs during the initial stage are outlined.

When constructing galaxy mock catalogs based on suits of dark matter halo catalogs generated with approximated, calibrated or machine-learning approaches, the assignment of intrinsic properties for such tracers is a step of paramount relevance, given that these can shape the abundance of mock galaxy cluster and the spatial distribution of mock galaxies. We explore the possibility to assign properties of dark matter halos within the context of calibrated/learning approaches, explicitly using clustering information. The goal is to retrieve the correct signal of primary and secondary large-effective bias as a function of properties reconstructed solely based on phase-space properties of the halo distribution and dark matter density field. The algorithm reconstructs a set halo properties (such as virial mass, maximum circular velocity, concentration and spin) constraint to reproduce both primary and secondary (or assembly) bias. The key ingredients of the algorithm are the implementation of individually-assigned large-scale effective bias, a multi-scale approach to account for halo exclusion and a hierarchical assignment of halo properties. The method facilitates the assignment of halo properties aiming at replicating the large-scale effective bias, both primary and secondary. This improves over over previous methods found in the literature, especially at the high mass population. We have designed an strategy to reconstruct the main properties of dark matter halos obtained by calibrated/learning algorithms, in a way that the expected one and two-point statistics (on large scales) replicates the signal from detailed N-body simulations. We encourage the application of this strategy (or the implementation of our algorithm) for the generation of mock catalogs of dark matter halos based on approximated methods.

With the James Webb Space Telescope (JWST) offering higher resolution data in space-based transmission spectroscopy, understanding the capabilities of our current atmospheric retrieval pipelines is essential. These new data cover wider wavelength ranges and at much higher spectral resolution than previous instruments have been able to offer. Therefore, it is often appealing to bin spectra to fewer points, better constrained in their transit depth, before using them as inputs for atmospheric retrievals. However, little quantitative analysis of the trade-off between spectral resolution and signal-to-noise ratio has been conducted thus far. As such, we produce a simulation replicating the observations of WASP-39b by the NIRSpec PRISM instrument on board JWST and assess the accuracy and consistency of retrievals while varying resolution and the average photometric error. While this probes a specific case we also plot `binning paths' in the resulting sensitivity maps to demonstrate the best attainable atmospheric parameter estimations starting from the position of the real JWST Early Release Science observation. We repeat this analysis on three different simulation setups where each includes an opaque cloud layer at a different height in the atmosphere. We find that a much greater resolution is needed in the case of a high cloud deck since features are already heavily muted by the presence of the clouds. In the other two cases, there are large `safe zones' in the parameter space. If these maps can be generalised, binning paths could inform future observations on how to achieve the most accurate retrieval results.

The Galactic disk is one of the main components of the Milky Way, which contributes most of the luminosity. Its structure is essential for understanding the formation and evolution of the Milky Way. Using 174,443 K-type dwarf stars observed by both LAMOST and Gaia DR3, we study the disk density profile in the local volume within 1,200 pc. In the azimuthal dimension, we find strong asymmetric signal of the thin disk. The surface density and the scale height of the southern disk significantly change versus the azimuthal angle at the same galactocentric distance $R$. Meanwhile, in the vertical dimension, the scale height of the northern disk has quite different trend than that of the southern one. The scale height of the southern disk shows a decreasing trend with $\phi\sim-2.5^\circ$, and change to an increasing one with $\phi\sim5.0^°$. Meanwhile, the scale height of the northern disk has a consistently smaller increase. Finally, we divide the entire sample into three subsamples based on metallicity and all three subsamples show significant non-axisymmetric and north-south asymmetric signals in the Galactic disk. Furthermore, we find that the scale height of the metal-poor ([Fe/H] $<$ -0.4 dex) subsample in the northern disk is greater than that of the metal-rich ([Fe/H] $>$ -0.1 dex) subsample. However, in the southern disk, the scale height exhibits varying relationships across different metallicity slices.

Rates of stellar tidal disruption events (TDEs) around supermassive black holes (SMBHs) have been extensively calculated using the loss cone theory, while theoretical work on TDE rates around intermediate-mass black holes (IMBHs) has been lacking. In this work, we aim to accurately calculate the IMBH TDE rates based on their black hole masses and the stellar profiles of their host galaxies obtained from the latest observations. We find that IMBH TDEs from the center of small galaxies have an overall rate comparable to SMBH TDEs, while off-nuclei IMBH TDEs from globular clusters have a much lower rate. Very interestingly, we show that the rate of IMBH TDE per galaxy generally increases with the black hole mass, which is opposite to the trend seen in SMBH TDEs. Furthermore, we report that IMBH TDEs typically occur in the pinhole regime, which means that deeply plunging events are more likely for IMBH TDEs compared to SMBH TDEs. We also calculate the volumetric TDE rates for IMBH and SMBH TDEs and compare with observed rates.

Linking atmospheric measurements to the bulk planetary composition and ultimately the planetary origin is a key objective in planetary science. In this work we identify the cases in which the atmospheric composition represents the bulk composition. We simulate the evolution of giant planets considering a wide range of planetary masses ($0.3-2~M_\mathrm{J}$), initial entropies ($8-11~k_\mathrm{B} m_\mathrm{u}^{-1}$), and primordial heavy-element profiles. We find that convective mixing is most efficient at early times (ages $\lesssim 10^7$ years) and that primordial composition gradients can be eroded. In several cases, however, the atmospheric composition can differ widely from the planetary bulk composition, with the exact outcome depending on the details. We show that the efficiency of convection is primarily controlled by the underlying entropy profile: for low primordial entropies of $8-9~k_\mathrm{B} m_\mathrm{u}^{-1}$ convective mixing can be inhibited and composition gradients can persist over billions of years. The scaling of mixing efficiency with mass is governed by the primordial entropy. For the same primordial entropy, low-mass planets mix more efficiently than high mass planets. If the primordial internal entropy would increase with mass, however, this trend could reverse. We also present a new analytical model that predicts convective mixing under the existence of composition (and entropy) gradients. Our results emphasize the complexity in the interpretation of atmospheric abundance measurements and show the great need to better understand the planetary formation process as it plays a key role in determining the planetary evolution and final structure.

The dust grain size distribution (GSD) likely varies significantly across different star-forming environments in the Universe, but the overall impact of this variation on star formation remains unclear. This ambiguity arises because the GSD interacts non-linearly with processes like heating/cooling, radiation, and chemistry, which have competing effects and different environmental dependencies. In this study, we investigate the effects of GSD variation on the thermochemistry and evolution of giant molecular clouds (GMCs). To achieve this, we conducted radiation-dust-magnetohydrodynamic simulations spanning a range of cloud masses and grain sizes, which explicitly incorporate the dynamics of dust grains within the full-physics framework of the STARFORGE project. We find that differences in grain size significantly alter the thermochemistry of GMCs. Specifically, we show that the leading-order effect is that larger grains, under fixed dust mass and dust-to-gas ratio conditions, result in lower dust opacities. This reduced opacity permits ISRF photons to penetrate more deeply and allows internal radiation field photons to permeate more extensively into the cloud, resulting in rapid gas heating and the inhibition of star formation. We find that star formation efficiency is highly sensitive to grain size, with an order of magnitude reduction in efficiency when grain size increases from 0.1 $\rm\mu m$ to 10 $\rm\mu m$. Additionally, we note that warmer gas suppresses the formation of low-mass stars. Moreover, as a consequence of the decreased opacities, we observe a greater proportion of gas residing in diffuse ionized structures.

The Synthetic Corona Outflow Model (SynCOM), an empirical model, simulates the solar corona's dynamics to match high-resolution observations, providing a useful resource for testing velocity measurement algorithms. SynCOM generates synthetic images depicting radial variability in polarized brightness and includes stochastic elements for plasma outflows and instrumental noise. It employs a predefined flow velocity probability distribution and an adjustable signal-to-noise ratio to evaluate different data analysis methods for coronal flows. By adjusting parameters to match specific coronal and instrumental conditions, SynCOM offers a platform to assess these methods for determining coronal velocity and acceleration. Validating these measurements would help to understand solar wind origins and support missions such as the Polarimeter to Unify the Corona and Heliosphere (PUNCH). In this study, we demonstrate how SynCOM can be employed to assess the precision and performance of two different flow tracking methods. By providing a ground-truth based on observational data, we highlight the importance of SynCOM in confirming observational standards for detecting coronal flows.

We examine Saturn's non-auroral (dayglow) emissions at Lyman-$\alpha$ observed by the {Cassini/UVIS} instrument from 2004 until 2016, to constrain meridional and seasonal trends in the upper atmosphere. We separate viewing geometry effects from trends driven by atmospheric properties, by applying a multi-variate regression to the observed emissions. The Lyman-$\alpha$ dayglow brightnesses depend on the incident solar flux, solar incidence angle, emission angle, and observed latitude. The emissions across latitudes and seasons show a strong dependence with solar incidence angle, typical of resonantly scattered solar flux and consistent with no significant internal source. We observe a bulge in Ly-$\alpha$ brightness that shifts with the summer season from the southern to the northern hemisphere. We estimate atomic hydrogen optical depths above the methane homopause level for dayside disk observations (2004-2016) by comparing observed Lyman-$\alpha$ emissions to a radiative transfer model. We model emissions from resonantly scattered solar flux and a smaller but significant contribution by scattered photons from the interplanetary hydrogen (IPH) background. During northern summer, inferred hydrogen optical depths steeply decrease with latitude towards the winter hemisphere from a northern hemisphere bulge, as predicted by a 2D seasonal photochemical model. The southern hemisphere mirrors this trend during its summer. However, inferred optical depths show substantially more temporal variation between 2004 and 2016 than predicted by the photochemical model.

Computation of a grid of self consistent 1D model atmospheres of cool stars, sub-stellar objects and exoplanets in the effective temperature range 300K to 3000K, including cloud formation, chemical non-equilibrium effects, and stellar irradiation. The models are called MSG, because they are based on an iterative coupling between three well tested codes, the MARCS stellar atmosphere code, the StaticWeather cloud formation code and the GGchem chemical equilibrium code. It includes up-to-date molecular and atomic opacities, cloud formation and advanced chemical equilibrium calculations, and involves new numerical methods at low temperatures to allow robust convergence. The coupling between the MARCS radiative transfer and GGchem chemical equilibrium computations has made it possibly effectively to reach convergence based on electron pressure for the warmer models and gas pressure for the cooler models, enabling self-consistent modelling of stellar, sub-stellar and exoplanetary objects in a very wide range of effective temperatures. Here we describe the basic details of the models, with illustrative examples of cloudy and irradiated models as well as models based on non-equilibrium chemistry. The qualitative changes in the relative abundances of TiO, H2O, CH4, NH3, and other molecules in our models follow the observationally defined M, L, T (and Y) sequences, but reveal more complex and depth dependent abundance changes, and therefore a spectral classification depending on more parameters. The self consistent coupling to Static-Weather cloud computations, allows detailed comparison between nucleation and observed relative dimming of different spectral bands, with advanced applications for new identification methods of potential exoplanetary biology.

The standard cosmological model, which assumes statistical isotropy and parity invariance, predicts the absence of correlations between even-parity and odd-parity observables of the cosmic microwave background (CMB). Contrary to these predictions, large-angle CMB temperature anomalies generically involve correlations between even-$\ell$ and odd-$\ell$ angular power spectrum $C_\ell$, while recent analyses of CMB polarization have revealed non-zero equal-$\ell$ $EB$ correlations. These findings challenge the conventional understanding, suggesting deviations from statistical isotropy, violations of parity, or both. Cosmic topology, which involves changing only the boundary conditions of space relative to standard cosmology, offers a compelling framework to potentially account for such parity-violating observations. Topology inherently breaks statistical isotropy, and can also break homogeneity and parity, providing a natural paradigm for explaining observations of parity-breaking observables without the need to add parity violation to the underlying microphysics. Our investigation delves into the harmonic space implications of topology for CMB correlations, using as an illustrative example $EB$ correlations generated by tensor perturbations under both parity-preserving and parity-violating scenarios. Consequently, these findings not only challenge the foundational assumptions of the standard cosmological model but also open new avenues for exploring the topological structure of the Universe through CMB observations.

The initial conditions on the anisotropies of the stochastic gravitational-wave background of cosmological origin (CGWB) largely depend on the mechanism that generates the gravitational waves. Since the CGWB is expected to be non-thermal, the computation of the initial conditions could be more challenging w.r.t. the Cosmic Microwave Background (CMB), whose interactions with other particles in the early Universe lead to a blackbody spectrum. In this paper, we show that the initial conditions for the cosmological background generated by quantum fluctuations of the metric during inflation deviate from adiabaticity. These primordial gravitational waves are indeed generated by quantum fluctuations of two independent degrees of freedom (the two polarization states of the gravitons). Furthermore, the CGWB plays a negligible role in the Einstein's equations, because its energy density is subdominant w.r.t. ordinary matter. Therefore, the only possible way to compute the initial conditions for inflationary gravitons is to perturb the energy-momentum tensor of the gravitational field defined in terms of the small-scale tensor perturbation of the metric. This new and self-consistent approach shows that a large, non-adiabatic initial condition is present even during the single-field inflation. Such a contribution enhances the total angular power spectrum of the CGWB compared to the standard adiabatic case, increasing also the sensitivity of the anisotropies to the presence of relativistic and decoupled particles in the early Universe. In this work we have also proved that our findings are quite general and apply to both single-field inflation and other scenarios in which the CGWB is generated by the quantum fluctuations of the metric, like the curvaton.

Context. Since the first publication of the Gaia catalogue a new view of our Galaxy has arrived. Its astrometric and photometric information has improved the precision of the physical parameters of open star clusters obtained from them. Aims. Using the Gaia DR3 catalogue, we aim to find physical stellar members including faint stars for 370 Galactic open clusters located within 1 kpc. We also estimate the age, metallicity, distance modulus and extinction of these clusters. Methods. We employ the HDBSCAN algorithm on both astrometric and photometric data to identify members in the open clusters. Subsequently, we refine the samples by eliminating outliers through the application of the Mahalanobis metric utilizing the chi-square distribution at a confidence level of 95%. Furthermore, we characterize the stellar parameters with the PARSEC isochrones. Results. We obtain reliable star members for 370 open clusters with an average parallax error of 0.16 mas. We identify about 40% more stars in these clusters compared to previous work using the Gaia DR2 catalogue, including faint stars as new members with G > 17. Before the clustering application we correct the parallax zero-point bias to avoid spatial distribution stretching that may affect clustering results. Our membership lists include merging stars identified by HDBSCAN with astrometry and photometry. We note that the use of photometry in clustering can recover up to 10% more stars in the fainter limit than clustering based on astrometry only, this combined with the selection of stars filtering them out by quality cuts significantly reduces the number of stars with huge parallax error. After clustering, we estimate age, Z, and AV from the photometry of the membership lists.

We present the results of the observational study of the blazar S5 0716+716 in the optical bands B, V, R, and I between March 2019 and August 2023 to investigate its variability on diverse time-scales. The blazar was followed up by the T60 robotic telescope in Turkey for 416 nights to obtain long-term variability during this period. In order to search for intraday variability of the object, we have carried out 21 nights of observations with the T100 telescope for at least 1 hour. The blazar showed a ~2.47 mag variation in the optical R-band during our monitoring period, the brightest state on 18.01.2020 (MJD 58866) as R=12.109$\pm$0.011 and the faintest state on 23.03.2019 (MJD 58565) as R=14.580$\pm$0.013. We employed the nested ANOVA test and the power enhanced F-test to quantify intraday variability which showed that the blazar was significantly variable in the R-band on 12 out of 21 nights. Correlation analysis of the light curves shows that the emission in the BVRI optical bands was strongly correlated both in the short and long term without any time lag. The blazar has likely quasi-periods of 186$\pm$30, 532$\pm$76 days in the optical R-band light curve according to the WWZ, and the LS periodogram. The IDV and LTV features are discussed within the frame of prospective scenarios.

We present a flow-based generative approach to emulate grids of stellar evolutionary models. By interpreting the input parameters and output properties of these models as multi-dimensional probability distributions, we train conditional normalizing flows to learn and predict the complex relationships between grid inputs and outputs in the form of conditional joint distributions. Leveraging the expressive power and versatility of these flows, we showcase their ability to emulate a variety of evolutionary tracks and isochrones across a continuous range of input parameters. In addition, we describe a simple Bayesian approach for estimating stellar parameters using these flows and demonstrate its application to asteroseismic datasets of red giants observed by the Kepler mission. By applying this approach to red giants in open clusters NGC 6791 and NGC 6819, we illustrate how large age uncertainties can arise when fitting only to global asteroseismic and spectroscopic parameters without prior information on initial helium abundances and mixing length parameter values. We also conduct inference using the flow at a large scale by determining revised estimates of masses and radii for 15,388 field red giants. These estimates show improved agreement with results from existing grid-based modelling, reveal distinct population-level features in the red clump, and suggest that the masses of Kepler red giants previously determined using the corrected asteroseismic scaling relations have been overestimated by 5-10%.

This manuscript revisits the phenomenological emergent dark energy model (PEDE) by confronting it with recent cosmological data from early and late times. In particular we analyze PEDE model by using the baryon acoustic oscillation (BAO) measurements coming from both Dark Energy Spectroscopy Instrument (DESI) data release 1 and Sloan Digital Sky Survey (SDSS). Additionally, the measurements from cosmic chronometers, supernovae type Ia (Pantheon+), quasars, hydrogen II galaxies and cosmic background radiation distance priors are considered. By performing a Bayesian analysis based on Monte Carlo Markov Chain, we find consistent results on the constraints when SDSS and DESI are considered. However, we find higher values on the Hubble constant than Supernova $H_0$ for the Equation of State (SH0ES) does although it is still in agreement, within $1\sigma$ confidence level, when BAO measurements are added. Furthermore, we estimate the age of the Universe younger $\sim3\%$ than the one predicted by the standard cosmology. Additionally, we report values of $q_0 = -0.771^{+0.007}_{-0.007}$, $z_T = 0.764^{+0.011}_{-0.011}$ for the deceleration parameter today and the deceleration-acceleration transition redshift, respectively.

In this report, we summarize the activities of the International Astronomical Consortium for High Energy Calibration (IACHEC) from the 15th IACHEC Workshop in Pelham, Germany. Sixty scientists directly involved in the calibration of operational and future high-energy missions gathered for 3.5 days to discuss the status of the cross-calibration between the current international complement of X-ray observatories and the possibilities to improve it. This summary consists of reports from the Working Groups with topics ranging across the identification and characterization of standard calibration sources, multi-observatory cross-calibration campaigns, appropriate and new statistical techniques, calibration of instruments and characterization of background, preservation of knowledge, and results for the benefit of the astronomical community.

We present 1.8 years of near-daily Swift monitoring of the bright, strongly variable Type 1 AGN Fairall 9. Totaling 575 successful visits, this is the largest such campaign reported to date. Variations within the UV/optical are well-correlated, with longer wavelengths lagging shorter wavelengths in the direction predicted by thin disk/lamp-post models. The correlations are improved by detrending; subtracting a second-order polynomial fit to the UV/optical light curves to remove long-term trends that are not of interest to this study. Extensive testing indicates detrending with higher-order polynomials removes too much intrinsic variability signal on reverberation timescales. These data provide the clearest detection to date of interband lags within the UV, indicating that neither emission from a large disk nor diffuse continuum emission from the broad-line region can independently explain the full observed lag spectrum. The observed X-ray flux variations are poorly correlated with those in the UV/optical. Further, subdivision of the data into four ~160 day light curves shows that the UV/optical lag spectrum is highly stable throughout the four periods, but the X-ray to UV lags are unstable, significantly changing magnitude and even direction from one period to the next. This indicates the X-ray to UV relationship is more complex than predicted by the simple reprocessing model often adopted for AGN. A bowl model (lamp-post irradiation and blackbody reprocessing on a disk with a steep rim) fit suggests the disk thickens at a distance (~10 lt-day) and temperature (~8000K) consistent with the inner edge of the BLR.

By taking into account the latest observations and theoretical constraints, we investigate the merger and post-merger of binary neutron stars (NSs) with numerical simulations employing hadronic and hybrid equations of state (EOSs). We name our hybrid stars Neutron-quark stars (NQS), because the transition from hadrons to quarks starts at a density lower than the central density of $\sim 1 M_{\odot}$ stars. The two scenarios of transition to quark matter, a strong first-order phase transition (1PT) or a crossover, feature either a drop to almost zero or a rapid increase (peak) in the square of the sound speed $c_s^2$, implying a softening or stiffening during the transition, respectively. Although the properties of NQSs in equilibrium may not be distinguishable from those of NSs, we find that the post-merger gravitational-wave (GW) main frequency $f_2$ for the crossover scenario is generally lower than that of hadronic models with the same tidal deformability, indicating that a crossover transition is in principle observable when both the inspiral and post-merger signals are detected. Since it is viable according to current multi-messenger constraints, we also consider an EOS with a 1PT taking place at 1.8 times the nuclear saturation density ($n_0$), with a stiff quark EOS ($c_s^2 = 2/3~c^2$) after the transition. It is the first time that such a binary merger is studied numerically in full general relativity. Although its $f_2$ is $\sim 300$ Hz higher than that of its baseline, the relation between $f_2$ and the tidal deformability of inspiralling stars is close to that for hadronic EOSs.

Coronal mass ejections (CMEs) stand as intense eruptions of magnetized plasma from the Sun, playing a pivotal role in driving significant changes of the heliospheric environment. Deducing the properties of CMEs from their progenitors in solar source regions is crucial for space weather forecasting. Deducing the properties of CMEs from their progenitors in solar source regions is crucial for space weather forecasting. The primary objective of this paper is to establish a connection between CMEs and their progenitors in solar source regions, enabling us to infer the magnetic structures of CMEs before their full development. To this end, we create a dataset comprising a magnetic flux rope series with varying projection shapes, sizes and toroidal fluxes, using the Regularized Biot-Savart Laws (RBSL). Thereafter, we simulate the propagation of these flux ropes from the solar surface to a distance of 25$R_{\odot}$ with our global coronal MHD model which is named COCONUT. Our parametric survey reveals significant impacts of source flux ropes on the consequent CMEs. We find that the projection shape can influence the magnetic structures of CMEs at 20$R_{\odot}$, albeit with minimal impacts on the propagation speed. However, these impacts diminish as source flux ropes become fat. In terms of toroidal flux, our simulation results demonstrate a pronounced correlation with the propagation speed of CMEs, as well as the successfulness in erupting. This work builds the bridge between the CMEs in the outer corona and their progenitors in solar source regions. Our parametric survey suggests that the projection shape, cross-section radius and toroidal flux of source flux ropes are crucial parameters in predicting magnetic structures and propagation speed of CMEs, providing valuable insights for space weather prediction.

We present a full analysis of galaxy major merger pair fractions, merger rates, and mass accretion rates, thus uncovering the role of mergers in galaxy formation at the earliest previously unexplored epoch of $4.5<z<11.5$. We target galaxies with masses $\log_{10}(\mathrm{M}_*/\mathrm{M}_\odot) = 8.0 - 10.0$, utilizing data from eight JWST Cycle-1 fields (CEERS, JADES GOODS-S, NEP-TDF, NGDEEP, GLASS, El-Gordo, SMACS-0723, MACS-0416), covering an unmasked area of 189.36 $\mathrm{arcmin}^2$. We develop a new probabilistic pair-counting methodology that integrates full photometric redshift posteriors and corrects for detection incompleteness to quantify close pairs with physical projected separations between 20 and 50 kpc. Our analysis reveals an increase in pair fractions up to $z = 8$, reaching $0.211 \pm 0.065$, followed by a statistically flat evolution to $z = 11.5$. We find that the galaxy merger rate increases from the local Universe up to $z = 6$ and then stabilizes at a value of $\sim 6$ Gyr$^{-1}$ up to $z = 11.5$. We fit both a power-law and a power-law + exponential model to our pair fraction and merger rate redshift evolution, finding that the latter model describes the trends more accurately, particularly at $z = 8.0 - 11.5$. In addition, we measure that the average galaxy increases its stellar mass due to mergers by a factor of $2.77 \pm 0.99$ from redshift $z = 10.5$ to $z = 5.0$. Lastly, we investigate the impact of mergers on galaxy stellar mass growth, revealing that mergers contribute $71 \pm 25\%$ as much to galaxy stellar mass increases as star formation from gas. This indicates that mergers drive about half of galaxy assembly at high redshift.