Papers for Tuesday, Jul 19 2022

We present a spectroscopic analysis of 44 low-luminosity host galaxies of Type Ia supernovae (SNe Ia) detected by the All-Sky Automated Survey for Supernovae (ASAS-SN), using the emission lines to measure metallicities and star formation rates. We find that although the star formation activity of our sample is representative of general galaxies, there is some evidence that the lowest-mass SN Ia host galaxies (log($M_\star/M_\odot$)$<8$) in our sample have high metallicities compared to general galaxies of similar masses. We also identify a subset of 5 galaxies with particularly high metallicities. This highlights the need for spectroscopic analysis of more low-luminosity, low-mass SN Ia host galaxies to test the robustness of these conclusions and their potential impact on our understanding of SN Ia progenitors.

Neutron star mergers have recently become a tool to study extreme gravity, nucleosynthesis, and the chemical composition of the Universe. To date, there has been one joint gravitational and electromagnetic observation of a binary neutron star merger, GW170817, as well as a solely gravitational observation, GW190425. In order to accurately identify and interpret electromagnetic signals of neutron star mergers, better models of the matter outflows generated by these mergers are required. We compare a series of ejecta models to see where they provide strong constraints on the amount of ejected mass expected from a system, and where systematic uncertainties in current models prevent us from reliably extracting information from observed events. We also examine 2396 neutron star equations of state compatible with GW170817 to see whether a given ejecta mass could be reasonably produced with a neutron star of said equation of state, and whether different ejecta models provide consistent predictions. We find that the difference between models is often comparable to or larger than the error generally assumed for these models, implying better constraints on the models are needed. We also note that the extrapolation of outflow models outside of their calibration window, while commonly needed to analyze gravitational wave events, is extremely unreliable and occasionally leads to completely unphysical results.

Modern grating manufacturing techniques suffer from inherent issues that limit their peak efficiencies. The anisotropic etching of silicon facilitates the creation of custom gratings that have sharp and atomically smooth facets, directly addressing these issues. We describe work to fabricate and characterize etched silicon echelles optimized for the far ultraviolet (FUV; 90 - 180 nm) bandpass. We fabricate two echelles that have similar parameters to the mechanically ruled grating flown on the CHESS sounding rocket. We demonstrate a 42% increase in peak order efficiency and an 83% decrease in interorder scatter using these gratings. We also present analysis on where the remaining efficiency resides. These demonstrated FUV echelle improvements benefit the faint source sensitivity and high-resolution performance of future UV observatories.

Reionization of helium is expected to occur at redshifts $z\sim3$ and have important consequences for quasar populations, galaxy formation, and the morphology of the intergalactic medium, but there is little known empirically about the process. Here we show that kinetic Sunyaev-Zeldovich (kSZ) tomography, based on the combination of CMB measurements and galaxy surveys, can be used to infer the primordial helium abundance as well as the time and duration of helium reionization. We find a high-significance detection at ${\sim10\sigma}$ can be expected from Vera Rubin Observatory and CMB-S4 in the near future. A more robust characterization of helium reionization will require next-generation experiments like MegaMapper (a proposed successor to DESI) and CMB-HD.

As one of the promising candidates of cold dark matter (DM), primordial black holes (PBHs) were formed due to the collapse of over-densed regions generated by the enhanced curvature perturbations during the radiation-dominated era. The enhanced curvature perturbations are expected to be non-Gaussian in some relevant inflation models and hence the higher-order loop corrections to the curvature power spectrum might be non-negligible as well as altering the abundance of PBHs. In this paper, we calculate the one-loop correction to the curvature power spectrum with local-type non-Gaussianities characterizing by $F_{\mathrm{NL}}$ and $G_{\mathrm{NL}}$ standing for the quadratic and cubic non-Gaussian parameters, respectively. Requiring that the one-loop correction be subdominant, we find a perturbativity condition, namely $|2cAF_{\mathrm{NL}}^2+6AG_{\mathrm{NL}}|\ll 1$, where $c$ is a constant coefficient which can be explicitly calculated in the given model and $A$ denotes the variance of Gaussian part of enhanced curvature perturbation, and such a perturbativity condition can provide a stringent constraint on the relevant inflation models for the formation of PBHs.

When studying transiting exoplanets it is common to assume a spherical planet shape. However short rotational periods can cause a planet to bulge at its equator, as is the case with Saturn whose equatorial radius is almost 10% larger than its polar radius. As a new generation of instruments comes online, it is important to continually assess the underlying assumptions of models to ensure robust and accurate inferences. We analyze bulk samples of known transiting planets and calculate their expected signal strength if they were to be oblate. We find that for noise levels below 100ppm, as many as 100 planets could have detectable oblateness. We also investigate the effects of fitting spherical planet models to synthetic oblate lightcurves. We find that this biases the retrieved parameters by several standard deviations for oblateness values > 0.1-0.2. When attempting to fit an oblateness model to both spherical and oblate lightcurves, we find that the sensitivity of such fits is correlated with both the SNR as well as the time sampling of the data, which can mask the oblateness signal. For typical values of these quantities for Kepler observations, it is difficult to rule out oblateness values less than ~0.25. This results in an accuracy wall of 10-15% for the density of planets which may be oblate. Finally, we find that a precessing oblate planet has the ability to mimic the signature of a long-period companion via transit timing variations, inducing offsets at the level of 10s of seconds.

We present a Bayesian hierarchical framework to analyze photometric galaxy survey data with stellar population synthesis (SPS) models. Our method couples robust modeling of spectral energy distributions with a population model and a noise model to characterize the statistical properties of the galaxy populations and real observations, respectively. By self-consistently inferring all model parameters, from high-level hyper-parameters to SPS parameters of individual galaxies, one can separate sources of bias and uncertainty in the data. We demonstrate the strengths and flexibility of this approach by deriving accurate photometric redshifts for a sample of spectroscopically-confirmed galaxies in the COSMOS field, achieving a performance competitive with publicly-released photometric redshift catalogs based on the same data. Prior to this work, this approach was computationally intractable in practice due to the heavy computational load of SPS model calls; we overcome this challenge using with neural emulators. We find that the largest photometric residuals are associated with poor calibration for emission line luminosities and thus build a framework to mitigate these effects. This combination of physics-based modeling accelerated with machine learning paves the path towards meeting the stringent requirements on the accuracy of photometric redshift estimation imposed by upcoming cosmological surveys. The approach also has the potential to create new links between cosmology and galaxy evolution through the analysis of photometric datasets.

We aim to hunt for massive binaries hosting a black hole companion (OB+BH) and establish the natal mass-ratio distribution of massive stars at the subsolar metallicity environment of the Large Magellanic Cloud (LMC). To this end, we use the shift-and-add grid disentangling technique to characterize the hidden companions in 51 SB1 O-type and evolved B-type binaries in the LMC monitored in the framework of the Tarantula Massive Binary Monitoring (TMBM). We find that, out of the 51 SB1 systems, 43 (84%) are found to have non-degenerate stellar companions, of which 28 are confident detections, and 15 are less certain (SB1: or SB2:). Of these 43 targets, one is found to be a triple (VFTS 64), and two are found to be quadruples (VFTS 120, 702). The remaining eight targets (16%) retain an SB1 classification. Aside from the unambiguous case of VFTS 243, analysed in detailed in a separate paper, we identify two additional OB+BH candidates: VFTS 514 and VFTS 779. Additional black holes may be present in the sample but at a lower probability. Our study firmly establishes a virtually flat natal mass-ratio distribution for O-type stars at LMC metallicity, covering the entire mass-ratio range (0.05 < q < 1) and periods in the range 0 < log P < 3 [d]. The nature of the OB+BH candidates should be verified through future monitoring, but the frequency of OB+BH candidates is generally in line with recent predictions at LMC metallicity.

Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray faint binary in the Large Magellanic Cloud. With an orbital period of 10.4-d, it comprises an O-type star of 25 solar masses and an unseen companion of at least nine solar masses. Our spectral analysis excludes a non-degenerate companion at a 5-sigma confidence level. The minimum companion mass implies that it is a black hole. No other X-ray quiet black hole is unambiguously known outside our Galaxy. The (near-)circular orbit and kinematics of VFTS 243 imply that the collapse of the progenitor into a black hole was associated with little or no ejected material or black-hole kick. Identifying such unique binaries substantially impacts the predicted rates of gravitational-wave detections and properties of core-collapse supernovae across the Cosmos.

A non-negligible source of systematic bias in cosmological analyses of galaxy surveys is the on-sky modulation caused by foregrounds and variable image characteristics such as observing conditions. Standard mitigation techniques perform a regression between the observed galaxy density field and sky maps of the potential contaminants. Such maps are ad-hoc, lossy summaries of the heterogeneous sets of co-added exposures that contribute to the survey. We present a methodology to address this limitation, and extract the spurious correlations between the observed distribution of galaxies and arbitrary stacks of single-epoch exposures. We study four types of galaxies (LRGs, ELGs, QSOs, LBGs) in the three regions of the DESI Legacy Surveys (North, South, DES), which results in twelve samples with varying levels and type of contamination. We find that the new technique outperforms the traditional ones in all cases, and is able to remove higher levels of contamination. This paves the way for new methods that extract more information from multi-epoch galaxy survey data and mitigate large-scale biases more effectively.

Resolved observations of debris discs can be used to derive radial profiles of Azimuthally-averaged Surface Density (ASD), which carries important information about the disc structure even in presence of non-axisymmetric features and has improved signal-to-noise characteristics. We develop a (semi-)analytical formalism allowing one to relate ASD to the underlying semi-major axis and eccentricity distributions of the debris particles in a straightforward manner. This approach does not involve the distribution of particle apsidal angles, thus simplifying calculations. It is a much faster, more flexible and effective way of calculating ASD than the Monte Carlo sampling of orbital parameters of debris particles. We present explicit analytical results based on this technique for a number of particle eccentricity distributions, including two cases of particular practical importance: a prescribed radial profile of eccentricity, and the Rayleigh distribution of eccentricities. We then show how our framework can be applied to observations of debris discs and rings for retrieving either the semi-major axis distribution or (in some cases) the eccentricity distribution of debris, thus providing direct information about the architecture and dynamical processes operating in debris discs. Our approach also provides a fast and efficient way of forward modeling observations. Applications of this technique to other astrophysical systems, e.g. the nuclear stellar disc in M31 or tenuous planetary rings, are also discussed.

If Type Ia supernovae (SNe~Ia) result from a white dwarf being ignited by Roche lobe overflow from a nondegenerate companion, then as the supernova explosion runs into the companion star its ejecta will be shocked, causing an early blue excess in the lightcurve. A handful of these excesses have been found in single-object studies, but inferences about the population of SNe~Ia as a whole have been limited because of the rarity of multiwavelength followup within days of explosion. Here we present a three-year investigation yielding an unbiased sample of nine nearby ($z<0.01$) SNe~Ia with exemplary early data. The data are truly multiwavelength, covering $UBVgri$ and Swift bandpasses, and also early, with an average first epoch 16.0 days before maximum light. Of the nine objects, three show early blue excesses. We do not find enough statistical evidence to reject the null hypothesis that SNe~Ia predominantly arise from Roche-lobe-overflowing single-degenerate systems ($p=0.94$). When looking at the objects' colors, we find the objects are almost uniformly near-UV-blue, in contrast to earlier literature samples which found that only a third of SNe~Ia are near-UV-blue, and we find a seemingly continuous range of $B-V$ colors in the days after explosion, again in contrast with earlier claims in the literature. This study highlights the importance of early, truly multiwavelength, high-cadence data in determining the progenitor systems of SNe~Ia and in revealing their diverse early behavior.

The bulk of the present-day stellar mass was formed in galaxies when the universe was less than half its current age (i.e., $1 \lesssim z \lesssim 3$). While this likely marks one of the most critical time periods for galaxy evolution, we currently do not have a clear picture on the radial extent and distribution of cold molecular gas and associated star formation within the disks of galaxies during this epoch. Such observations are essential to properly estimate the efficiency at which such galaxies convert their gas into stars, as well as to account for the various energetic processes that govern this efficiency. Long-wavelength (i.e., far-infrared--to--radio) observations are critical to penetrate the high-levels of extinction associated with dusty, infrared-bright galaxies that are driving the stellar mass assembly at such epochs. In this article we discuss how the next-generation Very Large Array will take a transformative step in our understanding of galaxy formation and evolution by delivering the ability to simultaneously study the relative distributions molecular gas and star formation on sub-kpc scales unbiased by dust for large populations of typical galaxies in the early universe detected by future far-infrared space missions.

