Papers for Tuesday, Apr 18 2023

Laser Interferometer Space Antenna (LISA) will observe gravitational waves from galactic binaries (GBs) of white dwarfs or neutron stars. Some of these objects are among the most magnetic astrophysical objects in the Universe. Magnetism, by secularly disrupting the orbit, can eventually affect the gravitational waves emission and could then be potentially detected and characterized after several years of observations by LISA. Currently, the data processing pipeline of the LISA Data Challenge (LDC) for GBs does not consider either magnetism or eccentricity. Recently, it was shown [Bourgoin et al. PRD 105, 124042 (2022)] that magnetism induces a shift on the gravitational wave frequencies. Additionally, it was argued that, if the binary's orbit is eccentric, the presence of magnetism could be detected by LISA. In this work, we explore the consequences of a future data analysis conducted on quasi-circular and magnetic GB systems using the current LDC tools. We first show that a single eccentric GB can be interpreted as several GBs and this can eventually bias population studies deduced from LISA's future catalog. Then, we confirm that for quasi-circular orbits, the secular magnetic energy of the system can be inferred if the signal-to-noise ratio of the second harmonic is high enough to be detected by traditional quasi-monochromatic source searching algorithms. LISA observations could therefore bring new insights on the nature and origin of magnetic fields in white dwarfs or neutron stars.

A growing number of far-infrared bright sources completely invisible in deep extragalactic optical surveys hint at an elusive population of z>4 dusty, star-forming galaxies. Cycle 1 JWST surveys are now detecting their rest-frame optical light, which provides key insight into their stellar properties and statistical constraints on the population as a whole. This work presents the JWST/NIRCam counterpart from the COSMOS-Web survey to a far-infrared SCUBA-2 and ALMA source, AzTECC71, which was previously undetected at wavelengths shorter than 850 microns. AzTECC71, amongst the reddest galaxies in COSMOS-Web with F277W - F444W~0.9, is undetected in NIRCam/F150W and F115W and fainter in F444W than other sub-millimeter galaxies identified in COSMOS-Web by 2-4 magnitudes. This is consistent with the system having both a lower stellar mass and higher redshift than the median dusty, star-forming galaxy. With deep ground- and space-based upper limits combined with detections in F277W, F444W and the far-IR including ALMA Band 6, we find a high probability (99%) that AzTECC71 is at z>4 with z_phot=5.7(+0.8,-0.7). This galaxy is massive (logM*/Msun~10.7) and IR-luminous (logLIR/Lsun~12.7), comparable to other optically-undetected but far-IR bright dusty, star-forming galaxies at z>4. This population of luminous, infrared galaxies at z>4 is largely unconstrained but comprises an important bridge between the most extreme dust-obscured galaxies and more typical high-redshift star-forming galaxies. If further far-IR-selected galaxies that drop out of the F150W filter in COSMOS-Web have redshifts z>4 like AzTECC71, then the volume density of such sources may be ~3-10x greater than previously estimated.

Over the last decade, observations have shown that the mean mass of ultra-high-energy cosmic rays (UHECRs) increases progressively toward the highest energies. However, the precise composition is still unknown, and several theoretical studies hint at the existence of a subdominant proton component up to the highest energies. Motivated by the exciting prospect of performing charged-particle astronomy with ultra-high-energy (UHE) protons we quantify the level of UHE-proton flux that is compatible with present multimessenger observations and the associated fluxes of neutral messengers produced in the interactions of the protons. We study this scenario with numerical simulations of two independent populations of extragalactic sources and perform a fit to the combined UHECR energy spectrum and composition observables, constrained by diffuse gamma-ray and neutrino observations. We find that up to of order $10\%$ of the cosmic rays at the highest energies can be UHE protons, although the result depends critically on the selected hadronic interaction model for the air showers. Depending on the maximum proton energy ($E_\text{max}^\text{p}$) and the redshift evolution of sources, the associated flux of cosmogenic neutrinos and UHE gamma rays can significantly exceed the multimessenger signal of the mixed-mass cosmic rays. Moreover, if $E_\text{max}^\text{p}$ is above the GZK limit, we predict a large flux of UHE neutrinos above EeV energies that is absent in alternate scenarios for the origin of UHECRs. We present the implications and opportunities afforded by these UHE proton, neutrino and photon fluxes for future multimessenger observations.

The {\it Gaia} eDR3 catalogue has recently been used to construct samples of nearby wide binaries to study the internal kinematics of these objects using relative velocities of the two component stars, $\Delta V$, total binary masses, $m_{B}$, and separations, $s$. For $s \gtrsim 0.035$ pc, these binaries probe the low acceleration $a<a_{0}$ regime over which the gravitational anomalies usually attributed to dark matter are observed in the flat rotation curves of spiral galaxies, where $a_{0}\approx 1.2\times 10^{10}$ is the acceleration scale of MOND. Such experiments test the degree of generality of these anomalies, by exploring the same acceleration regime using independent astronomical systems of vastly smaller mass and size. A signal above Newtonian expectations has been observed when $a<a_{0}$, alternatively interpreted as evidence of a modification in the relevant fundamental physics, or as being due to kinematic contaminants affecting the experiment; the presence of undetected stellar components, unbound encounters and spurious projection effects. Here I take advantage of the enhanced DR3 {\it Gaia} catalogue to perform a more rigorous and detailed study of the internal kinematics of wide binaries than what has previously been possible. Having internally determined accurate {\it Gaia} stellar masses and estimates of binary probabilities for each star using spectroscopic information, together with a larger sample of radial velocities, allows for a significant improvement in the analysis of wide binaries and careful exclusion of possible kinematic contaminants. Resulting $\Delta V$ vs. $s$ and $\Delta V$ vs. $m_{B}$ scalings accurately tracing Newtonian expectations for the high acceleration regime, but consistent with the distance and mass velocity scalings observed in spiral galaxies in the low acceleration one, are obtained.

We report strong and rapid X-ray variability found from the super-Eddington accreting quasar SDSS J081456.10+532533.5 at $z=0.1197$. It has a black-hole mass of $2.7\times10^{7}{M_{\odot}}$ and a dimensionless accretion rate of $\approx4$ measured from reverberation-mapping observations. It showed weak X-ray emission in the 2021 February Chandra observation, with the 2 keV flux density being $9.6^{+11.6}_{-4.6}$ times lower compared to an archival Swift observation. The 2 keV flux density is also $11.7^{+9.6}_{-6.3}$ times weaker compared to the expectation from its optical/UV emission. In a follow-up XMM-Newton observation 32 days later, the 2 keV flux density increased by a factor of $5.3^{+6.4}_{-2.4}$, and the spectra are best described by a power law modified with partial-covering absorption; the absorption-corrected intrinsic continuum is at a nominal flux level. Nearly simultaneous optical spectra reveal no variability, and there is only mild long-term optical/infrared variability from archival data (with a maximum variability amplitude of $\approx50\%$). We interpret the X-ray variability with an obscuration scenario, where the intrinsic X-ray continuum does not vary but the absorber has variable column density and covering factor along the line of sight. The absorber is likely the small-scale clumpy accretion wind that has been proposed to be responsible for similar X-ray variability in other super-Eddington accreting quasars.

We present multi-epoch spectroscopic follow-up of a sample of ellipsoidal variables selected from Gaia DR3 as candidates for hosting quiescent black holes (BHs) and neutron stars (NSs). Our targets were identified as BH/NS candidates because their optical light curves -- when interpreted with models that attribute variability to tidal distortion of a star by a companion that contributes negligible light -- suggest that the companions are compact objects. From the likely BH/NS candidates identified in recent work accompanying Gaia DR3, we select 14 of the most promising targets for follow-up. We obtained spectra for each object at 2-10 epochs, strategically observing near conjunction to best-constrain the radial velocity semi-amplitude. From the measured semi-amplitudes of the radial velocity curves, we derive minimum companion masses of $M_{2,\min} \leq 0.5 ~ M_{\odot}$ in all cases. Assuming random inclinations, the typical inferred companion mass is $M_2 \sim 0.15 ~ M_{\odot}$. This makes it unlikely that any of these systems contain a BH or NS, and we consider alternative explanations for the observed variability. We can best reproduce the observed light curves and radial velocities with models for unequal-mass contact binaries with starspots. Some of the objects in our sample may also be detached main-sequence binaries, or even single stars with pulsations or starspot variability masquerading as ellipsoidal variation. We provide recommendations for future spectroscopic efforts to further characterize this sample and more generally to search for compact object companions in close binaries.

Direct $N$-body simulations of a large number of particles, especially in the study of planetesimal dynamics and planet formation, have been computationally challenging even with modern machines. This work presents the combination of fully parallelized $N^2/2$ interactions and the incorporation of the GENGA code's close encounter pair grouping strategy to enable MIMD parallelization of the Symplectic Massive Body Algorithm (SyMBA) with OpenMP on multi-core CPUs in shared-memory environment. SyMBAp (SyMBA parallelized) preserves the symplectic nature of SyMBA and shows good scalability, with a speedup of 30.8 times with 56 cores in a simulation with 5,000 fully interactive particles.

Magnetic fields are predicted to have a crucial impact on the structure, evolution and chemistry of protoplanetary disks. However, a direct detection of the magnetic field towards these objects has yet to be achieved. In order to characterize protoplanetary disk magnetic fields, we investigate the impact of the Zeeman effect on the (polarized) radiative transfer of emission from paramagnetic molecules excited in protoplanetary disks. While the effects of the Zeeman effect are commonly studied in the circular polarization of spectral lines, we perform a comprehensive modeling also of the Zeeman-induced broadening of spectral lines and their linear polarization. We develop simplified radiative transfer models adapted to protoplanetary disks, which we compare to full three-dimensional polarized radiative transfer simulations. We find that the radiative transfer of circular polarization is heavily affected by the expected polarity-change of the magnetic field between opposite sides of the disk. In contrast, Zeeman broadening and linear polarization are relatively unaffected by this sign change due to their quadratic dependence on the magnetic field. We can match our simplified radiative transfer models to full polarization modeling with high fidelity, which in turn allows us to prescribe straight-forward methods to extract magnetic field information from Zeeman broadening and linear polarization observations. We find that Zeeman broadening and linear polarization observations are highly advantageous methods to characterize protoplanetary disk magnetic fields as they are both sensitive probes of the magnetic field and are marginally affected by any sign change of the disk magnetic field. Applying our results to existing circular polarization observations of protoplanetary disk spectral lines suggests that the current upper limits on the toroidal magnetic field strengths have to be raised.

We present analysis using a citizen science campaign to improve the cosmological measures from the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). The goal of HETDEX is to measure the Hubble expansion rate, $H(z)$, and angular diameter distance, $D_A(z)$, at $z =$ 2.4, each to percent-level accuracy. This accuracy is determined primarily from the total number of detected Lyman-$\alpha$ emitters (LAEs), the false positive rate due to noise, and the contamination due to [O II] emitting galaxies. This paper presents the citizen science project, Dark Energy Explorers, with the goal of increasing the number of LAEs, decreasing the number of false positives due to noise and the [O II] galaxies. Initial analysis shows that citizen science is an efficient and effective tool for classification most accurately done by the human eye, especially in combination with unsupervised machine learning. Three aspects from the citizen science campaign that have the most impact are 1) identifying individual problems with detections, 2) providing a clean sample with 100% visual identification above a signal-to-noise cut, and 3) providing labels for machine learning efforts. Since the end of 2022, Dark Energy Explorers has collected over three and a half million classifications by 11,000 volunteers in over 85 different countries around the world. By incorporating the results of the Dark Energy Explorers we expect to improve the accuracy on the $D_A(z)$ and $H(z)$ parameters at $z =$ 2.4 by 10 - 30%. While the primary goal is to improve on HETDEX, Dark Energy Explorers has already proven to be a uniquely powerful tool for science advancement and increasing accessibility to science worldwide.

