Papers for Thursday, Mar 09 2023

We present the first statistical study of spatially integrated non-Gaussian stellar kinematics spanning 7 Gyr in cosmic time. We use deep, rest-frame optical spectroscopy of massive galaxies (stellar mass $M_\star > 10^{10.5} {\rm M}_\odot$) at redshifts z = 0.05, 0.3 and 0.8 from the SAMI, MAGPI and LEGA-C surveys, to measure the excess kurtosis $h_4$ of the stellar velocity distribution, the latter parametrised as a Gauss-Hermite series. We find that at all redshifts where we have large enough samples, $h_4$ anti-correlates with the ratio between rotation and dispersion, highlighting the physical connection between these two kinematic observables. In addition, and independently from the anti-correlation with rotation-to-dispersion ratio, we also find a correlation between $h_4$ and $M_\star$, potentially connected to the assembly history of galaxies. In contrast, after controlling for mass, we find no evidence of independent correlation between $h_4$ and aperture velocity dispersion or galaxy size. These results hold for both star-forming and quiescent galaxies. For quiescent galaxies, $h_4$ also correlates with projected shape, even after controlling for the rotation-to-dispersion ratio. At any given redshift, star-forming galaxies have lower $h_4$ compared to quiescent galaxies, highlighting the link between kinematic structure and star-forming activity.

A number of neutron stars have been observed within the remnants of the core-collapse supernova explosions that created them. In contrast, black holes are not yet clearly associated with supernova remnants. Indeed, some observations suggest that black holes are ``born in the dark'', i.e. without a supernova explosion. Herein, we present a multi-wavelength analysis of the X-ray transient Swift J1728.9$-$3613, based on observations made with Chandra, ESO-VISTA, MeerKAT, NICER, NuSTAR, Swift, and XMM-Newton. Three independent diagnostics indicate that the system likely harbors a black hole primary. Infrared imaging signals a massive companion star that is broadly consistent with an A or B spectral type. Most importantly, the X-ray binary lies within the central region of the catalogued supernova remnant G351.9$-$0.9. Our deep MeerKAT image at 1.28~GHz signals that the remnant is in the Sedov phase; this fact and the non-detection of the soft X-ray emission expected from such a remnant argue that it lies at a distance that could coincide with the black hole. Utilizing a formal measurement of the distance to Swift J1728.9$-$3613 ($d = 8.4\pm 0.8$ kpc), a lower limit on the distance to G351.9$-$0.9 ($d \geq 7.5$ kpc), and the number and distribution of black holes and supernova remnants within the Milky Way, extensive simulations suggest that the probability of a chance superposition is $<1.7\%$ ($99.7\%$ credible interval). The discovery of a black hole within a supernova remnant would support numerical simulations that produce black holes and remnants, and thus provide clear observational evidence of distinct black hole formation channels. We discuss the robustness of our analysis and some challenges to this interpretation.

Most dynamical models of galaxies to date assume axisymmetry, which is not representative of a significant fraction of massive galaxies. We have built triaxial orbit-superposition Schwarzschild models of galaxies observed by the SAMI Galaxy Survey, in order to reconstruct their inner orbital structure and mass distribution. The sample consists of 153 passive galaxies with total stellar masses in the range $10^{9.5}$ to $10^{12} M_{\odot}$. We present an analysis of the internal structures and intrinsic properties of these galaxies as a function of their environment. We measure their environment using three proxies: central or satellite designation, halo mass and local $5^{th}$ nearest neighbour galaxy density. We find that although these intrinsic properties correlate most strongly with stellar mass, environment does play a secondary role: at fixed stellar mass, galaxies in the densest regions are more radially anisotropic. In addition, central galaxies, and galaxies in high local densities show lower values of edge-on spin parameter proxy \lam. We also find suggestions of a possible trend of the fractions of orbits with environment for lower-mass galaxies (between $10^{9.5}$ and $10^{11} M_{\odot}$) such that, at fixed stellar mass, galaxies in higher local densities and halo mass have higher fractions of hot orbits and lower fractions of warm orbits. Our results demonstrate that after stellar mass, environment does play a role in shaping present-day passive galaxies.

We study the initial development, structure and evolution of protoplanetary clumps formed in 3D resistive MHD simulations of self-gravitating disks. The magnetic field grows by means of the recently identified gravitational instability dynamo (Riols & Latter 2018; Deng et al. 2020). Clumps are identified and their evolution is tracked finely both backward and forward in time. Their properties and evolutionary path is compared to clumps in companion simulations without magnetic fields. We find that magnetic and rotational energy are important in the clumps' outer regions, while in the cores, despite appreciable magnetic field amplification, thermal pressure is most important in counteracting gravity. Turbulent kinetic energy is of a smaller scale than magnetic energy in the clumps. Compared to non-magnetized clumps, rotation is less prominent, which results in lower angular momentum in much better agreement with observations. In order to understand the very low sub-Jovian masses of clumps forming in MHD simulations, we revisit the perturbation theory of magnetized sheets finding support for a previously proposed magnetic destabilization in low-shear regions. This can help explaining why fragmentation ensues on a scale more than an order of magnitude smaller than that of the Toomre mass. The smaller fragmentation scale and the high magnetic pressure in clumps' envelopes explain why clumps in magnetized disks are typically in the super-Earth to Neptune mass regime rather than Super-Jupiters as in conventional disk instability. Our findings put forward a viable alternative to core accretion to explain widespread formation of intermediate-mass planets.

The origin and distribution of stellar-mass black hole spins are a rare window into the progenitor stars and supernova events that create them. Swift J1728.9-3613 is an X-ray binary, likely associated with the supernova remnant G351.9-0.9 (Balakrishnan et al. 2023). A NuSTAR X-ray spectrum of this source during its 2019 outburst reveals reflection from an accretion disk extending to the innermost stable circular orbit. Modeling of the relativistic Doppler shifts and gravitational redshifts imprinted on the spectrum measures a dimensionless spin parameter of $a=0.86\pm0.02$ ($1\sigma$ confidence), a small inclination angle of the inner accretion disk $\theta<10$ degrees, and a sub-solar iron abundance in the disk $A_{\rm Fe}<0.84$. This high spin value rules out a neutron star primary at the $5\;\sigma$ level of confidence. If the black hole is located in a still visible supernova remnant, it must be young. Therefore, we place a lower limit on the natal black hole spin of $a>0.82$, concluding that the black hole must have formed with a high spin. This demonstrates that black hole formation channels that leave a supernova remnant, and those that do not (e.g. Cyg X-1), can both lead to high natal spin with no requirement for subsequent accretion within the binary system. Emerging disparities between the population of high-spin black holes in X-ray binaries and the low-spin black holes that merge in gravitational wave events may therefore be explained in terms of different stellar conditions prior to collapse, rather than different environmental factors after formation.