We present a forecast study on the cross-correlation between cosmic shear tomography from the Chinese Survey Space Telescope (CSST), and CMB lensing from Ali CMB Polarization Telescope (AliCPT-1) in Tibet. We generate the correlated galaxy lensing and CMB lensing signals from the Gaussian realizations based on the inputted auto- and cross-spectra. As for the error budget, we consider the CMB lensing reconstruction noise based on the AliCPT-1 lensing reconstruction pipeline; the shape noise of the galaxy lensing measurement; CSST photo-$z$ error; photo-$z$ bias; intrinsic alignment effect. The AliCPT-1 CMB lensing mock data are generated according to two experimental stages, namely the "4 modules*yr" and "48 modules*yr" cases. We estimate the cross-spectra in 4 tomographic bins according to the CSST photo-$z$ distribution in the range of $z\in[0,4)$. After reconstructing the pseudo-cross-spectra from the realizations, we calculate the signal-to-noise ratio (SNR). By combining the 4 photo-z bins, the total cross-correlation SNR$\approx17$ (AliCPT-1 "4 modules*yr") and SNR$\approx26$ (AliCPT-1 "48 modules*yr"). Finally, we study the cosmological application of this cross-correlation signal. Due to the negative contribution to the galaxy lensing data, the exclusion of intrinsic alignment in the template fitting will lead to roughly a $0.6\sigma$ increasement in $\sigma_8$ but without changing the $S_8$ value. For AliCPT-1 first and second stages, the cross-correlation of CSST cosmic shear with CMB lensing give $\sigma_8=0.770\pm 0.034$ and $S_8=0.797\pm 0.028$ and $\sigma_8=0.801\pm 0.023$ and $S_8=0.813\pm 0.016$, respectively.

Designing compact instruments is the key for the scientific exploration by smaller spacecrafts such as cubesats or by deep space missions. Such missions require compact instrument designs to have minimal instrument mass. Here we present a proof of concept for miniaturization of the Global Oscillation Network Group GONG instrument. GONG instrument routinely obtains solar full disk Doppler and magnetic field maps of the solar photosphere using Ni 676 nm absorption line. A key concept for miniaturization of GONG optical design is to replace the bulky Lyot filter with a narrow-band interference filter and reduce the length of feed telescope. We present validation of the concept via numerical modeling as well as by proof of concept observations.

Aims. To search for these rare objects, we study 32 Galactic O-type stars that were reported as SB1s in the literature. In our sample we include Cyg X-1, which is known to host an accreting stellar-mass BH, and HD 74194, a supergiant fast X-ray transient, in order to validate our methodology. The final goal is to characterise the nature of the unseen companions to determine if they are MS stars, stripped helium stars, triples, or compact objects such as neutron stars or stellar-mass BHs. Methods. After measuring radial velocities and deriving orbital solutions for all the systems in our sample, we performed spectral disentangling to extract putative signatures of faint secondary companions from the composite spectra. We derived stellar parameters for the visible stars and estimated the mass ranges of the secondary stars using the binary mass function. Variability observed in the photometric TESS light curves was also searched for indications of the presence of putative companions, degenerate or not. Results. In 17 of the 32 systems reported as SB1s, we extract secondary signatures, down to mass ratios of ~0.15. For the 17 newly detected double-lined spectroscopic binaries (SB2s), we derive physical properties of the individual components and discuss why they have not been detected as such before. Among the remaining systems, we identify nine systems with possible NS or low-mass MS companions. For Cyg X-1 and HD 130298, we are not able to extract any signatures for the companions, and the minimum masses of their companions are estimated to be about 7Msun. Our simulations show that secondaries with such a mass should be detectable from our dataset, no matter their nature: MS stars, stripped helium stars or even triples. While this is expected for Cyg X-1, confirming our methodology, our simulations also strongly suggest that HD 130298 could be another candidate to host a stellar-mass BH.

The GammaTPC is an MeV-scale single-phase liquid argon time-projection-chamber gamma-ray telescope concept with a novel dual-scale pixel-based charge-readout system. It promises to enable a significant improvement in sensitivity to MeV-scale gamma-rays over previous telescopes. The novel pixel-based charge readout allows for imaging of the tracks of electrons scattered by Compton interactions of incident gamma-rays. The two primary contributors to the accuracy of a Compton telescope in reconstructing an incident gamma-ray's original direction are its energy and position resolution. In this work, we focus on using deep learning to optimize the reconstruction of the initial position and direction of electrons scattered in Compton interactions, including using probabilistic models to estimate predictive uncertainty. We show that the deep learning models are able to predict locations of Compton scatters of MeV-scale gamma-rays from simulated pixel-based data to better than 0.6 mm RMS error, and are sensitive to the initial direction of the scattered electron. We compare and contrast different deep learning uncertainty estimation algorithms for reconstruction applications. Additionally, we show that event-by-event estimates of the uncertainty of the locations of the Compton scatters can be used to select those events that were reconstructed most accurately, leading to improvement in locating the origin of gamma-ray sources on the sky.

We revisit the mechanism of vortex unpinning caused by the neutron-vortex scattering in the inner crust of a pulsar. The strain energy released by the crustquake is assumed to be absorbed in some part of the inner crust and causes pair-breaking quasi-neutron excitations from the existing free neutron superfluid in the bulk of the inner crust. The scattering of these quasi-neutrons with the vortex core normal neutrons unpins a large number of vortices from the thermally affected regions and results in pulsar glitches. We consider the geometry of a cylindrical shell of the affected pinning region to study the implications of the vortex unpinning in the context of pulsar glitches. We find that a pulsar can release about $\sim 10^{11} - 10^{13}$ vortices by this mechanism. These numbers are equivalent to the glitch size of orders $\sim 10^{-11} - 10^{-9}$ for Vela-like pulsars with the characteristic age $\tau \simeq 10^4$ years. For Crab-like younger pulsars, the glitch size is increased by one order of magnitude. We also suggest a possibility of a vortex avalanche triggered by the movement of the unpinned vortices. A rough estimate of the glitch size caused by an avalanche shows an encouraging result.

The present thesis is devoted to the non-parametric reconstruction of some cosmological parameters using diverse observational datasets. The Universe is assumed to be spatially homogeneous and isotropic, thus described by the FLRW metric. The first chapter provides a brief introduction to cosmology and focuses on the reconstruction methods. An assessment of the cosmic distance-duality relation is discussed in chapter 2. In chapter 3, a non-parametric reconstruction of the cosmological jerk parameter is carried out. Chapter 4 explores the possibility of a non-gravitational interaction in the cosmic dark sector. In chapter 5, attempts are made to revisit a non-parametric reconstruction of the cosmic deceleration parameter using various combinations of recently updated background datasets and the growth rate of structure measurements from the redshift-space distortions. Finally, chapter 6 contains the concluding remarks and relevant discussion regarding the overall work presented in the thesis.

We calculate the electron-positron pair production rate at the base of the jet of Cyg X-1 by collisions of photons from its hot accretion flow using the measurement of its average soft gamma-ray spectra by the Compton Gamma Ray Observatory and INTEGRAL satellites. We have found that this rate approximately equals the flow rate of the leptons emitting the observed synchrotron radio-to-IR spectrum of the jet, calculated using an extended jet model following that of Blandford \& Koenigl. This coincidence shows the jet composition is likely to be pair-dominated. The same coincidences were found before in the microquasar MAXI J1820+070 and in the radio galaxy 3C 120, which shows that the considered mechanism can be universal for at least some classes of relativistic jets. Furthermore, we recalculate the jet power of Cyg X-1. The presence of pairs can strongly reduces the power in the bulk motion of ions, which then limits the parameter space at which the jet can power the $\sim$5-pc nebular structure present in its vicinity.

The cessation of AGN activity in radio galaxies leads to a remnant phase during which jets are no longer sustained, but lobes can be detected for a period of time before they fade away due to radiative and dynamical energy losses. The time-scale of the remnant phase and AGN duty cycle are vital to understanding the evolution of radio galaxies. In this paper, we report new band-3 observations with the upgraded Giant Meterwave Radio Telescope (uGMRT) for five remnant radio galaxies. Our uGMRT observations reveal emission of low-surface-brightness in all five remnants with 400 MHz surface brightness in the range of 36$-$201 mJy arcmin$^{-2}$. With band-3 uGMRT observations, we discover wing-shaped radio morphology in one of our sample sources. Using radio observations at 150 MHz, 325 MHz, 400 MHz, and 1.5 GHz, we model the radio spectral energy distributions (SEDs) of our sample sources with the continuous injection-off model (CI$_{\rm OFF}$), that assumes an active phase with continuous injection followed by a remnant phase. We obtain total source ages ($t_{\rm s}$) in the range of 20.3 Myr to 41.4 Myr with $t_{\rm OFF}$/$t_{\rm s}$ distributed in the range of 0.16 to 0.63, which in turn suggests them to belong to different evolutionary phases. We note that, in comparison to the remnants reported in the literature, our sample sources tend to show lower spectral ages that can be explained by the combined effects of more dominant inverse Compton losses for our sources present at the relatively higher redshifts and possible rapid expansion of lobes in their less dense environments.

We present the design and performance of the TANSPEC, a medium-resolution $0.55-2.5~\mu$m cryogenic spectrometer and imager, now in operation at the 3.6-m Devasthal Optical Telescope (DOT), Nainital, India. The TANSPEC provides three modes of operation which include, photometry with broad- and narrow-band filters, spectroscopy with short slits of 20$^{\prime \prime}$ length and different widths (from 0.5$^{\prime \prime}$ to 4.0$^{\prime \prime}$) in cross-dispersed mode at a resolving power R of $\sim$2750, and spectroscopy with long slits of 60$^{\prime \prime}$ length and different widths (from 0.5$^{\prime \prime}$ to 4.0$^{\prime \prime}$) in prism mode at a resolving power R of $\sim$100-350. TANSPEC's imager mode provides a field of view of 60$^{\prime \prime} \times 60^{\prime \prime}$ with a plate scale of 0.245$^{\prime \prime}$/pixel on the 3.6-m DOT. The TANSPEC was successfully commissioned during April-May 2019 and the subsequent characterization and astronomical observations are presented here. The TANSPEC has been made available to the worldwide astronomical community for science observations from October 2020.

We present MUSE integral field stellar and ionized velocity maps for a sample of 14 barred galaxies. Most of these objects exhibit "S"-shape iso-velocities in the bar region indicative of the presence of streaming motions in the velocity fields. % By applying circular rotation models we observe that bars leave symmetric structures in the residual maps of the stellar velocity. %which demonstrates the capabilities of the MUSE instrument for detecting kinematic bar signatures. % We built non-circular rotation models using the \xs~tool to characterize the observed velocity fields. In particular we adopt bisymmetric models and a harmonic decomposition for a bar potential for describing the non-axisymmetric velocities. We find that both models reproduce the observed kinematic features. % The position angle of the oval distortion estimated from the bisymmetric model correlates with the photometric bar position angle $(\rho_{pearson} = 0.95)$, which suggest that non-circular velocities are caused by the bar. However because of the low amplitudes of the $s_3$ harmonic we can not rule out radial flows as possible source. % Because of the weak detection of \ha~in our objects we are not able to compare gas to stellar non-circular motions in our sample, although we show that when galaxies are gas rich the oval distortion is also observed but with larger amplitudes. % Finally, we do not find evidence that the amplitude of the non-circular motions is dependent on the bar size, stellar mass or the global SFR.

After the detection of GRB 170817A, the first unambiguous off-axis gamma-ray burst (GRB), several studies tried to understand the structure of GRB jets. The initial jet structure (directly produced by the central engine) can be partially preserved, or can be completely modified by the interaction with the environment. In this study, we perform three-dimensional, special relativistic hydrodynamics simulations of long GRB jets evolving through a massive progenitor star. Different jet scenarios were considered: Top-hat, Gaussian jets dominated by pressure or by kinetic energy, as well as a model of a supernova (SN) plus a jet both propagating through the progenitor. We found that, while propagating inside the progenitor star, jets with different initial structures are nearly indistinguishable. Kinetic dominated jets are faster and more collimated than pressure dominated jets. The dynamics of jets inside the progenitor star strongly depends on the presence of an associated SN, which can strongly decelerate the jet propagation. We show that the initial structure of GRB jets is preserved, or not, mainly depending on the jet duration ($t_j$) and on the jet break-out time ($t_{\rm bo}$) from the dense environment. The initial structure is preserved in long-lasting jets ($t_j\gg t_{\rm bo}$), while the interaction with the medium shapes short duration jets ($t_j\ll t_{\rm bo}$). Future observations of a large sample of off-axis GRB will constrain the density stratification of the environment, the initial structure of the jet, and the physics of the central engine itself.

Astrophysical disks that are sufficiently cold and dense are linearly unstable to the formation of axisymmetric rings as a result of the disk's gravity. In practice, spiral structures are formed, which may in turn produce bound fragments. We study a nonlinear dynamical path that can explain the development of spirals in a local model of a gaseous disk on the subcritical side of the gravitational instability bifurcation. Axisymmetric equilibria can be radially periodic or localized, in the form of standing solitary waves. The solitary solutions have an energy slightly larger than a smooth disk. They are further unstable to non-axisymmetric perturbations with a wide range of azimuthal wavenumbers. The solitary waves may act as a pathway to spirals and fragmentation.