In the near-UV and optical transmission spectrum of the hot Jupiter WASP-121b, recent observations have detected strong absorption features of Mg, Fe, Ca, and H$\alpha$, extending outside of the planet's Roche lobe. Studying these atomic signatures can directly trace the escaping atmosphere and constrain the energy balance of the upper atmosphere. To understand these features, we introduce a detailed forward model by expanding the capability of a one-dimensional model of the upper atmosphere and hydrodynamic escape to include important processes of atomic metal species. The hydrodynamic model is coupled to a Ly$\alpha$ Monte Carlo radiative transfer calculation to simulate the excited hydrogen population and associated heating/ionization effects. Using this model, we interpret the detected atomic features in the transmission spectrum of WASP-121b and explore the impact of metals and excited hydrogen on its upper atmosphere. We demonstrate the use of multiple absorption lines to impose stronger constraints on the properties of the upper atmosphere than the analysis of a single transmission feature can provide. In addition, the model shows that line broadening due to atmospheric outflow driven by the Roche lobe overflow is necessary to explain the observed line widths and highlights the importance of the high mass-loss rate caused by the Roche lobe overflow that requires careful consideration of the structure of the lower and middle atmosphere. We also show that metal species and excited state hydrogen can play an important role in the thermal and ionization balance of ultra-hot Jupiter thermospheres.

We have gathered optical-region spectra, derived model atmosphere parameters, and computed elemental abundances for 15 red giant stars in the open cluster NGC 7789. We focus on the light element group CNOLi that provides clues to evolutionary changes associated with internal fusion events and chemical mixing. We confirm and extend an early report that NGC 7789 stars 193 and 301 have anomalously large Li abundances, and that these values are apparently unconnected to any other elements' abundances in these stars. A companion study of He I lambda 10830 lines in both field stars and cluster members shows that star 301 has a strong He feature while star 193 does not. Possible explanations for the large Li abundances of these stars include helium flash-induced mixing events and binary interactions at some past or present times. In either case an internal eruption of energy could cause fresh synthesis of lithium via the Cameron-Fowler Beryllium transport mechanism. Rapid transport of lithium to the outer layers may have created significant chromospheric transient disturbances, producing enough helium ionization to allow for the strong lambda 10830 absorption in star 301.

We report on the disk-averaged absolute brightness temperatures of Venus measured at four microwave frequency bands with the Cosmology Large Angular Scale Surveyor (CLASS). We measure temperatures of 432.3 $\pm$ 2.8 K, 355.6 $\pm$ 1.3 K, 317.9 $\pm$ 1.7 K, and 294.7 $\pm$ 1.9 K for frequency bands centered at 38.8, 93.7, 147.9, and 217.5 GHz, respectively. We do not observe any dependence of the measured brightness temperatures on solar illumination for all four frequency bands. A joint analysis of our measurements with lower frequency Very Large Array (VLA) observations suggests relatively warmer ($\sim$ 7 K higher) mean atmospheric temperatures and lower abundances of microwave continuum absorbers than those inferred from prior radio occultation measurements.

Active galactic nuclei (AGN) in the early Universe are thought to be prominent sources of energy and ionizing photons that affected the growth of their host galaxy and their environment. However, it is still unclear how the supermassive black holes (SMBHs) that fuel these AGN grew to the observed high masses already at high redshifts. Observations of high-redshift SMBH progenitors or lower-mass AGN will thus help characterize the evolution of SMBHs and their impact on the surroundings. With the launch of the JWST, fainter objects at high redshifts can now be detected, including lower-mass AGN. We assess the observability of such low luminosity AGN, using the cosmological simulation code GIZMO to provide a realistic environment for black hole growth in the early Universe. Soon after the first stars are born in the simulation run, we insert stellar-remnant black hole seeds of various initial masses, between $300$ and $10^4 {\rm \ M}_{\odot}$, at the center of a dark matter halo and follow their growth until $z\sim6$. Such stellar black hole seeds placed in a typical high-$z$ environment do not significantly accrete and grow to reach masses that can be observed with the JWST under conditions of standard Bondi-Hoyle accretion, as energy input from stellar feedback and chaotic dynamics prevent efficient gas accretion onto the black holes. To be observed with the JWST, rarer but still physically feasible growth regimes, involving Eddington or super-Eddington accretion, would be required. Alternatively, AGN observability may be boosted under even rarer conditions of extreme gravitational lensing.

HD 53143 is a mature Sun-like star and host to a broad disk of dusty debris, including a cold outer ring of planetesimals near 90 AU. Unlike most other inclined debris disks imaged at visible wavelengths, the cold disk around HD 53143 appears as disconnected "arcs" of material, with no forward scattering side detected to date. We present new, deeper Hubble Space Telescope (HST) Space Telescope Imaging Spectrograph (STIS) coronagraphic observations of the HD 53143 debris disk and show that the forward scattering side of the disk remains undetected. By fitting our KLIP-reduced observations via forward modeling with an optically thin disk model, we show that fitting the visible wavelength images with an azimuthally symmetric disk with unconstrained orientation results in an unphysical edge-on orientation that is at odds with recent ALMA observations, while constraining the orientation to that observed by ALMA results in nearly isotropically scattering dust. We show that the HD53143 host star exhibits significant stellar variations due to spot rotation and revisit age estimates for this system.

We present measurements of the interferometrically-resolved binary star system 12 Com and the single giant star 31 Com in the cluster Coma Berenices. 12 Com is a double-lined spectroscopic binary system consisting of a G7 giant and an A3 dwarf at the cluster turnoff. Using an extensive radial velocity dataset and interferometric measurements from PTI and the CHARA array, we measured masses $M_1 =2.64 \pm 0.07 M_\odot$ and $M_2 =2.10 \pm 0.03 M_\odot$. Interferometry also allows us to resolve the giant, and measure its size as $R_1 = 9.12 \pm 0.12 \pm 0.01 R_\odot$. With the measured masses and radii, we find an age of $533 \pm 41 \pm 42$ Myr. For comparison, we measure the radius of 31 Com to be $8.36 \pm 0.15 R_\odot$. Based on the photometry and radius measurements, 12 Com A is likely the most evolved bright star in the cluster, large enough to be in the red giant phase, but too small to have core helium burning. Simultaneous knowledge of 12 Com A's mass and photometry puts strong constraints on convective core overshooting during the main sequence phase, which in turn reduces systematic uncertainties in the age. Increased precision in measuring this system also improves our knowledge of the progenitor of the cluster white dwarf WD1216+260.

After having left the heliosphere, Voyager 1 and Voyager 2 continue to travel through interstellar space. The Pioneer 10, Pioneer 11, and New Horizons spacecraft are also on paths to pass the heliopause. These spacecraft have communicated with the Deep Station Network (DSN) radio antennas in order to download scientific data and telemetry data. Outward transmissions from DSN travel to the spacecraft and beyond into interstellar space. These transmissions have encountered and will encounter other stars, introducing the possibility that intelligent life in other solar systems will encounter our terrestrial transmissions. We use the beamwidth of the transmissions between DSN and interstellar spacecraft to perform a search around the past and future positions of each spacecraft obtained from the JPL Horizons System. By performing this search over the Gaia Catalogue of Nearby Stars (GCNS), a catalogue of precisely mapped stars within 100 pc, we determine which stars the transmissions of these spacecraft will encounter. We highlight stars that are in the background of DSN transmissions and calculate the dates of these encounters to determine the time and place for potential intelligent extraterrestrial life to encounter terrestrial transmissions.

We place lower limits on the obliquities between debris disks and their host stars for 31 systems by comparing their disk and stellar inclinations. While previous studies did not find evidence for misalignment, we identify 6 systems with minimum obliquities falling between ~30{\deg}-60{\deg}, indicating that debris disks can be significantly misaligned with their stars. These high-obliquity systems span a wide range of stellar parameters with spectral types K through A. Previous works have argued that stars with masses below 1.2 $M_\odot$ (spectral types of ~F6) have magnetic fields strong enough to realign their rotation axes with the surrounding disk via magnetic warping; given that we observe high obliquities for relatively low-mass stars, magnetic warping alone is likely not responsible for the observed misalignment. Yet, chaotic accretion is expected to result in misalignments of ~20{\deg} at most and cannot explain the larger obliquities found in this work. While it remains unclear how primordial misalignment might occur and what role it plays in determining the spin-orbit alignment of planets, future work expanding this sample is critical towards understanding the mechanisms that shape these high-obliquity systems.

The soft X-ray excess in the spectra of active galactic nuclei is characterized by similar electron temperatures of 0.1 -- 0.3 keV and similar photon indices around 2.2 -- 3. It remains a puzzle why both values are not sensitive to the black hole mass nor accretion rate. Supposing that the scattering-dominated surface layer of an accretion disk can act as a warm corona, we construct a vertically one-zone model to understand what determines its temperature. By solving the equations of (1) the condition for the effective optical depth, (2) the energy balance, and (3) dominance of the Compton cooling over the bound-free cooling, we could reproduce the basic observational features of the soft excess, provided that anomalous heating takes place in the warm corona. The similar temperatures can be understood, since both of the anomalous heating and Compton cooling rates are proportional to the dissipation rate of the accretion energy, while similar photon indices are a natural consequence of the fact that observed photons are finally emitted from the layer of Compton $y\sim 1$. The warm corona solutions only exist at smaller radii, indicating the structure of a warm corona inside and a hot corona outside. The soft excess is not observed in black hole binaries, since disk temperatures are too high for the Compton scattering to work as cooling. The derived temperatures are somewhat underestimation, however. This may indicate a necessity of multi-zone corona structure. The stability of the warm corona and its consequences are briefly discussed.

Environmental seismic disturbances limit the sensitivity of LIGO gravitational wave detectors. Trains near the LIGO Livingston detector produce low frequency (0.5-10 Hz) ground noise that couples into the gravitational wave sensitive frequency band (10-100 Hz) through light reflected in mirrors and other surfaces. We investigate the effect of trains during the Advanced LIGO third observing run, and propose a method to search for narrow band seismic frequencies responsible for contributing to increases in scattered light. Through the use of the linear regression tool Lasso (least absolute shrinkage and selection operator) and glitch correlations, we identify the most common seismic frequencies that correlate with increases in detector noise as 0.6-0.8 Hz, 1.7-1.9 Hz, 1.8-2.0 Hz, and 2.3-2.5 Hz in the LIGO Livingston corner station.

