Papers for Thursday, Oct 26 2023

Coronal mass ejections (CMEs) are large-scale expulsions of plasma and magnetic fields from the Sun into the heliosphere and are the most important driver of space weather. The geo-effectiveness of a CME is primarily determined by its magnetic field strength and topology. Measurement of CME magnetic fields, both in the corona and heliosphere, is essential for improving space weather forecasting. Observations at radio wavelengths can provide several remote measurement tools for estimating both strength and topology of the CME magnetic fields. Among them, gyrosynchrotron (GS) emission produced by mildly-relativistic electrons trapped in CME magnetic fields is one of the promising methods to estimate magnetic field strength of CMEs at lower and middle coronal heights. However, GS emissions from some parts of the CME are much fainter than the quiet Sun emission and require high dynamic range (DR) imaging for their detection. This thesis presents a state-of-the-art calibration and imaging algorithm capable of routinely producing high DR spectropolarimetric snapshot solar radio images using data from a new technology radio telescope, the Murchison Widefield Array. This allows us to detect much fainter GS emissions from CME plasma at much higher coronal heights. For the first time, robust circular polarization measurements have been jointly used with total intensity measurements to constrain the GS model parameters, which has significantly improved the robustness of the estimated GS model parameters. A piece of observational evidence is also found that routinely used homogeneous and isotropic GS models may not always be sufficient to model the observations. In the future, with upcoming sensitive telescopes and physics-based forward models, it should be possible to relax some of these assumptions and make this method more robust for estimating CME plasma parameters at coronal heights.

Recent observations with JWST and ALMA have revealed extremely massive quiescent galaxies at redshifts of z=3 and higher, indicating both rapid onset and quenching of star formation. Using the cosmological simulation suite Magneticum Pathfinder we reproduce the observed number densities and stellar masses, with 36 quenched galaxies of stellar mass larger than 3e10Msun at z=3.42. We find that these galaxies are quenched through a rapid burst of star-formation and subsequent AGN feedback caused by a particularly isotropic collapse of surrounding gas, occurring on timescales of around 200Myr or shorter. The resulting quenched galaxies host stellar components which are kinematically fast rotating and alpha-enhanced, while exhibiting a steeper metallicity and flatter age gradient compared to galaxies of similar stellar mass. The gas of the galaxies has been metal enriched and ejected. We find that quenched galaxies do not inhabit the densest nodes, but rather sit in local underdensities. We analyze observable metrics to predict future quenching at high redshifts, finding that on shorter timescales <500Myr the ratio M_bh/M_* is the best predictor, followed by the burstiness of the preceding star-formation, t50-t90 (time to go from 50% to 90% stellar mass). On longer timescales, >1Gyr, the environment becomes the strongest predictor, followed by t50-t90, indicating that at high redshifts the consumption of old and lack of new gas are more relevant for long-term prevention of star-formation than the presence of a massive AGN. We predict that relics of such high-z quenched galaxies should best be characterized by a strong alpha enhancement.

Feedback from supermassive black holes (SMBHs) is believed to be a critical driver of the observed quenching of star formation and color bimodality of galaxies above the Milky Way mass scale. In recent years, various forms of SMBH feedback have been implemented as subgrid models in galaxy formation simulations, but most implementations have involved simplified prescriptions or coarse-grained models of the interstellar medium (ISM). We present the first set of FIRE-3 cosmological zoom-in simulations with AGN feedback evolved to $z\sim0$, examining a set of galaxies with halos in the mass range $10^{12}-10^{13}\,{\rm M_{\odot}}$. These high-resolution simulations combine detailed stellar and ISM physics with multi-channel AGN feedback including radiative feedback, mechanical outflows, and in some simulations, cosmic rays. We find that massive (>L*) galaxies in these simulations can match local scaling relations including the stellar mass-halo mass relation, the $M_{\rm BH}$-$\sigma$ relation, the size-mass relation, and the Faber-Jackson relation. Many of the massive galaxies in the simulations with AGN feedback have quenched star formation and elliptical morphologies, in qualitative agreement with observations. In contrast, simulations at the massive end without AGN feedback produce galaxies that are too massive and form stars at too high rates, are order-of-magnitude too compact, and have velocity dispersions well above Faber-Jackson. Despite these successes, the AGN physics models analyzed do not necessarily produce uniformly realistic galaxies across the full mass range studied when the feedback parameters are held constant, indicating that further refinements of the black hole modeling may be warranted.

Characterizing the physical properties of cool supergiants allows us to probe the final stages of a massive star's evolution before it undergoes core collapse. Despite their importance, the fundamental properties for these stars -- $T_{\rm eff}$ and $\log L/L_\odot$ -- are only known for a limited number of objects. The third data release of the Gaia mission contains precise photometry and low-resolution spectroscopy of hundreds of cool supergiants in the LMC with well-constrained properties. Using these data, we train a simple and easily-interpretable machine learning model to regress effective temperatures and luminosities with high accuracy and precision comparable to the training data. We then apply our model to 5000 cool supergiants, many of which have no previously-published $T_{\rm eff}$ or $L$ estimates. The resulting Hertzprung-Russell diagram is well-populated, allowing us to study the distribution of cool supergiants in great detail. Examining the luminosity functions of our sample, we find a notable flattening in the luminosity function of yellow supergiants above $\log L/L_\odot=5$, and a corresponding steepening of the red supergiant luminosity function. We place this finding in context with previous results, and present its implications for the infamous red supergiant problem.

A puzzling population of extremely massive quiescent galaxies at redshifts beyond z=3 has recently been revealed by JWST and ALMA, some of them with stellar ages that show their quenching times to be as high as z=6, while their stellar masses are already above 5e10Msun. These extremely massive yet quenched galaxies challenge our understanding of galaxy formation at the earliest stages. Using the hydrodynamical cosmological simulation suite Magneticum Pathfinder, we show that such massive quenched galaxies at high redshifts can be successfully reproduced with similar number densities as observed. The stellar masses, sizes, formation redshifts, and star formation histories of the simulated quenched galaxies match those determined with JWST. Following these quenched galaxies at z=3.4 forward in time, we find 20% to be accreted onto a more massive structure by z=2, and from the remaining 80% about 30% rejuvenate up to z=2, another 30% stay quenched, and the remaining 40% rejuvenated on a very low level of star formation. Stars formed through rejuvenation are mostly formed on the outer regions of the galaxies, not in the centres. Furthermore, we demonstrate that the massive quenched galaxies do not reside in the most massive nodes of the cosmic web, but rather live in side-nodes of approximately Milky-Way halo mass. Even at z=0, only about 10% end up in small-mass galaxy clusters, while most of the quenched galaxies at z=3.4 end up in group-mass halos, with about 20% actually not even reaching 1e13Msun in halo mass.

We present high cadence optical and ultraviolet observations of the Type II supernova (SN), SN 2022jox which exhibits early spectroscopic high ionization flash features of \ion{H}{1}, \ion{He}{2}, \ion{C}{4}, and \ion{N}{4} that disappear within the first few days after explosion. SN 2022jox was discovered by the Distance Less than 40 Mpc (DLT40) survey $\sim$0.75 days after explosion with followup spectra and UV photometry obtained within minutes of discovery. The SN reached a peak brightness of M$_V \sim$ $-$17.3 mag, and has an estimated $^{56}$Ni mass of 0.04 M$_{\odot}$, typical values for normal Type II SNe. The modeling of the early lightcurve and the strong flash signatures present in the optical spectra indicate interaction with circumstellar material (CSM) created from a progenitor with a mass loss rate of $\dot{M} \sim 10^{-3}-10^{-2}\ M_\odot\ \mathrm{yr}^{-1}$. There may also be some indication of late-time CSM interaction in the form of an emission line blueward of H$\alpha$ seen in spectra around 200 days. The mass-loss rate is much higher than the values typically associated with quiescent mass loss from red supergiants, the known progenitors of Type II SNe, but is comparable to inferred values from similar core collapse SNe with flash features, suggesting an eruptive event or a superwind in the progenitor in the months or years before explosion.

We use the Voids-within-Voids-within-Voids (VVV) simulations, a suite of successive nested N-body simulations with extremely high resolution (denoted, from low to high resolution, by L0 to L7), to test the Press-Schechter (PS), Sheth-Tormen (ST), and extended Press-Schechter (EPS) formulae for the halo abundance over the entire mass range, from mini-haloes of $10^{-6}\ \mathrm{M_\odot}$, to rich cluster haloes of $10^{15}\ \mathrm{M_\odot}$, at different redshifts, from $z=30$ to the present. We find that at $z=0$ and $z=2$, ST gives the best prediction for the results of L0, which has the mean cosmic density (i.e. overdensity $\delta=0$), at $10^{11-15} ~\mathrm{M_\odot}$. The higher resolution levels (L1-L7) are biased regions at various underdensities ($\delta<-0.6$). The EPS formalism takes this into account since it gives the mass function of a region conditioned, in this case, on having a given underdensity. The EPS gives good predictions for these higher levels, with deviations $\lesssim 20\%$, at $10^{-6-12.5} ~\mathrm{M_\odot}$. At $z \sim 7-15$, the ST predictions for L0 and the EPS predictions for L1-L7 show somewhat larger deviations from the simulation results. However, at even higher redshifts, $z \sim 30$, the EPS prediction fits the simulations well again. We further confirm our results by picking more subvolumes from the full L0 simulation, finding that our conclusions depend only weakly on the size and overdensity of the chosen region. Since at mean density the EPS reduces to the PS mass function, its good agreement with the higher-level simulations implies that the PS (or, even better, the ST) formula gives an accurate description of the total halo mass function in representative regions of the universe.