The detection with of the Atacama Large Millimeter Array (ALMA) of dust-rich high redshift galaxies whose cold dust emission is spatially disconnected from the ultraviolet emission bears a challenge for modelling their spectral energy distributions (SED) with codes based on an energy budget between the stellar and dust components. We test the validity of energy balance modelling on a nearby resolved galaxy with vastly different ultraviolet and infrared spatial distributions and infer what information can be reliably retrieved from the analysis of the full spectral energy distribution. We use 15 broadband images of the Antennae Galaxies ranging from far-ultraviolet to far-infrared and divide Arp 244 into 58 square ~1 kpc$^2$ regions. We fit the data with CIGALE to determine the star formation rate, stellar mass and dust attenuation of each region. We compare these quantities for the addition of the 58 regions to the ones obtained for Arp 244 as a whole and find that both estimates are consistent within one sigma. We present the spatial distribution of these physical parameters as well as the shape of the attenuation curve across the Antennae Galaxies . We also observe a flattening of the attenuation curves with increasing attenuation and dust surface density in agreement with the predictions of hydrodynamical simulations coupled with radiative transfer modelling.

Here we present the detailed investigation of AGNs in two Seyfert 1 galaxies, LEDA 3095839 and VII Zw 244. Both of them were observed within the photometric reverberation mapping project in Special Astrophysical Observatory of the Russian Academy of Sciences (SAO RAS), during which we earlier obtained the SMBHs masses. After that, both galaxies were observed in spectropolarimetric and polarimetric modes on the BTA 6 m telescope of the SAO RAS with the focal reducer SCORPIO-2. The linear polarization of the continuum and broad Balmer lines has been measured. It was found that (i) there were no signs of equatorial scattering in the LEDA 3095839 galaxy in the broad H-alpha line, and we were able to estimate the value of SMBH spin and the magnetic field strength in the disk from the level of continuum polarization; (ii) for the galaxy VII Zw 244, the presence of equatorial scattering was shown, due to which the mass of the SMBH was independently measured, the inclination angle of the system was obtained, and the value of the spin was estimated.

This work presents the first results of the Deep IFS View of Nuclei of Galaxies (DIVING$^\mathrm{3D}$) survey. We analysed the nuclear emission-line spectra of a sub-sample we call mini-DIVING$^\mathrm{3D}$, which includes all Southern galaxies with B < 11.2 and |b| > 15 degrees. We verified that $23\% \pm 4\%$ of the galaxies show nuclear emission-line properties characteristic of Low Ionization Nuclear Emission-Line Regions (LINERs). Diagnostic diagram analysis reveals an apparent dichotomy, not detected in previous studies, between objects classified as H II regions and as LINERs or Seyferts, with very few galaxies classified as transition objects. A possible explanation for this result is that at least part of the transition objects are composite systems, with a central LINER contaminated by the emission from circumnuclear H II regions. The higher spatial resolution of the DIVING$^\mathrm{3D}$ survey, in comparison with previous studies, allowed us to isolate the nuclear emission from circumnuclear contaminations, reducing the number of transition objects. We also propose an alternative scenario, in which the emission-line spectra of some transition objects are the result of shock heating by central outflows, together with photoionization by young stars. Clear evidence of active galactic nuclei (AGNs), in the optical and X-ray spectral bands, were detected in 69$\%$ of the LINERs in the mini-DIVING$^\mathrm{3D}$ sample. Considering the entire mini-DIVING$^\mathrm{3D}$ sample, evidence of AGNs were detected in 65$\%$ of the objects.

A discrete and exact algorithm for obtaining planetary systems is derived in a recent article (Eur. Phys. J. Plus 2022, 137:99). Here the algorithm is used to obtain planetary systems with forces different from the Newtonian inverse square gravitational forces. A Newtonian planetary system exhibits regular elliptical orbits, and here it is demonstrated that a planetary system with pure inverse forces also is stable and with regular orbits, whereas a planetary system with inverse cubic forces is unstable and without regular orbits. The regular orbits in a planetary system with inverse forces deviate, however, from the usual elliptical orbits by having revolving orbits with tendency to orbits with three or eight loops. Newton's Proposition 45 in $\textit{Principia}$ for the Moon's revolving orbits caused by an additional attraction to the gravitational attraction is confirmed, but whereas the additional inverse forces stabilize the planetary system, the additional inverse cubic forces can destabilize the planetary system at a sufficient strength.

Interactions of ultra-high energy cosmic rays (UHECRs) accelerated in specific astrophysical environments have been shown to shape the energy production rate of nuclei differently from that of the secondary neutrons escaping from the confinement zone. Here, we aim at testing a generic scenario of in-source interactions through a phenomenological modeling of the flux and composition of UHECRs. We fit a model in which nucleons and nuclei follow different particle energy distributions to the all-particle energy spectrum, proton spectrum below the ankle energy and distributions of maximum shower depths above this energy, as inferred at the Pierre Auger Observatory. We obtain that the data can be reproduced using a spatial distribution of sources that follows the density of extragalactic matter on both local and large scales, providing hence a realistic set of constraints for the emission mechanisms in cosmic accelerators, for their energetics and for the abundances of elements at escape from their environments. While the quasi mono-elemental increase of the cosmic-ray mass number observed on Earth from ${\simeq}\: 2\:$EeV up to the highest energies calls for nuclei accelerated with a hard spectral index, the inferred flux of protons down to ${\simeq}\: 0.6\:$ EeV is shown to require for this population a spectral index significantly softer than that generally obtained up to now. We demonstrate that modeling UHECR data across the ankle substantiate the conjecture of in-source interactions in a robust statistical framework, although pushing the mechanism to the extreme.

Context. The Galactic Center (GC) region has been studied in gamma rays in the past decades - the GC excess detected by Fermi is still not fully understood and the first detection of a PeVatron by H.E.S.S. indicates the existence of sources that can accelerate cosmic rays up to a PeV or higher. Aims. In this paper, we are investigating the origin of the PeVatron emission detected by H.E.S.S. by, for the first time, simulating cosmic rays in the GC in a realistic three-dimensional gas and photon field distribution and large-scale magnetic field. Methods. We solve the 3D transport equation with an anisotropic diffusion tensor using the approach of stochastic differential equations as implemented in the propagation software CRPropa 3.1 (Merten et al. 2017). We test five different source distributions for four different configurations of the diffusion tensor, i.e. with ratios of the perpendicular to parallel components $\epsilon = 0.001,\ 0.01,\ 0.1,\ 0.3$. Results. We find that the two-dimensional distribution of gamma rays as measured by H.E.S.S. is best fit by a model that considers three cosmic-ray sources, i.e. a central source, the SNR G0.9+0.1 and the source HESS J1746-285. The fit indicates that propagation is dominated by parallel diffusion with $ \espilon = 0.001 $. Conclusions. We find that the 3D propagation in the 3D gas and B-field configurations taken from Guenduez et al. (2020) can explain the general features of the data well. We predict that CTA should be able to identify emission from SgrB2 and the six dust ridge clouds that are included in our simulations and that should be detectable with a the expected resolution of CTA of 0.033$^\circ$.

Quasi-periodic fast-propagating (QFP) magnetosonic wave trains are commonly observed in the low corona at extreme ultraviolet wavelength bands. Here, we report the first white-light imaging observation of a QFP wave train propagating outwardly in the outer corona ranging from 2 to 4 solar Radii. The wave train was recorded by the Large Angle Spectroscopic Coronagraph on board the Solar and Heliospheric Observatory, and it was associated with a GOES M1.5 flare in NOAA active region AR12172 at the southwest limb of the solar disk. Measurements show that the speed and period of the wave train were about 218 km/s and 26 minutes, respectively. The extreme ultraviolet imaging observations taken by the Atmospheric Imaging Assembly on board the Solar Dynamic Observatory reveals that in the low corona the QFP wave train was associated with the failed eruption of a breakout magnetic system consisting of three low-lying closed loop systems enclosed by a high-lying large-scale one. Data analysis results show that the failed eruption of the breakout magnetic system was mainly because of the magnetic reconnection occurred between the two sided low-lying closed-loop systems. This reconnection enhances the confinement capacity of the magnetic breakout system because the upward-moving reconnected loops continuously feed new magnetic fluxes to the high-lying large-scale loop system. For the generation of the QFP wave train, we propose that it could be excited by the intermittent energy pulses released by the quasi-periodic generation, rapid stretching and expansion of the upward-moving, strongly bent reconnected loops.

We investigate the conditional colour-magnitude distribution (CCMD), namely the colour-magnitude distribution at fixed halo mass, of the central galaxies in semi-analytic galaxy formation model (SAM) and hydrodynamic simulations. We analyse the CCMD of central galaxies in each halo mass bin with the Gaussian mixture model and find that it can be decomposed into red and blue components nearly orthogonal to each other, a red component narrow in colour and extended in magnitude and a blue component narrow in magnitude and extended in colour. We focus on the SAM galaxies to explore the origin of the CCMD components by studying the relation between central galaxy colour and halo or galaxy properties. Central galaxy colour is correlated with halo assembly properties for low mass haloes and independent of them for high mass haloes. Galaxy properties such as central supermassive black hole mass, cold gas mass, and gas specific angular momentum can all impact central galaxy colour. These results are corroborated by an alternative machine learning analysis in which we attempt to predict central galaxy colour with halo and galaxy properties. We find that the prediction for colours of central galaxies can be significantly improved using both halo and galaxy properties as input compared to using halo properties alone. With the halo and galaxy properties considered here, we find that subtle discrepancies remain between predicted and original colour distribution for low mass haloes and that no significant determining properties are identified in massive haloes, suggesting modulations by additional stochastic processes in galaxy formation.

Homochirality is a generic and unique property of all biochemical life and is considered a universal and agnostic biosignature. Upon interaction with unpolarized light, homochirality induces fractional circular polarization in the light that is scattered from it, which can be sensed remotely. As such, it can be a prime candidate biosignature in the context of future life-detection missions and observatories. The linear polarizance of vegetation is also sometimes envisaged as a biosignature, although it does not share the molecular origin as circular polarization. It is known that the linear polarization of surfaces is strongly dependent on the phase angle. The relation between the phase angle and circular polarization stemming from macromolecular assemblies, such as in vegetation, however, remained unclear. We demonstrate in this study using the average of 27 different species that the circular polarization phase angle dependency of vegetation induces relatively small changes in spectral shape and mostly affects the signal magnitude. With these results we underline the use of circular spectropolarimetry as a promising agnostic biosignature complementary to the use of linear spectropolarimetry and scalar reflectance.

Earth is depleted in volatile elements relative to chondritic meteorites, its possible building blocks. The extent of this depletion increases with decreasing condensation temperature, and is approximated by a cumulative normal distribution, unlike that in any chondrite. However, moderately volatile elements, occupying the mid-range of the distribution, have chondritic isotope ratios, contrary to that expected from loss by partial vaporisation/condensation. Here we reconcile these observations by showing, using N-body simulations, that Earth accreted stochastically from many precursor bodies whose variable compositions reflect the temperatures at which they formed. Impact-induced atmospheric loss was efficient only when the proto-Earth was small, and elements that accreted thereafter retain near-chondritic isotope ratios. Earth's composition is reproduced when initial temperatures of planetesimal- to embryo-sized bodies are set by disk accretion rates of (1.08 $\pm$ 0.17) $\times$ 10$^{-7}$ solar masses/yr, although they may be perturbed by $^{26}$Al heating on bodies formed at different times. The model implies a heliocentric gradient in composition and rapid planetesimal formation within $\sim$ 1 Myr, in accord with radiometric volatile depletion ages of Earth.

The post-Keplerian(PK) parameters inferred from pulsar timing provide a convenient way to test Einstein's general theory of relativity. However, before obtaining a pure orbital decay $\dot{P}_b$ induced by gravitational wave radiation, which is one of the PK parameters, a number of factors need to be accounted for carefully. The effect of tidal dissipation on $\dot{P}_b$ has been thought of as negligible. Here, we investigate the data for possible effects of tidal dissipation on $\dot{P}_b$. The possibility of the tidal dissipation as a contributor to $\dot{P}_b$ in a large sample of millisecond pulsar binaries is investigated in detail. We collected a large sample of pulsar binaries with measured $\dot{P}_b$. All of the systems are millisecond pulsars. The residual $\dot{P}^{Res}_b$ of these systems was obtained by subtracting the three normal effects, that is to say the effect of Shklovskii, line-of-sight acceleration, and gravitational radiation. Assuming that tidal dissipation is responsible for such a residual $\dot{P}^{Res}_b$, the tidal parameters of these systems can be calculated and compared with the tidal models. The residual $\dot{P}^{Res}_b$ is distributed over the half positive and half negative. The dynamical tidal model can explain the residual $\dot{P}_b$ of millisecond pulsar-white dwarf binaries. And the Love number of the main-sequence companion of \object{PSR J1227-4853} can be derived as a reasonable value $k_2=0.177^{+0.098}_{-0.058}$ with the equilibrium tidal model. Those results are compatible with the scenario of tidal dissipation. Additionally, a weak correlation between the tidal parameter and orbital period is revealed, likely originating in the tidal process of the recycled stage which is worthy of further investigation.

