Papers for Thursday, Sep 15 2022

Measurements of weak gravitational lensing at low redshifts ($z\lesssim 0.5-1$), quantified by the parameter $S_8$, favor weaker matter clustering than that expected from the standard $\Lambda$CDM cosmological model with parameters determined by cosmic microwave background (CMB) measurements. However, the amplitude of matter clustering at higher redshifts, as probed by lensing of the CMB, is consistent with $\Lambda$CDM. This apparent paradox suggests a connection between the $S_8$ tension and the transition from matter to dark-energy domination. Here we show that the tension can be resolved by introducing a friction between dark matter and dark energy without altering the tightly constrained expansion history. The low-$S_8$ measurements favor (at $\gtrsim3\sigma$, in this one parameter model) a non-zero drag leading to a suppression of low-redshift power right around the transition from matter to dark-energy domination. We suggest ways to further probe the scenario.

We present a comprehensive census of the AGNs in the GOODS-S/HUDF region from the X-ray to the radio, covering both the obscured and unobscured populations. This work includes a robust analysis of the source optical-to-mid-IR SEDs featuring (semi-)empirical AGN and galaxy dust emission models and Baysian fitting techniques, ultra-deep VLA 3 and 6 GHz observations, and an integrated analysis of various AGN selection techniques, including X-ray properties, UV-to-MIR SED analysis, optical spectral features, mid-IR colors, radio loudness and spectral slope, and AGN variability. In total, we report $\sim$900 AGNs over the $\sim$170 arcmin$^2$ 3D-HST GOODS-S footprint, which has doubled the AGN number identified in the previous X-ray sample with $\sim$26\% of our sample undetected in the deepest Chandra image. With a summary of AGN demographics from different selection methods, we find that no one single band or technique comes close to selecting a complete AGN sample despite the great depth of the data in GOODS-S/HUDF. We estimate the yields of various approaches and explore the reasons for incompleteness. We characterize the statistical properties, such as source number density, obscuration fraction and luminosity function of the AGN sample in this field and discuss their immediate implications. We also provide some qualitative predictions of the AGN sample that might be discovered by the upcoming JWST surveys.

Among the multiple stellar populations in globular clusters (GCs) the very high-temperature H-burning regime, able to produce elements up to potassium, is still poorly explored. Here we present the first abundance analysis of K in 42 giants of NGC 6715 (M 54) with homogeneous abundances of light elements previously derived in our FLAMES survey. Owing to the large mass and low metallicity, a large excess of K could be expected in this GC, which is located in the nucleus of the Sagittarius dwarf galaxy. We actually found a spread in [K/Fe] spanning about 1 dex, with [K/Fe] presenting a significant anti-correlation with [O/Fe] ratios, regardless of the metallicity component in M 54. Evidence for a K-Mg anti-correlation also exists, but this is statistically marginal because of the lack of very Mg-poor stars in this GC. We found, however, a strong correlation between K and Ca. These observations clearl y show that the K enhancement in M 54 is probably due to the same network of nuclear reactions generating the phenomenon of multiple stellar populations, at work in a regime of very high temperature. The comparison with recent results in omega Cen is hampered by an unexplained trend with the temperatures for K abundances from optical spectroscopy, and somewhat by a limited sample size for infrared APOGEE data. There are few doubts, however, that the two most massive GCs in the Milky Way host a K-Mg anti-correlation.

We investigate impacts of long-wavelength gravitational waves (GWs) on nonlinear structure formation by utilizing the tidal separate universe simulations. Based on the equivalence of a long-wavelength GW to a uniform tidal field in a local frame, we provide a way to incorporate a long-wavelength GW into the tidal separate universe simulation as an effective anisotropic expansion. This methodology enables us to study effects of GWs on large-scale structure efficiently. We measure the anisotropic imprint in the local power spectrum from the tidal separate universe simulations with GWs, which corresponds to the scalar-scalar-tensor bispectrum in squeezed limit or the so-called power spectrum response to GWs. We also detect the halo tidal bias induced by GWs from the response of the halo-matter cross-power spectrum to GWs, as well as the linear shape bias (or the linear alignment coefficient) induced by GWs from the one-point function of the halo ellipticity. In contrast to the case of the tidal field induced by scalar perturbations, we discover that the wavenumber dependence of the temporal evolution of GWs naturally causes these biases to be scale-dependent. We also find that this scale dependence is well approximated by the second-order density induced by the coupling between scalar and tensor perturbation. This highlights that the structure formation, especially the process to determine the halo shape, is nonlocal in time. Our findings lay the foundation for predicting the impact of GWs on large-scale structure.

Blazars are active galactic nuclei that launch collimated, powerful jets of magnetized relativistic plasma. Their primary jet, whose emission typically spans from low-frequency radio to very high-energy ($\gtrsim0.1$ TeV) $\gamma$-rays (Blandford et al., 2019), is aligned towards our line of sight. Multiwavelength polarization is a crucial probe of the magnetic field structure and emission processes in such jets. Until now, sensitive polarization observations have been limited to the radio, infrared, and optical range, thereby leaving a gap in our knowledge of the physical conditions experienced by the most energetic particles. Here, we report the first-ever detection of X-ray polarization from the jet in an accreting supermassive black hole system, the blazar Markarian 501 (Mrk 501). The recently launched Imaging X-ray Polarimetry Explorer ($IXPE$, Weisskopf et al., 2022) measures a linear polarization degree ($\Pi$) over the 2-8 keV X-ray energy range of 10$\pm$2% with an electric vector position angle of 134$^\circ\pm$5$^\circ$, parallel to the radio jet. The X-ray $\Pi$ is more than a factor of 2 higher than the optical $\Pi$. We conclude that an energy-stratified relativistic electron population, i.e., an acceleration scenario where the higher energy particles emit from more magnetically ordered regions closer to the acceleration site, is the most likely explanation of the higher degree of polarization at X-ray energies. A second $IXPE$ observation conducted 16 days later yielded similar results, strengthening our conclusions.

We present comprehensive multi-component dynamical models of M54 (NGC6715), the nuclear star cluster of the Sagittarius dwarf galaxy (Sgr), which is undergoing a tidal disruption in the Milky Way halo. Previous papers in the series used a large MUSE mosaic data set to identify multiple stellar populations in the system and study their kinematic differences. Here we use Jeans-based dynamical models that fit the population properties (mean age and metallicity), spatial distributions, and kinematics simultaneously. They provide a solid physical explanation to our previous findings. The population-dynamical models deliver a comprehensive view of the whole system, and allow us to disentangle the different stellar populations. We explore their dynamical interplay and confirm our previous findings about the build-up of Sgr's nuclear cluster via contributions from globular cluster stars, Sgr inner field stars, and in-situ star formation. We explore various parameterisations of the gravitational potential and show the importance of a radially varying mass-to-light ratio for the proper treatment of the mass profile. We find a total dynamical mass within M54's tidal radius ($\sim75$ pc) of $1.60\pm0.07\times10^6$ Msun in excellent agreement with $N$-body simulations. The metal-poor globular cluster stars contribute about $65\%$ of the total mass or $1.04\pm0.05\times10^6$ Msun. The metal-rich stars can be further divided into young and intermediate age populations that contribute $0.32\pm0.02\times10^6$ Msun ($20\%$) and $0.24\pm0.02\times10^6$ Msun ($15\%$), respectively. Our population-dynamical models successfully distinguish the different stellar populations in Sgr's nucleus because of their different spatial distributions, ages, metallicities, and kinematic features.

The detection of a $0.2\,M_\odot$ extremely low-mass white dwarf (hereafter, EW) in a wide orbit ($P_{\rm orb}\approx450$ days) with a $1.1\,M_\odot$ main-sequence (MS) companion KIC 8145411 challenges our current understanding of how EWs form. The traditional channel for EW formation via mass transfer from the WD progenitor is expected to form EW binaries in tight orbits. Indeed, majority of known EWs are found in tight binaries with a median $P_{\rm orb}\approx 5.4$ hrs. Using numerical scattering experiments, we find that binary-binary strong encounters in star clusters can sufficiently widen the orbit of a typical EW binary to explain the observed wide orbit of the KIC 8145411 system. The $P_{\rm orb}$ distribution for EW binaries produced through binary-binary encounters is bimodal: one mode corresponds to the initial orbital period of the EW binary, while the other is near $P\sim$ few $10^2$ days, similar to the orbital period of the KIC 8145411 system. We find that the production of wide EW binaries that are also ejected from the cluster peaks at a star clusters mass of $\sim10^5\,M_\odot$ with a rate of $\sim10^{-3}\,\rm{Gyr^{-1}}$. Assuming that $50\%$ of all stars form in star clusters and an initial cluster mass function $\propto m^{-2}$, we estimate a galactic formation rate of $\sim3.64\times10^3\,\rm{Gyr^{-1}}$ for wide EW binaries.