We study the recent star formation histories (SFHs) of 575 intermediate-mass galaxies (IMGs, $10^{9} \leq M/M_{\odot} \leq 10^{10}$) in COSMOS at $0.3<z<0.4$ by comparing their H$\alpha$ and UV luminosities. These two measurements trace star formation rates (SFRs) on different timescales and together reveal fluctuations in recent activity. We compute $L_{{\rm H}\alpha}$ from Magellan IMACS spectroscopy while $L_{\rm UV}$ is derived from rest-frame 2800 $\text{\r{A}}$ photometry. Dust corrections are applied to each band independently. We compare the deviation of $L_{{\rm H}\alpha}$ and $L_{\rm UV}$ from their respective star forming sequences (i.e., $\Delta\log L_{{\rm H}\alpha}$ and $\Delta\log L_{\rm UV}$) and after accounting for observational uncertainties we find a small intrinsic scatter between the two quantities ($\sigma_{\delta} \lesssim 0.03$ dex). This crucial observational constraint precludes strong fluctuations in the recent SFHs of IMGs: simple linear SFH models indicate that a population of IMGs would be limited to only factors of $\lesssim 2$ change in SFR over $200$ Myr and $\lesssim 30\%$ on shorter timescales of $20$ Myr. No single characteristic SFH for IMGs, such as an exponentially rising/falling burst, can reproduce the individual and joint distribution of $\Delta\log L_{{\rm H}\alpha}$ and $\Delta\log L_{\rm UV}$. Instead, an ensemble of SFHs is preferred. Finally, we find that IMG SFHs predicted by recent hydrodynamic simulations, in which feedback drives rapid and strong SFR fluctuations, are inconsistent with our observations.

We consider the prospects for future ultrahigh energy cosmic ray and neutrino observations to constrain the evolution of sources producing a proton flux above 10 EeV (1 EeV = 10^18 eV). We find that strong constraints on the source evolution can be obtained by combining measurements of the cosmic ray proton fraction above 30 EeV with measurement of the neutrino flux at 1 EeV, if neutrinos are predominantly of cosmogenic origin. In the case that interactions in the source environment produce a significant astrophysical neutrino flux, constraints on the source evolution may require measurement of the observed proton fraction, as well as, the neutrino flux at multiple energies, such as 1 EeV and 10 EeV. Finally, we show that fits to current UHECR data favor models which result in a >30 EeV proton fraction and 1 EeV neutrino flux that could realistically be discovered by the next generation of experiments.

In local disk galaxies such as our Milky Way, older stars generally inhabit a thicker disk than their younger counterparts. Two competing models have attempted to explain this result: one in which stars first form in thin disks that gradually thicken with time through dynamical heating, and one in which stars form in thick disks at early times and in progressively thinner disks at later times. We use a direct measure of the thicknesses of stellar disks at high redshift to discriminate between these scenarios. Using legacy HST imaging from the CANDELS and GOODS surveys, we measure the rest-optical scale heights of 491 edge-on disk galaxies spanning 0.4 < z < 2.5. We measure a median intrinsic scale height for the full sample of 0.74 +/- 0.03 kpc, with little redshift evolution of both the population median and scatter. The median is consistent with the thick disk of the Milky Way today (0.6 - 1.1 kpc), but is smaller than the median scale height of local disks (~1.5 kpc) which are matched to our high-redshift sample by descendant mass. These findings indicate that (1) while disks as thick as the Milky Way's thick disk were in place at early times, (2) to explain the full disk galaxy population today, the stellar disks in galaxies need to on average physically thicken after formation.

The formation of molecular gas in interstellar clouds is a slow process, but is enhanced by gas compression. Magnetohydrodynamic (MHD) waves create compressed quasiperiodic linear structures, referred to as striations. Striations are observed at column densities where the atomic to molecular gas transition takes place. We explore the role of MHD waves in the CO chemistry in regions with striations within molecular clouds. We target a region with striations in the Polaris Flare cloud. We conduct a CO J=2-1 survey in order to probe the molecular gas properties. We use archival starlight polarization data and dust emission maps in order to probe the magnetic field properties and compare against the CO properties. We assess the interaction of compressible MHD wave modes with CO chemistry by comparing their characteristic timescales. The estimated magnetic field is 38 - 76 $\mu$G. In the CO integrated intensity map, we observe a dominant quasi-periodic intensity structure, which tends to be parallel to the magnetic field orientation and has a wavelength of one parsec approximately. The periodicity axis is $\sim$ 17 degrees off from the mean magnetic field orientation and is also observed in the dust intensity map. The contrast in the CO integrated intensity map is $\sim 2.4$ times larger than the contrast of the column density map, indicating that CO formation is enhanced locally. We suggest that a dominant slow magnetosonic mode with estimated period $2.1 - 3.4$ Myr, and propagation speed $0.30 - 0.45$ km~s$^{-1}$, is likely to have enhanced the formation of CO, hence created the observed periodic pattern. We also suggest that, within uncertainties, a fast magnetosonic mode with period 0.48 Myr and velocity $2.0$ km~s$^{-1}$ could have played some role in increasing the CO abundance. Quasiperiodic CO structures observed in striation regions may be the imprint of MHD wave modes.

Carbonaceous chondrites are some of the most primitive meteorites and derive from planetesimals that formed a few million years after the beginning of the solar system. Here, using new and previously published Cr, Ti, and Te isotopic data, we show that carbonaceous chondrites exhibit correlated isotopic variations that can be accounted for by mixing among three major constituents having distinct isotopic compositions, namely refractory inclusions, chondrules, and CI chondrite-like matrix. The abundances of refractory inclusions and chondrules are coupled and systematically decrease with increasing amount of matrix. We propose that these correlated abundance variations reflect trapping of chondrule precursors, including refractory inclusions, in a pressure maximum in the disk, which is likely related to the water ice line and the ultimate formation location of Jupiter. The variable abundance of refractory inclusions/chondrules relative to matrix is the result of their distinct aerodynamical properties resulting in differential delivery rates and their preferential incorporation into chondrite parent bodies during the streaming instability, consistent with the early formation of matrix-poor and the later accretion of matrix-rich carbonaceous chondrites. Our results suggest that chondrules formed locally from isotopically heterogeneous dust aggregates which themselves derive from a wide area of the disk, implying that dust enrichment in a pressure trap was an important step to facilitate the accretion of carbonaceous chondrite parent bodies or, more generally, planetesimals in the outer solar system.

Using the ${\rm F{\small LARES}}$: First Light And Reionisation Epoch Simulations suite we explore the consequences of a realistic star-dust model on the observed properties of galaxies. We find that the attenuation in the UV declines rapidly from the galactic centre, and more luminous galaxies have extended star formation that suffers less obscuration than their fainter counterparts. This gives rise to a non-linear relationship between the observed UV luminosity and the UV attenuation, that produces a double power-law shape to the UV luminosity function of the ${\rm F{\small LARES}}$ galaxies. Spatially distinct stellar populations within galaxies experience a wide range of dust attenuation due to variations in the dust optical depth along the line of sight, ranging from fully obscured to unobscured. The overall attenuation curve of the whole galaxy is then a complex combination of those individual lines of sight. We explore the manifestation of this effect to study the reliability of line ratios, in particular the Balmer decrement and the BPT diagram. We find the Balmer decrement predicted Balmer line attenuation to be very different from those expected from commonly used attenuation curves from literature, and the observed BPT line ratios shifted from their intrinsic dust-free values. Finally, we explore the variation in observed properties with viewing angle, finding average differences of $\sim0.3$ magnitudes in the UV attenuation.