Cyg X-1 is probably the most studied and best understood black-hole X-ray binary. Recently, its accretion geometry has been probed with the X-ray polarization. The position angle of the X-ray emitting flow was found to be aligned with the position angle of the radio jet in the plane of the sky. At the same time, the observed high polarization degree could be obtained only for a high inclination of the X-ray emitting flow, indicating a misalignment between the binary axis and the black hole spin. The jet, in turn, is believed to be directed by the spin axis, hence similar misalignment is expected between the jet and binary axes. We test this hypothesis using very long (up to about 26 years) multi-band radio observations. We find the misalignment of $20^\circ$-$30^\circ$; however, on the contrary to the earlier expectations, the misalignment lies in the plane of the sky and not along the line of sight. Furthermore, the presence of the misalignment questions our understanding of the evolution of this binary system.

The solar magnetic fields emerging from the photosphere into the chromosphere and corona are comprised of a combination of "closed" and "open" fields. The closed magnetic field lines are defined as those having both ends rooted in the solar surface, while the open field lines are those having one end extending out into interplanetary space and the other rooted at the Sun's surface. Since the early 2000's, the amount of total unsigned open magnetic flux estimated by coronal models have been in significant disagreement with in situ spacecraft observations, especially during solar maximum. Estimates of total open unsigned magnetic flux using coronal hole observations (e.g., using extreme ultraviolet (EUV) or Helium (He) I) are in general agreement with the coronal model results and thus show similar disagreements with in situ observations. While several possible sources producing these discrepancies have been postulated over the years, there is still no clear resolution to the problem. This paper provides a brief overview of the problem and summarizes some proposed explanations for the discrepancies. In addition, two different ways of estimating the total unsigned open magnetic flux are presented, utilizing the Wang-Sheeley-Arge (WSA) model, and one of the methods produce surprisingly good agreement with in situ observations. The findings presented here suggest that active regions residing near the boundaries of mid-latitude coronal holes are the probable source of the missing open flux. This explanation also brings in line many of the seemly contradictory facts that have made resolving this problem so difficult.

We analyzed Asteroid Terrestrial-impact Last Alert System (ATLAS), Zwicky Transient Facility (ZTF) and All-Sky Automated Survey for Supernovae (ASAS-SN) data of MASTER OT J055845.55+391533.4 and found that this object repeats superoutburst with a dip in the middle of the outburst followed by long and sometimes oscillating rebrightening, just like a WZ Sge-type dwarf nova or an AM CVn-type object. The mean supercycle was 298(8) d, too short for a WZ Sge star, but with only a few normal outbursts. We also observed the 2023 February-March superoutburst and established the superhump period of 0.05509(2) d. This period appears to exclude the possibility of an AM CVn star. Although the 2023 observations could not detect superhumps after the dip, the 2014, 2016 and 2021 data seem to suggest that low-amplitude superhumps were present during the rebrightening phase. We note that a dip during a superoutburst is a feature common to the unusual SU UMa-type dwarf nova MASTER OT J172758.09+380021.5 during some of its superoutbursts. These objects may comprise a new class of rebrightening phenomenon in SU UMa-type dwarf novae.

The Gaia-ESO Survey is a public spectroscopic survey that has targeted $\gtrsim10^5$ stars covering all major components of the Milky Way from the end of 2011 to 2018, delivering its public final release in May 2022. Unlike other spectroscopic surveys, Gaia-ESO is the only survey that observed stars across all spectral types with dedicated, specialised analyses: from O ($T_\mathrm{eff} \sim 30,000-52,000$~K) all the way to K-M ($\gtrsim$3,500~K). The physics throughout these stellar regimes varies significantly, which has previously prohibited any detailed comparisons between stars of significantly different type. In the final data release (internal data release 6) of the Gaia-ESO Survey, we provide the final database containing a large number of products such as radial velocities, stellar parameters and elemental abundances, rotational velocity, and also, e.g., activity and accretion indicators in young stars and membership probability in star clusters for more than 114,000 stars. The spectral analysis is coordinated by a number of Working Groups (WGs) within the Survey, which specialise in the various stellar samples. Common targets are analysed across WGs to allow for comparisons (and calibrations) amongst instrumental setups and spectral types. Here we describe the procedures employed to ensure all Survey results are placed on a common scale to arrive at a single set of recommended results for all Survey collaborators to use. We also present some general quality and consistency checks performed over all Survey results.

We investigate the stability of nearby disc galaxies and galaxies at redshift ($z$) equal to 4.5. We explore the connection between the stability parameter $(Q_{RW})$, star formation rate ($SFR$), gas fraction $(f^{Gas})$, and the time scale for growth of gravitational instabilities $(\tau)$. We find that, despite differences in morphology $91$ $\%$ of the nearby galaxies have a minimum value of stability parameter ($Q^{Min}_{RW}$) greater than $1$ indicating stability against the growth of axisymmetric instabilities. The spirals in our sample have higher median star formation rate, lower median $Q_{RW}$, a lower $f^{Gas}$ and small time scale for growth of gravitational instabilities than irregular galaxies. We find that the gravitational instabilities in spirals convert a large fraction of gas into stars quickly, depleting the gas reservoirs. On the other hand, star formation occurs more gradually over longer timescales in irregulars with a higher gas fraction. We then compare the stability of the nearby galaxies with galaxies at $z\,=\,4.5$. We find that net stability levels in the nearby galaxies and the galaxies at $z\,=\,4.5$ are primarily driven by the stellar disc suggesting the presence of an inherent mechanism that self-regulates the stability. Finally, upon removing the contribution of the dark matter to the total potential, the median $Q_{RW}$ for the nearby galaxies and galaxies at $z \,= \,4.5$ remains unchanged indicating that the baryons can self-regulate the stability levels, at least in a statistical sense.

Close galaxy flybys, interactions during which two galaxies inter-penetrate, are frequent and can significantly affect the evolution of individual galaxies. Equal-mass flybys are extremely rare and almost exclusively distant, while frequent flybys have mass ratios 0.1 or lower, with a secondary galaxy penetrating deep into the primary. This can result in comparable strengths of interaction S between the two classes of flybys and lead to essentially the same effects. To demonstrate this, emphasize and explore the role of the impact parameter $b$, we performed a series of N-body simulations of flybys with varying relative $b$ ranging from 0.114 to 0.272 of the virial radius of the primary. Two-armed spirals form during flybys, with radii of origin correlated with $b$ and strengths well approximated with an inverted S-curve. The impact parameter does not affect the shape of induced spirals, and the lifetimes of a distinguished spiral structure appear to be constant, 2 Gyr. Bars, with strengths anti-correlated with $b$, form after the encounter is over in simulations with interaction strengths $S\geq0.076$, but they are short-lived except for the stronger ones with $S\geq0.129$. We showcase an occurrence of double bar that survives for a long time in one of the simulations. Effects on the pre-existing bar instability are diverse. There is no uniform correlation between these effects and $b$, as they are secondary effects, happening later in a post-flyby stage. Bulges are resilient to flybys, while dark matter halos can significantly spin up in the amount anti-correlated with $b$. There is an offset angle between the angular momentum vector of the dark matter halo and that of a disc, and it correlates linearly with b. Flybys remain an important pathway for structural evolution within galaxies in the local Universe.

I present the first evidence of multiple populations in the globular cluster (GCs) 47Tucanae based on images collected with the near-infrared camera (NIRCam) on board the James Webb Space Telescope (JWST). While NIRCam photometry is poorly sensitive to multiple populations among stars brighter than the main-sequence (MS) knee, the M-dwarfs more-massive than 0.1 solar masses define a wide F115W-F322W2 color range due to multiple populations. The star-to-star color differences are mostly due to the different amounts of water vapor (hence oxygen) that affect the spectra of M-dwarfs. The chromosome map unveils an extended first population (1P) composed of M-dwarfs with different metallicities and three main groups of second-population (2P) stars that are depleted in oxygen with respect to the 1P. I present the discovery of an MS of very-low-mass stars (masses smaller than 0.1 solar masses) and tentatively associated it with a sequence composed of O-rich stars alone.

Three-dimensional (3D) kinetic maps of the Milky Way interstellar medium are an essential tool in studies of its structure and of star formation. We aim to assign radial velocities to Galactic interstellar clouds now spatially localized based on starlight extinction and star distances from Gaia and stellar surveys. We developed an automated search for coherent projections on the sky of clouds isolated in 3D extinction density maps on the one hand, and regions responsible for CO radio emissions at specific Doppler shifts on the other hand. The discrete dust structures were obtained by application of the Fellwalker algorithm to a recent 3D extinction density map. For each extinction cloud, a technique using a narrow sliding spectral window moved along the contour-bounded CO spectrum and geometrical criteria was used to select the most likely velocity interval. We applied the new contour-based technique to the 3D extinction density distribution within the volume encompassing the Taurus, Auriga, Perseus and California molecular complexes. From the 45 clouds issued from the decomposition, 42 were assigned a velocity. The remaining structures correspond to very weak CO emission or extinction. We used the non-automated assignments of radial velocities to clouds of the same region presented in paper I and based on KI absorption spectra as a validation test. The new fully automated determinations were found in good agreement with these previous measurements. Our results show that an automated search based on cloud contour morphology can be efficient and that this novel technique may be extended to wider regions of the Milky Way and at larger distance. We discuss its limitations and potential improvements after combination with other techniques.

Hydrodynamic atmospheric escape is considered an important process that shapes the evolution of sub-Jovian exoplanets, particularly those with short orbital periods. The metastable He line in the near-infrared at $1.083$ $\mu$m is a reliable tracer of atmospheric escape in hot exoplanets, with the advantage of being observable from the ground. However, observing escaping He in sub-Jovian planets has remained challenging due to the systematic effects and telluric contamination present in ground-based data. With the successful launch and operations of JWST, we now have access to extremely stable high-precision near-infrared spectrographs in space. Here we predict the observability of metastable He with JWST in two representative and previously well-studied warm Neptunes, GJ 436 b ($T_{\rm eq} = 687~{\rm K}$, $R_{\rm p} = 0.37~{\rm R_J}$) and GJ 1214 b ($T_{\rm eq} = 588~{\rm K}$, $R_{\rm p} = 0.25~{\rm R_J}$). Our simulated JWST observations for GJ 436 b demonstrate that a single transit with NIRSpec/G140H is sensitive to mass loss rates that are two orders of magnitude lower than what is detectable from the ground. Our exercise for GJ 1214 b show that the best configuration to observe the relatively weak outflows of warm Neptunes with JWST is with NIRSpec/G140H, and that NIRSpec/G140M and NIRISS/SOSS are less optimal. Since none of these instrument configurations can spectrally resolve the planetary absorption, we conclude that the 1D isothermal Parker-wind approximation may not be sufficient for interpreting such observations. More sophisticated models are critical for breaking the degeneracy between outflow temperature and mass-loss rate for JWST measurements of metastable He.