We present a formulation and numerical algorithm to extend the scheme for grey radiation magneto-hydrodynamics (MHD) developed by Jiang (2021) to include the frequency dependence via the multi-group approach. The entire frequency space can be divided into arbitrary number of groups in the lab frame, and we follow the time dependent evolution of frequency integrated specific intensities along discrete rays inside each group. Spatial transport of photons is done in the lab frame while all the coupling terms are solved in the fluid rest frame. Lorentz transformation is used to connect different frames. Radiation transport equation is solved fully implicitly in time while the MHD equations are evolved explicitly so that time step is not limited by the speed of light. A finite volume approach is used for transport in both spatial and frequency spaces to conserve radiation energy density and momentum. The algorithm includes photon absorption, electron scattering as well as Compton scattering, which is calculated by solving the Kompaneets equation. The algorithm is accurate for a wide range of optical depth conditions and can handle both radiation pressure and gas pressure dominated flows. It works for both Cartesian and curvilinear coordinate systems with adaptive mesh refinement. We provide a variety of test problems including radiating sphere, shadow test, absorption of a moving gas, Bondi type flows as well as a collection of test problems for thermal and bulk Compton scattering. We also discuss examples where frequency dependence can make a big difference compared with the grey approach.

We analyze HCN and HNC emission in the nearby starburst galaxy NGC 253 to investigate its effectiveness in tracing heating processes associated with star formation. This study uses multiple HCN and HNC rotational transitions observed using ALMA via the ALCHEMI Large Program. To understand the conditions and associated heating mechanisms within NGC 253's dense gas, we employ Bayesian nested sampling techniques applied to chemical and radiative transfer models which are constrained using our HCN and HNC measurements. We find that the volume density $n_{\text{H}_{2}}$ and cosmic ray ionization rate (CRIR) $\zeta$ are enhanced by about an order of magnitude in the galaxy's central regions as compared to those further from the nucleus. In NGC 253's central GMCs, where observed HCN/HNC abundance ratios are lowest, $n \sim 10^{5.5}$ cm$^{-3}$ and $\zeta \sim 10^{-12}$ s$^{-1}$ (greater than $10^4$ times the average Galactic rate). We find a positive correlation in the association of both density and CRIR with the number of star formation-related heating sources (supernova remnants, HII regions, and super hot cores) located in each GMC, as well as a correlation between CRIRs and supernova rates. Additionally, we see an anticorrelation between the HCN/HNC ratio and CRIR, indicating that this ratio will be lower in regions where $\zeta$ is higher. Though previous studies suggested HCN and HNC may reveal strong mechanical heating processes in NGC 253's CMZ, we find cosmic ray heating dominates the heating budget, and mechanical heating does not play a significant role in the HCN and HNC chemistry.

We are conducting a survey using twilight time on the Dark Energy Camera with the Blanco 4m telescope in Chile to look for objects interior to Earth's and Venus' orbits. To date we have discovered two rare Atira/Apohele asteroids, 2021 LJ4 and 2021 PH27, which have orbits completely interior to Earth's orbit. We also discovered one new Apollo type Near Earth Object (NEO) that crosses Earth's orbit, 2022 AP7. Two of the discoveries likely have diameters greater than 1 km. 2022 AP7 is likely the largest Potentially Hazardous Asteroid (PHA) discovered in about eight years. To date we have covered 624 square degrees of sky near to and interior to the orbit of Venus. The average images go to 21.3 mags in the r-band, with the best images near 22nd mag. Our new discovery 2021 PH27 has the smallest semi-major axis known for an asteroid, 0.4617 au, and the largest general relativistic effects (53 arcseconds/century) known for any body in the Solar System. The survey has detected about 15 percent of all known Atira NEOs. We put strong constraints on any stable population of Venus co-orbital resonance objects existing, as well as the Atira and Vatira asteroid classes. These interior asteroid populations are important to complete the census of asteroids near Earth, including some of the most likely Earth impactors that cannot easily be discovered in other surveys. Comparing the actual population of asteroids found interior to Earth and Venus with those predicted to exist by extrapolating from the known population exterior to Earth is important to better understand the origin, composition and structure of the NEO population.

We analyze a sample of 25 [Ne V] (${\lambda}$3426 $\r{A}$) emission-line galaxies at 1.4 < z < 2.3 using Hubble Space Telescope/Wide Field Camera 3 G102 and G141 grism observations from the CANDELS Lyman-${\alpha}$ Emission at Reionization (CLEAR) survey. [Ne V] emission probes extremely energetic photoionization (creation potential of 97.11 eV), and is often attributed to energetic radiation from active galactic nuclei (AGN), radiative supernova shocks, or an otherwise very hard ionizing spectrum from the stellar continuum. In this work, we use [Ne V] in conjunction with other rest-frame UV/optical emission lines ([O II] ${\lambda\lambda}$3726,3729 $\r{A}$, [Ne III] ${\lambda}$3869 $\r{A}$, [O III] ${\lambda\lambda}$4959,5007 $\r{A}$, H${\alpha}$+[N II] ${\lambda\lambda}$6548,6583 $\r{A}$, [S II] ${\lambda\lambda}$6716,6731 $\r{A}$), deep (2-7 Ms) X-ray observations (from Chandra), and mid-infrared imaging (from Spitzer) to study the origin of this emission and to place constraints on the nature of the ionizing engine. The majority of the [Ne V]-detected galaxies have properties consistent with ionization from AGN. However, for our [Ne V]-selected sample, the X-ray luminosities are consistent with local (z < 0.1) X-ray-selected Seyferts, but the [Ne V] luminosities are more consistent with those from z ~ 1 X-ray-selected QSOs. The excess [Ne V] emission requires either reduced hard X-rays, or a ~0.1 keV excess. We discuss possible origins of the apparent [Ne V] excess, which could be related to the "soft (X-ray) excess" observed in some QSOs and Seyferts, and/or be a consequence of a complex/anisotropic geometry for the narrow line region, combined with absorption from a warm, relativistic wind ejected from the accretion disk. We also consider implications for future studies of extreme high-ionization systems in the epoch of reionization (z > 6) with the James Webb Space Telescope.

Arbitrariness in the zero point of bolometric corrections is a nearly century-old paradigm leading to two more paradigms. "Bolometric corrections must always be negative," and "bolometric magnitude of a star ought to be brighter than its $V$ magnitude". Both were considered valid before IAU 2015 General Assembly Resolution B2, a revolutionary document that supersedes all three aforementioned paradigms. The purpose of this article is to initiate a new insight and a new understanding of the fundamental astrophysics and present new capabilities to obtain standard and more accurate stellar luminosities and gain more from accurate observations in the era after Gaia. The accuracy gained will aid in advancing stellar structure and evolution theories, and also Galactic and extragalactic research, observational cosmology and dark matter and dark energy searches.

Solar pores are efficient magnetic conduits for propagating magnetohydrodynamic wave energy into the outer regions of the solar atmosphere. Pore observations often contain isolated and/or unconnected structures, preventing the statistical examination of wave activity as a function of atmospheric height. Here, using high resolution observations acquired by the Dunn Solar Telescope, we examine photospheric and chromospheric wave signatures from a unique collection of magnetic pores originating from the same decaying sunspot. Wavelet analysis of high cadence photospheric imaging reveals the ubiquitous presence of slow sausage mode oscillations, coherent across all photospheric pores through comparisons of intensity and area fluctuations, producing statistically significant in-phase relationships. The universal nature of these waves allowed an investigation of whether the wave activity remained coherent as they propagate. Utilizing bi-sector Doppler velocity analysis of the Ca II 8542 {\AA} line, alongside comparisons of the modeled spectral response function, we find fine-scale 5 mHz power amplification as the waves propagate into the chromosphere. Phase angles approaching zero degrees between co-spatial bi-sectors spanning different line depths indicate standing sausage modes following reflection against the transition region boundary. Fourier analysis of chromospheric velocities between neighboring pores reveals the annihilation of the wave coherency observed in the photosphere, with examination of the intensity and velocity signals from individual pores indicating they behave as fractured wave guides, rather than monolithic structures. Importantly, this work highlights that wave morphology with atmospheric height is highly complex, with vast differences observed at chromospheric layers, despite equivalent wave modes being introduced into similar pores in the photosphere.

Massive stars are predominantly found in binaries and higher order multiples. While the period and eccentricity distributions of OB stars are now well established across different metallicity regimes, the determination of mass-ratios has been mostly limited to double-lined spectroscopic binaries. As a consequence, the mass-ratio distribution remains subject to significant uncertainties. Open questions include the shape and extent of the companion mass-function towards its low-mass end and the nature of undetected companions in single-lined spectroscopic binaries. In this contribution, we present the results of a large and systematic analysis of a sample of over 80 single-lined O-type spectroscopic binaries (SB1s) in the Milky Way and in the Large Magellanic Cloud (LMC). We report on the developed methodology, the constraints obtained on the nature of SB1 companions, the distribution of O star mass-ratios at LMC metallicity and the occurrence of quiescent OB+black hole binaries.