We investigate how the recovery of galaxy star formation rates (SFRs) using energy-balance spectral energy distribution (SED) fitting codes depends on their recent star formation histories (SFHs). We use the Magphys and Prospector codes to fit 6,706 synthetic spectral energy distributions (SEDs) of simulated massive galaxies at $1 < z < 8$ from the Feedback in Realistic Environments (FIRE) project. We identify a previously unknown systematic error in the Magphys results due to bursty star formation: the SFR estimates of individual galaxies can differ from the true values by as much as 1 dex, at large statistical significance ($>5\sigma$), depending on the details of their recent SFH. The SFRs inferred using Prospector do not exhibit this trend, likely because unlike Magphys, Prospector uses non-parametric SFHs. We urge caution when using Magphys, or other codes assuming parametric SFHs, to study galaxies where the average SFR may have changed significantly over the last $\sim$100 Myr, such as those which have recently quenched their star formation or those experiencing an ongoing burst. This concern is especially relevant, for example, when fitting JWST observations of very high-redshift galaxies.

Exoplanets have been detected around stars at various stages of their lives, ranging from young stars emerging from formation, to latter stages of evolution, including white dwarfs and neutron stars. Post main sequence stellar evolution can result in dramatic, and occasionally traumatic, alterations to the planetary system architecture, such as tidal disruption of planets and engulfment by the host star. The $\rho$ CrB system is a particularly interesting case of advanced main sequence evolution, due to the relative late age and brightness of the host star, its similarity to solar properties, and the harboring of four known planets. Here, we use stellar evolution models to estimate the expected trajectory of the $\rho$ CrB stellar properties, especially over the coming 1.0-1.5 billion years as it evolves off the main sequence. We show that the inner three planets (e, b, and c) are engulfed during the red giant phase and asymptotic giant branch, likely destroying those planets via either evaporation or tidal disruption at the fluid body Roche limit. The outer planet, planet d, is briefly engulfed by the star several times toward the end of the asymptotic giant branch, but the stellar mass loss and subsequent changing planetary orbit may allow the survival of the planet into the white dwarf phase of the stellar evolution. We discuss the implications of this outcome for similar systems, and describe the consequences for planets that may lie within the Habitable Zone of the system.

We perform a reanalysis of the BOSS CMASS DR12 galaxy dataset using a simulation-based emulator for the Wavelet Scattering Transform (WST) coefficients. Moving beyond our previous works, which laid the foundation for the first galaxy clustering application of this estimator, we construct a neural net-based emulator for the cosmological dependence of the WST coefficients and the 2-point correlation function multipoles, trained from the state-of-the-art suite of \textsc{AbacusSummit} simulations combined with a flexible Halo Occupation Distribution (HOD) galaxy model. In order to confirm the accuracy of our pipeline, we subject it to a series of thorough internal and external mock parameter recovery tests, before applying it to reanalyze the CMASS observations in the redshift range $0.46<z<0.57$. We find that a joint WST + 2-point correlation function likelihood analysis allows us to obtain marginalized 1$\sigma$ errors on the $\Lambda$CDM parameters that are tighter by a factor of $2.5-6$, compared to the 2-point correlation function, and by a factor of $1.4-2.5$ compared to the WST-only results. This corresponds to a competitive $0.9\%$, $2.3\%$ and $1\%$ level of determination for parameters $\omega_c$, $\sigma_8$ $\&$ $n_s$, respectively, and also to a $0.7\%$ $\&$ $2.5 \%$ constraint on derived parameters h and $f(z)\sigma_8(z)$, in agreement with the \textit{Planck} 2018 results. Our results reaffirm the constraining power of the WST and highlight the exciting prospect of employing higher-order statistics in order to fully exploit the power of upcoming Stage-IV spectroscopic observations.

Recent measurements of Jupiter's gravitational field (by Juno) and seismology of Saturn's rings (by Cassini) strongly suggest that both planets have a stably-stratified core that still possesses a primordial gradient in the concentration of heavy elements. The existence of such a "diffusely" stratified core has been a surprise as it was long expected that the Jovian planets should be fully convective and hence fully mixed. A vigorous zone of convection, driven by surface cooling, forms at the surface and deepens through entrainment of fluid from underneath. In fact, it was believed that this convection zone should grow so rapidly that the entire planet would be consumed in less than a million years. Here we suggest that two processes, acting in concert, present a solution to this puzzle. All of the giant planets are rapidly rotating and have a cooling rate that declines with time. Both of these effects reduce the rate of fluid entrainment into the convection zone. Through the use of an analytic prescription of entrainment in giant planets, we demonstrate that these two effects, rotation and dwindling surface cooling, result in a convection zone which initially grows but eventually stalls. The depth to which the convective interface asymptotes depends on the rotation rate and on the stratification of the stable interior. Conversely, in a nonrotating planet, or in a planet that maintains a higher level of cooling than current models suggest, the convection zone deepens forever, eventually spanning the entire planet.

Accreting binary white dwarf systems are among the sources expected to emanate gravitational waves that the Laser Interferometer Space Antenna (LISA) will detect. We investigate how accurately the binary parameters may be measured from LISA observations. We complement previous studies by performing our parameter estimation on binaries containing a low-mass donor with a thick, hydrogen-rich envelope. The evolution is followed from the early, pre-period minimum stage, in which the donor is non-degenerate, to a later, post-period minimum stage with a largely degenerate donor. We present expressions for the gravitational wave amplitude, frequency, and frequency derivative in terms of white dwarf parameters (masses, donor radius, etc.), where binary evolution is driven by gravitational wave radiation and accretion torques, and the donor radius and logarithmic change in radius ($\eta_{\rm d}$) due to mass loss are treated as model parameters. We then perform a Fisher analysis to reveal the accuracy of parameter measurements, using models from Modules for Experiments in Stellar Astrophysics (MESA) to estimate realistic fiducial values at which we evaluate the measurement errors. We find that the donor radius can be measured relatively well with LISA observations alone, while we can further measure the individual masses if we have an independent measurement of the luminosity distance from electromagnetic observations. When applied to the parameters of the recently-discovered white dwarf binary ZTF J0127+5258, our Fisher analysis suggests that we will be able to constrain the system's individual masses and donor radius using LISA's observations, given ZTF's measurement of the luminosity distance.

A stellar occultation of Gaia DR3 2646598228351156352 by the Centaur (2060) Chiron was observed from the South African Astronomical Observatory on 2018 November 28 UT. Here we present a positive detection of material surrounding Chiron from the 74-in telescope for this event. Additionally, a global atmosphere is ruled out at the tens of mircobar level for several possible atmospheric compositions. There are multiple 3-sigma drops in the 74-in light curve: three during immersion and two during emersion. Occulting material is located between 242-270 km from the center of the nucleus in the sky plane. Assuming the ring-plane orientation proposed for Chiron from the 2011 occultation, the flux drops are located at 352, 344, and 316 km (immersion), and 357, and 364 km (emersion) from the center, with normal optical depths of 0.26, 0.36, and 0.22 (immersion) and 0.26 and 0.18 (emersion), and equivalent widths between 0.7-1.3 km. This detection is similar to the previously proposed two-ring system and is located within the error bars of that ring-pole plane; however, the normal optical depths are less than half of the previous values, and three features are detected on immersion. These results suggest that the properties of the surrounding material have evolved between the 2011, 2018, and 2022 observations.

Long-period comet C/2020 S3 (Erasmus) reached perihelion at 0.398 au on UT 2020 December 12.67, making it a bright, near-Sun object. Images taken between mid-November and December 2020 using the HI-1 camera and COR2 coronagraph onboard STEREO-A, as well as the LASCO/C3 coronagraph onboard SoHO, show significant variations in the plasma tail position angles. To analyze these variations, a simple technique was developed to calculate the aberration angles. These angles are defined as the angle between the sun-comet line and the tail axis, measured in the orbital plane. The aberration angles were found to range from $1.2^\circ$ to $46.8^\circ$, with an average (median) value of approximately $20.3^\circ$ ($16.3^\circ$). By considering the aberration angles, the solar wind radial velocities during the observations were inferred to range from 73.9 km/s to 573.5 km/s, with mean (median) values of approximately 205.5 km/s (182.3 km/s). Throughout the observations, two periods were identified where the tails showed forward tilting, which cannot be explained by aberration alone. In one case, this anomalous position angle was sustained for at least 11 days and is possibly due to co-rotating interaction regions. In the other case, the tail exhibited dramatic excursions from 180$^\circ$ to 150$^\circ$ back to 210$^\circ$ over a limited period of around 34 hours. This behavior is tentatively explained as a consequence of the interaction with a halo Coronal Mass Ejection that was launched from NOAA 12786 and arrived at comet C/2020 S3 during the time when the tail displayed its wagging behavior.