Several redback and black widow millisecond pulsar binaries have episodes of flaring in X-rays and optical. We initially detected such behavior from the Fermi selected redback candidate 1FGL J0523.5$-$2529 during optical time-series monitoring. Triggered observations with the Neil Gehrels Swift Observatory over the next $\approx100$ days showed episodic flaring in X-rays with luminosity up to $8\times10^{33}$ erg s$^{-1}$ ($\sim100$ times the minimum), and a comparable luminosity in the optical/UV, with similar power-law spectra of $f_{\nu}\propto\nu^{-0.7}$. These are the most luminous flares seen in any non-accreting "spider" pulsar system, which may be related to the large size of the companion through the fraction of the pulsar wind that it or its ablated wind intercepts. Simultaneously with an optical flare, we see Balmer-line and He I emission, not previously known in this object, which is evidence of a stellar wind that may also inhibit detection of radio pulsations. The quiescent optical light curves, while dominated by ellipsoidal modulation, show evidence of variable non-uniform temperature that could be due either to large starspots or asymmetric heating of the companion by the pulsar. This may explain a previous measurement of unusual non-zero orbital eccentricity as, alternatively, distortion of the radial-velocity curve by the surface temperature distribution of the large companion.

In the standard picture, episodes of luminous quasar activity are directly related to supermassive black hole (SMBH) growth. The ionising radiation emitted over a quasar's lifetime alters the ionisation state of the surrounding intergalactic medium (IGM), enhancing the Ly$\alpha$ forest transmission -- so-called proximity effect -- which can be observed in absorption spectra of background sources. Owing to the finite speed of light, the transverse direction of the proximity effect is sensitive to the quasar's radiative history, resulting in `light echoes' that encode the growth history of the SMBH on Myr-timescales. In this paper, we introduce a new technique to photometrically map this quasar light echoes using Ly$\alpha$ forest tomography by using a carefully selected pair of narrow-band filters. A foreground narrow-band filter is used to measure Ly$\alpha$ forest transmission along background galaxies selected as Ly$\alpha$ emitters by a background narrow-band filter. This novel double narrow-band tomographic technique utilises the higher throughput and wider field of view of imaging over spectroscopy to efficiently reconstruct a two-dimensional map of Ly$\alpha$ forest transmission around a quasar. We present a fully Bayesian framework to measure the luminous quasar lifetime of a SMBH from photometric IGM tomography, and examine the observational requirements. This new technique provides an efficient strategy to map a large area of the sky with a modest observing time and to identify interesting regions to be examined by further deep 3D follow-up spectroscopic Ly$\alpha$ forest tomography.

Cataclysmic variables (CVs) exhibit a plethora of variable phenomena many of which require long, uninterrupted light curves to reveal themselves in detail. The month long datasets provided by TESS are well suited for this purpose. TESS has the additional advantage to have observed a huge number of stars, among them many CVs. Here, a search for periodic variations in a sample of CVs of the novalike and old novae subtypes is presented. In 10 of the 15 targets either previously unseen positive or negative superhumps or unusual features in known superhumps are identified. The TESS light curves demonstrate that the occurrence of superhumps in these types of CVs is not an exception but quite common. For 8 systems new or improved values for the orbital period are measured. In TV Col the long-sought optical manifestation of the white dwarf spin period is first seen in form of its orbital sideband. The mystery of multiple photometric periods observed in CP Pup in the past is explained by irregularly occurring anomalous states which are reflected in the light curve.

Recent astrochemical models and experiments have explained that complex organic molecules (COMs; molecules composed of six or more atoms) are produced on the dust grain mantles in cold and dense gas in prestellar cores. However, the detailed chemical processes and the roles of physical conditions on chemistry are still far from understood. To address these questions, we investigated twelve high-mass star-forming regions using the ALMA band 6 observations. They are associated with 44/95GHz class I and 6.7 GHz class II CH$_{3}$OH masers, indicative of undergoing active accretion. We found 28 hot cores with COMs emission among 68 continuum peaks at 1.3 mm and specified 10 hot cores associated with 6.7 GHz class II CH$_{3}$OH masers. Up to 19 COMs are identified including oxygen- and nitrogen-bearing molecules and their isotopologues in cores. The derived abundances show a good agreement with those from other low- and high-mass star-forming regions, implying that the COMs chemistry is predominantly set by the ice chemistry in the prestellar core stage. One clear trend is that the COMs detection rate steeply grows with the gas column density, which can be attributed to the efficient formation of COMs in dense cores. In addition, cores associated with a 6.7 GHz class II CH$_{3}$OH maser tend to be enriched with COMs. Finally, our results suggest that the enhanced abundances of several molecules in our hot cores could be originated by the active accretion as well as different physical conditions of cores.

We present mid-infrared properties of three submillimeter galaxies in the SMACS 0723 galaxy cluster field, combining the recently released images by JWST and archival data from HST and ALMA. The new JWST data allow resolving Sub-millimeter Galaxies (SMGs) for the first time and demonstrate that these have dust present both in the centers and in the disk, in contrast to the existing ALMA observations which only detect dust in the central region. For one of our sources, the best fitting surface brightness model has a disk with Sersic index of $\sim$ 0.9 and the residual image shows clear spiral arms, consistent with previous studies of SMG morphologies. The morphology results imply that the dust of SMGs that belong to the main sequence does not originate from major mergers. The dust attenuation in the central parts of galaxies will bias the half-light radius estimation, and we find the half light radius of the SMGs measured from the rest-frame NIR band could be $\sim$ 1.5 times smaller than the rest-frame V band, thus the previous stellar mass distribution based on the HST images alone may be affected by dust extinction. Finally, we are able to detect an SMG with $m_{\rm AB}\sim$ 27 in the mid-infrared, demonstrating the impressive sensitivity of JWST.

The Epoch of Reionization Spectrometer (EoR-Spec) is one of the instrument modules to be installed in the Prime-Cam receiver of the Fred Young Submillimeter Telescope (FYST). This six-meter aperture telescope will be built on Cerro Chajnantor in the Atacama Desert in Chile. EoR-Spec is designed to probe early star-forming regions by measuring the [CII] fine-structure lines between redshift z = 3.5 and z = 8 using the line intensity mapping technique. The module is equipped with a scanning Fabry-Perot interferometer (FPI) to achieve the spectral resolving power of about RP = 100. The FPI consists of two parallel and identical, highly reflective mirrors with a clear aperture of 14 cm, forming a resonating cavity called etalon. The mirrors are silicon-based and patterned with double-layer metamaterial anti-reflection coatings (ARC) on one side and metal mesh reflectors on the other. The double-layer ARCs ensure a low reflectance at one substrate surface and help tailor the reflectance profile over the FPI bandwidth. Here we present the design, fabrication processes, test setup, and characterization of silicon mirrors for the FPI.

For decades, physicists have analyzed various versions of a ``cosmic no-hair" conjecture, to understand under what conditions a spacetime that is initially spatially anisotropic and/or inhomogeneous will flow into an isotropic and homogeneous state. Wald's theorem, in particular, established that homogeneous but anisotropic spacetimes, if filled with a positive cosmological constant plus additional matter sources that satisfy specific energy conditions, will necessarily flow toward an (isotropic) de Sitter state at late times. In this paper we study the flow of homogeneous but anisotropic spacetimes toward isotropic states under conditions more general than those to which Wald's theorem applies. We construct an effective field theory (EFT) treatment for generic ``single-clock" systems in anisotropic spacetimes -- which are not limited to realizations compatible with scalar-field constructions -- and identify fixed points in the resulting phase space. We identify regions of this phase space that flow to isotropic fixed points -- including a de Sitter fixed point -- even in the absence of a bare cosmological constant, and for matter sources that do not obey the energy conditions required for Wald's theorem. Such flows into de Sitter reveal the emergence of an effective cosmological constant.

We present a simple criterion to predict the explodability of massive stars based on the density and entropy profiles before collapse. If a pronounced density jump is present near the Si/O interface, the star will likely explode. We develop a quantitative criterion by using $\sim 400$ 1D simulations where $\nu$-driven turbulence is included via mixing length theory. This criterion correctly identifies the outcome of the supernova $\sim 95 \%$ of the time. We also find no difference in how this criterion performs on two different sets of progenitors, evolved using two different stellar evolution codes: FRANEC and KEPLER. The explodability as a function of mass of the two sets of progenitors is very different, showing: (i) that uncertainties in the stellar evolution prescriptions influence the predictions of supernova explosions; (ii) the most important properties of the pre-collapse progenitor that influence the explodability are its density and entropy profiles. We highlight the importance that $\nu$-driven turbulence plays in the explosion by comparing our results to previous works.

We study spatially resolved scaling relations among stars, dust, and gas in ten nearby spiral galaxies. In a preceding paper Abdurro'uf et al. (2022), we have derived spatially resolved properties of the stellar population and dust by panchromatic spectral energy distribution (SED) fitting using piXedfit. Now, we investigate resolved star formation ($\Sigma_{\rm H_{2}}$--$\Sigma_{\rm SFR}$--$\Sigma_{*}$) and dust scaling relations. While the relations with all sub-galactic regions of the galaxies are reasonably tight ($\sigma \lesssim 0.3$ dex), we find that most of the scaling relations exhibit galaxy-to-galaxy variations in normalization and shape. Only two relations of $\Sigma_{\rm dust}$--$\Sigma_{\rm gas}$ and $\Sigma_{\rm dust}$--$\Sigma_{\rm H_{2}}$ do not show noticeable galaxy-to-galaxy variations among our sample galaxies. We further investigate correlations among the scaling relations. We find significant correlations among the normalization of the $\Sigma_{\rm H_{2}}$--$\Sigma_{\rm SFR}$--$\Sigma_{*}$ relations, which suggest that galaxies with higher levels of resolved $\text{H}_{2}$ fraction ($f_{\rm H_{2}}$) tend to have higher levels of resolved star formation efficiency (SFE) and specific star formation rate (sSFR). We also observe that galaxies with higher levels of resolved dust-to-stellar mass ratios tend to have higher levels of resolved sSFR, SFE, and $f_{\rm H_{2}}$. Moreover, we find that galaxies with higher global sSFR and less compact morphology tend to have higher levels of the resolved sSFR, SFE, and $f_{\rm H_{2}}$, which can explain the variations in the normalization of the $\Sigma_{\rm H_{2}}$--$\Sigma_{\rm SFR}$--$\Sigma_{*}$ relationships. Overall, we observe indications of the contributions of both global and local factors in governing the star formation process in galaxies.

We use the X-ray luminosity relation of radio-loud quasars (RLQs) to measure these luminosity distances as well as estimate cosmological parameters. We adopt four parametric models of X-ray luminosity to test luminosity correlation for RLQs and radio-intermediate quasars (RIQs) and give these cosmological distances. By Bayesian information criterion (BIC), the data suggest that the luminosity relation ${L_X} \propto L_{UV}^{{\gamma_{uv}}}L_{Radio}^{\gamma_{radio}'}$ for RLQs has better goodness of fit, relative to other models, which can be interpreted as this relation being preferred for RLQs. Meanwhile, we compare the results from flat-spectrum radio-loud quasars (FSRLQs) and steep-spectrum radio-loud quasars (SSRLQs), which indicate that their luminosity correlations are not exactly the same. We also consider dividing the RLQs sample into various redshift bins, which can be used to check if the X-ray luminosity relation depends on the redshift. Finally, we apply a combination of RLQs and SNla Pantheon to verify the nature of dark energy concerning whether or not its density deviates from the constant, and give the statistical results.

We present the optical photometric and spectroscopic observations of the nearby Type Ia supernova (SN) 2021hpr. The observations covered the phase of $-$14.37 to +63.68 days relative to its maximum luminosity in the $B$ band. The evolution of multiband light/color curves of SN 2021hpr is similar to that of normal Type Ia supernovae (SNe Ia) with the exception of some phases, especially a plateau phase that appeared in the $V-R$ color curve before peak luminosity, which resembles that of SN 2017cbv. The first spectrum we observed at t $\sim -$14.4 days shows a higher velocity for the Si II $\lambda$6355 feature ($\sim$ 21,000 km s$^{-1}$) than that of other normal Velocity (NV) SNe Ia at the same phase. Based on the Si II $\lambda$6355 velocity of $\sim$ 12,420 km s$^{-1}$ around the maximum light, we deduce that SN 2021hpr is a transitional object between high velocity (HV) and NV SNe Ia. Meanwhile, the Si II $\lambda$6355 feature shows a high velocity gradient (HVG) of about 800 km s$^{-1}$ day$^{-1}$ from roughly $-$14.37 to $-$4.31 days relative to the $B$-band maximum, which indicates that SN 2021hpr can also be classified as an HVG SN Ia. The evolution of SN 2021hpr is similar to that of SN 2011fe. Including SN 2021hpr, there have been six supernovae observed in the host galaxy NGC 3147, and the supernovae explosion rate in the last 50 yr is slightly higher for SNe Ia, while lower for SNe Ibc and SNe II it is lower than expected rate from the radio data. Inspecting the spectra, we find that SN 2021hpr has a metal-rich (12 + log(O/H) $\approx$ 8.648) circumstellar environment, where HV SNe tend to reside. Based on the decline rate of SN 2021hpr in the $B$ band, we determine the distance modulus of the host galaxy NGC 3147 using the Phillips relation to be 33.46 $\pm$ 0.21 mag, which is close to that found by previous works.

Protostellar outflows are one of the most outstanding features of star formation. The observational studies over several decades successfully demonstrated that outflows are ubiquitously associated with low- and high-mass protostars in the solar-metallicity Galactic condition. However, the environmental dependence of protostellar outflow properties is still poorly understood, particularly in the low-metallicity regime. Here we report the first detection of a molecular outflow in the Small Magellanic Cloud with 0.2 $Z_{\odot}$, using ALMA observations at a spatial resolution of 0.1 pc toward the massive protostar Y 191. The bipolar outflow is nicely illustrated by high-velocity wings of CO(3-2) emission with $\gtrsim$15 km s$^{-1}$. The evaluated properties of the outflow (momentum, mechanical force, etc.) are consistent with those of the Galactic counterparts. Our results suggest that the molecular outflows, i.e., the guidepost of the disk accretion at the small scale, might be universally associated with protostars across the metallicity range of$\sim$0.2-1 $Z_{\odot}$.

HARMONI is the first light visible and near-IR integral field spectrograph for the ELT covering a large spectral range from 450nm to 2450nm with resolving powers from 3500 to 18000 and spatial sampling from 60mas to 4mas. It can operate in two Adaptive Optics modes-SCAO and LTAO-or with no AO. The project is preparing for Final Design Reviews. The laser Tomographic AO (LTAO) system provides AO correction with very high sky-coverage thanks to two systems: the Laser Guide Star Sensors (LGSS) and the Natural Guide Star Sensors (NGSS). LGSS is dedicated to the analysis of the wavefront coming from 6 laser guide stars created by the ELT. It is made of 6 independent wavefront sensor (WFS) modules mounted on a rotator of 600mm diameter to stabilise the pupil onto the microlens array in front of the detector. The optical design accepts elongated spots of up to 16 arcsec with no truncation using a CMOS detector from SONY. We will present the final optical and mechanical design of the LGSS based on freeform lenses to minimize the numbers of optical components and to accommodate for the diversity of sodium layer configurations. We will focus on rotator design, illustrating how we will move 1 tons with 90" accuracy in restrictive environment. Finally, we will present the strategy to verify the system in HARMONI context. The main challenge for the verification being how to test an AO system without access to the deformable mirror, part of the ELT.

In the first few Myr the massive stars dynamically interact, produce runaways and affect the initial binary population. Observing and interpreting the dynamics of young massive clusters is key to our understanding of the star formation process and predicting the outcome of stellar evolution. We have studied NGC6611 in the Eagle Nebula (M16), a young massive cluster hosting 19 O stars. We used Gaia EDR3 data to determine the membership, age, cluster dynamics and the kinematics of the massive stars including runaways. The membership analysis yields 137 members located at a mean distance of 1706 $\pm$ 7 pc. The colour - absolute magnitude diagram reveals a blue and a red population of pre-main-sequence stars, consistent with two distinct populations of stars. In line with earlier studies, the youngest population has a mean extinction $A_V$ = 3.6 $\pm$ 0.1 mag and an age = 1.3 $\pm$ 0.2 Myr, while the older population of stars has a mean extinction $A_V$ = 2.0 $\pm$ 0.1 mag and an age = 7.5 $\pm$ 0.4 Myr. The latter population is more spatially extended than the younger generation of stars. We argue that most of the OB stars belong to the younger population. We identify 8 runaways originating from the center of NGC6611, consistent with the dynamical ejection scenario. We show that ~ 50% of the O stars have velocities comparable to or greater than the escape velocity. These O stars can be traced back to the center of NGC6611 with kinematic ages ranging from 0 to 2 Myr. This suggests that dynamical interactions played an important role in the early evolution of NGC6611, which is surprising considering the low current stellar density. Comparing this to simulations of young massive clusters, the required initial radius of 0.1-0.5 pc is not consistent with that of NGC6611. The O stars could have initially formed in wide binaries and possibly harden through dynamical interactions.

One of the main complications for the interpretation of reflectance spectra of airless planetary bodies is surface alteration by space weathering caused by irradiation by solar wind and micrometeoroid particles. We aim to evaluate the damage to the samples from H and laser irradiation and relate it to the observed alteration in the spectra. We used olivine (OL) and pyroxene (OPX) pellets irradiated by 5 keV H ions and individual fs laser pulses and measured their visible (VIS) and near-infrared (NIR) spectra. We observed the pellets with scanning and transmission electron microscopy. We studied structural, mineralogical, and chemical modifications in the samples and connected them to changes in the reflectance spectra. In both minerals, H irradiation induces partially amorphous sub-surface layers containing small vesicles. In OL pellets, these vesicles are more tightly packed than in OPX ones. Related spectral change is mainly in the VIS spectral slope. Changes due to laser irradiation are mostly dependent on material's melting temperature. Only the laser-irradiated OL contains nanophase Fe particles, which induce detectable spectral slope change throughout the measured spectral range. Our results suggest that spectral changes at VIS-NIR wavelengths are mainly dependent on thickness of (partially) amorphous sub-surface layers. Amorphisation smooths microroughness, increasing the contribution of volume scattering and absorption over surface scattering. Soon after exposure to the space environment, the appearance of partially amorphous sub-surface layers results in rapid changes in the VIS spectral slope. In later stages (onset of micrometeoroid bombardment), we expect an emergence of nanoparticles to also mildly affect the NIR spectral slope. An increase in dimensions of amorphous layers and vesicles in the more space-weathered material will only cause band-depth variation and darkening.

We present the luminosity functions and host galaxy properties of the Dark Energy Survey (DES) core-collapse supernova (CCSN) sample, consisting of 69 Type II and 50 Type Ibc spectroscopically and photometrically-confirmed supernovae over a redshift range $0.045<z<0.25$. We fit the observed DES $griz$ CCSN light-curves and K-correct to produce rest-frame $R$-band light curves. We compare the sample with lower-redshift CCSN samples from Zwicky Transient Facility (ZTF) and Lick Observatory Supernova Search (LOSS). Comparing luminosity functions, the DES and ZTF samples of SNe II are brighter than that of LOSS with significances of 3.0$\sigma$ and 2.5$\sigma$ respectively. While this difference could be caused by redshift evolution in the luminosity function, simpler explanations such as differing levels of host extinction remain a possibility. We find that the host galaxies of SNe II in DES are on average bluer than in ZTF, despite having consistent stellar mass distributions. We consider a number of possibilities to explain this -- including galaxy evolution with redshift, selection biases in either the DES or ZTF samples, and systematic differences due to the different photometric bands available -- but find that none can easily reconcile the differences in host colour between the two samples and thus its cause remains uncertain.

Spectropolarimetic campaigns have established that large-scale magnetic fields are present at the surfaces of approximately 10% of massive dwarf stars. However, there is a dearth of magnetic field measurements for their deep interiors. Asteroseismology of gravity-mode pulsations combined with rotating magneto-hydrodynamical calculations of the early-B main-sequence star HD 43317 constrain its magnetic field strength to be approximately $5\times10^{5}$ G just outside its convective core. This proof-of-concept study for magneto-asteroseismology opens a new window into the observational characterisation of magnetic fields inside massive stars.

In recent years, it has been discovered that massive stars commonly exhibit a non-coherent form of variability in their light curves referred to as stochastic low frequency (SLF) variability. Various physical mechanisms can produce SLF variability in such stars, including stochastic gravity waves excited at the interface of convective and radiative regions, dynamic turbulence generated in the near-surface layers, and clumpy winds. Gravity waves in particular are a promising candidate for explaining SLF variability as they can be ubiquitously generated in main sequence stars owing to the presence of a convective core, and because they provide the large-scale predominantly tangential velocity field required to explain macroturbulence in spectral line fitting. Here, I provide an overview of the methods and results of studying SLF variability in massive stars from time series photometry and spectroscopy.

We present new high-contrast images in near-infrared wavelengths ($\lambda$ = 1.04, 1.24, 1.62, 2.18 and 3.78$\mu$m) of the young variable star CQ Tau, aiming to constrain the presence of companions in the protoplanetary disc. We reached a Ks-band contrast of 14 magnitudes with SPHERE/IRDIS at separations greater than 0."4 from the star. Our mass sensitivity curve rules out giant planets above 4 M$_{\rm Jup}$ immediately outside the spiral arms at $\sim$60 au and above 2--3 M$_{\rm Jup}$ beyond 100 au to 5$\sigma$ confidence assuming hot-start models. We do, however, detect four spiral arms, a double-arc and evidence for shadows in scattered light cast by a misaligned inner disc. Our observations may be explained by an unseen close-in companion on an inclined and eccentric orbit. Such a hypothesis would also account for the disc CO cavity, disturbed kinematics.

Millimeter astronomy provides valuable information on the birthplaces of planetary systems. In order to compare theoretical models with observations, the dust component has to be carefully calculated. Here, we aim to study the effects of dust entrainment in photoevaporative winds and the ejection and drag of dust due to effects caused by radiation from the central star. We improved and extended the existing implementation of a two-population dust and pebble description in the global Bern/Heidelberg planet formation and evolution model. Modern prescriptions for photoevaporative winds were used and we account for settling and advection of dust when calculating entrainment rates. In order to prepare for future population studies with varying conditions, we explore a wide range of disk-, photoevaporation-, and dust-parameters. We find that if dust can grow to pebble sizes, that is, if they are resistant to fragmentation or turbulence is weak, drift dominates and the entrained mass is small but larger than under the assumption of no vertical advection of grains with the gas flow. For the case of fragile dust shattering at velocities of 1 m/s - as indicated in laboratory experiments -, an order of magnitude more dust is entrained which becomes the main dust removal process. Radiation pressure effects disperse massive, dusty disks on timescales of a few 100 Myr. These results highlight the importance of dust entrainment in winds as a solid mass removal process. Furthermore, this model extension lies the basis for future statistical studies of planet formation in their birth environment.

Supernova remnants are believed to be the main sites where Galactic cosmic rays originate. This scenario, however, fails to explain some of the features observed in the cosmic-ray spectrum. Microquasars have been proposed as additional candidates, because their non-thermal emission indicates the existence of efficient particle acceleration mechanisms in their jets. A promising scenario envisages the production of relativistic neutrons in the jets, that decay outside the system injecting relativistic protons to the surroundings. The first investigations of this scenario suggest that microquasars might be fairly alternative cosmic-ray sources. We aim at assessing the role played by the degree of collimation of the jet on the cosmic-ray energetics in the neutron-carrier scenario, as well as the properties of the emission region. Our goals are to explain the Galactic component of the observed proton cosmic-ray spectrum at energies higher than $\sim 10$ GeV and to relate the mentioned jet properties with the power and spectral index of the produced cosmic rays. We find that collimated jets, with compact acceleration regions close to the jet base, are very efficient sources that could deliver a fraction of up to $\sim 0.01$ of their relativistic proton luminosity into cosmic rays. Collimation is the most significant feature regarding efficiency; a well collimated jet might be $\sim 4$ orders of magnitude more efficient than a poorly collimated one. The main feature of the presented mechanism is the production of a spectrum with a steeper spectral index ($\sim 2.3$ at energies up to $\sim 10$ TeV) than in the supernova scenario, and closer to what is observed. The predictions of our model may be used to infer the total contribution of the population of Galactic microquasars to the cosmic ray population, and therefore to quantitatively assess their significance as cosmic-ray sources.

Rings around exoplanets (exorings) are one of the most expected discoveries in exoplanetary research. There is an increasing number of theoretical and observational efforts for detecting exorings, but none of them have succeeded yet. Most of those methods focus on the photometric signatures of exorings during transits, whereas less attention has been paid to light diffusely reflected: what we denote here as the bright side of the light curve. This is particularly important when we cannot detect the typical stellar flux drop produced by transiting exoplanets. Here, we endeavour to develop a general method to model the variations on the light curves of both ringed non-transiting and transiting exoplanets. Our model (dubbed as Pryngles) simulates the complex interaction of luminous, opaque, and semitransparent objects in planetary systems, discretizing their surface with small circular plane discs that resemble sequins or spangles. We perform several numerical experiments with this model, and show its incredible potential to describe the light curve of complex systems under various orbital, planetary, and observational configurations of planets, moons, rings, or discs. As our model uses a very general approach, we can capture effects like shadows or planetary/ring shine, and since the model is also modular we can easily integrate arbitrarily complex physics of planetary light scattering. A comparison against existing tools and analytical models of reflected light reveals that our model, despite its novel features, reliably reproduces light curves under common circumstances. Pryngles source code is written in PYTHON and made publicly available.

The neutron star low-mass X-ray binary SWIFT J1749.4-2807 is the only known eclipsing accreting millisecond X-ray pulsar. In this manuscript we perform a spectral characterization of the system throughout its 2021, two-week-long outburst, analyzing 11 NICER observations and quasi-simultaneous XMM-Newton and NuSTAR single observations at the outburst peak. The broadband spectrum is well-modeled with a black body component with a temperature of $\sim$0.6 keV, most likely consistent with a hot spot on the neutron star surface, and a Comptonisation spectrum with power-law index $\Gamma \sim 1.9$, arising from a hot corona at $\sim$12 keV. No direct emission from the disc was found, possibly due to it being too cool. A high truncation radius for the disc, i.e., at $\sim$20--30 R$_{G}$ , was obtained from the analysis of the broadened profile of the Fe line in the reflection component. The significant detection of a blue-shifted Fe XXVI absorption line at $\sim$7 keV indicates weakly relativistic X-ray disc winds, which are typically absent in the hard state of X-ray binaries. By comparing the low flux observed during the outburst and the one expected in a conservative mass-transfer, we conclude that mass-transfer in the system is highly non-conservative, as also suggested by the wind detection. Finally, using the Nicer spectra alone, we followed the system while it was fading to quiescence. During the outburst decay, as the spectral shape hardened, the hot spot on the neutron star surface cooled down and shrank, a trend which could be consistent with the pure power-law spectrum observed during quiescence.

We dedicate this thesis to the study of signatures coming from the primordial epochs of the universe. We will focus in particular on Primordial Black Holes (PBHs), which may be formed from perturbations generated during inflation and might comprise a fraction of the dark matter in the universe. In the first part of the thesis, we will address the PBH properties at the time of formation, that are their masses, spins and abundance, and investigate the generation of Gravitational Wave (GW) signals during their production. In the second part, we will describe the PBHs evolution across the cosmic history due to their assemble in binaries, phases of baryonic mass accretion and clustering effects. We will then discuss GW signatures coming from their coalescence, compare these predictions with present GW data detected by the LIGO/Virgo Collaboration (LVC) and assess the role of future GW experiments like 3G detectors and LISA in discovering these objects. Finally, in the third part, we will investigate some aspects of the interplay between black holes and fundamental physics in the early universe, focusing on the role of GWs to shed light on their properties.

Future actively cooled space-borne observatories for the far-infrared, loosely defined as a 1--10 THz band, can potentially reach a sensitivity limited only by background radiation from the Universe. This will result in an increase in observing speed of many orders of magnitude. A spectroscopic instrument on such an observatory requires large arrays of detectors with a sensitivity expressed as a noise equivalent power NEP = 3 $\times 10^{-20}$ $W\surd{Hz}$. We present the design, fabrication, and characterisation of microwave kinetic inductance detectors (MKIDs) for this frequency range reaching the required sensitivity. The devices are based on thin-film NbTiN resonators which use lens-antenna coupling to a submicron-width aluminium transmission line at the shorted end of the resonator where the radiation is absorbed. We optimised the MKID geometry for a low NEP by using a small aluminium volume of $\approx$ 1$\mu m^3$ and fabricating the aluminium section on a very thin (100 nm) SiN membrane. Both methods of optimisation also reduce the effect of excess noise by increasing the responsivity of the device, which is further increased by reducing the parasitic geometrical inductance of the resonator. We measure the sensitivity of eight MKIDs with respect to the power absorbed in the detector using a thermal calibration source filtered in a narrow band around 1.55 THz. We obtain a NEP$_{exp}(P_{abs})\:=\:3.1\pm0.9\times10^{-20}\:W\surd{Hz}$ at a modulation frequency of 200 Hz averaged over all measured MKIDs. The NEP is limited by quasiparticle trapping. The measured sensitivity is sufficient for spectroscopic observations from future, actively cooled space-based observatories. Moreover, the presented device design and assembly can be adapted for frequencies up to $\approx$ 10 THz and can be readily implemented in kilopixel arrays.

In this work, we carry out a suite of specially-designed numerical simulations to shed further light on dark matter (DM) subhalo survival at mass scales relevant for gamma-ray DM searches, a topic subject to intense debate nowadays. Specifically, we have developed and employed an improved version of DASH, a GPU $N$-body code, to study the evolution of low-mass subhaloes inside a Milky Way-like halo with unprecedented accuracy, reaching solar-mass and sub-parsec resolution in our simulations. We simulate subhaloes with varying mass, concentration, and orbital properties, and consider the effect of the gravitational potential of the Milky Way galaxy itself. More specifically, we analyze the evolution of both the bound mass fraction and annihilation luminosity of subhaloes, finding that most subhaloes survive until present time, even though in some cases they lose more than 99% of their mass at accretion. Baryons in the host induce a much more severe mass loss, especially when the subhalo orbit is more parallel to the galactic disk. Many of these subhaloes cross the solar galactocentric radius, thus making it easier to detect their annihilation fluxes from Earth. We find subhaloes orbiting a DM-only halo with a pericentre in the solar vicinity to lose 70-90% of their initial annihilation luminosity at redshift zero, which increases up to 99% when baryons are also included in the host. We find a strong relation between subhalo's mass loss and the effective tidal field at pericentre. Indeed, much of the dependence on concentration, orbital parameters, host potential and baryonic components can be explained through this single parameter. In addition to shedding light on the survival of low-mass galactic subhaloes, our results can provide detailed predictions that will aid current and future quests for the nature of DM.

The calibration of redshift distributions for photometric samples using spectroscopic surveys is plagued by the difficulty in modelling the selection functions of spectroscopic surveys. In this work, we analyse how these selection functions impact redshift inference and quantify the induced biases using local calibration tests in photometry space. The study is carried out using simulations that mimic the radial selection function of a spectroscopic survey and an accompanying mock catalog of a photometric galaxy survey catalog. We use a self-organizing map to partition the photometry space and perform a local $\chi^2$ test to study the probability calibration of redshift inferences that use the spectroscopic data for calibration. The goal of this work is to investigate the effect of uncorrected selection functions in the calibration data on redshift prediction accuracy and critically discuss mitigation methods. In particular we test culling-based bias correction techniques, that aim to remove redshift calibration biases by identifying regions in photometry with few spectroscopic calibration data, and propose avenues for future research. We found that removing regions in color-magnitude space that are underpopulated with spectroscopic calibration data does not remove all biases in redshift inference induced by the selection function.

Unknown class distributions in unlabelled astrophysical training data have previously been shown to detrimentally affect model performance due to dataset shift between training and validation sets. For radio galaxy classification, we demonstrate in this work that removing low angular extent sources from the unlabelled data before training produces qualitatively different training dynamics for a contrastive model. By applying the model on an unlabelled data-set with unknown class balance and sub-population distribution to generate a representation space of radio galaxies, we show that with an appropriate cut threshold we can find a representation with FRI/FRII class separation approaching that of a supervised baseline explicitly trained to separate radio galaxies into these two classes. Furthermore we show that an excessively conservative cut threshold blocks any increase in validation accuracy. We then use the learned representation for the downstream task of performing a similarity search on rare hybrid sources, finding that the contrastive model can reliably return semantically similar samples, with the added bonus of finding duplicates which remain after pre-processing.

Dark matter in the form of compact objects with mass $M_{\rm CO} \gtrsim 10 M_{\odot}$ can be constrained by its dynamical effects on wide binary stars. Motivated by the recent interest in Primordial Black Hole dark matter, we revisit these constraints. We improve on previous studies in several ways. Specifically, we i) implement a physically motivated model for the initial wide-binary semi-major axis distribution, ii) include unbound binaries, and iii) take into account the uncertainty in the relationship between semi-major axis and observed angular separation. These effects all tend to increase the predicted number of wide binaries (for a given compact object population). Therefore the constraints on the halo fraction in compact objects, $f_{\rm co}$, are significantly weakened. We find the fraction of halo dark matter in compact objects is $f_{\rm co} < 1$ for $M_{\rm co} \approx 300 \, M_{\odot}$, tightening with increasing $M_{\rm co}$ to $f_{\rm co} < 0.26$ for $M_{\rm co} \gtrsim 1000 \, M_{\odot}$.

Continued star formation over the lifetime of a galaxy suggests that low metalicity gas is steadily flowing in from the circumgalactic medium. Also, cosmological simulations of large-scale structure formation imply that gas is accreted onto galaxies from the halo inside which they formed. Direct observations are difficult, but in recent years observational indications of gas inflows from a circumgalactic medium were obtained. Here we suggest an indirect observational probe: looking for large-scale (exceeding few kpc) turbulence caused by the accretion. As a specific example we consider an accretion flow coplanar with the galaxy disk, and argue that Kelvin-Helmholtz turbulence will be generated. We employ a semi-analytic model of turbulence and derive the expected turbulence power spectrum. The latter turns out to be of a distinctive shape that can be compared with observational power spectra. As an illustrative example we use parameters of the Milky Way galaxy.

An overview is provided of the scientific goals of the Magellanic Cloud component of the STScI Directors Discretionary UV initiative ULLYSES, together with the complementary spectroscopic survey XShootU (VLT/Xshooter) and other ancillary datasets. Together, ULLYSES and XShootU permit the first comprehensive, homogeneous study of wind densities and velocities in metal-poor massive stars, plus UV/optical spectroscopic libraries for population synthesis models and a large number of interstellar sight-lines towards the Magellanic Clouds.

More than 20 years ago, Glendenning, Kettner and Weber proposed the existence of stable white dwarfs with a core of strange quark matter. More recently, by studying radial modes, Alford, Harris and Sachdeva concluded instead that those objects are unstable. We investigate the stability of these objects by looking again at their radial oscillations, while incorporating boundary conditions at the quark-hadron interface which correspond either to a rapid or to a slow conversion of hadrons into quarks. Our analysis shows that objects of this type are stable if the star is not strongly perturbed and ordinary matter cannot transform into strange quark matter because of the Coulomb barrier separating the two components. On the other hand, ordinary matter can be transformed into strange quark matter if the star undergoes a violent process, as in the preliminary stages of a type Ia supernova, and this causes the system to become unstable and to collapse into a strange quark star. In this way, accretion induced collapse of strange dwarfs can be facilitated and km-sized objects with subsolar masses can be produced.

The high-velocity, neutral hydrogen feature known as MI may be the result of a supernova that took place about 100,000 years ago at a distance of 163 pc. Low-velocity HI data show a clear cavity, a structure indicative of regions evacuated by old exploding stars, centered on the spatial coordinates of MI, (l,b) = (165o, 65.o5). The invisible companion of the yellow giant star, 56 Ursae Majoris, may be the remains of the supernova that evacuated the cavity and blasted MI itself outward at 120 km/s. The mass and energy of MI are easily in line with what is expected from a supernova. The X-rays seen by ROSAT are consistent with an origin in the resulting bow shock. Ironically, this scenario for MI only came together because we were exploring low-velocity gas in the direction of high-velocity clouds.

Future and on-going infrared and radio observatories such as JWST, METIS or ALMA will increase the amount of rest-frame IR spectroscopic data for galaxies by several orders of magnitude. While studies of the chemical composition of the ISM based on optical observations have been widely spread over decades for SFG and, more recently, for AGN, similar studies need to be performed using IR data. This regime can be especially useful in the case of AGN given that it is less affected by temperature and dust extinction, traces higher ionic species and can also provide robust estimations of the chemical abundance ratio N/O. We present a new tool based on a bayesian-like methodology to estimate chemical abundances from IR emission lines in AGN. We use a sample of 58 AGN with IR spectroscopic data retrieved from the literature to probe the validity of our method. The estimations of the chemical abundances based on IR lines in our sample are later compared with the corresponding abundances derived from the optical emission lines in the same objects. HII-CHI-Mistry-IR takes advantage of photoionization models, characterized by the chemical abundance ratios O/H and N/O and the ionization parameter $U$, to compare their predicted emission-line fluxes with a set of observed values. Instead of matching single emission lines, the code uses some specific emission-line ratios sensitive to the above free parameters. We report mainly solar and also subsolar abundances for O/H in the nuclear region for our sample of AGN, whereas N/O clusters around solar values. We find a discrepancy between the chemical abundances derived from IR and optical emission lines, being the latter higher than the former. This discrepancy, also reported by previous studies of the composition of the ISM in AGN from IR observations, is independent from the gas density or the incident radiation field to the gas.

We perform an expansion study of the Balmer dominated outer shock of the SNR 0519$-$69.0 in the LMC by using a combination of new HST WFC3 imagery obtained in 2020 and archival ACS images from 2010 and 2011. Thanks to the very long time baseline, our proper motion measurements are of unprecedented accuracy. We find a wide range of shock velocities, with the fastest shocks averaging 5280 km/s and the slowest grouping of shocks averaging just 1670 km/s. We compare the H_alpha images from HST with X-ray images from Chandra and mid-IR images from Spitzer, finding a clear anti-correlation between the brightness of the remnant in a particular location and the velocity of the blast wave at that location, supporting the idea that the bright knots of X-ray and IR emission result from an interaction with a dense inhomogeneous circumstellar medium. We find no evidence for X-ray emission, thermal or nonthermal, associated with the fastest shocks, as expected if the fastest velocities are the result of the blast wave encountering the lower density ambient medium of the LMC. We derive an age of the remnant of 670 +/- 70 yr, consistent with results derived from previous investigations.

We present photometric properties of 183 member galaxies in the Abell 426 cluster using the Sloan Digital Sky Survey (SDSS) imaging and spectroscopic observation. Detailed morphology based on visual classification followed by multi-component image decomposition of 179 galaxies is presented in the SDSS g, r, i-bands. More than 80 percent of the members are Early-type galaxies (ETGs), with elliptical, dwarf elliptical (dE), and lenticular morphology and follow the red-sequence in the color-magnitude diagram (CMD). With a few dEs and spirals in the blue cloud, the cluster CMD is nearly unimodal.The dEs are $\sim2$-mag fainter and follow a different Sersic index and central velocity dispersion distribution than their bright counterparts. Further, we establish the Kormendy relation (KR) and the Fundamental Plane relation (FPR) for 5 different samples of ETGs constructed based on derived physical parameters such as Sersic index, concentration, central velocity dispersion in g, r, i-bands. The mean r-band slope and zero-point of the KR are $3.02\pm0.1$ and $18.65\pm0.03$ in close agreement to other cluster ellipticals in the local and higher redshift. Kinematics-based ETG sample produces the least scatter in KR with zero-point getting brighter by $\sim1.3$-mag from g to i band. The dEs and other low-mass ETGs follow the KR with a similar slope but with $\sim1.3$-mag fainter zero-point and form a parallel KR. The bright ellipticals follow an FPR with $a=1.37\pm0.003$, $b=0.35\pm0.05$ and $c=- 9.37\pm0.02$ in the r-band; galaxies tend to deviate from this relation at the low-mass end. A catalog with morphology and 2D structural analysis is available online.

Theories of planet formation give contradicting results of how frequent close-in giant planets of intermediate mass stars (IMSs; $\rm 1.3\leq M_{\star}\leq 3.2\,M_{\rm \odot}$) are. Some theories predict a high rate of IMSs with close-in gas giants, while others predict a very low rate. Thus, determining the frequency of close-in giant planets of IMSs is an important test for theories of planet formation. We use the CoRoT survey to determine the absolute frequency of IMSs that harbour at least one close-in giant planet and compare it to that of solar-like stars. The CoRoT transit survey is ideal for this purpose, because of its completeness for gas-giant planets with orbital periods of less than 10 days and its large sample of main-sequence IMSs. We present a high precision radial velocity follow-up programme and conclude on 17 promising transit candidates of IMSs, observed with CoRoT. We report the detection of CoRoT-34b, a brown dwarf close to the hydrogen burning limit, orbiting a 1.1 Gyr A-type main-sequence star. We also confirm two inflated giant planets, CoRoT-35b, part of a possible planetary system around a metal-poor star, and CoRoT-36b on a misaligned orbit. We find that $0.12 \pm 0.10\,\%$ of IMSs between $1.3\leq M_{\star}\leq 1.6 M_{\rm \odot}$ observed by CoRoT do harbour at least one close-in giant planet. This is significantly lower than the frequency ($0.70 \pm 0.16\,\%$) for solar-mass stars, as well as the frequency of IMSs harbouring long-period planets ($\rm \sim 8\,\%$).

Context: By providing information about the location of scattering material along the line of sight (LoS) to pulsars, scintillation arcs are a powerful tool for exploring the distribution of ionized material in the interstellar medium. Here, we present observations that probe the ionized ISM on scales of $\sim$~0.001 -- 30~au. Aims: We have surveyed pulsars for scintillation arcs in a relatively unbiased sample with DM < 100 pc cm-3. We present multi-frequency observations of 22 low to moderate DM pulsars. Many of the 54 observations were also observed at another frequency within a few days. Methods: For all observations we present dynamic spectra, autocorrelation functions, and secondary spectra. We analyze these data products to obtain scintillation bandwidths, pulse broadening times, and arc curvatures. Results: We detect definite or probable scintillation arcs in 19 of the 22 pulsars and 34 of the 54 observations, showing that scintillation arcs are a prevalent phenomenon. The arcs are better defined in low DM pulsars. We show that well-defined arcs do not directly imply anisotropy of scattering. Only the presence of reverse arclets and a deep valley along the delay axis, which occurs in about 20\% of the pulsars in the sample, indicates substantial anisotropy of scattering. Conclusions: The survey demonstrates substantial patchiness of the ionized ISM on both au size scales transverse to the line of sight and on $\sim$~100~pc scales along it. We see little evidence for distributed scattering along most lines of sight in the survey.

Polarization phase curves of asteroids and other small airless bodies are influenced by the compositional and physical properties of their regolith. The mixing of minerals composing the regolith influences the negative polarization at small phase angles because it changes the multiple scattering properties of the medium. This work aims to demonstrate experimentally how the mixing effect influences the polarization phase curve at small phase angles for different mineralogies relevant for asteroids, and to determine how different aggregate sizes affect the negative polarization. We prepared a set of binary and ternary mixtures with different common minerals on asteroids and one set of the same mixture with different aggregate sizes. We measured their reflected light at 530 nm with full Stokes polarimetry at phase angles ranging from 0.8{\deg} to 30{\deg}. The mixing effect of the mixtures with both bright and dark minerals significantly changes the behavior of the phase curves in terms of minimum polarization, phase angle of the minimum, and inversion angle with respect to the mineral components that are mixed together. The changes in phase curve could explain the polarization observation of particular classes of asteroids (F and L class) and other asteroids with peculiar polarization curves or photometric properties. Furthermore, we demonstrate that the negative polarization is invariant to the presence of dust aggregates up to centimeter sizes.

The formation of galaxies can be understood in terms of the assembly patterns of each type of galactic component. To perform this kind of analysis, is necessary to define some criteria to separate those components. Decomposition methods based on dynamical properties are more physically motivated than photometry-based ones. We use the unsupervised Gaussian Mixture model of \texttt{galactic structure finder} to extract the components of a sub-sample of galaxies with Milky Way-like masses from the EAGLE simulations. A clustering in the space of first and second order dynamical moments of all identified substructures reveals five types of galaxy components: thin and thick disks, stellar halos, bulges and spheroids. We analyse the dynamical, morphological and stellar population properties of these five component types, exploring to what extent these properties correlate with each other, and how much they depend on the total galaxy stellar and dark matter halo masses. All galaxies contain a bulge, a stellar halo and a disk. 60% of objects host two disks (thin and thick), and 68% host also a spheroid. The dynamical disk-to-total ratio does not depend on stellar mass, but the median rotational velocities of the two disks do. Thin disks are well separated in stellar ages, [Fe/H] and $\alpha$-enhancement from the three dispersion-dominated components, while thick disks are in between. Except for thin disks, all components show correlations among their stellar population properties: older ages mean lower metallicities and larger $\alpha$-enhancement. Finally, we quantify the weak dependence of stellar population properties on each component's dynamics.

We present a first-look analysis of the JWST ERO data in the SMACS J0723.3-7327 cluster field. We begin by reporting 10 new spectroscopic redshifts from $\lambda_\mathrm{obs}=1.8-5.2\mu$m NIRSpec medium-resolution ($R=\lambda/\Delta\lambda = 1000$) data. These are all determined via multiple high-SNR emission line detections, with 5 objects at $1 < z < 3$ displaying multiple rest-frame near-infrared Hydrogen Paschen lines, and 5 objects at $5 < z < 9$ displaying rest-frame optical Oxygen and Hydrogen Balmer lines. For the 5 higher-redshift galaxies we extract fluxes in 6 NIRCam bands spanning $\lambda_\mathrm{obs}=0.8-5\mu$m and perform spectral energy distribution fitting using these data in combination with existing HST photometry. The $7 < z < 9$ objects exhibit a U-shaped pattern across the F277W, F356W and F444W bands, indicating the presence of a Balmer break seen in emission (Balmer jump) and high-equivalent-width [O\,\textsc{iii}] emission. This is indicative of an extremely young stellar population, with the bulk of the current mass having formed within the past 10 Myr. We report robust stellar masses and mean stellar ages from our spectral fitting, with the four $z > 6$ galaxies exhibiting low stellar masses from log$_{10}(M_*/$M$_\odot)=7.4-8.6$ and correspondingly young mean stellar ages of only a few Myr. This work highlights the critical importance of combining large upcoming NIRCam surveys with NIRSpec follow-up to measure the spectroscopic redshifts necessary to robustly constrain physical parameters.

Large fractions of metals are missing from the observable gas-phase in the interstellar medium (ISM) because they are incorporated into dust grains, a phenomenon called dust depletion. The study of dust depletion in the ISM is important to investigate the origin and evolution of metals and cosmic dust. Here we aim at characterizing the dust depletion of several metals from the Milky Way to distant galaxies. We collect ISM metal column densities from absorption-line spectroscopy in the literature, and in addition, we determine Ti and Ni column densities from a sample of 70 damped Lyman-$\alpha$ absorbers (DLAs) towards quasars, observed with UVES/VLT. We use ISM relative abundances to estimate the dust depletion of 18 metals (C, P, O, Cl, Kr, S, Ge, Mg, Si, Cu, Co, Mn, Cr, Ni, Al, Ti, Zn and Fe) for different environments (the Milky Way, the Magellanic Clouds (MCs), DLAs towards quasars and towards gamma-ray bursts). We observe linear relations between the depletion of each metal and the strength of dust depletion, which we trace with the observed [Zn/Fe]. In the neutral ISM of the MCs we find small deviations from linearity observed as an overabundance of the $\alpha$-elements Ti, Mg, S and an underabundance of Mn. The deviations disappear if we assume that all OB stars observed towards the MCs in our sample have an $\alpha$-element enhancement and Mn underabundance. This may imply that the MCs have been recently enriched in $\alpha$-elements, potentially due to recent bursts of star formation. The observed strong correlations of the depletion sequences of the metals all the way from low metallicity QSO-DLAs to the Milky Way suggest that cosmic dust has a common origin, independently of the star formation history, which varies significantly between these different galaxies. This supports the importance of grain growth in the ISM as a significant process of dust production.

We treat prospects for multimessenger astronomy with giant flares (GFs), a rare transient event featured by magnetars that can be as luminous as a hundred of the brightest supernovae ever observed. The beamed photons could correlate with an axion counterpart via resonant conversion in the magnetosphere. In a realistic parameter space, we find that the sensitivity limit to galactic GFs for currently viable experiments is $\mathrm{g}_{\phi \gamma}\!\gtrsim\!\mathrm{several}\!\times\!10^{-13}$ GeV$^{-1}$ \& $\mathrm{g}_{\phi e}\!\gtrsim\!\mathrm{few}\!\times\!10^{-12}$. We rule out the compatibility of axion flares with the recent XENON1T excess only due to the time persistence of the signal.

Superradiance in black holes is well-understood but a general treatment for superradiance in stars has until now been lacking. This is surprising given the ease with which we can observe isolated neutron stars and the array of signatures which would result from stellar superradiance. In this work, we present the first systematic pipeline for computing superradiance rates in rotating stars. Our method can be used with any Lagrangian describing the interaction between the superradiant field and the constituents of the star. Our scheme falls into two parts: firstly we show how field theory at finite density can be used to express the absorption of long wavelength modes into the star in terms of microphsyical scattering processes. This allows us to derive a damped equation of motion for the bosonic field. We then feed this into an effective theory for long wavelengths (the so-called worldline formalism) to describe the amplification of superradiant modes of arbitrary multipole moment for a rapidly rotating star. Our method places stellar superradiance on a firm theoretical footing and allows the calculation of the superradiance rate arising from any interaction between a bosonic field and stellar matter.

We consider the production of axion dark matter through the misalignment mechanism in the context of a nonstandard cosmological history involving early matter domination by a scalar field with a time-dependent decay rate. In cases where the temperature of the Universe experiences a temporary period of increase, Hubble friction can be restored in the evolution of the axion field, resulting in the possibility of up to three "crossings" of the axion mass and the Hubble expansion rate. This has the effect of dynamically resetting the misalignment mechanism to a new initial state for a second distinct phase of oscillation. The resultant axion mass required for the present dark matter relic density is never bigger than the standard-history window and can be smaller by more than three orders of magnitude, which can be probed by upcoming experiments such as ABRACADABRA, KLASH, ADMX, MADMAX, and ORGAN, targeting the axion-photon coupling. This highlights the possibility of exploring the cosmological history prior to Big Bang Nucleosynthesis through searches for axion dark matter beyond the standard window.

Extreme mass ratio inspirals (EMRIs) -- systems with a compact object orbiting a much more massive (e.g., galactic center) black hole -- are of interest both as a new probe of the environments of galactic nuclei, and their waveforms are a precision test of the Kerr metric. This work focuses on the effects of an external perturbation due to a third body around an EMRI system. This perturbation will affect the orbit most significantly when the inner body crosses a resonance with the outer body, and result in a change of the conserved quantities (energy, angular momentum, and Carter constant) or equivalently of the actions, which results in a subsequent phase shift of the waveform that builds up over time. We present a general method for calculating the changes in action during a resonance crossing, valid for generic orbits in the Kerr spacetime. We show that these changes are related to the gravitational waveforms emitted by the two bodies (quantified by the amplitudes of the Weyl scalar $\psi_4$ at the horizon and at $\infty$) at the frequency corresponding to the resonance. This allows us to compute changes in the action variables for each body, without directly computing the explicit metric perturbations, and therefore we can carry out the computation by calling an existing black hole perturbation theory code. We show that our calculation can probe resonant interactions in both the static and dynamical limit. We plan to use this technique for future investigations of third-body effects in EMRIs and their potential impact on waveforms for LISA.

Exotic objects such as the Ellis wormhole are expected to act as gravitational lenses. Much like their non-exotic counterparts, information about these lenses can be found by considering the strong and weak lensing fields they induce. In this work, we consider how weak gravitational lensing flexion can provide information beyond that of shear. We find that, in the case of an Ellis wormhole-type metric, directional flexion can distinguish between the case of a positive or negative convergence, while directional shear cannot. We also consider cosmic flexion, the flexion correlation function whose signal originates from the large-scale structure of the Universe, in the context of modified gravity. We find flexion to be a unique probe of parametric models of modified gravity, particularly in the case of scale-dependent phenomenological post-GR functions.

In this paper, we investigate the primordial perturbations of inflation model induced from the multi-field mimetic gravity, where there are two field during inflation, and thus both adiabatic and isocurvature perturbation modes are generated. We show that although it is true that the original adiabatic perturbation mode loses the kinetic term due to the constraint equation, by applying the curvaton mechanism where one of the field is viewed as curvaton field, the adiabatic perturbation can actually be transferred from the isocurvature one at the end of inflation. Detailed calculations are performed for both inflationary and the consequent matter-dominant epochs. Therefore, the so-called "non-propagating problem" of the adiabatic mode will actually do no harm to the multi-field mimetic inflation models.

We investigate the modifications in the neutrino flavor oscillations under the influence of a stochastic gravitational wave background (SGWB), in a scenario in which General Relativity is modified by an additional Chern-Simons (CS) term. Assuming that the dark matter halo is in the form of axions, the CS coupling modifies the pattern of the neutrino flavor oscillations at Earth up to a total suppression in some frequency range. At the same time, the SGWB in the halo could stimulate the axion decay into gravitons over a narrow frequency range, leading to a potentially detectable resonance peak in the enhanced SGWB strain. A consistent picture would require these features to potentially show up in neutrino detection from supernovae, gravitational wave detectors, and experiments aimed at the search for axions in the Milky Way halo.

Modified Newtonian equations for gravitational orbits in the expanding universe indicate that local gravitationally bounded systems like galaxies and planetary systems are unaffected by the expansion of the Universe. This result is derived under the assumption of the space expansion described by the standard FLRW metric. In this paper, an alternative metric is applied and the modified Newtonian equations are derived for the space expansion described by the conformal FLRW metric. As shown by Vavry\v{c}uk (Frontiers in Physics, 2022), this metric is advantageous, because it properly predicts the cosmic time dilation and fits the SNe Ia luminosity observations with no need to introduce dark energy. Surprisingly, the Newtonian equations based on the conformal FLRW metric behave quite differently than those based on the standard FLRW metric. In contrast to the common opinion that local systems resist the space expansion, the results for the conformal metric indicate that all local systems expand according to the Hubble flow. The evolution of the local systems with cosmic time is exemplified on numerical modelling of spiral galaxies. The size of the spiral galaxies grows consistently with observations and a typical spiral pattern is well reproduced. The theory predicts flat rotation curves without an assumption of dark matter surrounding the galaxy. The theory resolves challenges to the $\Lambda$CDM model such as the problem of faint satellite galaxies, baryonic Tully-Fisher relation or the radial acceleration relation. Furthermore, puzzles in the solar system are successfully explained such as the Pioneer anomaly, the Faint young Sun paradox, the Moon's and Titan's orbit anomalies or the presence of rivers on ancient Mars.

Many alternative theories of gravity screens a Yukawa-type potential. This article shows Keplerian-type parametrization as a solution of Yukawa type potential accurate equations of motion for two non-spinning compact objects moving in an eccentric orbit. A bound from the solar system is presented.

We present new bounds on the cosmic abundance of magnetic monopoles based on the survival of primordial magnetic fields during the reheating and radiation-dominated epochs. The new bounds can be stronger than the conventional Parker bound from galactic magnetic fields, as well as bounds from direct searches. We also apply our bounds to monopoles produced by the primordial magnetic fields themselves through the Schwinger effect, and derive additional conditions for the survival of the primordial fields.

The problem of finding a vacuum definition for a single quantum field in curved space-times is discussed under a new geometrical perspective. The phase space dynamics of the quantum field modes are mapped to curves in a 2-dimensional hyperbolic metric space, in which distances between neighbor points are shown to be proportional to the Bogoliubov coefficients associated with their corresponding mode solutions in phase space. The vacuum state for each mode is then defined as the unique trajectory from which all mapped phase space solutions move within thin annular regions around it. This property implies the stability of the vacuum state: solutions evolved from a point in this trajectory stay close to it as both evolve, and the particle creation is therefore minimized. The new approach is applied to the well-known cases of the time-independent dynamics, where the solutions draw circles around this curve, and when the adiabatic approximation is valid. The analysis is then extended to time-dependent cases in which the adiabatic approximation is not applicable, in the super-Hubble or low-frequency regimes. It is shown that stability trajectories can also be found in these situations, and stable quantum vacua can be obtained. This new formalism is applied to two situations: de Sitter space, where the Bunch-Davies vacuum is obtained in a completely different manner through an analysis in the super-Hubble regime, and in the context of cosmological bouncing models, in which the contracting phase is dominated by a cosmological constant in the asymptotic past. A new vacuum state for cosmological perturbations is proposed in this situation.

Cosmology is most typically analyzed using standard co-moving coordinates, in which the galaxies are (on average, up to presumably small peculiar velocities) "at rest", while "space" is expanding. But this is merely a specific coordinate choice; and it is important to realise that for certain purposes other, (sometimes radically different) coordinate choices might also prove useful and informative, but without changing the underlying physics. Specifically, herein we shall consider the k=0 spatially flat FLRW cosmology but in Painleve-Gullstrand coordinates -- these coordinates are very explicitly not co-moving: "space" is now no longer expanding, although the distance between galaxies is still certainly increasing. Working in these Painleve-Gullstrand coordinates provides an alternate viewpoint on standard cosmology, and the symmetries thereof, and also makes it somewhat easier to handle cosmological horizons. With a longer view, we hope that investigating these Painleve-Gullstrand coordinates might eventually provide a better framework for understanding large deviations from idealized FLRW spacetimes. We illustrate these issues with a careful look at the Kottler and McVittie spacetimes.

In the context of the modified teleparallel $f(T, B)$-theory of gravity, we consider a homogeneous and anisotropic background geometry described by the Kantowski-Sachs line element. We derive the field equations and investigate the existence of exact solutions. Furthermore, the evolution of the trajectories for the field equations is studied by deriving the stationary points at the finite and infinite regimes. For the $f(T,B)=T+F\left( B\right) $ theory, we prove that for a specific limit of the function $F\left( B\right) $, the anisotropic Universe has the expanding and isotropic Universe as an attractor with zero spatial curvature. We remark that there are no future attractors where the asymptotic solution describes a Universe with nonzero spatial curvature.

We study the asymptotic dynamics of $f(T, B)$-theory in an anisotropic Bianchi III background geometry. We show that an attractor always exists for the field equations, which depends on a free parameter provided by the specific $f(T, B)$ functional form. The attractor is an accelerated spatially flat FLRW or non-accelerated LRS Bianchi III geometry. Consequently, the $f(T, B)$-theory provides a spatially flat and isotropic accelerated Universe.

The Noether symmetry analysis is applied for the analysis of the field equations in an anisotropic background in $f(T,B)$-theory. We consider the $f\left( T,B\right) =T+F\left( B\right) $ which describes a small deviation from TEGR introduced by the boundary scalar $B$. For the Bianchi\ I, Bianchi III and Kantowski-Sachs geometries there exists a minisuperspace description and Noether's theorems are applied. We investigate the existence of invariant point transformations. We find that for the Bianchi I spacetime the gravitational field equations are Liouville integrable for the $F\left( B\right) =-\frac{B}{\lambda}\ln B$ theory. The analytic solution is derived and the application of Noether symmetries to the Wheeler-DeWitt equation of quantum cosmology is discussed.

In the standard model of cosmology, the background evolution of the Universe can in general be adequately described by general relativity and a uniform and isotropic metric minimally coupled with a collection of perfect fluids. These fluids are usually described by their energy-momentum tensor, which can be derived from the fluid's Lagrangian density. Under general relativity, the Lagrangian density is only relevant to the extent that it results in the correct energy-momentum tensor for a specific perfect fluid. This is not the case in theories that feature a nonminimal coupling (NMC) between the matter fields and gravity. In such cases, the on-shell Lagrangian density of the matter fields appears explicitly in the equations of motion, in addition to their energy-momentum tensor. The determination of the correct on-shell Lagrangian density for a particular fluid is therefore of paramount importance in order to provide an accurate description of the corresponding cosmological implications. In essence, this is the problem tackled in this thesis. We have aimed at addressing three key points. We covered some of the results in the literature regarding the Lagrangian density of cosmic fluids, and cleared up some misunderstandings regarding the freedom of choice (or lack thereof) of its on-shell form, both in general relativity and in theories featuring an NMC. In addition, we derived the correct Lagrangian density for fluids composed of solitonic particles with fixed rest mass and structure. Secondly, we studied the thermodynamic behaviour of perfect fluids of this type in the context of theories featuring an NMC between gravity and the matter fields. Finally, we used these results to derive novel cosmological constraints on specific NMC gravity models, using data from cosmic microwave background, big-bang nucleosynthesis, type Ia supernovae and baryon acoustic oscillations observations.

A recently proposed extension of the geodesic equations of motion, where the worldline traced by a test particle now depends on the scalar curvature, is used to study the formation of galaxies and galactic rotation curves. This extension is applied to the motion of a fluid in a spherical geometry, resulting in a set of evolution equations for the fluid in the nonrelativistic and weak gravity limits. Focusing on the stationary solutions of these equations and choosing a specific class of angular momenta for the fluid in this limit, we show that dynamics under this extension can result in the formation of galaxies with rotational velocity curves (RVC) that are consistent with the Universal Rotation Curve (URC), and through previous work on the URC, the observed rotational velocity profiles of 1100 spiral galaxies. In particular, a spectrum of RVCs can form under this extension, and we find that the two extreme velocity curves predicted by it brackets the ensemble of the URCs constructed from these 1100 velocity profiles. We also find that the asymptotic behavior of the URC is consistent with that of the most probable asymptotic behavior of the RVCs predicted by the extension. A stability analysis of these stationary solutions is also done, and we find them to be stable in the galactic disk, while in the galactic hub they are stable if the period of oscillations of perturbations is longer than $0.91_{\pm0.31}$ to $1.58_{\pm 0.46}$ billion years.

In this work we analize how the inclusion of extra mesonic degrees of freedom affect the finite density solitons crystals of the Skyrme model. In particular, the first analytic examples of hadronic crystals at finite baryon density in both the Skyrme $\omega$-mesons model as well as for the Skyrme $\rho$-mesons theory are constructed. These configurations have arbitrary topological charge and describe crystals of baryonic tubes surrounded by a cloud of vector-mesons. In the $\omega$-mesons case, it is possible to reduce consistently the complete set of seven coupled non-linear field equations to just two integrable differential equations; one ODE for the Skyrmion profile and one PDE for the $\omega$-mesons field. This analytical construction allows to show explicitly how the inclusion of $\omega$-mesons in the Skyrme model reduces the repulsive interaction energy between baryons. In the Skyrme $\rho$-mesons case, it is possible to construct analytical solutions using a meron-type ansatz and fixing one of the couplings of the $\rho$-mesons action in terms of the others. We show that, quite remarkably, the values obtained for the coupling constants by requiring the consistency of our ansatz are very close to the values used in the literature to reduce nuclei binding energies of the Skyrme model without vector-mesons. Moreover, our analytical results are in qualitative agreement with the available results on the nuclear spaghetti phase.