Efficient radiation at second and/or higher harmonics of Wce has been suggested to circumvent the escaping difficulty of the electron cyclotron maser emission mechanism when it is applied to solar radio bursts, such as spikes. In our earlier study, we developed a three-step numerical scheme to connect the dynamics of energetic electrons within a large-scale coronal loop structure with the microscale kinetic instability energized by the obtained nonthermal velocity distribution and found that direct and efficient harmonic X-mode (X2 for short) emission can be achieved due to the strip-like features of the distribution. That study only considered the radiation from the loop top at a specific time. Here we present the emission properties along the loop at different locations and timings. We found that, in accordance with our earlier results, few to several strip-like features can appear in all cases, and the first two strips play the major role in exciting X2 and Z (i.e., the slow extraordinary mode) that propagate quasi-perpendicularly. For the four sections along the loop, significant excitation of X2 is observed from the upper two sections, and the strongest emission is from the top section. In addition, significant excitation of Z is observed for all loop sections, while there is no significant emission of the fundamental X mode. The study provides new insight into coherent maser emission along the coronal loop structure during solar flares.

We present the photometry of 16 91T/99aa-like Type Ia Supernovae (SNe Ia) observed by the Las Cumbres Observatory. We also use an additional set of 21 91T/99aa-like SNe Ia and 87 normal SNe Ia from the literature for an analysis of the standardizability of the luminosity of 91T/99aa-like SNe. We find that 91T/99aa-like SNe are 0.2 mag brighter than normal SNe Ia, even when fully corrected by the light curve shapes and colors. The weighted root-mean-square of 91T/99aa-like SNe (with $z_{CMB}>0.01$) Hubble residuals is $0.25\pm0.03$ mag, suggesting that 91T/99aa-like SNe are also excellent relative distance indicators to $\pm$12%. We compare the Hubble residuals with the pseudo-equivalent width (pEW) of Si II $\lambda\lambda$6355 around the date of maximum brightness. We find that there is a broken linear correlation in between those two measurements for our sample including both 91T/99aa-like and normal SNe Ia. As the $pEW_{max}$(Si II $\lambda\lambda$6355) increasing, the Hubble residual increases when $pEW_{max}$(Si II $\lambda\lambda$6355)$<55.6$ \r{A}. However, the Hubble residual stays constant beyond this. Given that 91T/99aa-like SNe possess shallower Si II lines than normal SNe Ia, the linear correlation at $pEW_{max}$(Si II $\lambda\lambda$6355)$<55.6$ \r{A} can account for the overall discrepancy of Hubble residuals derived from the two subgroups. Such a systematic effect needs to be taken into account when using SNe Ia to measure luminosity distances.

Mass loss is a key property to understand stellar evolution and in particular for low-metallicity environments. Our knowledge has improved dramatically over the last decades both for single and binary evolutionary models. However, episodic mass loss although definitely present observationally, is not included in the models, while its role is currently undetermined. A major hindrance is the lack of large enough samples of classified stars. We attempted to address this by applying an ensemble machine-learning approach using color indices (from IR/Spitzer and optical/Pan-STARRS photometry) as features and combining the probabilities from three different algorithms. We trained on M31 and M33 sources with known spectral classification, which we grouped into Blue/Yellow/Red/B[e] Supergiants, Luminous Blue Variables, classical Wolf-Rayet and background galaxies/AGNs. We then applied the classifier to about one million Spitzer point sources from 25 nearby galaxies, spanning a range of metallicites ($1/15$ to $\sim3~Z_{\odot}$). Equipped with spectral classifications we investigated the occurrence of these populations with metallicity.

We analyze the eccentric response of a low mass coplanar circumbinary disc to secular tidal forcing by a Keplerian eccentric orbit central binary. The disc acquires a forced eccentricity whose magnitude depends on the properties of the binary and disc. The largest eccentricities occur when there is a global apsidal resonance in the disc. The driving frequency by the binary is its apsidal frequency that is equal to zero. A global resonance occurs when the disc properties permit the existence of a zero apsidal frequency free eccentric mode. Resonances occur for different free eccentric modes that differ in the number of radial nodes. For a disc not at resonance, the eccentricity distribution has somewhat similar form to the eccentricity distributions in discs at resonance that have the closest matching disc aspect ratios. For higher disc aspect ratios, the forced eccentricity distribution in a 2D disc is similar to that of the fundamental free mode. The forced eccentricity distribution in a 3D disc is similar to that of higher order free modes, not the fundamental mode, unless the disc is very cool. For parameters close to resonance, large phase shifts occur between the disc and binary eccentricities that are locked in phase. Forced eccentricity may play an important role in the evolution of circumbinary discs and their central binaries.

Strong solar flares occur in $\delta $-spots characterized by the opposite-polarity magnetic fluxes in a single penumbra. Sunspot formation via flux emergence from the convection zone to the photosphere can be strongly affected by convective turbulent flows. It has not yet been shown how crucial convective flows are for the formation of $\delta$-spots. The aim of this study is to reveal the impact of convective flows in the convection zone on the formation and evolution of sunspot magnetic fields. We simulated the emergence and transport of magnetic flux tubes in the convection zone using radiative magnetohydrodynamics code R2D2. We carried out 93 simulations by allocating the twisted flux tubes to different positions in the convection zone. As a result, both $\delta$-type and $\beta$-type magnetic distributions were reproduced only by the differences in the convective flows surrounding the flux tubes. The $\delta $-spots were formed by the collision of positive and negative magnetic fluxes on the photosphere. The unipolar and bipolar rotations of the $\delta$-spots were driven by magnetic twist and writhe, transporting magnetic helicity from the convection zone to the corona. We detected a strong correlation between the distribution of the nonpotential magnetic field in the photosphere and the position of the downflow plume in the convection zone. The correlation could be detected $20$-$30$ h before the flux emergence. The results suggest that high free energy regions in the photosphere can be predicted even before the magnetic flux appears in the photosphere by detecting the downflow profile in the convection zone.

We present an experimental study of photon statistics for high-contrast imaging with the Microwave Kinetic Inductance Detector (MKID) Exoplanet Camera (MEC) located behind the Subaru Coronagraphic Extreme Adaptive Optics System (SCExAO) at the Subaru Telescope. We show that MEC measures the expected distributions for both on-axis companion intensity and off-axis intensity which manifests as quasi-static speckles in the image plane and currently limits high-contrast imaging performance. These statistics can be probed by any MEC observation due to the photon-counting capabilities of MKID detectors. Photon arrival time statistics can also be used to directly distinguish companions from speckles using a post-processing technique called Stochastic Speckle Discrimination (SSD). Here, we we give an overview of the SSD technique and highlight the first demonstration of SSD on an extended source -- the protoplanetary disk AB Aurigae. We then present simulations that provide an in-depth exploration as to the current limitations of an extension of the SSD technique called Photon-Counting SSD (PCSSD) to provide a path forward for transitioning PCSSD from simulations to on-sky results. We end with a discussion of how to further improve the efficacy of such arrival time based post-processing techniques applicable to both MKIDs, as well as other high speed astronomical cameras.

Decayless kink oscillations of solar coronal loops are studied in terms of a low-dimensional model based on a randomly driven Rayleigh oscillator with coefficients experiencing random fluctuations. The model considers kink oscillations as natural modes of coronal loops, decaying by linear resonant absorption. The damping is counteracted by random motions of the loop footpoints and the interaction of the loop with external quasi-steady flows with random fluctuations. In other words, the model combines the self-oscillatory and randomly driven mechanisms for the decayless behaviour. The random signals are taken to be of the stationary red noise nature. In the noiseless case, the model has an asymptotically stationary oscillatory solution, i.e., a kink self-oscillation. It is established that the kink oscillation period is practically independent of noise. This finding justifies the seismological estimations of the kink and Alfv\'en speeds and the magnetic field in an oscillating loop by kink oscillations, based on the observed oscillation period. The oscillatory patterns are found to be almost harmonic. Noisy fluctuations of external flows modulate the amplitude of the almost monochromatic oscillatory pattern symmetrically, while random motions of the loop footpoints cause antisymmetric amplitude modulation. Such modulations are also consistent with the observed behaviour.

Magnetic fields can drastically change predictions of evolutionary models of massive stars via mass-loss quenching, magnetic braking, and efficient angular momentum transport, which we aim to quantify in this work. We use the MESA software instrument to compute an extensive main-sequence grid of stellar structure and evolution models, as well as isochrones, accounting for the effects attributed to a surface fossil magnetic field. The grid is densely populated in initial mass (3-60 M$_\odot$), surface equatorial magnetic field strength (0-50 kG), and metallicity (representative of the Solar neighbourhood and the Magellanic Clouds). We use two magnetic braking and two chemical mixing schemes and compare the model predictions for slowly-rotating, nitrogen-enriched ("Group 2") stars with observations in the Large Magellanic Cloud. We quantify a range of initial field strengths that allow for producing Group 2 stars and find that typical values (up to a few kG) lead to solutions. Between the subgrids, we find notable departures in surface abundances and evolutionary paths. In our magnetic models, chemical mixing is always less efficient compared to non-magnetic models due to the rapid spin-down. We identify that quasi-chemically homogeneous main sequence evolution by efficient mixing could be prevented by fossil magnetic fields. We recommend comparing this grid of evolutionary models with spectropolarimetric and spectroscopic observations with the goals of i) revisiting the derived stellar parameters of known magnetic stars, and ii) observationally constraining the uncertain magnetic braking and chemical mixing schemes.