The variability of active galactic nuclei (AGNs) is ubiquitous but has not yet been understood. Measuring the optical variation properties of AGNs, such as variation timescale and amplitude, and then correlating them with their fundamental physical parameters, have long served as a critical way of exploring the origin of AGN variability and the associated physics of the accretion process in AGNs. Obtaining accurate variation properties of AGNs is thus essential. It has been found that the damped random walk (DRW) process can well describe the AGN optical variation, however, there is a controversy over how long a minimal monitoring baseline is required to obtain unbiased variation properties. In this work, we settle the controversy by exhaustively scrutinizing the complex combination of assumed priors, adopted best-fit values, ensemble averaging methods, and fitting methods. Then, the newly proposed is an optimized solution where unbiased variation properties of an AGN sample possessing the same variation timescale can be obtained with a minimal baseline of about 10 times their variation timescale. Finally, the new optimized solution is used to demonstrate the positive role of time domain surveys to be conducted by the Wide Field Survey Telescope in improving constraints on AGN variation properties.

Recent studies of Galactic evolution revealed that the dynamics of the stellar component might be one of the key factors when considering galactic habitability. We run an N-body simulation model of the Milky Way, which we evolve for 10 Gyr, to study the secular evolution of stellar orbits and the resulting galactic habitability-related properties, i.e., the density of the stellar component and close stellar encounters. The results indicate that radial migrations are not negligible, even in a simple axisymmetric model with mild levels of dynamical heating, and that the net outward diffusion of the stellar component can populate galactic outskirts with habitable systems. Habitable environment is also likely even at sub-Solar galactocentric radii, because the rate of close encounters should not significantly degrade habitability. Stars that evolve from non-circular to stable nearly-circular orbits typically migrate outwards, settling down in a broad Solar neighborhood. The region between $R \approx 3$ kpc and $R \approx 12$ kpc represents the zone of radial mixing, which can blur the boundaries of the Galactic Habitable Zone, as it has been conventionally understood. The present-day stable population of the stars in the Solar neighborhood originates from this radial mixing zone, with most of the stars coming from the inner regions. The Solar system can be considered as a typical Milky Way habitable system because it migrated outwards from the metal-rich inner regions of the Disk and has a circular orbit in the present epoch. We conclude that the boundaries of the Galactic Habitable Zone cannot be sharply confined for a given epoch because of the mixing caused by the stellar migrations and secular evolution of stellar orbits.

We present a comprehensive analysis of the variable stars projected on the field of the Galactic bulge globular cluster NGC 6522, offering valuable insights into their characteristics. Using proper motions from Gaia DR3, we aim to distinguish between field stars and true cluster members. For an accurate color-magnitude diagram of the member stars, we produced a differential reddening map. We detect and discuss the peculiarities of variable stars of the types RR Lyrae, type II Cepheids, Long Period Variables (LPV) and eclipsing binaries, whose light curves are available through the OGLE III and IV databases. Notably, we explore the variable V24, which shows a prominent phase modulation resulting from period changes in a time scale of ~1100 days. The variable stars among the cluster members serve as indicators of the cluster metallicity and distance; these determinations are based on their light curves. With the Fourier light curve decomposition of three RRc stars, we have derived the following cluster parameters: the metallicity in the spectroscopic scale [Fe/H]$_{\rm UVES}$=--1.16 $\pm$0.09; and the mean distance $D=8.77 \pm 0.16$ kpc.

In cosmology, the quest for primordial $B$-modes in cosmic microwave background (CMB) observations has highlighted the critical need for a refined model of the Galactic dust foreground. We investigate diffusion-based modeling of the dust foreground and its interest for component separation. Under the assumption of a Gaussian CMB with known cosmology (or covariance matrix), we show that diffusion models can be trained on examples of dust emission maps such that their sampling process directly coincides with posterior sampling in the context of component separation. We illustrate this on simulated mixtures of dust emission and CMB. We show that common summary statistics (power spectrum, Minkowski functionals) of the components are well recovered by this process. We also introduce a model conditioned by the CMB cosmology that outperforms models trained using a single cosmology on component separation. Such a model will be used in future work for diffusion-based cosmological inference.

We present the multiwavelength flaring activity of the blazar AO 0235+164 during its recent active period from 2013 to 2019. From a discrete correlation function (DCF) analysis, we find a significant (>95%) correlation between radio and $\gamma$-ray light curves with flares at longer wavelengths following flares at shorter wavelengths. We identify a new jet component in 43 GHz VLBA data that was ejected from the radio core on MJD $57246^{+26}_{-30}$ (2015 August 12), during the peak of the 2015 radio flare. From the analysis of the jet component, we derived a Doppler factor of $\delta_{\rm var}=28.5\pm8.4$, a bulk Lorentz factor of $\Gamma=16.8^{+3.6}_{-3.1}$, and an intrinsic viewing angle of $\theta_{\rm v}=1.42^{+1.07}_{-0.52}\textrm{ degrees}$. Investigation of the quasi-simultaneous radio data revealed a partially absorbed spectrum with the turnover frequency varying in the range of $10-70$ GHz and the peak flux density varying in the range of $0.7-4$ Jy. We find the synchrotron self-absorption magnetic field strength to be $B_{\rm SSA}=15.3^{+12.6}_{-14.0}$ mG at the peak of the 2015 radio flare, which is comparable to the equipartition magnetic field strength of $B_{\rm EQ}=43.6^{+10.6}_{-10.4}$ mG calculated for the same epoch. Additional analysis of the radio emission region in the relativistic jet of AO 0235+164 suggests that it did not significantly deviate from equipartition during its recent flaring activity.

Photometric redshift estimation, an essential process in astronomy for distance estimation, obtains the redshift of celestial structures by utilizing the magnitude of objects in varying wavelength filters. This research capitalized on a dataset of 50,000 objects from the Sloan Digital Sky Survey, comprising 5 bands of magnitudes and their corresponding redshift labels. Typically, studies use spectral distribution templates (SED) for redshift prediction. However, these templates are expensive and hard to obtain, especially with larger datasets. The paper explores approaches for Data-Driven methodology instead of template based prediction. Adopting both a decision tree regression model and a Fully Connected Neural Network (FCN) for analysis, the FCN significantly outperformed the decision tree regressor, achieving an impressive root mean square error (RMSE) of 0.009 compared to the decision tree's RMSE above 0.16. The strong performance of the FCN highlights its ability to capture intricate relationships in astronomical data, holding the potential for data-driven redshift estimation, which will help advance next generation surveys.

Fast radio bursts (FRBs) are transient radio signals with millisecond-duration, large dispersion measure (DM) and extremely high brightness temperature. Among them, FRB 20180916B has been found to have a 16-day periodic activity. However, the physical origin of the periodicity is still a mystery. Here, we utilize the comprehensive observational data to diagnose the periodic models. We find that the ultra-long rotation model is the most probable one for the periodic activity. However, this model cannot reproduce the observed rotation measure (RM) variations. We propose a self-consistent model, i.e., a massive binary containing a slowly rotational neutron star and a massive star with large mass loss, which can naturally accommodate the wealth of observational features for FRB 20180916B. In this model, the RM variation is periodic, which can be tested by future observations.

We report the discovery of SDSS~J022932.28+713002.7, a nascent extremely low-mass (ELM) white dwarf (WD) orbiting a massive ($> 1\,M_\odot$ at 2$\sigma$ confidence) companion with a period of 36 hours. We use a combination of spectroscopy, including data from the ongoing SDSS-V survey, and photometry to measure the stellar parameters for the primary pre-ELM white dwarf. The lightcurve of the primary WD exhibits ellipsoidal variation, which we combine with radial velocity data and $\tt{PHOEBE}$ binary simulations to estimate the mass of the invisible companion. We find that the primary WD has mass $M_1$ = $0.18^{+0.02}_{-0.02}$ M$_\odot$ and the unseen secondary has mass $M_2$ = $1.19^{+0.21}_{-0.14}$ M$_\odot$. The mass of the companion suggests that it is most likely a near-Chandrasekhar mass white dwarf or a neutron star. It is likely that the system recently went through a Roche lobe overflow from the visible primary onto the invisible secondary. The dynamical configuration of the binary is consistent with the theoretical evolutionary tracks for such objects, and the primary is currently in its contraction phase. The measured orbital period puts this system on a stable evolutionary path which, within a few Gyrs, will lead to a contracted ELM white dwarf orbiting a massive compact companion.