The energy injected from dark matter annihilation and decay processes potentially raises the ionisation of the intergalactic medium and leaves visible footprints on the anisotropy maps of the cosmic microwave background (CMB). Galactic foregrounds emission in the microwave bands contaminate the CMB measurement and may affect the search for dark matter's signature. In this paper, we construct a full CMB data and foreground simulation based on the design of the next-generation ground-based CMB experiments. The foreground residual after the components separation on maps is fully considered in our data analysis, accounting for various contamination from the emission of synchrotron, thermal dust, free-free and spinning dust. We analyse the corresponding sensitivity on dark matter parameters from the temperature and polarization maps, and we find that the CMB foregrounds leave a non-zero yet controllable impact on the sensitivity. Comparing with statistics-only analysis, the CMB foreground residual leads to a factor of 7%-23% weakening on energy-injection constraints, depending on the specific dark matter process and experimental configuration. Strong limits on dark matter annihilation rate and decay lifetime can be expected after foreground subtraction.

Using the $\delta N$ formalism we calculate the one-loop correction to the large-scale power spectrum of the curvature perturbation in the standard scenario where primordial black holes are formed in the early universe thanks to a phase of ultra-slow-roll in single-field inflation. We explicitly show that one-loop corrections are negligible when the transition from the ultra-slow-roll to the slow-roll phase is smooth. We conclude that the PBH formation scenario through a ultra-slow-roll phase is viable.

The late phases of the orbital evolution of an Earth-like planet around a Sun-like star are revisited considering the effect of the density fluctuations associated with convective motions inside the star. Such fluctuations produce a random perturbation of the stellar outer gravitational field that excites a small residual eccentricity in the orbit of the planet counteracting the effects of tides that tend to circularize the orbit. We compute the power spectrum of the outer gravitational field fluctuations of the star in the quadrupole approximation and study their effects on the orbit of the planet using a perturbative approach. The residual eccentricity is found to be a stochastic variable showing a Gaussian distribution. Adopting a model of the stellar evolution of our Sun computed with MESA, we find that the Earth will be engulfed close to the tip of the red giant branch evolution phase. We find a maximum mean value of the residual eccentricity of about 0.026 immediately before the engulfment. Considering an Earth-mass planet with an initial orbital semimajor axis sufficiently large to escape engulfment, we find that the mean value of the residual eccentricity is greater than 0.01 for an initial separation up to about 1.4 au. The engulfment of the Earth by the red giant Sun is found to be a stochastic process, in contrast to the deterministic character assumed in previous studies. If an Earth-like planet escapes engulfment, its orbit around its remnant white dwarf star will be moderately eccentric. Such a residual eccentricity on the order of a few hundredths can play a relevant role in sustaining the pollution of the white dwarf atmosphere by asteroids and comets as observed in several objects.

Earth's water, intrinsic oxidation state, and metal core density are fundamental chemical features of our planet. Studies of exoplanets provide a useful context for elucidating the source of these chemical traits. Planet formation and evolution models demonstrate that rocky exoplanets commonly formed with hydrogen-rich envelopes that were lost over time. These findings suggest that Earth may also have formed from bodies with H$_2$-rich primary atmospheres. Here we use a self-consistent thermodynamic model to show that Earth's water, core density, and overall oxidation state can all be sourced to equilibrium between H$_2$-rich primary atmospheres and underlying magma oceans in its progenitor planetary embryos. Water is produced from dry starting materials resembling enstatite chondrites as oxygen from magma oceans reacts with hydrogen. Hydrogen derived from the atmosphere enters the magma ocean and eventually the metal core at equilibrium, causing metal density deficits matching that of Earth. Oxidation of the silicate rocks from solar-like to Earth-like oxygen fugacities also ensues as Si, along with H and O, alloys with Fe in the cores. Reaction with hydrogen atmospheres and metal-silicate equilibrium thus provides a simple explanation for fundamental features of Earth's geochemistry that is consistent with rocky planet formation across the galaxy.

The Australia Telescope Large Area Survey (ATLAS) and the VLA survey in the XMM-LSS/VIDEO deep field provide deep ($\approx 15$ ${\mu}$Jybeam$^{-1}$) and high-resolution ($\approx$ 4.5--8 arcsec) radio coverage of the three XMM-SERVS fields (W-CDF-S, ELAIS-S1, and XMM-LSS). These data cover a total sky area of 11.3 deg$^2$ and contain $\approx 11000$ radio components. Furthermore, about 3~deg$^2$ of the XMM-LSS field also has deeper MIGHTEE data that achieve a median RMS of 5.6 ${\mu}$Jy beam$^{-1}$ and detect more than 20000 radio sources. We analyze all these radio data and find source counterparts at other wavebands utilizing deep optical and IR surveys. The nature of these radio sources is studied using radio-band properties (spectral slope and morphology), and the IR-radio correlation. %and spectral energy distribution. Radio AGNs are selected and compared with those selected using other methods (e.g. X-ray). We found 1765 new AGNs that were not selected using X-ray and/or MIR methods. We constrain the FIR-to-UV SEDs of radio AGNs using {\sc cigale} and investigate the dependence of radio AGN fraction upon galaxy stellar mass and star-formation rate.

By characterizing the contribution of stray light to large datasets from the NuSTAR X-ray observatory collected over 2012--2017, we report a measurement of the cosmic X-ray background in the 3--20 keV energy range. These data represent $\sim20\%$ sky coverage while avoiding Galactic Ridge X-ray emission and are less weighted by deep, survey fields than previous measurements with NuSTAR. Images in narrow energy bands are stacked in detector space and spatially fit with a model representing the stray light and uniform pattern expected from the cosmic X-ray background and the instrumental background, respectively. We establish baseline flux values from Earth-occulted data and validate the fitting method on stray light observations of the Crab, which further serve to calibrate the resulting spectra. We present independent spectra of the cosmic X-ray background with the FPMA and FPMB detector arrays, which are in excellent agreement with the canonical characterization by HEAO 1 and are $10\%$ lower than most subsequent measurements; $F_{\rm{3-20~keV}}^{FPMA} = 2.63 \times 10^{-11}~\rm{erg~s^{-1}~cm^{-2}~deg^{-2}}$ and $F_{\rm{3-20~keV}}^{FPMB} = 2.58 \times 10^{-11}~\rm{erg~s^{-1}~cm^{-2}~deg^{-2}}$. We discuss these results in light of previous measurements of the cosmic X-ray background and consider the impact of systematic uncertainties on our spectra.

Swirls are ubiquitous in the solar atmosphere. They are believed to be related to the excitation of different modes of magnetohydrodynamic waves and pulses, as well as spicules. However, statistical studies of their collective behaviour are rare. In this paper, we aim to study the collective, as well as the behaviour of individual photospheric and chromospheric swirls detected by the automated swirl detection algorithm (ASDA) from observations obtained by the Swedish 1-m Solar Telescope and the Hinode satellite. Detailed analysis of six different parameters of photospheric and chromospheric swirls is performed employing the wavelet analysis. Two clusters of periods with significant wavelet power, one from $3-8$ minutes and the other from $10-14$ minutes, have been found. The former coincides with the dominant period of the global $p$-mode spectrum. Wavelet and Fast Fourier Transform (FFT) analysis of example swirls also reveals similar periods. These results suggest that global $p$-modes might be important for triggering photospheric and thus chromospheric swirls. A novel scenario of global $p$-modes providing energy and mass fluxes to the upper solar atmosphere via generating swirls, Alfv\'en pulses and spicules is then proposed.

We use the Athena++ Monte Carlo (MC) radiation transfer module to post-process simulation snapshots from non-relativistic Athena++ radiation magnetohydrodynamic (RMHD) simulations. These simulations were run using a gray (frequency integrated) approach but were also restarted and ran with a multi-group approach that accounts for Compton scattering with a Kompaneets operator. These simulations produced moderately super-Eddington accretion rates onto a 6.62 $M_\odot$ black hole. Since we only achieve inflow equilibrium out to 20-25 gravitational radii, we focus on the hard X-ray emission. We provide a comparison between the MC and RMHD simulations showing that the treatment of Compton scattering in the gray RMHD simulations underestimates the gas temperature in the regions above and below the accretion disk. In contrast, the restarted multi-group snapshots provides a treatment for the radiation field that is more consistent with the MC calculations, and result in post-processed spectra with harder X-ray emission compared to their gray snapshot counterparts. We characterize these MC post-processed spectra using commonly employed phenomenological models used for spectral fitting. We also attempt to fit our MC spectra directly to observations of the ultraluminous X-ray source (ULX) NGC 1313 X-1, finding best fit values that are competitive to phenomenological model fits, indicating that first principle models of super-Eddington accretion may adequately explain the observed hard X-ray spectra in some ULX sources.

We create mock X-ray observations of hot gas in galaxy clusters with a new extension of L-Galaxies semi-analytic model of galaxy formation, which includes the radial distribution of hot gas in each halo. Based on the model outputs, we first build some mock light cones, then generate mock spectra with SOXS package and derive the mock images in the light cones. Using the mock data, we simulate the mock X-ray spectra for ROSAT all-sky survey, and compare the mock spectra with the observational results. Then, we consider the design parameters of HUBS mission and simulate the observation of the halo hot gas for HUBS as an important application of our mock work. We find: (1) Our mock data match the observations by current X-ray telescopes. (2) The survey of hot baryons in resolved clusters by HUBS is effective below redshift 0.5, and the observations of the emission lines in point-like sources at z>0.5 by HUBS help us understand the hot baryons in the early universe. (3) By taking the advantage of the large simulation box and flexibility in semi-analytic models, our mock X-ray observations provide the opportunity to make target selection and observation strategies for forthcoming X-ray facilities.

The quasi-periodic oscillations (QPOs) found in active galactic nuclei (AGNs) are a very interesting observational phenomenon implying an unknown physical mechanism around supermassive black holes. Several AGNs have been found to have QPO phenomena in the X-ray energy band. Long-duration X-ray observations were collected and reduced for six AGNs with a suspected QPO. The Gaussian process (GP) model celerite was used to fit the light curves and to search for the quasi-periodicity behavior. The power spectral density and parameter posterior distributions of each light curve were calculated with the optimal model. Of the six AGNs, only RE J1034+396 was found to have a QPO of about 3600 s. The other five sources do not show QPO modulation behavior. We propose that a hot spot on the accretion disk is a possible physical mechanism resulting in this quasi-periodic behavior of AGNs.

In this work, we present Stellar Spectra Factory (SSF), a tool to generate empirical-based stellar spectra from arbitrary stellar atmospheric parameters. The relative flux-calibrated empirical spectra can be predicted by SSF given arbitrary effective temperature, surface gravity, and metallicity. SSF constructs the interpolation approach based on the SLAM, using ATLAS-A library as the training dataset. SSF is composed of 4 data-driven sub-models to predict empirical stellar spectra. SSF-N can generate spectra from A to K type and some M giant stars, covering 3700 < Teff < 8700 K, 0 < logg < 6 dex, and -1.5 < [M/H] < 0.5 dex. SSF-gM is mainly used to predict M giant spectra with 3520 < Teff < 4000K and -1.5 < [M/H] < 0.4 dex. SSF-dM is for generating M dwarf spectra with 3295 < Teff < 4040K, -1.0 < [M/H] < 0.1 dex. And SSF-B can predict B-type spectra with 9000 < Teff < 24000K and -5.2< MG < 1.5 mag. The accuracy of the predicted spectra is validated by comparing the flux of predicted spectra to those with same stellar parameters selected from the known spectral libraries, MILES and MaStar. The averaged difference of flux over optical wavelength between the predicted spectra and the corresponding ones in MILES and MaStar is less than 5%. More verification is conducted between the magnitudes calculated from the integration of the predicted spectra and the observations in PS1 and APASS bands with the same stellar parameters. No significant systematic difference is found between the predicted spectra and the photomatric observations. The uncertainty is 0.08mag in r band for SSF-gM when comparing with the stars with the same stellar parameters selected from PS1. And the uncertainty becomes 0.31mag in i band for SSF-dM when comparing with the stars with the same stellar parameters selected from APASS.