Classical Be stars, regardless of spectral subtype, display multiperiodic light modulations in the frequency range 0.1-12 d$^{-1}$, when observed with high-cadence and long duration. This behaviour is attributed to non-radial pulsations and/or rotation of the Be star. The main goal of this work is to investigate the fast photometric variability of the optical counterparts to Be/X-ray binaries and compare the general patterns of such variability with the Galactic population of classical Be stars. We analyzed 21 sources with TESS. High-cadence photometry with two ground-based telescopes was also performed for 15 sources. Standard Fourier analysis and least-squares fitting methods were employed in the frequency analysis. All sources exhibit intra-night light variations with intensity variations of 0.01-0.06 mag in the ground-based observations and up to 5% in flux in TESS observations. This variability manifests itself as multi-periodic signals in the frequency range 0.2-12 d$^{-1}$. We find that the patterns of variability of the Be stars in Be/X-ray binaries agree with that of classical early-type Be stars in terms of the general shape of the periodograms. Based on the general shape and number of peaks in the periodograms, Be/X-ray binaries can be classified into different types. The most common case is the presence of groups of closely-spaced frequencies (67%), followed by sources that exhibit isolated signals (18%). The remaining type of sources displays frequency spectra characterized by a mixed pattern of stochastic variability and high-frequency peaks. This study reveals that short-term optical photometric variability is a very common, if not ubiquitous, feature intrinsic to the Be optical companions in Be/X-ray binaries. This variability is mainly attributed to pulsations that originate in the stellar interior.

H_2O submillimeter emission is a powerful diagnostic of the molecular interstellar medium in a variety of sources, including low- and high-mass star forming regions of the Milky Way, and from local to high redshift galaxies. However, the excitation mechanism of these lines in galaxies has been debated, preventing a basic consensus on the physical information that H_2O provides. Both radiative pumping due to H_2O absorption of far-infrared photons emitted by dust and collisional excitation in dense shocked gas have been proposed to explain the H_2O emission. Here we propose two basic diagnostics to distinguish between the two mechanisms: 1) in shock excited regions, the ortho-H_2O 3_{21}-2_{12} 75um and the para-H_2O 2_{20}-1_{11} 101um rotational lines are expected to be in emission while, if radiative pumping dominates, both far-infrared lines are expected to be in absorption; 2) based on statistical equilibrium of H_2O level populations, the radiative pumping scenario predicts that the apparent isotropic net rate of far-infrared absorption in the 3_{21}-2_{12} (75um) and 2_{20}-1_{11} (101um) lines should be higher than or equal to the apparent isotropic net rate of submillimeter emission in the 3_{21}-3_{12} (1163 GHz) and 2_20-2_{11} (1229 GHz) lines, respectively. Applying both criteria to all 16 galaxies and several galactic high-mass star-forming regions where the H_2O 75um and submillimeter lines have been observed with Herschel/PACS and SPIRE, we show that in most (extra)galactic sources the H_2O submillimeter line excitation is dominated by far-infrared pumping, with collisional excitation of the low-excitation levels in some of them. Based on this finding, we revisit the interpretation of the correlation between the luminosity of the H_2O 988 GHz line and the source luminosity in the combined galactic and extragalactic sample.

One of the possible mechanisms for heating the solar atmosphere is the magnetic reconnection occurring at different spatio-temporal scales. The discovery of fast bursty nanojets due to reconnection in the coronal loops has been linked to nanoflares and considered as possible mechanism for coronal heating. The occurrence of these jets mostly in the direction inwards to the loop were observed in the past. In this study, we report ten reconnection nanojets, four with directions inward while six moving outward to the loop, in observations from High-resolution Coronal Imager 2.1 (Hi-C 2.1) and Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO). We determined the maximum length, spire width, speed, and lifetimes of these jets and studied their correlations. We found that outward moving jets with higher speeds are longer in length and duration while the inward moving jets show opposite behaviour. Average duration of the outward jets is $\approx$42 s and inwards jets is $\approx$24s. We identified jets with subsonic speeds below 100 km s$^{-1}$ to high-speed over 150 km s$^{-1}$. These jets can be identified in multiple passbands of AIA extending from upper transition region to the corona suggesting their multi-thermal nature.

In this paper, we report an interesting kinematic phenomenon around the halos' edge related to the `splashback' radius. After the shell-crossing, cosmic flow exhibits various rotational morphologies via stream-mixing. Vorticity is generated in a particular way that coincides with the large-scale structure. Notably, one specific flow morphology concentrates around halos. A detailed examination reveals a sharp change in the logarithmic derivative of its volume fraction at the location of the splashback radius defined as the outermost caustic structure. Such a feature encodes valuable phase space information and provides a new perspective on understanding the dynamical evolution of halos. As a volume-weighted quantity, the profile of flow morphology is purely kinematic. Unlike other related studies, the rotational flow morphologies capture the anisotropic phase structure in the multi-stream region.

Aims. Coronal Mass Ejections (CMEs) are the most fascinating explosion in the solar system; however, their formation is still not fully understood. Methods. Here, we investigate a well-observed CME on 2021 May 07 that showed a typical three-component structure and was continuously observed from 0 to 3 Rsun by a combination of SDO/AIA (0--1.3 Rsun), PROBA2/SWAP (0--1.7 Rsun) and MLSO/K-Cor (1.05--3 Rsun). Furthermore, we compare the morphological discrepancy between the CME white-light bright core and EUV blob. In the end, we explore the origin of various radio bursts closely related to the interaction of the CME overexpansion with nearby streamer. Results. An interesting finding is that the height increases of both the CME leading front and bright core are dominated by the overexpansion during the CME formation. The aspect ratios of the CME bubble and bright core, quantifying the overexpansion, are found to decrease as the SO/STIX 4--10 keV and GOES 1--8 A soft X-ray flux of the associated flare increases near the peaks, indicating an important role of the flare reconnection in the first overexpansion. The CME bubble even takes place a second overexpansion although relatively weak, which is closely related to the compression with a nearby streamer and likely arises from an ideal MHD process. Moreover, the CME EUV blob is found to be relatively lower and wider than the CME white-light bright core, may correspond to the bottom part of the growing CME flux rope. The interaction between the CME and the streamer leads to two type II radio bursts, one normally drifting and one stationary, which are speculated to be induced at two different sources of the CME-driven shock front. The bidirectional electrons evidenced by series of "C-shaped" type III bursts suggest that the interchange reconnection be also involved during the interaction of the CME and streamer.

Cosmic dust models are key ingredients in advancing our understanding of astronomical environments as diverse as interstellar clouds in galaxies, circumstellar envelopes around evolved and young stars, and protoplanetary disks. Such models consist of several dust populations, each with different compositions and size distributions. They may also consider different grain shapes, although most models assume spherical grains. All include a component of silicate dust. The absorption and emission properties of these dust components are calculated from the optical constants of each dust material which have various experimental, phenomenological, and theoretical origins depending on the model. We aim to provide the community with new sets of optical constants for amorphous silicate dust analogues at low temperatures. The analogues consist of four Mg-rich silicate samples of stoichiometry ranging from enstatite to olivine, and of eight samples of Mg and Fe rich silicates with a pyroxene stoichiometry and differing magnesium and iron content. We calculated the optical constants from transmission measurements using the Kramers-Kronig relations, assuming that the grains are small compared to the wavelength and prolate in shape with axis ratios of 1.5 and 2 for the Mg and Fe rich samples, respectively. New optical constants for silicate dust analogues were calculated over the wavelength range from 5 to 800-1000 microns, depending on the sample, and at temperatures of 10, 30, 100, 200, and 300 K.We determined the uncertainties on the derived optical constants based on the assumptions used to calculate them. To facilitate the use of these data in cosmic dust models, we provide optical constants extrapolated outside the measured spectral range into the UV-NIR and mm-cm wavelength ranges, as well as formulae that can be used to interpolate them at any temperature in the range 10 - 300 K.