The first JWST spectroscopy of the luminous galaxy GN-z11 simultaneously both established its redshift at $z=10.6$ and revealed a rest-ultraviolet spectrum dominated by signatures of highly-ionized nitrogen, which has so far defied clear interpretation. Here we present a reappraisal of this spectrum in the context of both detailed nebular modeling and nearby metal-poor reference galaxies. The N IV] emission enables the first nebular density measurement in a star-forming galaxy at $z>10$, and reveals evidence for extremely high densities $n_e\gtrsim 10^5$ $\mathrm{cm^{-3}}$. We definitively establish with a suite of photoionization models that regardless of ionization mechanism and accounting for depletion and this density enhancement, an ISM substantially enriched in nitrogen ($[\mathrm{N/O}]=+0.52$) is required to reproduce the observed lines. A search of local UV databases confirms that nearby metal-poor galaxies power N IV] emission, but that this emission is uniformly associated with lower densities than implied in GN-z11. We compare to a unique nearby galaxy, Mrk~996, where a high concentration of Wolf-Rayet stars and their CNO-processed wind ejecta produce a UV spectrum remarkably similar to that of both GN-z11 and the Lyc-leaking super star cluster in the Sunburst Arc. Collating this evidence in the context of Galactic stellar abundances, we suggest that the peculiar nitrogenic features prominent in GN-z11 may be a unique signature of intense and densely clustered star formation in the evolutionary chain of the present-day globular clusters, consistent with in-situ early enrichment with nuclear-processed stellar ejecta on a massive scale. Combined with insight from local galaxies, these and future JWST data open a powerful new window onto the physical conditions of star formation and chemical enrichment at the highest redshifts.

We report on the use of a kinetic-inductance traveling-wave parametric amplifier (KITWPA) as the first amplifier in the readout chain of a microwave superconducting quantum interference device (SQUID) multiplexer (umux). This umux is designed to multiplex signals from arrays of low temperature detectors such as superconducting transition-edge sensor microcalorimeters. When modulated with a periodic flux-ramp to linearize the SQUID response, the flux noise improves, on average, from $1.6$ $\mu\Phi_0/\sqrt{\mathrm{Hz}}$ with the KITWPA off, to $0.77$ $\mu\Phi_0/\sqrt{\mathrm{Hz}}$ with the KITWPA on. When statically biasing the umux to the maximally flux-sensitive point, the flux noise drops from $0.45$ $\mu\Phi_0/\sqrt{\mathrm{Hz}}$ to $0.2$ $\mu\Phi_0/\sqrt{\mathrm{Hz}}$. We validate this new readout scheme by coupling a transition-edge sensor microcalorimeter to the umux and detecting background radiation. The combination of umux and KITWPA provides a variety of new capabilities including improved detector sensitivity and more efficient bandwidth utilization.

We report a peculiar variable blue star in the globular cluster NGC 6397, using Hubble Space Telescope optical imaging. Its position in the colour-magnitude diagrams, and its spectrum, are consistent with this star being a helium core white dwarf (He WD) in a binary system. The optical light curve shows a periodicity at 18.5 hours. We argue that this periodicity is due to the rotation of the WD and possibly due to magnetic spots on the surface of the WD. This would make this object the first candidate magnetic He WD in any globular cluster (GC), and the first candidate magnetic WD in a detached binary system in any GC and one of the few He WDs with a known rotation period and of magnetic nature. Another possibility is that this system is a He WD in a binary system with another WD or another degenerate object, which would make this object one of the few candidate non-accreting double degenerate binaries in any GC.

With the start of a new Great Observatories era, there is renewed concern that the demand for these forefront facilities, through proposal pressure, will exceed conventional peer-review management's capacity for ensuring an unbiased and efficient selection. There is need for new methods, strategies, and tools to facilitate those reviews. Here, we describe PACMan2, an updated tool for proposal review management that utilizes machine learning models and techniques to topically categorize proposals and reviewers, to match proposals to reviewers, and to facilitate proposal assignments, mitigating some conflicts of interest. We find that the classifier has cross-validation accuracy of $80.0\pm2.2\%$ on proposals for time on the Hubble Space Telescope and the James Webb Space Telescope.

The nebular He II {\lambda}1640 emission line is observed in star forming galaxies out to large distances and can be used to constrain the properties of sources of He+ ionizing photons. For this purpose, it is crucial to understand which are the main stellar sources of these photons. In some nearby metal poor starburst galaxies, nebular He II {\lambda}4686 (optical equivalent) is accompanied by a broad underlying component, which is generally attributed to formation in the winds of classical (He burning) Wolf Rayet stars, primarily of the WN subtype. In such cases, the origin of the nebular component has been proposed to be the escape of He+ ionizing photons from the winds of the WN stars, at least partially. We use archival long slit observations obtained with Focal Reducer Low Dispersion Spectrograph (FORS1) on the Very Large Telescope to look for nebular He II {\lambda}4686 emission south of the WN6h + WN6(7) close binary in HD 5980. We only find broad He II {\lambda}4686 emission, as far as aprox. 7.6 pc from the binary. A comparison with observations obtained with Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope, at a similar orbital phase, shows that the FORS1 broad He II emission is likely contamination from the multiple star system HD 5980. We use models to show that no significant He+ ionizing flux is expected from the WN stars in HD 5980 and that when similar stars are present in a coeval stellar population, the O stars can be far greater emitters of He+ ionizing radiation.

V906 Carinae was one of the best observed novae of recent times. It was a prolific dust producer and harboured shocks in the early evolving ejecta outflow. Here, we take a close look at the consequences of these early interactions through study of high-resolution UVES spectroscopy of the nebular stage and extrapolate backwards to investigate how the final structure may have formed. A study of ejecta geometry and shaping history of the structure of the shell is undertaken following a spectral line SHAPE model fit. A search for spectral tracers of shocks in the nova ejecta is undertaken and an analysis of the ionised environment. Temperature, density and abundance analyses of the evolving nova shell are presented.

As some of the most compact stellar objects in the universe, neutron stars are unique cosmic laboratories. The study of neutron stars provides an ideal theoretical testbed for investigating both physics at supra-nuclear densities as well as fundamental physics. Their global astrophysical properties however depend strongly on the star's internal structure, which is currently unknown due to uncertainties in the equation of state. In recent years, a lot of work has revealed the existence of universal relations between stellar quantities that are insensitive to the equation of state. At the same time, the fields of multimessenger astronomy and machine learning have both advanced significantly. As such, there has been a confluence of research into their combination and the field is growing. In this paper, we develop universal relations for rapidly rotating neutron stars, by using supervised machine learning methods, thus proposing a new way of discovering and validating such relations. The analysis is performed for tabulated hadronic, hyperonic, and hybrid EoS-ensembles that obey the multimessenger constraints and cover a wide range of stiffnesses. The relations discussed could provide an accurate tool to constrain the equation of state of nuclear matter when measurements of the relevant observables become available.