GREX-PLUS (Galaxy Reionization EXplorer and PLanetary Universe Spectrometer) is a mission candidate for a JAXA's strategic L-class mission to be launched in the 2030s. Its primary sciences are two-fold: galaxy formation and evolution and planetary system formation and evolution. The GREX-PLUS spacecraft will carry a 1.2 m primary mirror aperture telescope cooled down to 50 K. The two science instruments will be onboard: a wide-field camera in the 2-8 $\mu$m wavelength band and a high resolution spectrometer with a wavelength resolution of 30,000 in the 10-18 $\mu$m band. The GREX-PLUS wide-field camera aims to detect the first generation of galaxies at redshift $z>15$. The GREX-PLUS high resolution spectrometer aims to identify the location of the water ``snow line'' in proto-planetary disks. Both instruments will provide unique data sets for a broad range of scientific topics including galaxy mass assembly, origin of supermassive blackholes, infrared background radiation, molecular spectroscopy in the interstellar medium, transit spectroscopy for exoplanet atmosphere, planetary atmosphere in the Solar system, and so on.

(abridged) Far-infrared (FIR) line emission provides key information about the gas cooling and heating due to shocks and UV radiation associated with the early stages of star formation. Gas cooling via FIR lines might, however, depend on metallicity. We aim to quantify the FIR line emission and determine the spatial distribution of the CO rotational temperature, ultraviolet (UV) radiation field, and H2 number density toward the embedded cluster Gy 3-7 in the CMa-l224 star-forming region, whose metallicity is expected to be intermediate between that of the LMC and the Solar neighborhood. By comparing the total luminosities of CO and [O I] toward Gy 3-7 with values found for low- and high-mass protostars extending over a broad range of metallicities, we also aim to identify the possible effects of metallicity on the FIR line cooling within our Galaxy. We studied SOFIA/FIFI-LS spectra of Gy 3-7 covering several FIR lines. The spatial extent of CO high-J (J>14) emission resembles that of the elongated 160 um continuum emission detected with Herschel. The CO transitions from J=14-13 to J=16-15 are detected throughout the cluster and show a median rotational temperature of 170+/-30 K on Boltzmann diagrams. Comparisons to other protostars observed with Herschel show a good agreement with intermediate-mass sources in the inner Galaxy. Assuming an origin of the [O I] and high-J CO emission in UV-irradiated C-shocks, we obtained pre-shock H2 number densities of 10^4-5 cm-3 and UV radiation field strengths of 0.1-10 Habing fields. Far-IR line observations reveal ongoing star formation in Gy 3-7, dominated by intermediate-mass Class 0/I young stellar objects. The ratio of molecular-to-atomic far-IR line emission shows a decreasing trend with bolometric luminosities of the protostars. However, it does not indicate that the low-metallicity has an impact on the line cooling in Gy 3-7.

Recently, the characterisation of binary systems of neutron stars has become central in various fields such as gravitational waves, gamma-ray bursts (GRBs), and the chemical evolution of galaxies. In this work, we explore possible observational proxies that can be used to infer some characteristics of the delay time distribution (DTD) of neutron star mergers (NSMs). We construct a sample of model galaxies that fulfils the observed galaxy stellar mass function, star formation rate versus mass relation, and the cosmic star formation rate density. The star formation history of galaxies is described with a log-normal function characterised by two parameters: the position of the maximum and the width of the distribution. We assume a theoretical DTD that mainly depends on the lower limit and the slope of the distribution of the separations of the binary neutron stars systems at birth. We find that the current rate of NSMs ($\mathcal{R}=320^{+490}_{-240}$ Gpc$^{-3}$yr$^{-1}$) requires that $\sim0.3$ per cent of neutron star progenitors lives in binary systems with the right characteristics to lead to a NSM within a Hubble time. We explore the expected relations between the rate of NSMs and the properties of the host galaxy. We find that the most effective proxy for the shape of the DTD of NSMs is the current star formation activity of the typical host. At present, the fraction of short-GRBs observed in star-forming galaxies favours DTDs with at least $\sim40\%$ of mergers within $100$ Myr. This conclusion will be put on a stronger basis with larger samples of short-GRBs with host association (e.g. $600$ events at $z \leq 1$)

We study the acceleration of charged particles by ultra-relativistic shocks using test-particle Monte-Carlo simulations. Two field configurations are considered: (i) shocks with uniform upstream magnetic field in the plane of the shock, and (ii) shocks in which the upstream magnetic field has a cylindrical geometry. Particles are assumed to diffuse in angle due to frequent non-resonant scattering on small-scale fields. The steady-state distribution of particles' Lorentz factors is shown to approximately satisfy $dN/d\gamma \propto \gamma^{-2.2}$ provided the particle motion is scattering dominated on at least one side of the shock. For scattering dominated transport, the acceleration rate scales as $t_{\rm acc}\propto t^{1/2}$, though recovers Bohm scaling $t_{\rm acc}\propto t$ if particles become magnetised on one side of the shock. For uniform field configurations, a limiting energy is reached when particles are magnetised on both sides of the shock. For the cylindrical field configuration, this limit does not apply, and particles of one sign of charge will experience a curvature drift that redirects particles upstream. For the non-resonant scattering model considered, these particles preferentially escape only when they reach the confinement limit determined by the finite system size, and the distribution approaches the escapeless limit $dN/d\gamma \propto \gamma^{-1}$. The cylindrical field configuration resembles that expected for jets launched by the Blandford $\&$ Znajek mechanism, the luminous jets of AGN and GRBs thus provide favourable sites for the production of ultra-high energy cosmic rays.

The Dark Energy Camera (DECam) is a sensitive, wide field instrument mounted at the prime focus of the 4 m V. Blanco Telescope in Chile. Beside its main objectives, i.e. understanding the growth and evolution of structures in the Universe, the camera offers the opportunity to observe a 3 deg2 field of view in one single pointing and, with an adequate cadence, to identify the variable sources contained. In this paper, we present the result of a DECam observational campaign toward the LMC and give a catalogue of the observed variable sources. We considered all the available DECam observations of the LMC, acquired during 32 nights over a period of two years (from February 2018 to January 2020), and set up a specific pipeline for detecting and characterizing variable sources in the observed fields. Here, we report on the first 15 deg2 in and around the LMC as observed by DECam, testing the capabilities of our pipeline. Since many of the observed fields cover a rather crowded region of the sky, we adopted the ISIS subtraction package which, even in these conditions, can detect variables at a very low signal to noise ratio. All the potentially identified variable sources were then analyzed and each light curve tested for periodicity by using the Lomb-Scargle and Schwarzenberg-Czerny algorithms. Furthermore, we classified the identified sources by using the UPSILoN neural network. This analysis allowed us to find 70 981 variable stars, 1266 of which were previously unknown. We estimated the period of the variables and compared it with the available values in the catalogues. Moreover, for the 1266 newly detected objects, an attempted classification based on light curve analysis is presented.

We have used the ``dynamical clock'' to measure the level of dynamical evolution reached by three Galactic globular clusters (namely, NGC 3201, NGC 6316 and NGC 6440). This is an empirical method that quantifies the level of central segregation of blue stragglers stars (BSSs) within the cluster half-mass radius by means of the $A^+_{rh}$ parameter, defined as the area enclosed between the cumulative radial distribution of BSSs and that of a lighter population. The total sample with homogeneous determinations of $A^+_{rh}$ now counts a gran-total of 59 clusters: 52 old GCs in the Milky Way (including the three investigated here), 5 old clusters in the Large Magellanic Cloud, and 2 young systems in the Small Magellanic Cloud. The three objects studied here nicely nest into the correlation between $A^+_{rh}$ and the central relaxation time defined by the previous sample, thus proving and consolidating the use of the dynamical clock as an excellent tracer of the stage of star cluster dynamical evolution in different galactic environments. Finally, we discuss the advantages of using the dynamical clock as an indicator of star cluster dynamical ages, compared to the present-day central relaxation time.

We investigate a scenario where primordial black holes (PBHs) can be the progenitors of supermassive black holes (SMBHs) observed at $z\sim6$. To this end, we carried out clustering analysis using a sample of 81 quasars at $5.88 <z<6.49$, which is constructed in Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, and 11 quasars in the same redshift range selected from the literature. The resulting angular auto-correlation function (ACF) can be fitted to a power-law form of $\omega_\theta = 0.045^{+0.114}_{-0.106}~\theta^{-0.8}$ over a scale of $0.2\!-\!10$ degrees. We compare the ACF of the quasars to that predicted for the PBH model at $z\sim 6$ and found that such a scenario is excluded for a broad range of parameter space, from which we can conclude that a scenario with PBHs as SMBHs is not viable. We also discuss a model in which SMBHs at $z \sim 6$ originate from the direct collapse of PBH clumps and argue that the observed ACF excludes such a scenario in the context of our PBH model.

We present a novel approach for estimating cosmological parameters, $\Omega_m$, $\sigma_8$, $w_0$, and one derived parameter, $S_8$, from 3D lightcone data of dark matter halos in redshift space covering a sky area of $40^\circ \times 40^\circ$ and redshift range of $0.3 < z < 0.8$, binned to $64^3$ voxels. Using two deep learning algorithms, Convolutional Neural Network (CNN) and Vision Transformer (ViT), we compare their performance with the standard two-point correlation (2pcf) function. Our results indicate that CNN yields the best performance, while ViT also demonstrates significant potential in predicting cosmological parameters. By combining the outcomes of Vision Transformer, Convolution Neural Network, and 2pcf, we achieved a substantial reduction in error compared to the 2pcf alone. To better understand the inner workings of the machine learning algorithms, we employed the Grad-CAM method to investigate the sources of essential information in activation maps of the CNN and ViT. Our findings suggest that the algorithms focus on different parts of the density field and redshift depending on which parameter they are predicting. This proof-of-concept work paves the way for incorporating deep learning methods to estimate cosmological parameters from large-scale structures, potentially leading to tighter constraints and improved understanding of the Universe.

M31 and M33 are the closest spiral galaxies and the largest members (together with the Milky Way) of the Local group, which makes them interesting targets for indirect dark matter searches. In this paper we present studies of the expected sensitivity of the Cherenkov Telescope Array (CTA) to an annihilation signal from weakly interacting massive particles from M31 and M33. We show that a 100 h long observation campaign will allow CTA to probe annihilation cross-sections up to $\langle\sigma\upsilon\rangle\approx 5\cdot10^{-25}~\mathrm{cm^{3}s^{-1}}$ for the $\tau^{+}\tau^{-}$ annihilation channel (for M31, at a DM mass of 0.3 TeV), improving the current limits derived by HAWC by up to an order of magnitude. We present an estimate of the expected CTA sensitivity, by also taking into account the contributions of the astrophysical background and other possible sources of systematic uncertainty. We also show that CTA might be able to detect the extended emission from the bulge of M31, detected at lower energies by the Fermi/LAT.