The classification of gamma-ray-detected blazar candidates of uncertain type (BCU) is a relevant problem in extragalactic gamma-ray astronomy. Here we report the optical spectroscopic characterization, using two 3-4~m class telescopes, Telescopio Nazionale Galileo and Devasthal Optical Telescope, of 27 BCUs detected with the Fermi Large Area Telescope. Since the identification of emission lines is easier in broad-line blazars, which usually exhibit low frequency peaked (synchrotron peak frequency $\leqslant10^{14}$ Hz) spectral energy distribution, we primarily target such BCUs. We found that 8 out of 27 sources exhibit broad emission lines in their optical spectra, 3 of them have redshifts $>$1 and the farthest one is at $z=2.55$. The optical spectra of 2 of the 19 remaining objects are dominated by the absorption spectra of the host galaxy, and there is a tentative detection of the Lyman-$\alpha$ absorption feature in one source. The spectra of the remaining 16 objects, on the other hand, are found to be featureless.

We investigate accretion onto a central star, with the size, rotation rate, and magnetic dipole of a young stellar object, to study the flow pattern (velocity and density) of the fluid within and outside of the disc. We perform resistive MHD simulations of thin $\alpha$-discs, varying the parameters such as the stellar rotation rate and magnetic field, and (anomalous) coefficients of viscosity and resistivity in the disc. To provide a benchmark for the results and to compare with known analytic results, we also perform purely hydrodynamic simulations (HD) for the same problem. Although obtained for a different situation with differing inner boundary condition, the disc structure in the HD simulations closely follows the analytic solution of Klu\'zniak and Kita (2000) -- in particular a region of "midplane" backflow exists in the right range of radii, depending on the viscosity parameter. In the MHD solutions, whenever the magnetic Prandtl number does not exceed a certain critical value, the midplane backflow exists throughout the accretion disc, extending all the way down to the inner transition zone where the disc transitions to a magnetic funnel flow. For values of the magnetic Prandtl number close to the critical value the backflow and the inner disc undergo a quasiperiodic radial oscillation, otherwise the backflow is steady, as is the disc solution. From our results, supplemented by our reading of the literature, we conclude that midplane backflow is a real feature of at least some accretion discs, whether HD $\alpha$-discs or MHD discs, including ones driven by MRI turbulence.

Cosmic messengers (gamma rays, cosmic rays, neutrinos and gravitational waves) provide a powerful complementary way to search for Lorentz invariance violating effects to laboratory-based experiments. The long baselines and high energies involved make Cherenkov telescopes, air-shower arrays, neutrino telescopes and gravitational wave detectors unique tools to probe the expected tiny effects that the breaking of Lorentz invariance would cause in the propagation of these messengers, in comparison with the standard scenario. In this chapter we explain the expected effects that the mentioned detectors can measure and summarize current results of searches for Lorentz violation.

As the James Webb Space Telescope (JWST) has become fully operational, early-release data are now available to begin building the tools and calibrations for precision point-source photometry and astrometry in crowded cluster environments. Here, we present our independent reduction of NIRCam imaging of the metal-poor globular cluster M92, which were collected under Director's Discretionary Early Release Science programme ERS-1334. We derived empirical models of the Point Spread Function (PSF) for filters F090W, F150W, F277W, and F444W, and find that these PSFs: (i) are generally under-sampled (FWHM~2 pixel) in F150W and F444W and severely under-sampled (FWHM~1 pixel) in F090W and F277W; (ii) have significant variation across the field of view, up to ~15-20 %; and (iii) have temporal variations of ~3-4 % across multi-epoch exposures. We deployed our PSFs to determine the photometric precision of NIRCam for stars in the crowded, central regions of M92, measured to be at the ~0.01 mag level. We use these data to construct the first JWST colour-magnitude diagrams of a globular cluster. Employing existing stellar models, we find that the data reach almost the bottom of the M92 main sequence (~0.1 M$_{\odot}$), and reveal 24 white dwarf candidate members of M92 in the brightest portion of the white dwarf cooling sequence. The latter are confirmed through a cross-match with archival HST UV and optical data. We also detect the presence of multiple stellar populations along the low-mass main sequence of M92.

A critical consideration in the design of next generation gravitational wave detectors is isolation from seismic vibrations that introduces various coherent and incoherent noises to the interferometers at different frequencies. We present the results of a detailed low frequency ambient seismic noise characterization (0.1--10~Hz) at Gingin High Optical Power Facility in Western Australia using a seismic array. The dominant noise sources below 1~Hz is microseism (0.06--1~Hz), strongly correlated with swell and sea heights measured by nearby buoy stations. Above 1~Hz, the seismic spectrum is dominated by wind induced seismic noise with a diurnal variation that prevents characterizing the background anthropogenic noise sources based on their daily power variations. We use f-k beamforming to distinguish between coherent and incoherent wind induced seismic noise. This allows the separation of some anthropogenic noise from wind induced noise based on the temporal variation of spatio-spectral properties. We show that the seismic coherency is reduced by wind induced seismic noise for wind speeds above 6~m/s. Furthermore, there are several spectral peaks between 4--9~Hz associated with the interaction of wind with a 40~m tall tower among which one at 4.2~Hz is strongest and coherent. By comparing our results with the properties of seismic noise at Virgo, we demonstrate that while the secondary microseism noise level is two orders of magnitude higher in Gingin (0.2~Hz), the anthropogenic noise level is three orders of magnitude lower between 2 and 4~Hz due to the absence of nearby road traffic. It is also at least one order of magnitude lower between 4 and 10~Hz due to the sparse population in Gingin.

It is well known that major solar eruptions are often produced by active regions with continual photospheric shearing and converging motions. Here, through high accuracy magnetohydrodynamics simulation, we show how solar eruption is initiated in a single bipolar configuration as driven by first shearing and then converging motions at the bottom surface. Different from many previous simulations, we applied the converging motion without magnetic diffusion, thus it only increases the magnetic gradient across the polarity inversion line but without magnetic flux cancellation. The converging motion at the footpoints of the sheared arcade creates a current sheet in a quasi-static way, and the eruption is triggered by magnetic reconnection of the current sheet, which supports the same scenario as shown in our previous simulation with only shearing motion. With the converging motion, the current sheet is formed at a lower height and has a higher current density than with shearing motion alone, which makes reconnection more effective and eruption stronger. Moreover, the converging motion renders a fast decay rate of the overlying field with height and thus favorable for an eruption. This demonstrate that the converging flow is more efficient to create the current sheet and more favorable for eruption than by solely the shearing flow.

As Setti & Woltjer noted back in 1973, quasars could be used to construct the Hubble diagram but the actual application was not that straightforward. It took years to implement the idea successfully. Most of the ways to use quasars for cosmology now require an advanced understanding of their structure, step by step. We briefly review this progress, with unavoidable personal bias, and concentrate on bright unobscured sources. We will mention the problem of the gas flow character close to the innermost stable circular orbit close to the black hole, discussed 50 years ago, which later led to the development of the slim disk scenario, but was recently revived in the context of Magnetically Arrested Disks (MAD) and Standard and Normal Evolution (SANE) disk models. We also discuss the hot/warm corona issue, which is still under discussion and complicates the analysis of the X-ray reflection. We present the scenario of the formation of the low ionization part of the Broad Line Region as a mostly failed wind powered by radiation pressure acting on dust (FRADO - Failed Radiatively Driven Dusty Outflow model). Next, we discuss the cosmological constraints that are currently achievable with quasars, mostly concentrating on light echo methods (continuum time delays and spectral line time delays with respect to the continuum) which are (or should be) incorporating the progress mentioned above. Finally, we briefly mention future prospects in this direction.

We examine astronomical observations that would be achievable over a future timeline corresponding to the documented history of human civilization so far, $\sim 10^4$ years. We examine implications for measurements of the redshift drift, evolution of the CMB, and cosmic parallax. A number of events that are rare on the scale of centuries will become easily observable on a timescale $\sim 10^4$ years. Implications for several measurements related to gravity are discussed.

We present a model of compact stars with a dark matter core. The hadronic equation of state is based on the parity doublet model and does not present a phase transition to quark matter. Instead, a strong first-order phase transition to dark matter described by a constant speed of sound model leads to the scenario of compact star mass twins. Compact star structural properties which obey state-of-the-art measurements and constraints are presented.

Blue large-amplitude pulsators (BLAPs) represent a new and rare class of hot pulsating stars with unusually large amplitudes and short periods. The evolutionary path that could give rise to such kinds of stellar configurations is unclear. Here we report on a comprehensive study of the peculiar BLAP discovered by the Tsinghua University - Ma Huateng Telescopes for Survey (TMTS), TMTS J035143.63+584504.2 (TMTS-BLAP-1). This new BLAP has an 18.9 min pulsation period and is similar to the BLAPs with a low surface gravity and an extended helium-enriched envelope, suggesting that it is a low-gravity BLAP at the shortest-period end. In particular, the long-term monitoring data reveal that this pulsating star has an unusually large rate of period change, P_dot/P=2.2e-6/yr. Such a significant and positive value challenges its origins from both helium-core pre-white-dwarfs and core helium-burning subdwarfs, but is consistent with that derived from shell helium-burning subdwarfs. The particular pulsation period and unusual rate of period change indicate that TMTS-BLAP-1 is at a short-lived (~10^6 yr) phase of shell-helium ignition before the stable shell-helium burning; in other words, TMTS-BLAP-1 is going through a "Hertzsprung gap" of hot subdwarfs.