Photoionization and its inverse, electron-ion recombination, are key processes that influence many astrophysical plasmas (and gasses), and the diagnostics that we use to analyse the plasmas. In this review we provide a brief overview of the importance of photoionization and recombination in astrophysics. We highlight how the data needed for spectral analyses, and the required accuracy, varies considerably in different astrophysical environments. We then discuss photoionization processes, highlighting resonances in their cross-sections. Next we discuss radiative recombination, and low and high temperature dielectronic recombination. The possible suppression of low temperature dielectronic recombination (LTDR) and high temperature dielectronic recombination (HTDR) due to the radiation field and high densities is discussed. Finally we discuss a few astrophysical examples to highlight photoionization and recombination processes.

We present the first results from a new survey for high-redshift $(z\gtrsim5)$ gravitationally lensed quasars and close quasar pairs. We carry out candidate selection based on the colors and shapes of objects in public imaging surveys, then conduct follow-up observations to confirm the nature of high-priority candidates. In this paper, we report the discoveries of J0025--0145 ($z=5.07$) which we identify as an {intermediately-lensed quasar, and J2329--0522 ($z=4.85$) which is a kpc-scale close quasar pair. The {\em Hubble Space Telescope (HST)} image of J0025--0145 shows a foreground lensing galaxy located $0\farcs6$ away from the quasar. However, J0025--0145 does not exhibit multiple lensed images of the quasar, and we identify J0025--0145 as an intermediate lensing system (a lensing system that is not multiply imaged but has a significant magnification). The spectrum of J0025--0145 implies an extreme Eddington ratio if the quasar luminosity is intrinsic, which could be explained by a large lensing magnification. The {\em HST} image of J0025--0145 also indicates a tentative detection of the quasar host galaxy in rest-frame UV, illustrating the power of lensing magnification and distortion in studies of high-redshift quasar host galaxies. J2329--0522 consists of two resolved components with significantly different spectral properties, and a lack of lensing galaxy detection under sub-arcsecond seeing. We identify it as a close quasar pair, which is the highest confirmed kpc-scale quasar pair to date. We also report four lensed quasars and quasar pairs at $2<z<4$, and discuss possible improvements to our survey strategy.

Here we present our current update of CLOUDY on gas-phase chemical reactions for the formation and destruction of the SiS molecule, its energy levels, and collisional rate coefficients with H$_2$, H, and He over a wide range of temperatures. As a result, henceforth the spectral synthesis code CLOUDY predicts SiS line intensities and column densities for various astrophysical environments.

We present observations from X-ray to mid-infrared wavelengths of the most energetic non-quasar transient ever observed, AT2021lwx. Our data show a single optical brightening by a factor $>100$ to a luminosity of $7\times10^{45}$ erg s$^{-1}$, and a total radiated energy of $1.5\times10^{53}$ erg, both greater than any known optical transient. The decline is smooth and exponential and the ultra-violet - optical spectral energy distribution resembles a black body with temperature $1.2\times10^4$ K. Tentative X-ray detections indicate a secondary mode of emission, while a delayed mid-infrared flare points to the presence of dust surrounding the transient. The spectra are similar to recently discovered optical flares in known active galactic nuclei but lack some characteristic features. The lack of emission for the previous seven years is inconsistent with the short-term, stochastic variability observed in quasars, while the extreme luminosity and long timescale of the transient disfavour the disruption of a single solar-mass star. The luminosity could be generated by the disruption of a much more massive star, but the likelihood of such an event occurring is small. A plausible scenario is the accretion of a giant molecular cloud by a dormant black hole of $10^8 - 10^9$ solar masses. AT2021lwx thus represents an extreme extension of the known scenarios of black hole accretion.

In this work, we studied the Galactic Center supermassive black hole (SMBH), Sagittarius A* (Sgr A*), with the KVN and VERA Array (KaVA)/East Asian VLBI Network (EAVN) monitoring observations. Especially on 13 May 2019, Sgr A* experienced an unprecedented bright near infra-red (NIR) flare; so, we find a possible counterpart at 43 GHz (7 mm). As a result, a large temporal variation of the flux density at the level 15.4%, with the highest flux density of 2.04 Jy, is found on 11 May 2019. Interestingly, the intrinsic sizes are also variable, and the area and major-axis size show a marginal correlation with flux density with >2{\sigma}. Thus, we interpret that the emission region at 43 GHz follows the larger-when-brighter relation in 2019. The possible origins are discussed with an emergence of a weak jet/outflow component and the position angle change of the rotation axis of the accretion disk in time.

We study the formation and early evolution of young stellar objects (YSOs) using three-dimensional non-ideal magnetohydrodynamic (MHD) simulations to investigate the effect of cosmic ray ionization rate and dust fraction (or amount of dust grains) on circumstellar disk formation. Our simulations show that a higher cosmic ray ionization rate and a lower dust fraction lead to (i) a smaller magnetic resistivity of ambipolar diffusion, (ii) a smaller disk size and mass, and (iii) an earlier timing of outflow formation and a greater angular momentum of the outflow. In particular, at a high cosmic ray ionization rate, the disks formed early in the simulation are dispersed by magnetic braking on a time scale of about 104 years. Our results suggest that the cosmic ray ionization rate has a particularly large impact on the formation and evolution of disks, while the impact of the dust fraction is not significant.

When the Galactic Cosmic Rays (GCRs) entering the heliosphere, they encounter the solar wind plasma, and their intensity is reduced, so-called solar modulation. The modulation is caused by the combination of a few factors, such as particle energies, solar activity and solar disturbance. In this work, a 2D numerical method is adopted to simulate the propagation of GCRs in the heliosphere with SOLARPROP, and to overcome the time-consuming issue, the machine learning technique is also applied. With the obtained proton local interstellar spectra (LIS) based on the observation from Voyager 1 and AMS-02, the solar modulation parameters during the solar maximum activity of cycle 24 have been found. It shows the normalization and index of the diffusion coefficient indeed reach a maximal value in February 2014. However, after taking into account the travel time of particles with different energies, the peak time was found postponed to November 2014 as expected. The nine-month late is so-called time lag.

A long-standing issue in astrobiology is whether planets orbiting the most abundant type of stars, M-dwarfs, can support liquid water and eventually life. A new study shows that subglacial melting may provide an answer, significantly extending the habitability region, in particular around M-dwarf stars, which are also the most promising for biosignature detection with the present and near-future technology.

The cosmic origin of the elements, the fundamental chemical building blocks of the Universe, is still uncertain. Binary interactions play a key role in the evolution of many massive stars, yet their impact on chemical yields is poorly understood. Using the MESA stellar evolution code we predict the chemical yields ejected in wind mass loss and the supernovae of single and binary-stripped stars. We do this with a large 162 isotope nuclear network at solar-metallicity. We find that binary-stripped stars are more effective producers of the elements than single stars, due to their increased mass loss and an increased chance to eject their envelopes during a supernova. This increased production by binaries varies across the periodic table, with Fluorine and Potassium being more significantly produced by binary-stripped stars than single stars. We find that the C12/C13 could be used as an indicator of the conservativeness of mass transfer, as C13 is preferentially ejected during mass transfer while C12 is preferentially ejected during wind mass loss. We identify a number of gamma-ray emitting radioactive isotopes that may be used to help constrain progenitor and explosion models of core-collapse supernovae with next-generation gamma-ray detectors. For single stars we find V44 and Mn52 are strong probes of the explosion model, while for binary-stripped stars it is Cr48. Our findings highlight that binary-stripped stars are not equivalent to two single stars and that detailed stellar modelling is needed to predict their final nucleosynthetic yields.

(Abridged) Context. Most massive stars are located in multiple stellar systems. Magnetic fields are believed to be essential in the accretion and ejection processes around single massive protostars. Aims. Our aim is to unveil the influence of magnetic fields in the formation of multiple massive stars, in particular on the fragmentation modes and properties of the multiple protostellar system. Methods. Using RAMSES, we follow the collapse of a massive pre-stellar core with (non-ideal) radiation-(magneto-)hydrodynamics. We choose a setup which promotes multiple stellar system formation. Results. In the purely hydrodynamical models, we always obtain (at least) binary systems. When more than two stars are present, their gravitational interaction triggers mergers until there are two stars left. The following gas accretion increases their orbital separation and hierarchical fragmentation occurs so that both stars host a comparable disk and stellar system which then form similar disks as well. We identify several modes of fragmentation: Toomre-unstable disk fragmentation, arm-arm collision and arm-filament collision. Disks grow in size until they fragment and become truncated as the newly-formed companion gains mass. When including magnetic fields, the picture evolves: the primary disk produces less fragments, arm-filament collision is absent. Magnetic fields reduce the initial orbital separation but do not affect its further evolution, which is mainly driven by gas accretion. With magnetic fields, the growth of individual disks is regulated even in the absence of fragmentation or truncation. Conclusions. Hierarchical fragmentation is seen in unmagnetized and magnetized models. Magnetic fields, including non-ideal effects, are important because they remove certain fragmentation modes and limit the growth of disks, which is otherwise only limited through fragmentation.

The studies of FUors, EXors and other young eruptive stars are very important for the understanding of the earliest stages of pre-main-sequence evolution. We describe the current situation in this field. This is the short version of the review, presented in ``Non-stationary processes in the protoplanetary disks and their observational manifestations'' conference, Crimean Astrophysical Observatory

Common-mode choke inductors are useful tools for resolving grounding issues in large detector systems. Using inductive components on cryogenic pre-amplifier boards has so far been prevented by the poor low-temperature performance of common ferrite materialssuch as NiZn and MnZn. Recentlydeveloped nanocrystallineand amorphous ferrite materials promise improved performance up to the point where using magnetics at liquid mitrogen temperatures becomes feasible. This research applies the work of Yin et al. on characterizing ferrite materials by constructing and testing a common mode choke inductor for use on detector pre-amplifiers for the ELT first generation instruments.

Context. Extreme precision radial velocity (RV) surveys seeking to detect planets at RV semi-amplitudes of 10 cm/s are facing numerous challenges. One of those challenges is convective blueshift caused by stellar granulation and its suppression through magnetic activity which plays a significant role in hiding planetary signals in stellar jitter. Aims. Previously we found that for main sequence stars, convective blueshift as an observational proxy for the strength of convection near the stellar surface strongly depends on effective temperature. In this work we investigate 242 post main sequence stars, covering the subgiant, red giant, and asymptotic giant phases and empirically determine the changes in convective blueshift with advancing stellar evolution. Methods. We used the third signature scaling approach to fit a solar model for the convective blueshift to absorption-line shift measurements from a sample of coadded HARPS spectra, ranging in temperature from 3750 K to 6150 K. We compare the results to main sequence stars of comparable temperatures but with a higher surface gravity. Results. We show that convective blueshift becomes significantly stronger for evolved stars compared to main sequence stars of a similar temperature. The difference increases as the star becomes more evolved, reaching a 5x increase below 4300 K for the most evolved stars. The large number of stars in the sample, for the first time, allowed for us to empirically show that convective blueshift remains almost constant among the entire evolved star sample at roughly solar convection strength with a slight increase from the red giant phase onward. We discover that the convective blueshift shows a local minimum for subgiant stars, presenting a sweet spot for exoplanet searches around higher mass stars, by taking advantage of their spin-down during the subgiant transition.

We present a new analysis of the rest-frame UV and optical spectra of a sample of three $z>8$ galaxies discovered behind the gravitational lensing cluster RX J2129.4+0009. We combine these observations with those of a sample of $z>7.5$ galaxies from the literature, for which similar measurements are available. As already pointed out in other studies, the high [OIII]$\lambda$5007/[OII]$\lambda$3727 ratios ($O_{32}$) and steep UV continuum slopes ( $\beta$ ) are consistent with the values observed for low redshift Lyman continuum emitters, suggesting that such galaxies contribute to the ionizing budget of the intergalactic medium. We construct a logistic regression model to estimate the probability of a galaxy being a Lyman continuum emitter based on the measured $M_{UV}$, $\beta$, and $O_{32}$ values. Using this probability together with the UV luminosity function, we construct an empirical model that estimates the contribution of high redshift galaxies to reionization based on these observable quantities. Our analysis shows that at $z\sim8$, the average escape fraction of the galaxy population (i.e., including both LyC emitters and non-emitters) varies with $M_{UV}$, with brighter galaxies having larger $f_{esc}$. For $M_{UV}$ $< -$19 we find an average escape fraction of 20$\%$, decreasing to almost zero for $M_{UV}$$>-16$. Galaxies with intermediate UV luminosity ($-19 <$ $M_{UV}$ $< -16$) contribute half of the ionizing photons. The relative contribution of faint versus bright galaxies depends on redshift, with UV bright galaxies ($-23 <$ $M_{UV}$ $< -19$) becoming more important over time and reaching $\approx 40\%$ at the end of reionization around $z=6$.

We report the results of a search for long-period ($100<P<600$ days) periodic variability in SDSS Stripe 82 standards catalog. The SDSS coverage of Stripe 82 enables such a search because there are on average 20 observations per band in $ugriz$ bands for about 1 million sources, collected over about 6 years, with a faint limit of $r\sim22$ mag and precisely calibrated 1-2% photometry. We calculated the periods of candidate variable sources in this sample using the Lomb-Scargle periodogram and considered the three highest periodogram peaks in each of the $gri$ filters as relevant. Only those sources with $gri$ periods consistent within 0.1% were later studied. We use the Kuiper statistic to ensure uniform distribution of data points in phased light curves. We present 5 sources with the spectra consistent with quasar spectra and plausible periodic variability. This SDSS-based search bodes well for future sensitive large-area surveys, such as the Rubin Observatory Legacy Survey of Space and Time, which, due to its larger sky coverage (about a factor of 60) and improved sensitivity ($\sim2$ mag), will be more powerful for finding such sources.

Recent detections of extremely short-timescale microlensing events imply the existence of a large population of Earth- to Neptune-mass planets that appear to have no host stars. However, it is currently unknown whether these objects are truly free-floating planets or whether they are in wide orbits around a distant host star. Here, we present the analysis of high-resolution imaging observations of six free-floating planet candidates collected with the Keck telescope. If these candidates were actually wide-orbit planets, the light of the host would appear at a separation of 40-60 mas from the microlensing source star. No such stars are detected. We carry out injection and recovery simulations to estimate the sensitivity to putative host stars at different separations. Depending on the object, the presented observations rule out 10-35% of potential hosts assuming that the probability of hosting a planet does not depend on the host mass. If we assume that more massive stars are more likely to host a planet, we rule out 18-63% of potential hosts. We argue that deeper observations, for example with JWST, are needed to confidently confirm or refute the free-floating planet hypothesis.

When considering the migration of Jupiter and Saturn, a classical result is to find the planets migrating outwards and locked in the 3:2 mean motion resonance (MMR). These results were obtained in the framework of viscously accreting discs, in which the observed stellar accretion rates constrained the viscosity values. However, it has recently been shown observationally and theoretically that discs are probably less viscous than previously thought. Therefore, in this paper, we explore the dynamics of pairs of giant planets in low-viscosity discs. We performed two-dimensional hydrodynamical simulations using the grid-based code FARGOCA. In contrast to classical viscous discs, we find that the outer planet never crosses the 2:1 resonance and the pair does not migrate outwards. After a wide parameter exploration, including the mass of the outer planet, we find that the planets are primarily locked in the 2:1 MMR and in some cases in the 5:2 MMR. We explain semi-analytically why it is not possible for the outer planet to cross the 2:1 MMR in a low-viscosity disc. We find that pairs of giant planets migrate inwards in low-viscosity discs. Although, in some cases, having a pair of giant planets can slow down the migration speed with respect to a single planet. Such pairs of slowly migrating planets may be located, at the end of the disc phase, in the population of exoplanets of 'warm Jupiters'. However, the planets never migrate outwards. These results could have strong implications on the Solar System's formation scenarios if the Sun's protoplanetary disc had a low viscosity.

Atmospheric muon neutrinos are produced by meson decays in cosmic-ray-induced air showers. The flux depends on meteorological quantities such as the air temperature, which affects the density of air. Competition between decay and re-interaction of those mesons in the first particle production generations gives rise to a higher neutrino flux when the air density in the stratosphere is lower, corresponding to a higher temperature. A measurement of a temperature dependence of the atmospheric $\nu_{\mu}$ flux provides a novel method for constraining hadro\-nic interaction models of air showers. It is particularly sensitive to the production of kaons. Studying this temperature dependence for the first time requires a large sample of high-energy neutrinos as well as a detailed understanding of atmospheric properties. We report the significant ($> 10 \sigma$) observation of a correlation between the rate of more than 260,000 neutrinos, detected by IceCube between 2012 and 2018, and atmospheric temperatures of the stratosphere, measured by the Atmospheric Infrared Sounder (AIRS) instrument aboard NASA's AQUA satellite. For the observed 10$\%$ seasonal change of effective atmospheric temperature we measure a 3.5(3)$\%$ change in the muon neutrino flux. This observed correlation deviates by about 2-3 standard deviations from the expected correlation of 4.3$\%$ as obtained from theoretical predictions under the assumption of various hadronic interaction models

Obscuration events in type I active galactic nuclei (AGN) have been detected more frequently in recent years. The strong flux decrease in the soft X-ray band between observations has been caused by clouds with large column densities transiting our line-of-sight (LOS) and covering the central AGN. Another event has been captured in NGC 3227 at the end of 2019. We aim to determine the nature of the observed spectral variability in 2019 obscuration event. We split the two XMM-Newton observations from 2019 into timing bins of length $\sim$ 10 ks. We used the SPEX code to analyse the 0.35-10 keV EPIC-PN spectra of each timing bin. In the first observation (Obs 1), there is a strong anti-correlation between the column density ($N_H$) of the obscurer and the continuum normalisations of the X-ray power-law and soft Comptonisation components ($N_{pow}$ and $N_{comt}$, respectively). The powerlaw continuum models the hard X-rays produced by the corona, and the Comptonisation component models the soft X-ray excess and emission from the accretion disk. Through further testing we conclude that the continuum is likely to drive the observed variability, but we cannot rule out a possible contribution from NH of the obscurer if it fully transverses across the ionising source within our LOS during the observation. The ionisation parameter ($\xi$) of the obscurer is not easily constrained, and therefore it is not clear whether it varies in response to changes in ionising continuum. The second observation (Obs 2) displays a significantly lower count rate due to the combination of a high NH and covering fraction of the obscurer, and a lower continuum flux. The observed variability seen during the obscuration event of NGC 3227 in 2019 is likely driven by the continuum, but the obscurer varies at the same time, making it difficult to distinguish between the two possibilities with full certainty.

A full understanding of how unusually large "Ultra Diffuse Galaxies" (UDGs) fit into our conventional theory of dwarf galaxy formation remains elusive, despite the large number of objects identified locally. A natural extension of UDG research is the study of similar galaxies at higher redshift to establish how their properties may evolve over time. However, this has been a challenging task given how severely cosmological surface brightness dimming inhibits our ability to detect low-surface brightness galaxies at high-$z$. Here, we present a sample of low stellar surface density galaxies (LDGs) at moderate redshift, likely the progenitors of local UDGs, identified in deep near-IR {\it JWST} observations of the El Gordo cluster at $z=0.87$. By stacking 8 NIRCAM filters, reach an apparent surface brightness sensitivity of 24.59~mag~arcsec$^{-2}$, faint enough to be complete to the bright end of the LDG population. Our analysis identifies significant differences between this population and UDGs observed locally, including differences in their color and size distributions, which suggest that the UDG progenitors at high-$z$ are bluer and more extended than UDGs at $z=0$. This suggests that multiple mechanisms are responsible for UDG formation and that prolonged transformation of cluster dwarfs is not a primary UDG formation mechanism at high-$z$. Furthermore, we find a slight overabundance of LDGs in El Gordo, and, in contrast to findings in local clusters, our analysis does not show a deficit of LDGs in the center of El Gordo, implying that tidal destruction of LDGs is significant between $z=0.87$ and $z=0$.

The Large Interferometer For Exoplanets (LIFE) initiative aims to develop a space based mid-infrared (MIR) nulling interferometer to measure the thermal emission spectra of temperate terrestrial exoplanets. We investigate how well LIFE could characterize a cloudy Venus-twin exoplanet to: (1) test our retrieval routine on a realistic non-Earth-like MIR spectrum of a known planet, (2) investigate how clouds impact retrievals, (3) refine the LIFE requirements derived in previous Earth-centered studies. We run retrievals for simulated LIFE observations of a Venus-twin exoplanet orbiting a Sun-like star located 10 pc from the observer. By assuming different models (cloudy and cloud-free) we analyze the performance as a function of the quality of the LIFE observation. This allows us to determine how well atmosphere and clouds are characterizable depending on the quality of the spectrum. Our study shows that the current minimal resolution ($R=50$) and signal-to-noise ($S/N=10$ at $11.2\mu$m) requirements for LIFE suffice to characterize the structure and composition of a Venus-like atmosphere above the cloud deck if an adequate model is chosen. However, we cannot infer cloud properties. The accuracy of the retrieved planet radius ($R_{pl}$), equilibrium temperature ($T_{eq}$), and Bond albedo ($A_B$) depend on the choice of model. Generally, a cloud-free model performs best and thus the presence of clouds cannot be inferred. This model dependence of retrieval results emphasizes the importance of developing a community-wide best-practice for atmospheric retrieval studies. If we consider higher quality spectra (especially $S/N=20$), we can infer the presence of clouds and pose first constraints on their structure.

The magnetic activity of a star -- which modulates the stellar wind outflow -- shapes the immediate environments of orbiting planets and induces atmospheric loss thereby impacting their habitability. We perform a detailed parameter space study using three dimensional magnetohydrodynamic simulations to understand the effect of changing stellar wind magnetic field and planetary magnetic field strengths on planetary magnetospheric topology and atmospheric losses. It is observed that the relative strengths of stellar and planetary magnetic fields play a significant role in determining the steady state magnetospheric configuration and atmospheric erosion. When the stellar field is strengthened or the planetary field is weakened, stellar magnetic field accumulation occurs at the day-side of the planet which forces the magnetopause to shift closer to its surface. The magnetotail opens up leading to the formation of Alfv\'{e}n wings in the night-side wake region. We demonstrate how reconnection processes and wind conditions lead to the bifurcation of the magnetotail current sheet. With increasing stellar wind magnetic field strength, the day-side reconnection point approaches the planet thereby enhancing mass loss. We establish an analytic equation which successfully captures the modeled mass-loss rate variations of planets with changing magnetic field strengths. Our results are relevant for understanding how the interplay of stellar and planetary magnetism influence (exo)planetary environments and their habitability in star-planet systems with differing relative magnetic field strengths, or in a single star-planet system over the course of their evolution with age.

Stellar standard candles provide the absolute calibration for measuring the Universe's local expansion rate, H0, which differs by ~8% from the value inferred using the Cosmic Microwave Background assuming the concordance cosmological model, LambdaCDM. This Hubble tension indicates a need for important revisions of fundamental physics. However, the accuracy of the H0 measurement based on classical Cepheids has been challenged by a measurement based on the Tip of the Red Giant Branch (TRGB) method. A resolution of the Cepheids vs. TRGB dispute is needed to demonstrate well understood systematics and to corroborate the need for new physics. Here, we present an unprecedented 1.39% absolute calibration of the TRGB distance scale based on small-amplitude red giant stars (OSARGs). Our results improve by 20% upon previous calibrations and are limited by the accuracy of the distance to the Large Magellanic Cloud. This precision gain is enabled by the realization that virtually all stars at the TRGB are variable - a fact not previously exploited for TRGB calibration. Using observations and extensive simulations, we demonstrate that OSARGs yield intrinsically precise and accurate TRGB measurements thanks to the shape of their luminosity function. Inputting our calibration to the Carnegie Chicago Hubble Program's H0 analysis yields a value of H0 = 71.8 +/- 1.5 km/s/Mpc, in < 1 sigma agreement with the Cepheids-based H0 value and in 2.8 sigma tension with the early-Universe value.

In this work we consider a class of interacting vacuum corresponding to a generalised Chaplygin gas (gCg) cosmology. In particular we analyse two different scenarios at perturbation level for the same background interaction characterised by the parameter $\alpha$: (i) matter that follows geodesics, corresponding to homogeneous vacuum, and (ii) a covariant ansatz for vacuum density perturbations. In the latter case, we show that the vacuum perturbations are very tiny as compared to matter perturbations on sub-horizon scales. In spite of that, depending on the value of the Chaplygin gas parameter $\alpha$, vacuum perturbations suppress or enhance the matter growth rate as compared to the case (i). We use Cosmic Microwave Background (CMB), type Ia supernovae (SNe) and Redshift Space Distortion (RSD) measurements to test the observational viability of the model. We found that the mean value of our joint analysis clearly favours a positive interaction, i.e., an energy flux from dark matter to dark energy, with $\alpha \approx 0.143$ in both cases, while the cosmological standard model, recovered for $\alpha$=0, is ruled out by 3$\sigma$ confidence level. Noteworthy, the positive value of interaction can alleviate both the $H_0$ and $S_8$ tension for the dataset considered here.

A recent observational study found that the projected spatial distributions of the satellites of bright, isolated host galaxies tend to be lopsided with respect to the locations of the hosts. Here, we examine the spatial distributions of the satellites of a large number of bright, isolated host galaxies that were obtained from mock redshift surveys of a LCDM simulation. Host galaxies and their satellites were identified using selection criteria that are identical to those used in the observational study, allowing a direct comparison of the results for the simulated and observed systems. To characterize the spatial distribution of the satellites, we adopt two statistics: [1] the pairwise clustering of the satellites and [2] the Mean Resultant Length. In agreement with the observational study, we find a strong tendency for satellites in the simulation to be located on the same side of their host, and the signal is most pronounced for the satellites of blue hosts. These lopsided satellite distributions are not solely attributable to incompleteness of the observed satellite catalog or the presence of objects that have been falsely identified as satellites. In addition, satellites that joined their hosts' halos in the distant past (> 8 Gyr) show a pronounced lopsidedness in their spatial distributions and, therefore, the lopsidedness is not solely attributable to late-time accretion of satellites.

Modern sciences and astrophysics in particular study objects and phenomena not visible in physical terms, that is they cannot be investigated with the eyes or analogous optical systems. Nevertheless, they make intensive use of visual representations, showing data in a figurative way, using lights, colors, and shapes familiar to the user and aesthetically pleasant. Besides being inaccessible for Blind and Visually Impaired (BVI) users, such figurative visual representation can lead to misunderstandings about the real nature of the represented object if the code of representation is not declared. We argue that multi-sensory representations clearly arbitrary, i.e., that do not aim to imitate reality, are a valid choice for an effective meaning-making process of astronomical science for all. In an equity perspective, multi-sensory representations also create an effective common ground for inclusion among people with diverse abilities, skills, and learning styles, in the framework of Universal Design for Learning. In order to investigate our hypothesis we designed two mono-sensory representation (one only haptic and the other only acoustic) and tested them in individual and group workshops with both sighted and BVI users. We then used our results to guide the design of a multi-sensory representation of non-visible astronomical data including visual, acoustic, and haptic stimuli. We tested this representation as well, in order to refine and propose it to the public. The result is the exhibit Sense the Universe, to be used for outreach and education. Sense the Universe was part of a museum exhibition attended both by sighted and BVI users. Our findings suggest the validity of multi-sensory representations for a truly and effective engagement in scientific learning, both in terms of intelligibility and persistence of scientific contents and of a more equal access to scientific culture.

We obtain the non-equilibrium condensate of the Chern Simons density induced by a misaligned homogeneous coherent axion field in linear response. The Chern-Simons dynamical susceptibility is simply related to the axion self-energy, a result that is valid to leading order in the axion coupling but to all orders in the couplings of the gauge fields to other fields within or beyond the standard model except the axion. The induced Chern-Simons density requires renormalization which is achieved by vacuum subtraction. For ultralight axions of mass $m_a$ coupled to electromagnetic fields with coupling $g$, the renormalized high temperature contribution post-recombination is $\langle \vec{E}\cdot\vec{B}\rangle(t) = -\frac{g\,\pi^2\,T^4}{15} \,\overline{a}(t)+ \frac{g \,m^2_a\,T}{16\,\pi}\,\dot{\overline{a}}(t) $ with $\overline{a}(t)$ the dynamical homogeneous axion condensate. We conjecture that emergent axion-like quasiparticle excitations in condensed matter systems may be harnessed to probe cosmological axions and the Chern-Simons condensate. Furthermore, it is argued that a misaligned axion can also induce a non-abelian Chern-Simons condensate of similar qualitative form, and can also ``seed'' chiral symmetry breaking and induce a neutral pion condensate after the QCD phase transition.

We construct a viable model of the vector coherent oscillation dark matter. The vector boson is coupled to the inflaton through the kinetic function so that the effective Hubble mass term is cancelled out. In order to avoid strong constraints from isocurvature perturbation and statistically anisotropic curvature perturbation, the inflaton is arranged so that it does not contribute to the observed large scale curvature perturbation and we introduce a curvaton. We found viable vector coherent oscillation dark matter scenario for the wide vector mass range from $10^{-21}\,{\rm eV}$ to $1\,{\rm eV}$.

The Laser Interferometer Space Antenna (LISA) will provide us with a unique opportunity to observe the early inspiral phase of supermassive binary black holes (SMBBHs) in the mass range of $10^5-10^6\,M_{\odot}$, that lasts for several years. It will also detect the merger and ringdown phases of these sources. Therefore, such sources are extremely useful for multiparameter tests of general relativity (GR), where parametrized deviations from GR at multiple post-Newtonian orders are simultaneously measured, thus allowing for a rigorous test of GR. However, the correlations of the deviation parameters with the intrinsic parameters of the system make multiparameter tests extremely challenging to perform. We demonstrate the use of principal component analysis (PCA) to obtain a new set of deviation parameters, which are best-measured orthogonal linear combinations of the original deviation parameters. With the observation of an SMBBH of total redshifted mass, $\sim\mathrm{7\times10^5\,M_{\odot}}$ at a luminosity distance of 3 Gpc, we can estimate the five most dominant PCA parameters, with 1-$\sigma$ statistical uncertainty of $\lesssim 0.2$. The two most dominant PCA parameters can be bounded to $\sim \mathcal{O}(10^{-4})$, while the third and fourth-dominant ones to $\sim \mathcal{O}(10^{-3})$. Measurement of the PCA parameters with such unprecedented precision with LISA makes them an excellent probe to test the overall PN structure of the GW phase evolution.

The QCD axion bubbles can be formed due to an extra Peccei-Quinn (PQ) symmetry breaking in the early Universe. In this paper, we investigate the QCD axion bubbles formation from the PQ symmetry broken by hidden $SU(N)_H$ gauge interactions after inflation, which leads to the multiple vacua. The axion acquires a light mass and then settles down into different vacua. The QCD axion bubbles are formed when the conventional QCD axion arises during the QCD phase transition. In our scenario, the QCD axions that start to oscillate at the large values $\sim2\pi/3$ can lead to the high density axion bubbles with $N=2$. The cosmological implications of the QCD axion bubbles are also discussed, such as the primordial black holes (PBHs) and the axion miniclusters. We find that the PBH mass is lager than $\sim\mathcal{O}(5\times10^5)M_\odot$ for the axion scale $f_a\sim\mathcal{O}(10^{16})\, \rm GeV$.

In this work, we study general features of a regime where gauge fields produced during inflation cause a strong backreaction on the background evolution and its impact on the spectrum and the correlation length of gauge fields. With this aim, the gradient-expansion formalism previously proposed for the description of inflationary magnetogenesis in purely kinetic or purely axial coupling models, is extended to the case when both types of coupling are present. As it is formulated in position space, this method allows us to self-consistently take into account the backreaction of generated gauge fields on the inflationary background because it captures the nonlinear evolution of all physically relevant gauge-field modes at once. Using this extended gradient-expansion formalism, suitable for a wide range of inflationary magnetogenesis models, we study the gauge-field production in a specific generalization of the Starobinsky $R^2$-model with a nonminimal coupling of gauge fields to gravity. In the Einstein frame, this model implies, in addition to an asymptotically flat inflaton potential, also a nontrivial form of kinetic and axial coupling functions which decrease in time and, thus, are potentially suitable for the generation of gauge fields with a scale-invariant or even red-tilted power spectrum. The numerical analysis shows, however, that backreaction, which unavoidably occurs in this model for the interesting range of parameters, strongly alters the behavior of the spectrum and does not allow to obtain a sufficiently large correlation length for the magnetic field. The oscillatory behavior of the generated field, caused by the retarded response of the gauge field to changes of the inflaton velocity, was revealed.

We study the prospects for detection of continuous gravitational signals from "normal" Galactic neutron stars. We use a synthetic population generated by evolving stellar remnants in time, according to several models. We also briefly treat the case of recycled neutron stars. We consider the most recent constraints set by all-sky searches for continuous gravitational waves and use them for our detectability criteria. We discuss detection prospects for the current and the next generation of gravitational wave detectors. We find that neutron stars whose ellipticity is solely caused by magnetic deformations cannot produce any detectable signal, not even by 3rd-generation detectors. Currently detectable sources all have $B<10^{12}$G and show a strong correlation between magnetic field and ellipticity: the larger the magnetic field, the larger the ellipticity. According to our models, 3rd-generation detectors as the Einstein Telescope and Cosmic Explorer will be able to detect up to $\approx$ 250 more sources than the current detectors. Continuous waves from recycled neutron stars will likely remain elusive to detection by current detectors but should be detectable with the next generation of detectors.

Previous studies on high-energy gamma-ray burst neutrinos from IceCube suggest a neutrino speed variation at the Lorentz violation~(LV) scale of $\sim 6.4\times 10^{17}$~GeV, with opposite velocity variances between neutrinos and antineutrinos. Within a space-time foam model, inspired by string theory, we develop an approach to describe the suggested neutrino/antineutrino propagation properties with both Lorentz invariance and CPT symmetry breaking. A threshold analysis on the bremsstrahlung of electron-positron pair~($\nu\rightarrow\nu ee^{+}$) for the superluminal~(anti)neutrino is performed. We find that, due to the energy violation caused by the quantum foam, such reaction may be restricted to occur at sufficient high energies and could even be kinematically forbidden. Constraints on neutrino LV from vacuum $ee^{+}$ pair emission are naturally avoided. Future experiments are appealed to test further the CPT violation of cosmic neutrinos and/or neutrino superluminality.