When gravitational waves (GWs) from a spinning neutron star arrive from behind the Sun, they are subjected to gravitational lensing that imprints a frequency-dependent modulation on the waveform. This modulation traces the projected solar density and gravitational potential along the path as the Sun passes in front of the neutron star. We calculate how accurately the solar density profile can be extracted from the lensed GWs using a Fisher analysis. For this purpose, we selected three promising candidates (the highly spinning pulsars J1022+1001, J1730-2304, and J1745-23) from the pulsar catalog of the Australia Telescope National Facility. The lensing signature can be measured with $3 \sigma$ confidence when the signal-to-noise ratio (SNR) of the GW detection reaches $100 \, (f/300 {\rm Hz})^{-1}$ over a one-year observation period (where $f$ is the GW frequency). The solar density profile can be plotted as a function of radius when the SNR improves to $\gtrsim 10^4$.

To tackle the still unsolved and fundamental problem of the role of Active Galactic Nuclei (AGN) feedback in shaping galaxies, in this work we implement a new physical treatment of AGN-driven winds into our semi-analytic model of galaxy formation. To each galaxy in our model, we associate solutions for the outflow expansion and the mass outflow rates in different directions, depending on the AGN luminosity, on the circular velocity of the host halo, and on gas content of the considered galaxy. To each galaxy we also assign an effective radius derived from energy conservation during merger events, and a stellar velocity dispersion self-consistently computed via Jeans modelling. We derive all the main scaling relations between Black hole (BH) mass and total/bulge stellar mass, velocity dispersion, host halo dark matter mass, and star formation efficiency. We find that our improved AGN feedback mostly controls the dispersion around the relations but plays a subdominant role in shaping slopes and/or normalizations of the scaling relations. Including possible limited-resolution selection biases in the model provides better agreement with the available data. The model does not point to any more fundamental galactic property linked to BH mass, with velocity dispersion playing a similar role with respect to stellar mass, in tension with present data. In line with other independent studies carried out on comprehensive semi-analytic and hydrodynamic galaxy-BH evolution models, our current results signal either an inadequacy of present cosmological models of galaxy formation in fully reproducing the local scaling relations, in terms of both shape and residuals, and/or point to an incompleteness issue affecting the local sample of dynamically-measured BHs.

Over the last few years, many studies have found an empirical relationship between the abundance of a star and its age. Here we estimate spectroscopic stellar ages for 178 825 red-giant stars observed by the APOGEE survey with a median statistical uncertainty of 17%. To this end, we use the supervised machine learning technique XGBoost, trained on a high-quality dataset of 3 060 red-giant and red-clump stars with asteroseismic ages observed by both APOGEE and Kepler. After verifying the obtained age estimates with independent catalogues, we investigate some of the classical chemical, positional, and kinematic relationships of the stars as a function of their age. We find a very clear imprint of the outer-disc flare in the age maps and confirm the recently found split in the local age-metallicity relation. We present new and precise measurements of the Galactic radial metallicity gradient in small age bins between 0.5 and 12 Gyr, confirming a steeper metallicity gradient for 2-5 Gyr old populations and a subsequent flattening for older populations mostly produced by radial migration. In addition, we analyse the dispersion about the abundance gradient as a function of age. We find a clear power-law trend (with an exponent $\beta\approx0.15$) for this relation, indicating a smooth radial migration history in the Galactic disc over the past 7-9 Gyr. Departures from this power law are detected at ages of 8 Gyr (possibly related to the Gaia Sausage/Enceladus merger) and 2.75 Gyr (possibly related to an enhancement of the star-formation rate in the Galactic disc). Finally, we confirm previous measurements showing a steepening in the age-velocity dispersion relation at around 9 Gyr, but now extending it over a large extent of the Galactic disc (5 kpc < RGal < 13 kpc). [Abridged]

Although double-peaked narrow emission-line galaxies have been studied extensively in the past years, only a few are reported with the green pea galaxies (GPs). Here we present our discovery of five GPs with double-peaked narrow [OIII] emission lines, referred to as DPGPs, selected from the LAMOST and SDSS spectroscopic surveys. We find that these five DPGPs have blueshifted narrow components more prominent than the redshifted components, with velocity offsets of [OIII]$\lambda$5007 lines ranging from 306 to 518 $\rm km\, s^{-1}$ and full widths at half maximums (FWHMs) of individual components ranging from 263 to 441 $\rm km\, s^{-1}$. By analyzing the spectra and the spectral energy distributions (SEDs), we find that they have larger metallicities and stellar masses compared with other GPs. The H$\alpha$ line width, emission-line diagnostic, mid-infrared color, radio emission, and SED fitting provide evidence of the AGN activities in these DPGPs. They have the same spectral properties of Type 2 quasars. Furthermore, we discuss the possible nature of the double-peaked narrow emission-line profiles of these DPGPs and find that they are more likely to be dual AGN. These DPGP galaxies are ideal laboratories for exploring the growth mode of AGN in the extremely luminous emission-line galaxies, the co-evolution between AGN and host galaxies, and the evolution of high-redshift galaxies in the early Universe.

We present a long-term overview of the atmospheric phase stability at the Atacama Large Millimeter/submillimeter Array (ALMA) site, using >5 years of data, that acts as the successor to the studies summarized two decades ago by Evans et al 2003. Importantly, we explore the atmospheric variations, the `phase RMS', and associated metadata of over 17000 accrued ALMA observations taken since Cycle 3 (2015) by using the Bandpass calibrator source scans. We indicate the temporal phase RMS trends for average baseline lengths of 500, 1000, 5000, and 10000m, in contrast to the old stability studies that used a single 300m baseline phase monitor system. At the ALMA site, on the Chajnantor plateau, we report the diurnal variations and monthly changes in the phase RMS on ALMA relevant baselines lengths, measured directly from data, and we reaffirm such trends in atmospheric transmission (via Precipitable Water Vapour - PWV). We confirm that day observations have respectively higher phase RMS and PWV in contrast to night, while the monthly variations show Chilean winter (June - August) providing the best, high-frequency and long-baseline observing conditions - low (stable) phase RMS and low PWV. Yet, not all good phase stability condition occur when the PWV is low. Measurements of the phase RMS as a function of short timescales, 30 to 240s, that tie with typical target source scan times, and as a function of baseline length indicate that phase variations are smaller for short timescales and baselines and larger for longer timescales and baselines. We illustrate that fast-switching phase-referencing techniques, that allow short target scan times, could work well in reducing the phase RMS to suitable levels specifically for high-frequencies (Band 8, 9 and 10), long-baselines, and the two combined.

We present a gas selection for Xe-based gas electron multiplier (GEM) detectors, Gas Multiplier Counters (GMCs) onboard the CubeSat X-ray observatory NinjaSat. To achieve an energy bandpass of 2-50 keV, we decided to use a Xe-based gas mixture at a pressure of 1.2 atm that is sensitive to high-energy X-rays. In addition, an effective gain of over 300 is required for a single GEM so that the 2 keV X-ray signal can be sufficiently larger than the electrical noise. At first, we measured the effective gains of GEM in nine Xe-based gas mixtures (combinations of Xe, Ar, CO2, CH4, and dimethyl ether; DME) at 1.0 atm. The highest gains were obtained with Xe/Ar/DME mixtures, while relatively lower gains were obtained with Xe/Ar/CO2, Xe/Ar/CH4, and Xe+quencher mixtures. Based on these results, we selected the Xe/Ar/DME (75%/24%/1%) mixture at 1.2 atm as the sealed gas for GMC. Then we investigated the dependence of an effective gain on the electric fields in the drift and induction gaps ranging from 100-650 V cm$^{-1}$ and 500-5000 V cm$^{-1}$, respectively, in the selected gas mixture. The effective gain weakly depended on the drift field while it was almost linearly proportional to the induction field: 2.4 times higher at 5000 V cm$^{-1}$ than at 1000 V cm$^{-1}$. With the optimal induction and drift fields, the flight model GMC achieves an effective gain of 460 with an applied GEM voltage of 590 V.

A source lying near hyperbolic umbilic (HU) leads to a ring-like image formation, constituting four images with high magnification factors and lying in a small region of the lens plane. Since (based on our earlier work) the observed number of HU image formations in cluster lenses is expected to increase in future, it is timely to investigate them in more detail. Like fold and cusp, HU also satisfies the magnification relation, i.e., the signed magnification sum of the four images equals zero. This work presents a detailed study of HU magnification relation ($R_{\rm hu}$) considering the elliptical Navarro-Frenk-White (eNFW) lens profile suitable for cluster scale dark matter halos. Our results show that for an isolated eNFW lens, $R_{\rm hu}$ is more sensitive to ellipticity than its mass or concentration parameter. An ellipticity greater than 0.3 results in $R_{\rm hu}$ lying close to zero with a small scatter around it. A substructure near the HU image formation causes the average $R_{\rm hu}$ value to deviate from zero and increases the scatter, with the amount of deviation depending on the image type near which the substructure lies. However, a population of substructures in the lens plane (equivalent to the galaxy lenses inside the cluster) does not significantly shift the average $R_{\rm hu}$ value from zero but increases the scatter around it. We find that $R_{\rm hu} \simeq 0$ for HU image formation in the Abell 1703 cluster. Repeating this test in other clusters where HU formations are discovered can be a useful indicator of substructure in cluster halos.

OB associations are important probes of recent star formation and Galactic structure. In this study, we focus on the Auriga constellation, an important region of star formation due to its numerous young stars, star-forming regions and open clusters. We show using \textit{Gaia} data that its two previously documented OB associations, Aur OB1 and OB2, are too extended in proper motion and distance to be genuine associations, encouraging us to revisit the census of OB associations in Auriga with modern techniques. We identify 5617 candidate OB stars across the region using photometry, astrometry and our SED fitting code, grouping these into 5 high-confidence OB associations using HDBSCAN. Three of these are replacements to the historical pair of associations - Aur OB2 is divided between a foreground and a background association - while the other two associations are completely new. We connect these OB associations to the surrounding open clusters and star-forming regions, analyse them physically and kinematically, constraining their ages through a combination of 3D kinematic traceback, the position of their members in the HR diagram and their connection to clusters of known age. Four of these OB associations are expanding, with kinematic ages up to a few tens of Myr. Finally, we identify an age gradient in the region spanning several associations that coincides with the motion of the Perseus spiral arm over the last $\sim$20 Myr across the field of view.

Some hydrogen-poor supernovae (SNe) are found to undergo interaction with dense circumstellar matter (CSM) that may originate from mass eruption(s) just prior to core-collapse. We model the interaction between the remaining star and the bound part of the erupted CSM that eventually fall back to the star. We find that while fallback initially results in a continuous CSM down to the star, feedback processes from the star can push the CSM to large radii of $\gtrsim 10^{15}$ cm from several years after the eruption. In the latter case, a tenuous bubble surrounded by a dense and detached CSM extending to $\gtrsim 10^{16}$ cm is expected. Our model offers a natural unifying explanation for the diverse CSM structures seen in hydrogen-poor SNe, such as Type Ibn/Icn SNe that show CSM signatures soon after explosion, and the recently discovered Type Ic SNe 2021ocs and 2022xxf ("the Bactrian") with CSM signatures seen only at late times.

We present mid-infrared spectroscopy of Jupiter's mid-to-high latitudes using Gemini-North/TEXES (Texas Echelon Cross Echelle Spectrograph) on March 17-19, 2017. These observations capture Jupiter's hydrocarbon auroral emissions before, during and after the arrival of a solar wind compression on March 18th, which highlights the coupling between the polar stratosphere and external space environment. In comparing observations on March 17th and 19th, we observe a brightening of the CH$_4$, C$_2$H$_2$ and C$_2$H$_4$ emissions in regions spatially coincident with the northern, duskside main auroral emission (henceforth, MAE). In inverting the spectra to derive atmospheric information, we determine that the duskside brightening results from an upper stratospheric (p < 0.1 mbar/z > 200 km) heating (e.g. $\Delta T$ = 9.1 $\pm$ 2.1 K at 9 $\mu$bar at 67.5$^\circ$N, 162.5$^\circ$W) with negligible heating at deeper pressures. Our interpretation is that the arrival of the solar wind enhancement drove magnetospheric dynamics through compression and/or viscous interactions on the flank. These dynamics accelerated currents and/or generated higher Poynting fluxes, which ultimately warmed the atmosphere through Joule heating and ion-neutral collisions. Poleward of the southern MAE, temperature retrievals demonstrate that auroral-related heating penetrates as deep as the 10-mbar level, in contrast to poleward of the northern MAE, where heating is only observed as deep as $\sim$3 mbar. We suggest this results from the south having higher Pedersen conductivities, and therefore stronger currents and acceleration of the neutrals, as well as the poleward heating overlapping with the apex of Jupiter's circulation thereby inhibiting efficient horizontal mixing/advection.

Context. The Optical Gravitational Lensing Experiment (OGLE) observed around 450,000 eclipsing binaries (EBs) towards the Galactic Bulge. Decade-long photometric observations such as these provide an exceptional opportunity to thoroughly examine the targets. However, observing dense stellar fields such as the Bulge may result in blends and contamination by close objects. Aims. We searched for periodic variations in the residual light curves of EBs in OGLE-IV and created a new catalogue for the EBs that contain `background' signals after the investigation of the source of the signal. Methods. From the about half a million EB systems, we selected those that contain more than 4000 data points. We fitted the EB signal with a simple model and subtracted it. To identify periodical signals in the residuals, we used a GPU-based phase dispersion minimisation python algorithm called cuvarbase and validated the found periods with Lomb-Scargle periodograms. We tested the reliability of our method with artificial light curves. Results. We identified 354 systems where short-period background variation was significant. In these cases, we determined whether it is a new variable or just the result of contamination by an already catalogued nearby one. We classified 292 newly found variables into EB, $\delta$ Scuti, or RR Lyrae categories, or their sub-classes, and collected them in a catalogue. We also discovered four new doubly eclipsing systems and one eclipsing multiple system with a $\delta$ Scuti variable, and modelled the outer orbits of the components.