The Double Asteroid Redirection Test (DART) spacecraft will impact into the asteroid Dimorphos on September 26, 2022 as a test of the kinetic impactor technique for planetary defense. The efficiency of the deflection following a kinetic impactor can be represented using the momentum enhancement factor, Beta, which is dependent on factors such as impact geometry and the specific target material properties. Currently, very little is known about Dimorphos and its material properties that introduces uncertainty in the results of the deflection efficiency observables, including crater formation, ejecta distribution, and Beta. The DART Impact Modeling Working Group (IWG) is responsible for using impact simulations to better understand the results of the DART impact. Pre-impact simulation studies also provide considerable insight into how different properties and impact scenarios affect momentum enhancement following a kinetic impact. This insight provides a basis for predicting the effects of the DART impact and the first understanding of how to interpret results following the encounter. Following the DART impact, the knowledge gained from these studies will inform the initial simulations that will recreate the impact conditions, including providing estimates for potential material properties of Dimorphos and Beta resulting from DARTs impact. This paper summarizes, at a high level, what has been learned from the IWG simulations and experiments in preparation for the DART impact. While unknown, estimates for reasonable potential material properties of Dimorphos provide predictions for Beta of 1-5, depending on end-member cases in the strength regime.

We investigated the quasi-periodic oscillation (QPO) features in the accretion-powered X-ray pulsar Cen X-3 observed by Insight-HXMT. For two observations in 2020 when Cen X-3 was in an extremely soft state, the power density spectrum revealed the presence of obvious QPO features at $\sim$40 mHz with an averaged fractional rms amplitude of $\sim9\%$. We study the mHz QPO frequency and rms amplitude over orbital phases, and find that the QPO frequency is $\sim$33-39 mHz at the orbital phase of 0.1-0.4, increasing to $\sim$37-43 mHz in the orbital phase of 0.4-0.8, but has no strong dependence on X-ray intensity. We also carried out an energy-dependent QPO analysis, the rms amplitude of the mHz QPOs have a decreasing trend as the energy increases from 2 to 20 keV. In addition, the QPO time-lag analysis shows that the time delay is $\sim 20$ ms (a hard lag) in the range of $\sim$5-10 keV, and becomes negative (time lag of $-(20-70)$ ms) above $\sim 10$ keV. The different QPO theoretical models are summarized and discussed. In the end, we suggest that these energy-dependent timing features as well as the origin of mHz QPOs in Cen X-3 may be ascribed to an instability when the accretion disk is truncated near the corotation radius.

The dynamics and the structure of the solar corona are determined by its magnetic field. Measuring coronal magnetic fields is, however, extremely hard. The polarization of the low-frequency radio emissions is one of the few observational probes of magnetic fields in the mid and high corona. Polarimetric calibration and imaging of the Sun at these frequencies is challenging. The brightness temperature and degree of polarization of the low-frequency solar radio emissions can vary by several orders of magnitude. These emissions also show dramatic spectral and temporal variations. Hence, to study these radio emissions, one needs high dynamic range spectro-polarimetric snapshot imaging. The Murchison Widefield Array (MWA), a Square Kilometre Array (SKA) precursor, is exceptionally well-suited for this purpose. Calibration and imaging of solar data to extract this information are, however, significant challenges in themselves - requiring a deep understanding of the instrument, capable sophisticated algorithms, and their reliable implementation. To meet these challenges we have developed an unsupervised and robust polarization calibration and imaging software pipeline. Here we present the architecture and some implementation details of this pipeline. It delivers high-fidelity and high-dynamic-range full polarimetric solar radio images at high spectro-temporal resolutions. We expect this pipeline to enable exciting new science with instruments like the MWA. We also hope that by not requiring a significant prior background in radio interferometric imaging, this pipeline will encourage wider use of radio imaging data in the larger solar physics community. The algorithm implemented here can easily be adapted for future arrays like the SKA.

Neutrino physics in the early Universe is key to our understanding of later cosmological stages, such as primordial nucleosynthesis (BBN) or the formation of large-scale structures. The coming decade promises new experimental results to explore and constrain cosmological models even more precisely - which requires robust theoretical predictions. This PhD thesis presents a study of the evolution of neutrinos in the first seconds after the Big Bang, more precisely when the temperature of the Universe is of the order of one mega-electronvolt. This evolution is obtained numerically by solving kinetic equations for which we propose a new derivation. A first application is the calculation of the so-called "standard" decoupling in order to calculate the cosmological parameter quantifying the energy density of the primordial relativistic species, $N_\mathrm{eff}$, to a precision of a few ten-thousandths. This study has highlighted the possibility of effectively describing the phenomenon of flavour oscillations, taking advantage of the large separation of time scales involved. Such an approximation is then adapted and validated in the case of non-zero asymmetries between neutrinos and antineutrinos. Finally, we study semi-analytically the consequences of incomplete neutrino decoupling on BBN, in order to understand how the primordial abundances of helium and deuterium are affected by this physics.

The largest and most close-in exoplanets would reflect enough star light to enable its ground-based photometric detection under the condition of a high to moderate albedo. We present the results of an observing campaign of secondary eclipse light curves of three of the most suitable exoplanet targets, WASP-43b, WASP-103b, and TrES-3b. The observations were conducted with meter-sized telescopes in the blue optical broadband filters Johnson B and Johnson V. We do not detect a photometric dimming at the moment of the eclipse, and derive a best-fit eclipse depth by an injection-recovery test. These depth values are then used to infer low geometric albedos ranging from zero to 0.18 with an uncertainty of 0.12 or better in most cases. This work illustrates the potential of ground-based telescopes to provide wavelength-resolved reflection properties of selected exoplanets even at short optical wavelengths, which otherwise are only accessible by the Hubble Space Telescope.

Vorticity is central to the nature of, and dynamical processes in turbulence, including turbulence in astrophysical fluids. The results of \cite{Raymond20a,Raymond20b} on vorticity in the post-shock fluid of the Cygnus Loop supernova remnant are therefore of great interest. We consider the degree to which spectroscopic measurements of an optically-thin line, the most common type of astronomical velocimetry, can yield unambiguous measurements of the vorticity in a fluid. We consider an ideal case of observations in the plane of a flow which may or may not contain vorticity. In one case, the flow possesses vorticity in a direction perpendicular to the plane of observations. In the other case, the flow is irrotational (zero vorticity) by construction. The observationally-deduced vorticity (referred to as the {\em pseudovorticity}) is inferred from spatial differences in the line-of-sight component of velocity, and assumptions of symmetry. My principal result is that in the case of the vortical flow, the pseudovorticity is a reasonable match for the true vorticity. However, and importantly, the pseudovorticity in the case of the irrotational flow field is also nonzero, and comparable in magnitude to that for a vortical flow. The conclusion of this paper is that while astronomical spectroscopic observations may yield a good estimate of the vorticity in a remote fluid, the robustness of such an inference cannot be insured.

Many prominences are supported by magnetic flux ropes. One important question is how we can determine whether the flux rope is weakly-twisted or strongly-twisted. In this paper, we attempted to check whether prominences supported by weakly-twisted and strongly-twisted flux ropes can manifest different features so that we might distinguish the two types of magnetic structures by their appearance. We performed pseudo three-dimensional simulations of two magnetic flux ropes with different twists. We found that the resulting two prominences differ in many aspects. The prominence supported by a weakly-twisted flux rope is composed mainly of transient threads, forming high-speed flows inside the prominence. Its horns are evident. Conversely, the one supported by a highly-twisted flux rope consists mainly of stable quasi-stationary threads, including longer independently trapped threads and shorter magnetically connected threads. It is also revealed that the prominence spine deviates from the flux rope axis in the vertical direction and from the photospheric polarity inversion line projected on the solar surface, especially for the weakly-twisted magnetic flux rope. The two types of prominences differ significantly in appearance. It is also suggested that a piling-up of short threads in highly-twisted flux ropes might account for the vertical-like threads in some prominences.

The acceleration of LMXB 4U 1820-30 that derived from its orbital-period derivative $\dot P_{\rm b}$ was supposed to be the evidence for an Intermediate-mass Black Hole (IMBH) in the Galactic globular cluster (GC) NGC 6624. However, we find that the anomalous $\dot P_{\rm b}$ is mainly due to the gravitational wave emission, rather than the acceleration in cluster potential. Using the standard structure models of GCs, we simulate acceleration distributions for pulsars in the central region of the cluster. By fitting the acceleration of J1823-3021A with the simulated distribution profiles (maximum values), it is suggested that an IMBH with mass $M\gtrsim 950^{+550}_{-350}~M_{\odot}$ may reside in the cluster center. We further show that the second period derivative $\ddot P$ of J1823-3021A is probably due to the gravitational perturbation of a nearby star.