The astronomy community has witnessed an explosive growth in the use of deep-learning techniques based on neural networks since the mid-2010s. The widespread adoption of these nature-inspired technologies has helped astronomers tackle previously insurmountable problems and provided an unprecedented opportunity for new discoveries. However, one of the primary tools of today's optical astronomy is neither natural nor efficient: their photo-sensing devices. Specifically, the modern CCD camera - like that of the cutting-edge Rubin Observatory - requires an internal clock to regularly expose the sensor to light, consumes a large amount of energy and information bandwidth, and has a limited dynamic range. On the contrary, biological eyes lack an internal clock and a shutter, have much higher pixel density but consume significantly less energy and bandwidth, and can adapt to bright and low light conditions. Inspired by the nature of the eyes, M. Mahowald and C. Mead introduced the revolutionary concept of a silicon retina sensor in 1991. Also known as event-based cameras (EBCs), these types of devices operate in a vastly different way compared to conventional CCD-based imaging sensors. EBCs mimic the operating principles of optic nerves and continuously produce a stream of events, with each event generated only when a pixel detects a change in light intensity. EBCs do not have fixed exposure times, have high dynamic range, require low power for operation, and can capture high-speed phenomena. These properties are important requirements for Cherenkov telescopes as well as other high-speed optical astronomy. This work presents the opportunities and challenges of using EBCs in those cases, and proposes a low-cost approach to experimentally assess the feasibility of this innovative technique.

Previous numerical simulations of double-peaked SNe IIb light curves have demonstrated that the radius and mass of the hydrogen-rich envelope of the progenitor star can significantly influence the brightness and timescale of the early-time light curve around the first peak. In this study, we investigate how Thomson scattering and chemical mixing in the SN ejecta affect the optical light curves during the early stages of the SNe IIb using radiation hydrodynamics simulations. By comparing the results from two different numerical codes (i.e., \stella{} and \snec{}), we find that the optical brightness of the first peak can be reduced by more than a factor of 3 due to the effect of Thomson scattering that causes the thermalization depth to be located below the Rosseland-mean photosphere, compared to the corresponding case where this effect is ignored. We also observe a short-lived plateau-like feature lasting for a few days in the early-time optical light curves of our models, in contrast to typical observed SNe IIb that show a quasi-linear decrease in optical magnitudes after the first peak. A significant degree of chemical mixing between the hydrogen-rich envelope and the helium core in SN ejecta is required to reconcile this discrepancy between the model prediction and observation. Meanwhile, to properly reproduce the first peak, a significant mixing of \nifs{} into the hydrogen-rich outermost layers should be restricted. Our findings indicate that inferring the SN IIb progenitor structure from a simplified approach that ignores these two factors may introduce substantial uncertainty.

The Herschel Gould Belt Survey showed that stars form in dense filaments in nearby molecular clouds. Recent studies suggest that massive filaments are bound by the slow shocks caused by accretion flows onto the filaments. The slow shock is known to be unstable to the corrugation deformation of the shock front. The corrugation instability could convert the accretion flow's ram pressure into turbulent pressure that influences the width of the filament, which, according to theory, determines the self-gravitational fragmentation scale and core mass. In spite of its importance, the effect of slow shock instability on star-forming filaments has not been investigated. In addition, the linear dispersion relation obtained from the ideal magnetohydrodynamics (MHD) analysis shows that the most unstable wavelength of shock corrugation is infinitesimally small (or mean free path). In the scale of dense filaments, the effect of ambipolar diffusion can suppress the instability at small scales. This study investigates the influence of ambipolar diffusion on the instability of the slow shock. We perform two-dimensional MHD simulations to examine the linear growth of the slow shock instability, considering the effect of ambipolar diffusion. The results demonstrate that the most unstable scale of slow shock instability is approximately five times the length scale of ambipolar diffusion l_AD calculated using post-shock variables, where, l_AD corresponds to the scale where the magnetic Reynolds number for ambipolar diffusivity is unity.