Aims. We introduce a new deep learning tool that estimates stellar parameters (such as effective temperature, surface gravity, and extinction) of young low-mass stars by coupling the Phoenix stellar atmosphere model with a conditional invertible neural network (cINN). Our networks allow us to infer the posterior distribution of each stellar parameter from the optical spectrum. Methods. We discuss cINNs trained on three different Phoenix grids: Settl, NextGen, and Dusty. We evaluate the performance of these cINNs on unlearned Phoenix synthetic spectra and on the spectra of 36 Class III template stars with well-characterised stellar parameters. Results. We confirm that the cINNs estimate the considered stellar parameters almost perfectly when tested on unlearned Phoenix synthetic spectra. Applying our networks to Class III stars, we find good agreement with deviations of at most 5--10 per cent. The cINNs perform slightly better for earlier-type stars than for later-type stars like late M-type stars, but we conclude that estimations of effective temperature and surface gravity are reliable for all spectral types within the network's training range. Conclusions. Our networks are time-efficient tools applicable to large amounts of observations. Among the three networks, we recommend using the cINN trained on the Settl library (Settl-Net), as it provides the best performance across the largest range of temperature and gravity.

Observations collected with the Atacama Large Millimeter/submillimeter Array (ALMA) and analysis of broadband X-ray spectra have recently suggested the presence of a central compact object (CCO) in SN 1987A. However, no direct evidence of the CCO has been found yet. Here we analyze Chandra X-ray observations of SN 1987A collected in 2007 and 2018, and synthesize the 2027 Chandra and 2037 Lynx spectra of the faint inner region of SN 1987A. We estimate the temporal evolution of the upper limits of the intrinsic luminosity of the putative CCO in three epochs (2018, 2027 and 2037). We find that these upper limits are higher for higher neutron star (NS) kick velocities due to the increased absorption from the surrounding cold ejecta. We compare NS cooling models with both the intrinsic luminosity limits obtained from the X-ray spectra, and the ALMA constraints with the assumption that the observed blob of SN 1987A is primarily heated by thermal emission. We find that the synthetic Lynx spectra are crucial to constrain physical properties of the CCO, which will be confirmed by future observations in the 2040s. We draw our conclusions based on two scenarios, namely the non-detection and detection of NS by Lynx. If the NS is not detected, its kick velocity should be ~700 km/s. Furthermore, the non-detection of the NS would suggest rapid cooling processes around the age of 40 years, implying strong crust superfluidity. Conversely, in the case of NS detection, the mass of the NS envelope must be high.

We present the first detection of the baryon acoustic oscillations (BAO) signal obtained using unblinded data collected during the initial two months of operations of the Stage-IV ground-based Dark Energy Spectroscopic Instrument (DESI). From a selected sample of 261,291 Luminous Red Galaxies spanning the redshift interval 0.4 < z < 1.1 and covering 1651 square degrees with a 57.9% completeness level, we report a ~5 sigma level BAO detection and the measurement of the BAO location at a precision of 1.7%. Using a Bright Galaxy Sample of 109,523 galaxies in the redshift range 0.1 < z < 0.5, over 3677 square degrees with a 50.0% completeness, we also detect the BAO feature at ~3 sigma significance with a 2.6% precision. These first BAO measurements represent an important milestone, acting as a quality control on the optimal performance of the complex robotically-actuated, fiber-fed DESI spectrograph, as well as an early validation of the DESI spectroscopic pipeline and data management system. Based on these first promising results, we forecast that DESI is on target to achieve a high-significance BAO detection at sub-percent precision with the completed 5-year survey data, meeting the top-level science requirements on BAO measurements. This exquisite level of precision will set new standards in cosmology and confirm DESI as the most competitive BAO experiment for the remainder of this decade.

Line intensity mapping (LIM) can provide a powerful means to constrain the theory of gravity and the nature of dark energy at low and high redshifts by mapping the large-scale structure (LSS) over many redshift epochs. In this paper, we investigate the potential of the next generation ground-based millimeter-wavelength LIM surveys in constraining several models beyond $\Lambda$CDM, involving either a dynamic dark energy component or modifications of the theory of gravity. Limiting ourselves to two-point clustering statistics, we consider the measurements of auto-spectra of several CO rotational lines (from J=2-1 to J=6-5) and the [CII] fine structure line in the redshift range of $0.25<z<12$. We consider different models beyond $\Lambda$CDM, each one with different signatures and peculiarities. Among them, we focus on Jordan-Brans-Dicke and axion-driven early dark energy models as examples of well-studied scalar-tensor theories acting at late and early times respectively. Additionally, we consider three phenomenological models based on an effective description of gravity at cosmological scales. We show that LIM surveys deployable within a decade (with $\sim 10^8$ spectrometer hours) have the potential to improve upon the current bounds on all considered models significantly. The level of improvements range from a factor of a few to an order of magnitude.

The detection of terrestrial planets by radial velocity and photometry is hindered by the presence of stellar signals. Those are often modeled as stationary Gaussian processes, whose kernels are based on qualitative considerations, which do not fully leverage the existing physical understanding of stars. Our aim is to build a formalism which allows to transfer the knowledge of stellar activity into practical data analysis methods. In particular, we aim at obtaining kernels with physical parameters. This has two purposes: better modelling signals of stellar origin to find smaller exoplanets, and extracting information about the star from the statistical properties of the data. We consider several observational channels such as photometry, radial velocity, activity indicators, and build a model called FENRIR to represent their stochastic variations due to stellar surface inhomogeneities. We compute analytically the covariance of this multi-channel stochastic process, and implement it in the S+LEAF framework to reduce the cost of likelihood evaluations from $O(N^3)$ to $O(N)$. We also compute analytically higher order cumulants of our FENRIR model, which quantify its non-Gaussianity. We obtain a fast Gaussian process framework with physical parameters, which we apply to the HARPS-N and SORCE observations of the Sun, and constrain a solar inclination compatible with the viewing geometry. We then discuss the application of our formalism to granulation. We exhibit non-Gaussianity in solar HARPS radial velocities, and argue that information is lost when stellar activity signals are assumed to be Gaussian. We finally discuss the origin of phase shifts between RVs and indicators, and how to build relevant activity indicators. We provide an open-source implementation of the FENRIR Gaussian process model with a Python interface.

In order to predict the cosmological abundance of dark matter, an estimation of particle rates in an expanding thermal environment is needed. For thermal dark matter, the non-relativistic regime sets the stage for the freeze-out of the dark matter energy density. We compute transition widths and annihilation, bound-state formation, and dissociation cross sections of dark matter fermion pairs in the unifying framework of non-relativistic effective field theories at finite temperature, with the thermal bath modeling the thermodynamical behaviour of the early universe. We reproduce and extend some known results for the paradigmatic case of a dark fermion species coupled to dark gauge bosons. The effective field theory framework allows to highlight their range of validity and consistency, and to identify some possible improvements.

The annihilation of TeV-scale leptophilic dark matter into electron-positron pairs (hereafter $e^+e^-$) will produce a sharp cutoff in the local cosmic-ray $e^+e^-$ spectrum at an energy matching the dark matter mass. At these high energies, $e^+e^-$ cool quickly due to synchrotron interactions with magnetic fields and inverse-Compton scattering with the interstellar radiation field. These energy losses are typically modelled as a continuous process. However, inverse-Compton scattering is a stochastic energy-loss process where interactions are rare but catastrophic. We show that when inverse-Compton scattering is modelled as a stochastic process, the expected $e^+e^-$ flux from dark matter annihilation is about a factor of $\sim$2 larger near the dark matter mass than in the continuous model. This greatly enhances the detectability of heavy dark matter.

The freeze-in mechanism has been shown to allow the simultaneous generation of cosmic dark matter and a viable matter-antimatter asymmetry in the universe. When the underlying interactions are described by higher-dimensional, non-renormalizable operators, the relevant freeze-in processes take place close to the highest considered cosmic temperatures. In this paper we study how the presence of a fluid that temporarily dominates the energy content of the early universe affects the predictions of this ``Ultraviolet Freeze-In Baryogenesis'' scenario. We find that this additional cosmic component has a significant impact on the predictions of concrete microscopic models, allowing for reheating temperatures which are much lower than those required in the simplest cosmological scenario. Moreover, we show that inflationary observables can constrain the parameter space of such models, once the latter are examined in conjunction with concrete models of inflation.