Type Ia supernovae (SNae Ia), standardisable candles that allow tracing the expansion history of the Universe, are instrumental in constraining cosmological parameters, particularly dark energy. State-of-the-art likelihood-based analyses scale poorly to future large datasets, are limited to simplified probabilistic descriptions, and must explicitly sample a high-dimensional latent posterior to infer the few parameters of interest, which makes them inefficient. Marginal likelihood-free inference, on the other hand, is based on forward simulations of data, and thus can fully account for complicated redshift uncertainties, contamination from non-SN Ia sources, selection effects, and a realistic instrumental model. All latent parameters, including instrumental and survey-related ones, per-object and population-level properties, are implicitly marginalised, while the cosmological parameters of interest are inferred directly. As a proof of concept, we apply truncated marginal neural ratio estimation (TMNRE), a form of marginal likelihood-free inference, to BAHAMAS, a Bayesian hierarchical model for SALT parameters. We verify that TMNRE produces unbiased and precise posteriors for cosmological parameters from up to 100 000 SNae Ia. With minimal additional effort, we train a network to infer simultaneously the O(100 000) latent parameters of the supernovae (e.g. absolute brightnesses). In addition, we describe and apply a procedure that utilises local amortisation of the inference to convert the approximate Bayesian posteriors into frequentist confidence regions with exact coverage. Finally, we discuss the planned improvements to the model that are enabled by using a likelihood-free inference framework, like selection effects and non-Ia contamination.

We present James Webb Space Telescope (JWST) Mid-InfraRed Instrument (MIRI) integral-field spectroscopy of the nearby merging, luminous infrared galaxy, NGC 7469. This galaxy hosts a Seyfert type-1.5 nucleus, a highly ionized outflow, and a bright, circumnuclear star-forming ring, making it an ideal target to study AGN feedback in the local Universe. We take advantage of the high spatial/spectral resolution of JWST/MIRI to isolate the star-forming regions surrounding the central active nucleus and study the properties of the dust and warm molecular gas on ~100 pc scales. The starburst ring exhibits prominent Polycyclic Aromatic Hydrocarbon (PAH) emission, with grain sizes and ionization states varying by only ~30%, and a total star formation rate of $\rm 10 - 30 \ M_\odot$/yr derived from fine structure and recombination emission lines. Using pure rotational lines of H2, we detect 1.2$\times$10$^{7} \rm \ M_\odot$ of warm molecular gas at a temperature higher than 200 K in the ring. All PAH bands get significantly weaker towards the central source, where larger and possibly more ionized grains dominate the emission. However, the bulk of the dust and molecular gas in the ring appears unaffected by the ionizing radiation or the outflowing wind from the AGN. These observations highlight the power of JWST to probe the inner regions of dusty, rapidly evolving galaxies for signatures of feedback and inform models that seek to explain the co-evolution of supermassive black holes and their hosts.

We present a segment-level wavefront stability error budget for space telescopes essential for exoplanet detection. We use a detailed finite element model to relate the temperature gradient at the location of the primary mirror to wavefront variations on each of the segment. We apply the PASTIS sensitivity model forward approach to allocate static tolerances in physical units for each segment, and transfer these tolerances to the temporal domain via a model of the WFS&C architecture in combination with a Zernike phase sensor and science camera. We finally estimate the close-loop variance and limiting contrast for the segments' thermo-mechanical modes.

We present Fabry-Perot observations in the H$\alpha$ and [S II] lines to study the kinematics of the Magellanic-type dwarf irregular galaxy NGC 1569, these observations allowed us to computed the H$\alpha$ velocity field of this galaxy. Doing a detailed analysis of the velocity along the line-of-sight and H$\alpha$ velocity profiles, we identified the origin of most of the motions in the innermost parts of the galaxy and discarded the possibility of deriving a rotation curve that traces the gravitational well of the galaxy. We analysed the kinematics of the ionised gas around 31 supernova remnants previously detected in NGC 1569 by other authors, in optical and radio emission. We found that the H$\alpha$ velocity profiles of the supernova remnants are complex indicating the presence of shocks. Fitting these profiles with several Gaussian functions, we computed their expansion velocities which rank from 87 to 188 km s$^{-1}$ confirming they are supernova remnants. Also, we determined the physical properties such as electron density, mechanical energy, and kinematic age for 30 of the 31 supernova remnants and found they are in the radiative phase with an energy range from 1 to 39$\times$10$^{50}$ erg s$^{-1}$ and an age from 2.3 to 8.9$\times$10$^4$ yr. Finally, we estimated the Surface Brightness - Diameter ($\Sigma$-D) Relation for NGC 1569 and obtained a slope $\beta$ = 1.26$\pm$0.2, comparable with the $\beta$ value obtained for supernova remnants in galaxies M31 and M33.

In this work we present the first results from a new ray-tracing tool to calculate cosmological distances in the context of fully nonlinear general relativity. We use this tool to study the ability of the general cosmographic representation of luminosity distance, as truncated at third order in redshift, to accurately capture anisotropies in the "true" luminosity distance. We use numerical relativity simulations of cosmological large-scale structure formation which are free from common simplifying assumptions in cosmology. We find the general, third-order cosmography is accurate to within 1% for redshifts to z\approx 0.034 when sampling scales strictly above 100 Mpc/h, which is in agreement with an earlier prediction. We find the inclusion of small-scale structure generally spoils the ability of the third-order cosmography to accurately reproduce the full luminosity distance for wide redshift intervals, as might be expected. For a simulation sampling small-scale structures, we find a +/- 5% variance in the monopole of the ray-traced luminosity distance at z \approx 0.02. Further, all 25 observers we study here see a 9--20% variance in the luminosity distance across their sky at z \approx 0.03, which reduces to 2--5% by z \approx 0.1. These calculations are based on simulations and ray tracing which adopt fully nonlinear general relativity, and highlight the potential importance of fair sky-sampling in low-redshift isotropic cosmological analysis.

VLA 1623 West is an ambiguous source that has been described as a shocked cloudlet as well as a protostellar disk. We use deep ALMA 1.3 and 0.87 millimeter observations to constrain its shape and structure to determine its origins better. We use a series of geometric models to fit the uv visibilities at both wavelengths with GALARIO. Although the Real visibilities show structures similar to what has been identified as gaps and rings in protoplanetary disks, we find that a modified Flat-Topped Gaussian at high inclination provides the best fit to the observations. This fit agrees well with expectations for an optically thick, highly inclined disk. Nevertheless, we find that the geometric models consistently yield positive residuals at the four corners of the disk at both wavelengths. We interpret these residuals as evidence that the disk is flared in the millimeter dust. We use a simple toy model for an edge-on flared disk and find that the residuals best match a disk with flaring that is mainly restricted to the outer disk at $R \gtrsim 30$ au. Thus, VLA 1623W may represent a young protostellar disk where the large dust grains have not yet had enough time to settle into the mid-plane. This result may have implications for how disk evolution and vertical dust settling impact the initial conditions leading to planet formation.

The streaming instability (SI) is one of the most promising candidates for triggering planetesimal formation by producing dense dust clumps that undergo gravitational collapse. Understanding how the SI operates in realistic protoplanetary disks (PPDs) is therefore crucial to assess the efficiency of planetesimal formation. Modern models of PPDs show that large-scale magnetic torques or winds can drive laminar gas accretion near the disk midplane. In a previous study, we identified a new linear dust-gas instability, the azimuthal drift SI (AdSI), applicable to such accreting disks and is powered by the relative azimuthal motion between dust and gas that results from the gas being torqued. In this work, we present the first nonlinear simulations of the AdSI. We show that it can destabilize an accreting, dusty disk even in the absence of a global radial pressure gradient, which is unlike the classic SI. We find the AdSI drives turbulence and the formation of vertically-extended dust filaments that undergo merging. In dust-rich disks, merged AdSI filaments reach maximum dust-to-gas ratios exceeding 100. Moreover, we find that even in dust-poor disks the AdSI can increase local dust densities by two orders of magnitude. We discuss the possible role of the AdSI in planetesimal formation, especially in regions of an accreting PPD with vanishing radial pressure gradients.

S-type planets, which orbit one component of multiple-star systems, place strong constraints on the planet formation and evolution models. A notable case study is Kepler-444, a triple-star system whose primary is orbited by five planets smaller than Venus in a compact configuration, and for which the stellar binary companion revolves around the primary on a highly eccentric orbit. Having access to the most precise up-to-date masses and orbital parameters is highly valuable to understand formation and evolution processes. We provide the first full dynamical exploration of this system, with the goal to refine those parameters. The planetary system does not appear in any of low-order two or three-planet mean-motion resonances (MMR). We provide the most precise up-to-date dynamical parameters for the planets and the stellar binary companion, using an approach that makes use of the Numerical Analysis of Fundamental Frequencies (NAFF) fast chaos indicator. The orbit of the latter is constrained by new observations from HIRES and Gaia, and also by the stability analysis. This update further challenges the planets formation processes. We also test the dynamical plausibility of a sixth planet in the system, following hints observed in the Hubble Space Telescope (HST) data. We find that this putative planet could exist over a broad range of masses, and with an orbital period roughly comprised between 12 and 20 days. We note an overall good agreement of the system with short-term orbital stability. This suggests that a diverse range of planetary system architectures could be found in multiple-star systems, potentially further challenging the planet formation models.