We investigate the nature of correlations in the small-scale flux statistics of the Lyman-$\alpha$ (Ly$\alpha$) forest across redshift bins. Understanding these correlations is important for unbiased cosmological and astrophysical parameter inference using the Ly$\alpha$ forest. We focus on the 1-dimensional flux power spectrum (FPS) and mean flux ($\bar F$) simulated using the semi-numerical lognormal model we developed in earlier work. The lognormal model can capture the effects of long wavelength modes with relative ease as compared to full smoothed particle hydrodynamical (SPH) simulations that are limited by box volume. For a single redshift bin of size $\Delta z\simeq 0.1$, we show that the lognormal model predicts positive cross-correlations between $k$-bins in the FPS, and a negative correlation for $\bar F\times$ FPS, in qualitative agreement with SPH simulations and theoretical expectations. For measurements across two neighbouring redshift bins of width $\Delta z$ each (obtained by 'splitting' skewers of length $2\Delta z$ in half), the lognormal model predicts an anti-correlation for FPS $\times$ FPS and a positive correlation for $\bar F\times$ FPS, caused by long wavelength modes. This is in contrast to SPH simulations which predict a negligible magnitude for cross-redshift correlations derived from such `split' skewers, and we discuss possible reasons for this difference. Finally, we perform a preliminary test of the impact of neglecting long wavelength modes on parameter inference, finding that whereas the correlation structure of neighbouring redshift bins has relatively little impact, the absence of long wavelength modes in the model can lead to $\gtrsim2-\sigma$ biases in the inference of astrophysical parameters. Our results motivate a more careful treatment of long wavelength modes in analyses that rely on the small scale Ly$\alpha$ forest for parameter inference.

Pseudodisks are dense structures formed perpendicular to the direction of the magnetic field during the gravitational collapse of a molecular cloud core. Numerical simulations of the formation of pseudodisks are usually computationally expensive with conventional CPU codes. To demonstrate the proof-of-concept of a fast computing method for this numerically costly problem, we explore the GPU-powered MHD code Astaroth, a 6th-order finite difference method with low adjustable finite resistivity implemented with sink particles. The formation of pseudodisks is physically and numerically robust and can be achieved with a simple and clean setup for this newly adopted numerical approach for science verification. The method's potential is illustrated by evidencing the dependence on the initial magnetic field strength of specific physical features accompanying the formation of pseudodisks, e.g. the occurrence of infall shocks and the variable behavior of the mass and magnetic flux accreted on the central object. As a performance test, we measure both weak and strong scaling of our implementation to find most efficient way to use the code on a multi-GPU system. Once suitable physics and problem-specific implementations are realized, the GPU-accelerated code is an efficient option for 3-D magnetized collapse problems.

RX J0440.9+4431 is an accreting X-ray pulsar (XRP) that remained relatively unexplored until recently, when major X-ray outburst activity enabled more in-depth studies. Here, we report on the discovery of ${\sim}0.2$ Hz quasi-periodic oscillations (QPOs) from this source observed with $Fermi$-GBM. The appearance of QPOs in RX J0440.9+4431 is thricely transient, that is, QPOs appear only above a certain luminosity, only at certain pulse phases (namely corresponding to the peak of its sine-like pulse profile), and only for a few oscillations at time. We argue that this newly discovered phenomenon (appearance of thricely transient QPOs -- or ATTO) occurs if QPOs are fed through an accretion disk whose inner region viscosity is unstable to mass accretion rate and temperature variations. Such variations are triggered when the source switches to the super-critical accretion regime and the emission pattern changes. We also argue that the emission region configuration is likely responsible for the observed QPOs spin-phase dependence.

Recently, several major pulsar timing array (PTA) collaborations have assembled strong evidence for the existence of a gravitational-wave background at frequencies around the nanohertz regime. Assuming that the PTA signal is attributed to scalar-induced gravitational waves, we jointly employ the PTA data from the NANOGrav 15-year data set, PPTA DR3, and EPTA DR2 to probe the conditions of the early Universe. Specifically, we explore the equation of state parameter ($w$), the sound speed ($c_s$), and the reheating temperature ($T_\mathrm{rh}$), finding $w=0.60^{+0.32}_{-0.39}$, $c_s\gtrsim 0.09$, and $T_\mathrm{rh}\lesssim 0.2\,\mathrm{GeV}$ for a lognormal power spectrum of the curvature perturbation. Furthermore, we compute Bayes factors to compare different models against the radiation domination model ($c_s^2 = w = 1/3$), effectively excluding the pressure-less fluid domination model. Our study underscores the significance of scalar-induced gravitational waves as a powerful tool to explore the nature of the early Universe.

In this paper we present the European Low Frequency Survey (ELFS), a project that will enable foregrounds-free measurements of primordial $B$-mode polarization to a level 10$^{-3}$ by measuring the Galactic and extra-Galactic emissions in the 5--120\,GHz frequency window. Indeed, the main difficulty in measuring the B-mode polarization comes not just from its sheer faintness, but from the fact that many other objects in the Universe also emit polarized microwaves, which mask the faint CMB signal. The first stage of this project will be carried out in synergy with the Simons Array (SA) collaboration, installing a 5.5--11 GHz coherent receiver at the focus of one of the three 3.5\,m SA telescopes in Atacama, Chile ("ELFS on SA"). The receiver will be equipped with a fully digital back-end based on the latest Xilinx RF System-on-Chip devices that will provide frequency resolution of 1\,MHz across the whole observing band, allowing us to clean the scientific signal from unwanted radio frequency interference, particularly from low-Earth orbit satellite mega-constellations. This paper reviews the scientific motivation for ELFS and its instrumental characteristics, and provides an update on the development of ELFS on SA.

In this paper, using the latest Pantheon+ sample of Type Ia supernovae (SNe Ia), Baryon Acoustic Oscillation (BAO) measurements, and observational Hubble data (OHD), we carry out a joint constraint on the Hubble constant $H_0$, the spatial curvature $\Omega_{\rm K}$, and the sound horizon at the end of drag epoch $r_{\rm d}$. To be model-independent, four cosmography models, i.e., the Taylor series in terms of redshift $y_1=z/(1+z)$, $y_2=\arctan(z)$, $y_3=\ln(1+z)$, and the Pad\'e approximants, are used without the assumption of flat Universe. The results show that the $H_0$ is anti-correlated with $\Omega_{\rm K}$ and $r_{\rm d}$, indicating smaller $\Omega_{\rm K}$ or $r_{\rm d}$ would be helpful in alleviating the Hubble tension. And the values of $H_0$ and $r_{\rm d}$ are consistent with the estimate derived from the Planck Cosmic Microwave Background (CMB) data based on the flat $\Lambda$CDM model, but $H_0$ is in 2.3$\sim$3.0$\sigma$ tension with that obtained by \cite{Riess2022} in all these cosmographic approaches. Meanwhile, a flat Universe is preferred by the present observations under all approximations except the third order of $y_1$ and $y_2$ of the Taylor series. Furthermore, according to the values of the Bayesian evidence, we found that the flat $\Lambda$CDM remains to be the most favored model by the joint datasets, and the Pad\'e approximant of order (2,2), the third order of $y_3$ and $y_1$ are the top three cosmographic expansions that fit the datasets best, while the Taylor series in terms of $y_2$ are essentially ruled out.

To understand the evolution of dust properties in molecular clouds in the course of the star formation process, we constrain the changes in the dust emissivity index from star-forming filaments to prestellar and protostellar cores to T Tauri stars. Using the NIKA2 continuum camera on the IRAM 30~m telescope, we observed the Taurus B211/B213 filament at 1.2\,mm and 2\,mm with unprecedented sensitivity and used the resulting maps to derive the dust emissivity index $\beta$. Our sample of 105 objects detected in the $\beta$ map of the B211/B213 filament indicates that, overall, $\beta$ decreases from filament and prestellar cores ($\beta \sim 2\pm0.5$) to protostellar cores ($\beta \sim 1.2 \pm 0.2$) to T-Tauri protoplanetary disk ($\beta < 1$). The averaged dust emissivity index $\beta$ across the B211/B213 filament exhibits a flat ($\beta \sim 2\pm0.3$) profile. This may imply that dust grain sizes are rather homogeneous in the filament, start to grow significantly in size only after the onset of the gravitational contraction/collapse of prestellar cores to protostars, reaching big sizes in T Tauri protoplanetary disks. This evolution from the parent filament to T-Tauri disks happens on a timescale of about 1-2~Myr.

We present new and archival Australia Telescope Compact Array and Atacama Large Millimeter/submillimeter Array data of the Small Magellanic Cloud supernova remnant 1E 0102.2-7219 at 2100, 5500, 9000, and 108000 MHz; as well as Hi data provided by the Australian Square Kilometre Array Pathfinder. The remnant shows a ring-like morphology with a mean radius of 6.2 pc. The 5500 MHz image reveals a bridge-like structure, seen for the first time in a radio image. This structure is also visible in both optical and X-ray images. In the 9000 MHz image we detect a central feature that has a flux density of 4.3 mJy but rule out a pulsar wind nebula origin, due to the lack of significant polarisation towards the central feature with an upper limit of 4 per cent. The mean fractional polarisation for 1E 0102.2-7219 is 7 +- 1 and 12 +- 2 per cent for 5500 and 9000 MHz, respectively. The spectral index for the entire remnant is -0.61 +- 0.01. We estimate the line-of-sight magnetic field strength in the direction of 1E 0102.2-7219 of ~44 microG with an equipartition field of 65 +- 5 microG. This latter model, uses the minimum energy of the sum of the magnetic field and cosmic ray electrons only. We detect an Hi cloud towards this remnant at the velocity range of ~160-180 km s-1 and a cavity-like structure at the velocity of 163.7-167.6 km s-1. We do not detect CO emission towards 1E 0102.2-7219.

We consider the finite velocity of the ejecta of a type Ia supernova (SN Ia) in the double detonation (DDet) scenario with a white dwarf (WD) mass-donor companion, and find that the runaway velocity of the surviving (mass donor) WD is lower than its pre-explosion orbital velocity by about 8-11 %. This implies that the fastest runaway WDs in the Galaxy, if come from the DDet scenario, require even more massive WDs than what a simple calculation that neglects the finite ejecta velocity gives. This extreme set of initial conditions makes such binaries less common. We also tentatively find that the inner ejecta largely deviate from spherical symmetry, unlike observed SN Ia remnants. We argue that these findings support the claim that the DDet scenario leads mostly to peculiar SNe Ia but not to normal SNe Ia.

We present the first test of coasting cosmological models with gravitational-wave standard sirens observed in the first three observing runs of the LIGO-Virgo-KAGRA detector network. We apply the statistical galaxy catalog method adapted to coasting cosmologies and infer constraints on the $H_0$ Hubble constant for the three fixed values of the curvature parameter $k=\left\{ -1,0,+1 \right\}$ in $H_0^2 c^{-2}$ units. The maximum posteriors and $68.3\%$ highest density intervals we obtained from a combined analysis of $46$ dark siren detections and a single bright siren detection are $H_0=\left\{68.1^{+8.5}_{-5.6},67.5^{+8.3}_{-5.2},67.1^{+6.6}_{-5.8} \right\}~\mathrm{km\ s^{-1}\ Mpc^{-1}}$, respectively. All our constraints on $H_0$ are consistent within one sigma with the latest measurement of $H_0$ applying the differential age method, which provides a constraint on $H_0$ in coasting cosmologies independently from $k$. Our results constrain all cosmological models with $a(t)\propto t$ linear expansion in the luminosity distance and redshift range of the $47$ LIGO-Virgo detections, i.e. $d_\mathrm{L}\lesssim 5~\mathrm{Gpc}$ and $z\lesssim 0.8$, which practically include all (both strictly linear and quasi-linear) models in the coasting model family. As we have found, the coasting models and the $\Lambda$CDM model fit equally well to the applied set of gravitational-wave detections.

Measuring the electron diffusion coefficient is the most straightforward task in the study of gamma-ray pulsar halos. The updated measurements of the spatial morphology and spectrum of the Geminga halo by the HAWC experiment enable us to constrain parameters beyond the diffusion coefficient, including the size of the slow-diffusion zone and the electron injection spectrum from the pulsar wind nebulae (PWN). Based on the two-zone diffusion model, we find that the slow-diffusion zone size ($r_*$) around Geminga is within the range of $30-70$pc. The lower boundary of this range is determined by the goodness of fit of the model to the one-dimensional morphology of the Geminga halo. The upper limit is derived from fitting the gamma-ray spectrum of the Geminga halo, along with the expectations for the power-law index of the injection spectrum based on simulations and PWN observations, i.e., $p\gtrsim1$. With $r_*$ set at its lower limit of $30$~pc, we obtain the maximum $p$ permitted by the HAWC spectrum measurement, with an upper limit of $2.17$ at a $3\sigma$ significance. Moreover, we find that when $r_*=30$pc and $p=2.17$, the predicted positron spectrum generated by Geminga at Earth coincides with the AMS-02 measurement in the $50-500$GeV range.

We aim to compare and contrast the typical shapes of synthetic Ca II 854.2 nm spectra found in Bifrost simulations having different magnetic activity with the spectral shapes found in a quiet Sun observation from the Swedish 1-m Solar Telescope (SST). We use clustering techniques to extract the typical Ca II 854.2 nm profile shapes synthesized from Bifrost simulations with varying amounts of magnetic activity. We degrade the synthetic profiles to observational conditions and repeat the clustering, and we compare our synthetic results with actual observations. While the mean spectra for our high resolution simulations compare reasonably well with the observations, we find that there are considerable differences between the clusters of observed and synthetic intensity profiles, even after the synthetic profiles have been degraded. The typical absorption profiles from the simulations are both narrower and display a steeper transition from the inner wings to the line core. Furthermore, even in our most quiescent simulation we find a far larger fraction of profiles with local emission around the core, or other exotic profile shapes, than in the observations. Looking into the atmospheric structure for a selected set of synthetic clusters, we find distinct differences in the temperature stratification for the clusters most and least similar to the observations. The narrow and steep profiles are associated with either weak gradients in temperature, or temperatures rising to a local maximum in the line wing forming region before sinking to a minimum in the line core forming region. The profiles that display less steep transitions show extended temperature gradients that are steeper in the range $-3 \lesssim \log \tau_{5000} \lesssim -1$.

We present the results of the investigation of the nova V1674 Her (2021), recognised as the swiftest classical nova, with $t_2 \sim 0.90$ days. The distance to the nova is estimated to be 4.97 kpc. The mass and radius of the WD are calculated to be $\sim~1.36~M_\odot$ and $\sim 0.15~R_\oplus$, respectively. Over the course of one month following the outburst, V1674 Her traversed distinct phases -- pre-maxima, early decline, nebular, and coronal -- displaying a remarkably swift transformation. The nebular lines emerged on day 10.00, making it the classical nova with the earliest observed commencement to date. We modelled the observed optical spectrum using the photoionization code \textsc{cloudy}. From the best-fitting model we deduced different physical and chemical parameters associated withe the system. The temperature and luminosity of the central ionizing sources are found in the range of $1.99 - 2.34~\times 10^5$ K and $1.26 - 3.16~ \times 10^{38}$ \ergs, respectively. Elements such as He, O, N, and Ne are found to be overabundant compared to solar abundance in both the nebular and coronal phases. According to the model, Fe II abundance diminishes while Ne abundance increases, potentially elucidating the rare hybrid transition between Fe and He/N nova classes. The ejected mass across all epochs spanned from $3.42 - 7.04~ \times 10^{-5}~M_\odot$. Morpho-kinematic modelling utilising \textsc{shape} revealed that the nova V1674 Her possesses a bipolar structure with an equatorial ring at the centre and an inclination angle of i = 67$\pm$ 1.5$^{\circ}$.

At the beginning of the year 2020, MAGIC reported a very-high-energy (VHE) flaring activity from the FSRQ QSO B1420+326. It is now the fourth known most distant blazar (z=0.682) with an observed VHE gamma-ray emission. In this work, we investigate the effect of photon-axion-like particle (ALP) oscillations in the gamma-ray spectra measured by Fermi-LAT and MAGIC around the flaring state. We set 95% confidence level (C.L.) upper limit on the ALP parameters and obtain a constraint on the photon-ALP coupling constant $g_{a\gamma} < 2\times10^{-11}$ GeV$^{-1}$ for ALP masses $m_{a} \sim 10^{-10}-10^{-9}$ eV. Assuming the hadronic origin of VHE photons, we also estimate the expected neutrino flux from this source and the contribution to diffuse neutrino flux from QSO B1420+326-like FSRQs at sub-PeV energies. Furthermore, we study the implications of photon-ALP oscillations on the counterpart $\gamma$-rays of the sub-PeV neutrinos. Finally, we investigate a viable scenario of invisible neutrino decay to ALPs on the gamma-ray spectra and diffuse $\gamma$-ray flux at sub-PeV energies. Interestingly, we find that for the choice of neutrino lifetime $\tau_{2}/m_{2} = 10^3$ s eV$^{-1}$, the $\gamma$-ray flux has a good observational sensitivity towards LHAASO-KM2A.

Most high-energy constructions that realise a phase of cosmic inflation contain many degrees of freedom. Yet, cosmological observations are all consistent with single-field embeddings. We show how volume selection effects explain this apparent paradox. Due to quantum diffusion, different regions of space inflate by different amounts. In regions that inflate most, and eventually dominate the volume of the universe, a generic mechanism is unveiled that diverts the inflationary dynamics towards single-field attractors. The formalism of constrained stochastic inflation is developed to this end.

We present gas-phase metallicity measurements for 583 emission line galaxies at $0.3<z<0.85$, including 388 dwarf galaxies with $log(M_{\star}/M_{\odot}) < 9.5$, and explore the dependence of the metallicity on the stellar mass and star formation properties of the galaxies. Metallicities are determined through the measurement of emission lines in very deep ($\sim$7 hr exposure) Keck/DEIMOS spectra taken primarily from the HALO7D survey. We measure metallicity with three strong-line calibrations (O3H$\beta$, R23, and O3O2) for the overall sample, as well as with the faint [Ne III]$\lambda$3869 and [O III]$\lambda$4363 emission lines for 112 and 17 galaxies where robust detections were possible. We construct mass-metallicity relations (MZR) for each calibration method, finding MZRs consistent with other strong-line results at comparable redshift, as well as with $z\sim0$ galaxies. We quantify the intrinsic scatter in the MZR as a function of mass, finding it increases with lower stellar mass. We also measure a weak but significant correlation between increased MZR scatter and higher specific star formation rate. We find a weak influence of SFR in the fundamental metallicity relation as well, with an SFR coefficient of $\alpha=0.21$. Finally, we use the flux ratios of the [O II]$\lambda\lambda$3727,3729 doublet to calculate gas electron density in $\sim$1000 galaxies with $log(M_{\star}/M_{\odot}) < 10.5$ as a function of redshift. We measure low electron densities ($n_e\sim25$ cm$^{-3}$) for $z<1$ galaxies, again consistent with $z\approx0$ conditions, but measure higher densities ($n_e\sim100$ cm$^{-3}$) at $z>1$. These results all suggest that there is little evolution in star-forming interstellar medium conditions from $z\sim1$ to $z=0$, confirmed with a more complete sample of low-mass galaxies than has previously been available in this redshift range.

In this paper we report the detection of five strong Gamma-Ray Bursts (GRBs) by the High-Energy Particle Detector (HEPD-01) mounted on board the China Seismo-Electromagnetic Satellite (CSES-01), operational since 2018 on a Sun-synchronous polar orbit at a $\sim$ 507 km altitude and 97$^\circ$ inclination. HEPD-01 was designed to detect high-energy electrons in the energy range 3 - 100 MeV, protons in the range 30 - 300 MeV, and light nuclei in the range 30 - 300 MeV/n. Nonetheless, Monte Carlo simulations have shown HEPD-01 is sensitive to gamma-ray photons in the energy range 300 keV - 50 MeV, even if with a moderate effective area above $\sim$ 5 MeV. A dedicated time correlation analysis between GRBs reported in literature and signals from a set of HEPD-01 trigger configuration masks has confirmed the anticipated detector sensitivity to high-energy photons. A comparison between the simultaneous time profiles of HEPD-01 electron fluxes and photons from GRB190114C, GRB190305A, GRB190928A, GRB200826B and GRB211211A has shown a remarkable similarity, in spite of the different energy ranges. The high-energy response, with peak sensitivity at about 2 MeV, and moderate effective area of the detector in the actual flight configuration explain why these five GRBs, characterised by a fluence above $\sim$ 3 $\times$ 10$^{-5}$ erg cm$^{-2}$ in the energy interval 300 keV - 50 MeV, have been detected.

Since 2018, the potential for a high-energy neutrino telescope, named the Pacific Ocean Neutrino Experiment (P-ONE), has been thoroughly examined by two pathfinder missions, STRAW and STRAW-b, short for short for Strings for Absorption Length in Water. The P-ONE project seeks to install a neutrino detector with a one cubic kilometer volume in the Cascadia Basin's deep marine surroundings, situated near the western shores of Vancouver Island, Canada. To assess the environmental conditions and feasibility of constructing a neutrino detector of that scale, the pathfinder missions, STRAW and STRAW-b, have been deployed at a depth of 2.7 km within the designated site for P-ONE and were connected to the NEPTUNE observatory, operated by Ocean Networks Canada (ONC). While STRAW focused on analyzing the optical properties of water in the Cascadia Basin, \ac{strawb} employed cameras and spectrometers to investigate the characteristics of bioluminescence in the deep-sea environment. This report introduces the STRAW-b concept, covering its scientific objectives and the instrumentation used. Furthermore, it discusses the design considerations implemented to guarantee a secure and dependable deployment process of STRAW-b. Additionally, it showcases the data collected by battery-powered loggers, which monitored the mechanical stress on the equipment throughout the deployment. The report also offers an overview of STRAW-b's operation, with a specific emphasis on the notable advancements achieved in the data acquisition (DAQ) system and its successful integration with the server infrastructure of ONC.

We use a carefully selected subsample of 1053 confirmed exoplanets from the NASA Exoplanet Archive to construct empirical power-law exoplanet mass-radius-temperature ($M$-$R$-$T$) relations. Using orthogonal distance regression to account for errors in both mass and radius, we allow the data to decide: 1) the number of distinct planetary regimes; 2) whether the boundaries of these regimes are best described by broken power laws joined at mass break points, or by discontinuous power laws motivated by changes in equations of state and temperature. We find strong support from the data for three distinct planetary $M$-$R$ regimes and for those regimes to be discontinuous. Our most successful model involves an $M$-$R$-$T$ relation in which ice/rock (rocky) and ice-giant (neptunian) planets are segregated by a pure-ice equation of state, whilst neptunes and gas giant (jovian) planets are segregated by a mass break at $M_{\rm br} = 115\pm19~M_{\oplus}$. The rocky planet regime is shown to follow $M \propto R^{0.34\pm0.01}$, whilst neptunes have $M\propto R^{0.55\pm0.02}$. Planets in both regimes are seen to extend to similar maximum masses. In the jovian regime, we find that $M \propto R^{0.00\pm0.01}T^{0.35\pm 0.02}$, where $T$ is the planet equilibrium temperature. This implies that, for jovian planets detected so far, equilibrium temperature alone provides a robust estimator of mass.

High-resolution solar spectroscopy provides a wealth of information from photospheric and chromospheric spectral lines. However, the volume of data easily exceeds hundreds of millions of spectra on a single observation day. Therefore, methods are needed to identify spectral signatures of interest in multidimensional datasets. Background-subtracted activity maps (BaSAMs) have previously been used to locate features of solar activity in time series of images and filtergrams. This research note shows how this method can be extended and adapted to spectral data.

Inflation is a highly favoured theory for the early Universe. It is compatible with current observations of the cosmic microwave background and large scale structure and is a driver in the quest to detect primordial gravitational waves. It is also, given the current quality of the data, highly under-determined with a large number of candidate implementations. We use a new method in symbolic regression to generate all possible simple scalar field potentials for one of two possible basis sets of operators. Treating these as single-field, slow-roll inflationary models we then score them with an information-theoretic metric ("minimum description length") that quantifies their efficiency in compressing the information in the Planck data. We explore two possible priors on the parameter space of potentials, one related to the functions' structural complexity and one that uses a Katz back-off language model to prefer functions that may be theoretically motivated. This enables us to identify the inflaton potentials that optimally balance simplicity with accuracy at explaining the Planck data, which may subsequently find theoretical motivation. Our exploratory study opens the door to extraction of fundamental physics directly from data, and may be augmented with more refined theoretical priors in the quest for a complete understanding of the early Universe.

In this paper, we re-assess the ability of the acoustic early dark energy (ADE) model to resolve the Hubble tension in light of the new Pantheon+ and S$H_0$ES data on the one hand, and the BOSS LRG and eBOSS QSO data, analysed under the effective field theory of large-scale structures (ETFofLSS) on the other hand. We find that the Pantheon+ data, which favor a larger $\Omega_m$ value than the Pantheon data, have a strong constraining power on the ADE model, while the EFTofLSS analysis of the BOSS and eBOSS data only slightly increases the constraints. We establish that the ADE model is now ruled out as a solution to the Hubble tension, with a remaining tension of $3.6\sigma$. In addition, we find that the axion-like early dark energy model performs better when confronted to the same datasets, with a residual tension of $2.5\sigma$. This work shows that the Pantheon+ data can have a decisive impact on models which aim to resolve the Hubble tension.

In this work, we examine the list of 63 candidates to T Tauri star (TTS) in the TAMC identified by their ultraviolet (UV) and infrared colours (IR) measured from data obtained by the Galaxy Evolution Explorer all sky survey (GALEX-AIS) and the Two Microns All Sky Survey (2MASS), respectively. The objective of this work is twofold: evaluate whether they are pre-main sequence (PMS) stars and evaluate the goodness of the UV-IR colour-colour diagram to detect PMS stars in wide-fields. The astrometric properties of these sources have been retrieved from the Gaia DR3 catalogue and used to evaluate their membership probability. Several classification algorithms have been tested to search for the kinematical groups but the final classification has been made with k-means++ algorithms. Membership probability has been evaluated by applying Logistic Regression. In addition, spectroscopic information available in the archive of the Large Sky Area Multi Object Fiber Spectroscopic Telescope has been used to ascertain their PMS nature when available. About 20% of the candidates share the kinematics of the TAMC members. Among them, HD 281691 is a G8-type field star located in front of the cloud and HO Aur is likely a halo star given the very low metallicity provided by Gaia. The rest are three known PMS stars (HD 30171, V600 Aur and J04590305+3003004), two previously unknown accreting M-type stars (J04510713+1708468 and J05240794+2542438) and, five additional sources, which are very likely PMS stars. Most of these new sources are concentrated at low galactic latitudes over the Auriga-Perseus region.

Rapidly rotating and strongly magnetized protoneutron stars (PNSs) created in core-collapse supernovae can drive relativistic magnetized winds. Ions and neutrons can be co-accelerated while they remain coupled through elastic collisions. We investigate the nucleosynthesis and subsequent nuclear disintegration, and find that relativistic neutrons can be generated in such magnetized winds. Upon eventual decoupling, resulting inelastic collisions with ejecta lead to pion production, resulting in $0.1-10\,{\rm GeV}$ neutrinos. Following this scenario presented in Murase, Dasgupta & Thompson, Phys. Rev. D, 89, 043012 (2014), we numerically calculate the spectra and light curves of quasithermal neutrino emission and find that power-law tails are formed without cosmic-ray acceleration. In the event of a Galactic supernova, $\sim 10-1000$ neutrino events could be detected with Hyper-Kamiokande, KM3Net-ORCA and IceCube-Upgrade for PNSs with surface magnetic field $B_{\rm dip}\sim 10^{13-15}\,{\rm G}$ and initial spin period $P_i \sim 1-30\,{\rm ms}$. Successful detection will enable us to study supernovae as multienergy neutrino sources and may provide clues to the roles of PNSs in diverse classes of transients.

In recent decades, the use of optical detection systems for meteor studies has increased dramatically, resulting in huge amounts of data being analyzed. Automated meteor detection tools are essential for studying the continuous meteoroid incoming flux, recovering fresh meteorites, and achieving a better understanding of our Solar System. Concerning meteor detection, distinguishing false positives between meteor and non-meteor images has traditionally been performed by hand, which is significantly time-consuming. To address this issue, we developed a fully automated pipeline that uses Convolutional Neural Networks (CNNs) to classify candidate meteor detections. Our new method is able to detect meteors even in images that contain static elements such as clouds, the Moon, and buildings. To accurately locate the meteor within each frame, we employ the Gradient-weighted Class Activation Mapping (Grad-CAM) technique. This method facilitates the identification of the region of interest by multiplying the activations from the last convolutional layer with the average of the gradients across the feature map of that layer. By combining these findings with the activation map derived from the first convolutional layer, we effectively pinpoint the most probable pixel location of the meteor. We trained and evaluated our model on a large dataset collected by the Spanish Meteor Network (SPMN) and achieved a precision of 98\%. Our new methodology presented here has the potential to reduce the workload of meteor scientists and station operators and improve the accuracy of meteor tracking and classification.

We examine the cosmological consequences of the heavy quarks in KSVZ-type axion models. We find that their presence often causes an early matter domination phase, altering the evolution of the Universe. This extends the axion mass into the region where standard cosmology leads to overproduction, and allows for a greater number of axion models with non-renormalizable terms to be viable. Quantitatively, we find that decays proceeding through effective terms of up to dimension 9 ($d=9$) remain consistent with cosmological constraints, in contrast with the result $d\leq5$ previously found in the literature. As a consequence, the heavy quarks can be much heavier and the axion mass window with the correct relic density for dark matter is extended by orders of magnitude, down to $m_a\approx 6\times 10^{-9} \,{\rm eV}$. This is achieved without resorting to fine-tuning of the initial misalignment angle, bolstering the motivation for many future axion haloscope experiments. Additionally, we explore how these models can be probed through measurements of the number of relativistic degrees of freedom at recombination.

The first generation of exascale systems will include a variety of machine architectures, featuring GPUs from multiple vendors. As a result, many developers are interested in adopting portable programming models to avoid maintaining multiple versions of their code. It is necessary to document experiences with such programming models to assist developers in understanding the advantages and disadvantages of different approaches. To this end, this paper evaluates the performance portability of a SYCL implementation of a large-scale cosmology application (CRK-HACC) running on GPUs from three different vendors: AMD, Intel, and NVIDIA. We detail the process of migrating the original code from CUDA to SYCL and show that specializing kernels for specific targets can greatly improve performance portability without significantly impacting programmer productivity. The SYCL version of CRK-HACC achieves a performance portability of 0.96 with a code divergence of almost 0, demonstrating that SYCL is a viable programming model for performance-portable applications.

AcubeSAT is an open-source CubeSat mission aiming to explore the effects of microgravity and radiation on eukaryotic cells using a compact microfluidic LoC platform. It is developed by SpaceDot, a volunteer, interdisciplinary student team at the Aristotle University of Thessaloniki and supported by the "Fly Your Satellite! 3" program of the ESA Education Office. The scientific data of the mission is comprised of microscope images captured through the on-board integrated camera setup. As the total size of the payload data is expected to be close to 2GB over 12 months, a fast and efficient downlink fulfilling the restrictive power, cost and complexity budgets is required. Currently, there is no open-source communications system design which fully supports these specific constraints, so we opted to develop our own solutions. The antenna system underwent multiple iterations as the design matured, a process highly aided by the feedback received from the ESA experts. The final communications system configuration consists of an S-band microstrip antenna operating at 2.4GHz and a UHF deployable antenna, for the payload data and TM&TC respectively, both in-house designed. In this paper, we will present AcubeSAT's antenna system iterations that span over 3 years, as well as the rationale and analysis results behind each. The development decisions will be highlighted throughout the paper in an effort to aid in the future development of such a low-cost CubeSat mission communications system.

A very simple non-singular inflationary model is presented where the unique matter content is a radiation fluid. The model slowly contracts from a very large, almost empty and flat spacetime and realizes a bounce. It is then launched to a quasi-de Sitter inflationary expansion with more than sixty e-folds, which smoothly changes to the usual classical, decelerated radiation-dominated expansion before nucleosynthesis. The initial contracting and final expanding phases are classical, but the intermediate bounce and inflationary phases are induced by quantum cosmological effects emerging from a Gaussian wave function quickly moving in configuration space. During this quantum era, a huge number of photons is created. The scale factor describing all this rich evolution is a surprisingly simple analytic function of conformal time. The cosmological scalar perturbations arising from quantum vacuum fluctuations in the far past of the model present an almost scale invariant spectrum with an amplitude compatible with observations for reasonable values of the free parameters of the model.

Right-handed neutrinos ($\nu_{R}$) offer an intriguing portal to new physics in hidden sectors where dark matter (DM) may reside. In this work, we delve into the simplest hidden sector involving only a real scalar exclusively coupled to $\nu_{R}$, referred to as the $\nu_{R}$-philic scalar. We investigate the viability of the $\nu_{R}$-philic scalar to serve as a DM candidate, under the constraint that the coupling of $\nu_{R}$ to the standard model is determined by the seesaw relation and is responsible for the observed DM abundance. By analyzing the DM decay channels and solving Boltzmann equations, we identify the viable parameter space. In particular, our study reveals a lower bound ($\sim10^{4}$ GeV) on the mass of $\nu_{R}$ for the $\nu_{R}$-philic scalar to be DM. The DM mass may vary from sub-keV to sub-GeV. Within the viable parameter space, monochromatic neutrino lines from DM decay can be an important signal for DM indirect detection.

We consider the Brans-Dicke theory in non-metricity gravity, which belongs to the family of symmetric teleparallel scalar-tensor theories. Our focus lies in exploring the implications of the conformal transformation, as we derive the conformal equivalent theory in the Einstein frame, distinct from the minimally coupled scalar field theory. The fundamental principle of the conformal transformation suggests the mathematical equivalence of the related theories. However, to thoroughly analyze the impact on physical variables, we investigate the spatially flat Friedmann--Lema\^{\i}tre--Robertson--Walker geometry, defining the connection in the non-coincidence gauge. We construct exact solutions for the cosmological model in one frame and compare the physical properties in the conformal related frame. Surprisingly, we find that the general physical properties of the exact solutions remain invariant under the conformal transformation. Finally, we construct, for the first time, an analytic solution for the symmetric teleparallel scalar-tensor cosmology.

We consider the formation of primordial black holes (PBHs), during the radiation-dominated Universe, generated from the collapse of super-horizon curvature fluctuations that are overlapped with others on larger scales. Using a set of different curvature profiles, we show that the threshold for PBH formation (defined as the critical peak of the compaction function) can be decreased by several percentages, thanks to the overlapping between the fluctuations. In the opposite case, when the fluctuations are sufficiently decoupled the threshold values behave as having the fluctuations isolated (isolated peaks). We find that the analytical estimates of arXiv:1907.13311 can be used accurately when applied to the corresponding peak that is leading to the gravitational collapse. We also study in detail the dynamics and estimate the final PBH mass for different initial configurations, showing that the profile dependence has a significant effect on that.

This paper presents Dandelion, a new dish antenna experiment searching for dark photons (DPs) with masses around the meV that will start acquiring data by the end of 2023. A spherical mirror acts as a conversion surface between DPs and standard photons that converge to a matrix of 418 Kinetic Inductance Detectors cooled down to 150 mK. A tilt of the mirror at 1 Hz moves the expected signal over the pixels thus enabling a continuous background measurement. The expected signal has two modulations: a spatial modulation providing a directional signature for the unambiguous discovery of a DP, and an intensity modulation allowing the determination of the polarization of the DP. For masses near the meV, the inflationary production of longitudinal and transverse DPs are mutually excluded, thus the polarization determination by Dandelion could shed a new light on the inflation phase of the early universe. A first Dandelion prototype operating for 30 days would improve by more than one order of magnitude the current exclusion limits on DPs at the meV mass scale and would probe this region with an unprecedented discovery potential based on directional detection.

In this paper we propose to apply copula entropy (CE) to photometric redshifts. CE is used to measure the correlations between photometric measurements and redshifts and then the measurements associated with high CEs are selected for predicting redshifts. We verified the proposed method on the SDSS quasar data. Experimental results show that the accuracy of photometric redshifts is improved with the selected measurements compared to the results with all the measurements used in the experiments, especially for the samples with high redshifts. The measurements selected with CE include luminosity magnitude, the brightness in ultraviolet band with standard deviation, and the brightness of the other four bands. Since CE is a rigorously defined mathematical concept, the models such derived is interpretable.

We perform a matched asymptotic expansion to find an analytic formula for the trajectory of a light ray in a Schwarzschild metric, in a power series expansion in the deviation of the impact parameter from its critical value. We present results valid to second sub leading order in this expansion. We use these results to find an analytic expansion for the angular location of the $n^{th}$ Einstein Ring (at large $n$) resulting from a star that lies directly behind a black hole but not necessarily far from it. The small parameter for this expansion is $e^{-\pi (2n+1) }$: our formulae are accurate to third order in this parameter.

Quickly localizing the identified white dwarf (WD) binaries is the basic requirement for the space-based gravitational wave (GW) detection. In fact, the amplitude of GW signals are modulated by the periodic motion of GW detectors on the solar orbit. The intensity of the observed signals is enhanced according to the observation time beyond a year to enhance a high signal to noise ratio (SNR). As data gap exists, the completeness of the data observed for a long time depends on filling gaps in the data. Actually, in a year period, the GW sources have a best observation orbit position of GW detectors, where the detector response intensity of GW is maximum. Thus, the best positions, where the direction of GW source is perpendicular to the detection arms, can be searched for the verified GW sources of the sky map to enhance SNR too. For the three arms response intensity of the GW signals changing more clearly with the location of the GW sources relative to the detector, the noises and the suppression of noise by time delay interferometer are ignored. In the four chosen sources, the two verification WD binaries: J0806 and V407 Vul are observed at the best orbit positions by TAIJI for the short time of 2 and 3 days respectively. The intensities of those GWs are above the values of the TAIJI sensitivity curve, significantly. Compared with a single detector, the network of two detectors does not significantly improve the accuracy of location of the verification binaries. The reason of that result is that one GW source can not be perpendicular to both detectors of TAIJI and LISA. These results imply that the searching of GW signals and parameter estimation of GW sources from the experimental data of the space-based mission do not ignore the orbit positions relevant to GW sources.

There may exist stellar-mass binary black holes (BBH) which merge while orbiting nearby a supermassive black hole (SMBH). In such a triple system, the SMBH will modulate the gravitational waveform of the BBH through orbital Doppler shift and de Sitter precession of the angular momentum. Future space-based GW observatories focused on the milli- and decihertz band will be uniquely poised to observe these waveform modulations, as the GW frequency from stellar-mass BBHs varies slowly in this band while modulation effects accumulate. In this work, we apply the Fisher information matrix formalism to estimate how well space-borne GW detectors can measure properties of BBH+SMBH hierarchical triples using the GW from orbiting BBH. We extend previous work by considering the more realistic case of an eccentric orbit around the SMBH, and notably include the effects of orbital pericenter precession. We find that for detector concepts such as LISA, B-DECIGO, and TianGO, we can extract the SMBH mass and semimajor axis of the orbit with a fractional uncertainty below the 0.1% level over a wide range of triple system parameters. Furthermore, we find that the effects of pericenter precession and orbital eccentricity significantly improve our ability to measure this system. We also find that while LISA could measure these systems, the decihertz detector concepts B-DECIGO and TianGO would enable better sensitivity to the triple's parameters.