By directly inverting the observational data of several neutron star observables in the three dimensional parameter space of the constant speed of sound (CSS) model while fixing all hadronic Equation of State parameters at their currently known most probable values, we constrain the three parameters of the CSS model and their correlations. Using two lower radius limits of $R_{2.01}=11.41$ km and $R_{2.01}=12.2$ km for PSR J0740+6620 obtained from two independent analyses using different approaches by the Neutron Star Interior Composition Explorer (NICER) Collaboration, the speed of sound squared $c_{\rm QM}^2$ in quark matter is found to have a lower limit of $0.35$ and $0.43$ in unit of $c^2$, respectively, above its conformal limit of $c_{\rm QM}^2<1/3$. Moreover, an approximately linear correlation between the first-order hadron-quark transition density $\rho_t$ and its strength $\Delta\varepsilon$ is found.

We investigate the mass range and the corresponding free-streaming length scale of dark matter produced non-thermally from decay of heavy objects which can be either dominant or sub-dominant at the moment of decay. We show that the resulting dark matter could be very light well below keV scale with a free-streaming length satisfying the Lyman-{\alpha} constraints. We demonstrate two explicit examples for such light cold dark matter.

We study the cross-correlation between the stochastic gravitational-wave background (SGWB) generated by binary black hole (BBH) mergers across the universe and the distribution of galaxies across the sky. We use the anisotropic SGWB measurement obtained using data from the third observing run (O3) of Advanced LIGO detectors and galaxy over-density obtained from the Sloan Digital Sky Survey (SDSS) spectroscopic catalog. We compute, for the first time, the angular power spectrum of their cross-correlation. Instead of integrating the SGWB across frequencies, we analyze the cross-correlation in 10 Hz wide SGWB frequency bands to study the frequency dependence of the cross-correlation angular power spectrum. Finally, we compare the observed cross-correlation to the spectra predicted by astrophysical models. We apply a Bayesian formalism to explore the parameter space of the theoretical models, and we set constraints on a set of (effective) astrophysical parameters describing the galactic process of gravitational wave (GW) emission. Parameterizing with a Gaussian function the astrophysical kernel describing the local process of GW emission at galactic scales, we find the 95\% upper limit on kernel amplitude to be $2.7 \times 10^{-32}$ erg cm$^{-3}$s$^{-1/3}$ when ignoring the shot noise in the GW emission process, and $2.16 \times 10^{-32}$ erg cm$^{-3}$s$^{-1/3}$ when the shot noise is included in the analysis. As the sensitivity of the LIGO-Virgo-KAGRA network improves, we expect to be able to set more stringent bounds on this kernel function and constrain its parameters.

We construct four equation of state (EoS) tables, tabulated over a range of temperatures, densities, and charge fractions, relevant for neutron star applications such as simulations of neutron star mergers. The EoS are computed from a relativistic mean-field theory constrained by the pure neutron matter EoS from chiral effective field theory, inferred properties of isospin-symmetric nuclear matter, and astrophysical observations of neutron star structure. To model nuclear matter at low densities, we attach an EoS that models inhomogeneous nuclear matter at arbitrary temperatures and charge fractions. The four EoS tables we develop are available from the CompOSE EoS repository https://compose.obspm.fr/eos/297 and https://gitlab.com/ahaber/qmc-rmf-tables.

One of the three testaments in favor of the big bang theory is the prediction of the primordial elemental abundances in the big-bang nucleosynthesis (BBN). The Standard BBN is a parameter-free theory due to the precise knowledge of the baryon-to-photon ratio of the Universe obtained from studies of the anisotropies of cosmic microwave background radiation. Although the computed abundances of light elements during primordial nucleosynthesis and those determined from observations are in good agreement throughout a range of nine orders of magnitude, there is still a disparity of $^7$Li abundance overestimated by a factor of $\sim 2.5$ when calculated theoretically. The number of light neutrino flavors, the neutron lifetime and the baryon-to-photon ratio in addition to the astrophysical nuclear reaction rates determine the primordial abundances. We previously looked into the impact of updating baryon-to-photon ratio and neutron lifetime and changing quite a few reaction rates on the yields of light element abundances in BBN. In this work, calculations are performed using new reaction rates for $^3$H(p,$\gamma$)$^4$He, $^6$Li(p,$\gamma$)$^7$Be, $^7$Be(p,$\gamma$)$^8$B, $^{13}$N(p,$\gamma$)$^{14}$O, $^7$Li(n,$\gamma$)$^8$Li and $^{11}$B(n,$\gamma$)$^{12}$B along with the latest measured value of neutron lifetime. We observe from theoretical calculations that these changes result in marginal improvement over a sizable twelve percent reduction in the abundance of $^7$Li achieved earlier.

The nocturnal mesopause region of the Earth's atmosphere radiates chemiluminescent emission from various roto-vibrational bands of hydroxyl (OH), which is therefore a good tracer of the chemistry and dynamics at the emission altitudes. Intensity variations can, e.g., be caused by the general circulation, gravity waves, tides, planetary waves, and the solar activity. While the basic OH response to the different dynamical influences has been studied quite frequently, detailed comparisons of the various individual lines are still rare. Such studies can improve our understanding of the OH-related variations as each line shows a different emission profile. We have therefore used about 90,000 spectra of the X-shooter spectrograph of the Very Large Telescope at Cerro Paranal in Chile in order to study 10 years of variations of 298 OH lines. The analysis focuses on climatologies of intensity, solar cycle effect, and residual variability (especially with respect to time scales of hours and about 2 days) for day of year and local time. For a better understanding of the resulting variability patterns and the line-specific differences, we applied decomposition techniques, studied the variability depending on time scale, and calculated correlations. As a result, the mixing of thermalized and nonthermalized OH level populations clearly influences the amplitude of the variations. Moreover, the local times of the variability features shift depending on the effective line emission height, which can mainly be explained by the propagation of the migrating diurnal tide. This behavior also contributes to remarkable differences in the effective solar cycle effect.

The concept of quarkyonic matter presents a promising alternative to the conventional models used to describe high-density matter and provides a more nuanced and detailed understanding of the properties of matter under extreme conditions that exist in astrophysical bodies. The aim of this study is to showcase the effectiveness of utilizing the quarkyonic model, in combination with the relativistic mean-field formalism, to parameterize the equation of state at high densities. Through this approach, we intend to investigate and gain insights into various fundamental properties of a static neutron star, such as its compositional ingredients, speed of sound, mass-radius profile, and tidal deformability. The obtained results revealed that the quarkyonic matter equation of state (EOS) is capable of producing a heavy neutron star with the mass range of $\sim$ $2.8 M_\odot$. The results of our inquiry have demonstrated that the EOS for quarkyonic matter not only yields a neutron star with a significantly high mass but also showcases a remarkable degree of coherence with the conformal limit of the speed of sound originating from deconfined QCD matter. Furthermore, we have observed that the tidal deformability of the neutron star, corresponding to the EOSs of quarkyonic matter, is in excellent agreement with the observational constraints derived from the GW170817 and GW190425 events. This finding implies that the quarkyonic model is capable of forecasting the behavior of neutron stars associated with binary merger systems. This aspect has been meticulously scrutinized in terms of merger time, gravitational wave signatures, and collapse times using numerical relativity simulations.

Interplanetary (IP) shocks are disturbances commonly observed in the solar wind. IP shock impacts can cause a myriad of space weather effects in the Earth's magnetopause, inner magnetosphere, ionosphere, thermosphere, and ground magnetic field. The shock impact angle, measured as the angle the shock normal vector performs with the Sun-Earth line, has been shown to be a very important parameter that controls shock geoeffectivess. An extensive review provided by Oliveira and Samsonov (2018) summarized all the work known at the time with respect to shock impact angles and geomagnetic activity; however, this topic has had some progress since Oliveira and Samsonov (2018) and the main goal of this mini review is to summarize all achievements to date in the topic to the knowledge of the author. Finally, this mini review also brings a few suggestions and ideas for future research in the area of IP shock impact angle geoeffectiveness.

Within its Voyage 2050 planning cycle, the European Space Agency (ESA) is considering long-term large class science mission themes. Gravitational-wave astronomy is among the topics under study. This paper presents "LISAmax", a gravitational-wave interferometer concept consisting of three spacecraft located close to the Sun-Earth libration points L3, L4 and L5, forming a triangular constellation with an arm length of 259 million kilometers (to be compared to LISA's 2.5 million kilometer arms). This is the largest triangular formation that can be reached from Earth without a major leap in mission complexity and cost. The sensitivity curve of such a detector is at least two orders of magnitude lower in amplitude than that of LISA. Depending on the choice of other instrument parameters, this makes the detector sensitive to gravitational waves in the micro-Hertz range and opens a new window for gravitational-wave astronomy, not covered by any other planned detector concept. We analyze in detail the constellation stability for a 10-year mission in the full numerical model and compute the orbit transfers using a European launcher and chemical propulsion. The payload design parameters are assessed, and the expected sensitivity curve is compared with a number of potential gravitational-wave sources. No show stoppers are identified at this point of the analysis.

Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated signals from strong lensing by 1) performing targeted searches for subthreshold signals, 2) calculating the degree of overlap amongst the intrinsic parameters and sky location of pairs of signals, 3) comparing the similarities of the spectrograms amongst pairs of signals, and 4) performing dual-signal Bayesian analysis that takes into account selection effects and astrophysical knowledge. We also search for distortions to the gravitational waveform caused by 1) frequency-independent phase shifts in strongly lensed images, and 2) frequency-dependent modulation of the amplitude and phase due to point masses. None of these searches yields significant evidence for lensing. Finally, we use the non-detection of gravitational-wave lensing to constrain the lensing rate based on the latest merger-rate estimates and the fraction of dark matter composed of compact objects.

A variety of supergravity and string models involve hidden sectors where the hidden sectors may couple feebly with the visible sectors via a variety of portals. While the coupling of the hidden sector to the visible sector is feeble its coupling to the inflaton is largely unknown. It could couple feebly or with the same strength as the visible sector which would result in either a cold or a hot hidden sector at the end of reheating. These two possibilities could lead to significantly different outcomes for observables. We investigate the thermal evolution of the two sectors in a cosmologically consistent hidden sector dark matter model where the hidden sector and the visible sector are thermally coupled and their thermal evolution occurs without the assumption of separate entropy conservation for each sector. Within this framework we analyze several phenomena to illustrate their dependence on the initial conditions. These include the allowed parameter space of models, dark matter relic density, proton-dark matter cross section, effective massless neutrino species at BBN time, self-interacting dark matter cross-section, where self-interaction occurs via exchange of dark photon, and Sommerfeld enhancement. Finally fits to the velocity dependence of dark matter cross sections from galaxy scales to the scale of galaxy clusters is given. The analysis indicates significant effects of the initial conditions on the observables listed above. The analysis is carried out within the framework where dark matter is constituted of dark fermions and the mediation between the visible and the hidden sector occurs via the exchange of dark photons. The techniques discussed here may have applications for a wider class of hidden sector models using different mediations between the visible and the hidden sectors to explore the impact of Big Bang initial conditions on observable physics.