Searches for dark matter decaying into photons constrain its lifetime to be many orders of magnitude larger than the age of the Universe. A corollary statement is that the abundance of any particle that can decay into photons over cosmological timescales is constrained to be much smaller than the cold dark-matter density. We show that an $\textit{irreducible}$ freeze-in contribution to the relic density of axions is in violation of that statement in a large portion of the parameter space. This allows us to set stringent constraints on axions in the mass range $100\rm \;eV - 100\; MeV$. At $10\rm \; keV$ our constraint on a photophilic axion is $g_{a\gamma \gamma} \lesssim 8.1 \times 10^{-14}~{\rm GeV}^{-1}$, almost three orders of magnitude stronger than the bounds established using horizontal branch stars; at $100~{\rm keV}$ our constraint on a photophobic axion coupled to electrons is $g_{aee} \lesssim 8.0 \times 10^{-15}$, almost four orders of magnitude stronger than present results. Although we focus on axions, our argument is more general and can be extended to, for instance, sterile neutrinos.

We investigate the gravitational lensing by spinning Proca stars and the shadows and lensing by Kerr black holes (BHs) with synchronised Proca hair, discussing both theoretical aspects and observational constraints from the Event Horizon Telescope (EHT) M87* and Sgr A* data. On the theoretical side, this family of BHs interpolates between Kerr-like solutions -- exhibiting a similar optical appearance to that of Kerr BHs -- to very non-Kerr like solutions, exhibiting exotic features such as cuspy shadows, egg-like shadows and ghost shadows. We interpret these features in terms of the structure of the fundamental photon orbits, for which different branches exist, containing both stable and unstable orbits, with some of the latter not being shadow related. On the observational side, we show that current EHT constraints are compatible with all such BHs that could form from the growth of the superradiant instability of Kerr BHs. Unexpectedly, given the (roughly) 10% error bars in the EHT data -- and in contrast to their scalar cousin model --, some of the BHs with up to 40% of their energy in their Proca hair are compatible with the current data. We estimate the necessary resolution of future observations to better constrain this model.

We apply the recently developed analysis of 16 years of INTEGRAL/SPI data including a dark matter spatial template to derive bounds on dark matter candidates lighter than WIMPs (like sterile neutrinos or axion-like particles) decaying into line or continuum electromagnetic final state channels. The bounds obtained are the strongest to date for dark matter masses between $\sim $60 keV and $\sim$16 MeV experiencing two-body decays producing photon lines.

The flux of high energy neutrinos and photons produced in a blazar could get attenuated when they propagate through the dark matter spike around the central black hole and the halo of the host galaxy. Using the observation by IceCube of a few high-energy neutrino events from TXS 0506+056, and their coincident gamma ray events, we obtain new constraints on the dark matter-neutrino and dark matter-photon scattering cross sections. Our constraints are orders of magnitude more stringent than those derived from considering the attenuation through the intergalactic medium and the Milky Way dark matter halo. When the cross-section increases with energy, our constraints are also stronger than those derived from the CMB and large-scale structure.

Developing an effective automatic classifier to separate genuine sources from artifacts is essential for transient follow-ups in wide-field optical surveys. The identification of transient detections from the subtraction artifacts after the image differencing process is a key step in such classifiers, known as real-bogus classification problem. We apply a self-supervised machine learning model, the deep-embedded self-organizing map (DESOM) to this "real-bogus" classification problem. DESOM combines an autoencoder and a self-organizing map to perform clustering in order to distinguish between real and bogus detections, based on their dimensionality-reduced representations. We use 32x32 normalized detection thumbnails as the input of DESOM. We demonstrate different model training approaches, and find that our best DESOM classifier shows a missed detection rate of 6.6% with a false positive rate of 1.5%. DESOM offers a more nuanced way to fine-tune the decision boundary identifying likely real detections when used in combination with other types of classifiers, for example built on neural networks or decision trees. We also discuss other potential usages of DESOM and its limitations.

The propagation of gravitational waves can be described in terms of null geodesics by using the geometrical optics approximation. However, at large but finite frequencies the propagation is affected by the spin-orbit coupling corrections to geometrical optics, known as the gravitational spin Hall effect. Consequently, gravitational waves follow slightly different frequency- and polarization-dependent trajectories. We study the potential for detection of the gravitational spin Hall effect in hierarchical triple black hole systems, consisting of an emitting binary orbiting a more massive body, acting as a gravitational lens. We calculate the difference in time of arrival with respect to the geodesic propagation and find that it follows a simple power law dependence on frequency with a fixed exponent. We calculate the gravitational spin Hall-corrected waveform and its mismatch with respect to the original waveform. The waveform carries a measurable imprint of the strong gravitational field if the source, lens and observer are sufficiently aligned or for generic observers if the source is close enough to the lens. We demonstrate that the gravitational spin Hall effect can be detected, providing an interesting avenue to probe general relativity and the environments of compact binary systems.

We have performed the first direct measurement of two resonances of the $^7$Be($\alpha,\gamma$)$^{11}$C reaction with unknown strengths using an intense radioactive $^7$Be beam and the DRAGON recoil separator. We report on the first measurement of the 1155 and 1110 keV resonance strengths of $1.73 \pm 0.25(stat.) \pm 0.40(syst.)$ eV and $125 ^{+27}_{-25}(stat.) \pm 15(syst.)$ meV, respectively. The present results have reduced the uncertainty in the $^7$Be($\alpha,\gamma$)$^{11}$C reaction rate to $\sim$ 9.4-10.7% over T = 1.5-3 GK, which is relevant for nucleosynthesis in the neutrino-driven outflows of core-collapse supernovae ($\nu p$-process). We find no effect of the new, constrained reaction rate on $\nu p$-process nucleosynthesis.

A possible mechanism to explain the origin of the light $p$-nuclei in the Galaxy is the nucleosynthesis in the proton-rich neutrino-driven wind ejecta of core-collapse supernovae via the $\nu p$-process. However this production scenario is very sensitive to the underlying supernova dynamics and the nuclear physics input. As far as the nuclear uncertainties are concerned, the breakout from the $pp$-chains via the $^7$Be($\alpha,\gamma$)$^{11}$C reaction has been identified as an important link which can influence the nuclear flow and therefore the efficiency of the $\nu p$-process. However its reaction rate is poorly known over the relevant temperature range, T = 1.5-3 GK. We report on the first direct measurement of two resonances of the $^7$Be($\alpha,\gamma$)$^{11}$C reaction with previously unknown strengths using an intense radioactive $^7$Be beam from the ISAC facility and the DRAGON recoil separator in inverse kinematics. We have decreased the $^7$Be($\alpha,\gamma$)$^{11}$C reaction rate uncertainty to $\sim$ 9.4-10.7% over the relevant temperature region.

We discuss the cosmological inflation in the context of $f(Q)$ non-metric gravity, where $Q$ is the non-metric scalar. After introducing conformal transformations for $f(Q)$ gravity, we first focus on the potential-slow-roll inflation by studying the corresponding potentials for different forms of the function $f(Q)$ in the Einstein frame. Secondly, we investigate the Hubble-slow-roll inflation for three classes of inflationary potentials considered for the specific form $f(Q)\propto Q^{2}$, in the Jordan frame. We compare results in both approaches with observations coming from Planck and BICEP2/Keck array satellites. Observational constraints on the parameters space of the models are obtained as well.

Experimental investigations of the fine plasma structure of interplanetary shocks are extremely difficult to conduct due to their small thickness and high speed relative to the spacecraft. We study the variations in the parameters of twice-ionized helium ions (4He++ ions or {\alpha}-particles) in the solar wind plasma during the passage of interplanetary shocks and Earth's bow shock. We use data with high time resolution gathered by the BMSW (Bright Monitor of Solar Wind) instrument installed on the SPEKTR-R satellite, which operated between August 2011 and 2019. The MHD parameters of He++ ions (the bulk velocity Va, temperature Ta, absolute density Na, and helium abundance Na/Np) are analyzed for 20 interplanetary shocks and compared with similar parameters for 25 Earth's bow shock crossings. Measurements from the WIND, Cluster and THEMIS satellites were also analyzed. The correlations in the changes in helium abundance Na/Np with the parameters {\beta}_{i}, {\theta}_{Bn} and M_{ms} were investigated. A correlation between Na/Np and the angle {\theta}_{Bn} was found: the lower the value of {\theta}_{Bn}, the greater the drop in helium abundance Na/Np falls behind the IP shock front. For Earth's bow shock crossings, we found a significant increase in the helium abundance Na/Np in quasi-perpendicular events.