Papers for Monday, Jun 21 2021

The FeII/MgII line flux ratio has been used as an indicator of the Fe/Mg abundance ratio in the broad line region (BLR) of active galactic nuclei (AGNs). On the basis of archival rest-frame UV spectra obtained via the Hubble Space Telescope and the Sloan Digital Sky Survey, we investigate the FeII/MgII ratios of type 1 AGNs at z < 2. Over wide dynamic ranges of AGN properties (i.e., black hole mass, AGN luminosity, and Eddington ratio), we confirm that the FeII/MgII ratio strongly correlates with Eddington ratio but not with black hole mass, AGN luminosity, or redshift. Our results suggest that the metallicity in the BLR are physically related to the accretion activity of AGNs, but not to the global properties of galaxies (i.e., galaxy mass and luminosity). With regard to the relation between the BLR metallicity and the accretion rate of AGNs, we discuss that metal cooling may play an important role in enhancing the gas inflow into the central region of host galaxies, resulting in the high accretion rate of AGNs.

We present a simple regulator-type framework designed specifically for modelling formation of dwarf galaxies. We explore sensitivity of model predictions for the stellar mass--halo mass and stellar mass--metallicity relations to different modelling choices and parameter values. Despite its simplicity, when coupled with realistic mass accretion histories of haloes from simulations and reasonable choices for model parameter values, the framework can reproduce a remarkably broad range of observed properties of dwarf galaxies over seven orders of magnitude in stellar mass. In particular, we show that the model can simultaneously match observational constraints on the stellar mass-halo mass relation, as well as observed relations between stellar mass and gas phase and stellar metallicities, gas mass, size, and star formation rate, as well as general form and diversity of star formation histories (SFHs) of observed dwarf galaxies. The model can thus be used to predict photometric properties of dwarf galaxies hosted by dark matter haloes in $N$-body simulations, such as colors, surface brightnesses, and mass-to-light ratios and to forward model observations of dwarf galaxies. We present examples of such modelling and show that colors and surface brightness distributions of model galaxies are in good agreement with observed distributions for dwarfs in recent observational surveys. We also show that in contrast with the common assumption, the absolute magnitude-halo mass relation is generally predicted to have a non-power law form in the dwarf regime, and that the fraction of haloes that host detectable ultrafaint galaxies is sensitive to reionization redshift (zrei) and is predicted to be consistent with observations for zrei<~9.

Since planet occurrence and primordial atmospheric retention probability increase with period, the occurrence-weighted median planets discovered by transit surveys may bear little resemblance to the low-occurrence, short-period planets sculpted by atmospheric escape ordinarily used to calibrate mass--radius relations and planet formation models. An occurrence-weighted mass--radius relation for the low-mass planets discovered so far by transit surveys orbiting solar-type stars requires both occurrence-weighted median Earth-mass and Neptune-mass planets to have a few percent of their masses in hydrogen/helium (H/He) atmospheres. Unlike the Earth that finished forming long after the protosolar nebula was dissipated, these occurrence-weighted median Earth-mass planets must have formed early in their systems' histories. The existence of significant H/He atmospheres around Earth-mass planets confirms an important prediction of the core-accretion model of planet formation. It also implies core masses $M_{\text{c}}$ in the range $2~M_{\oplus}\lesssim M_{\text{c}}\lesssim 8~M_{\oplus}$ that can retain their primordial atmospheres. If atmospheric escape is driven by photoevaporation due to extreme ultraviolet (EUV) flux, then our observation requires a reduction in the fraction of incident EUV flux converted into work usually assumed in photoevaporation models. If atmospheric escape is core driven, then the occurrence-weighted median Earth-mass planets must have large Bond albedos. In contrast to Uranus and Neptune that have at least 10% of their masses in H/He atmospheres, these occurrence-weighted median Neptune-mass planets are H/He poor. The implication is that they experienced collisions or formed in much shorter-lived and/or hotter parts of their parent protoplanetary disks than Uranus and Neptune's formation location in the protosolar nebula.

We investigate whether the recently suggested rotation and crystallization driven dynamo can explain the apparent increase of magnetism in old metal polluted white dwarfs. We find that the effective temperature distribution of polluted magnetic white dwarfs is in agreement with most/all of them having a crystallizing core and increased rotational velocities are expected due to accretion of planetary material which is evidenced by the metal absorption lines. We conclude that a rotation and crystallization driven dynamo offers not only an explanation for the different occurrence rates of strongly magnetic white dwarfs in close binaries, but also for the high incidence of weaker magnetic fields in old metal polluted white dwarfs.

We introduce and apply a new approach to probe the response of galactic stellar haloes to the interplay between cosmological merger histories and galaxy formation physics. We perform dark-matter-only, zoomed simulations of two Milky Way-mass hosts and make targeted, controlled changes to their cosmological histories using the genetic modification technique. Populating each history's stellar halo with a semi-empirical, particle-tagging approach then enables a controlled study, with all instances converging to the same large-scale structure, dynamical and stellar mass at $z=0$ as their reference. These related merger scenarios alone generate an extended spread in stellar halo mass fractions (1.5 dex) comparable to the observed population. Largest scatter is achieved by growing late ($z\leq1$) major mergers that spread out existing stars to create massive, in-situ dominated stellar haloes. Increasing a last major merger at $z\sim2$ brings more accreted stars into the inner regions, resulting in smaller scatter in the outskirts which are predominantly built by subsequent minor events. Exploiting the flexibility of our semi-empirical approach, we show that the diversity of stellar halo masses across scenarios is reduced by allowing shallower slopes in the stellar mass--halo mass relation for dwarf galaxies, while it remains conserved when central stars are born with hotter kinematics across cosmic time. The merger-dependent diversity of stellar haloes thus responds distinctly to assumptions in modelling the central and dwarf galaxies respectively, opening exciting prospects to constrain star formation and feedback at different galactic mass-scales with the coming generation of deep, photometric observatories.

Observations of the ICM in the outskirts of the Virgo cluster with Suzaku have found the gas mass fraction in the northern direction to be significantly above the expected level, indicating that there may be a very high level of gas clumping on small scales in this direction. Here, we explore the XMM-Newton data in the outskirts of Virgo, dividing it into a Voronoi tessellation to separate the bulk ICM component from the clumped ICM component. As the nearest galaxy cluster, Virgo's large angular extent allows the spatial scale of the tessellation to be much smaller than has been achieved using the same technique on intermediate redshift clusters, allowing us to probe gas clumping on the scales down to 5$\times$5 kpc. We find that the level of gas clumping in the outskirts to the north is relatively mild, ($\sqrt{C} < 1.1$), suggesting that our point-source detection procedure may have excluded a significant fraction of clumps. While correcting for clumping brings the gas mass fraction at $r_{200}$ into agreement with the universal gas mass fraction, the values outside $r_{200}$ remain significantly above it. This may suggest that non-thermal pressure support in the outskirts to the north is significant, and we find that a non-thermal pressure support at the level of 20 per cent of the total pressure outside $r_{200}$ can explain the high gas mass fraction to the north, which is in agreement with the range expected from simulations.

Supernova (SN) rates serve as an important probe of star-formation models and initial mass functions. Near-infrared seeing-limited ground-based surveys typically discover a factor of 3-10 fewer SNe than predicted from far-infrared (FIR) luminosities owing to sensitivity limitations arising from both a variable point-spread function (PSF) and high dust extinction in the nuclear regions of star-forming galaxies. This inconsistency has potential implications for our understanding of star-formation rates and massive-star evolution, particularly at higher redshifts, where star-forming galaxies are more common. To resolve this inconsistency, a successful SN survey in the local universe must be conducted at longer wavelengths and with a space-based telescope, which has a stable PSF to reduce the necessity for any subtraction algorithms and thus residuals. Here we report on a two-year Spitzer/IRAC 3.6 um survey for dust-extinguished SNe in the nuclear regions of forty luminous infrared galaxies (LIRGs) within 200 Mpc. The asymmetric Spitzer PSF results in worse than expected subtraction residuals when implementing standard template subtraction. Forward-modeling techniques improve our sensitivity by ~1.5 magnitudes. We report the detection of 9 SNe, five of which were not discovered by optical surveys. After adjusting our predicted rates to account for the sensitivity of our survey, we find that the number of detections is consistent with the models. While this search is nonetheless hampered by a difficult-to-model PSF and the relatively poor resolution of Spitzer, it will benefit from future missions, such as Roman Space Telescope and JWST, with higher resolution and more symmetric PSFs.

Gravitational waves (GWs) provide unobscured insight into the birthplace of neutron stars (NSs) and black holes in core-collapse supernovae (CCSNe). The nuclear equation of state (EOS) describing these dense environments is yet uncertain, and variations in its prescription affect the proto-neutron star (PNS) and the post-bounce dynamics in CCSNe simulations, subsequently impacting the GW emission. We perform axisymmetric simulations of CCSNe with Skyrme-type EOSs to study how the GW signal and PNS convection zone are impacted by two experimentally accessible EOS parameters, (1) the effective mass of nucleons, $m^\star$, which is crucial in setting the thermal dependence of the EOS, and (2) the isoscalar incompressibility modulus, $K_{\rm{sat}}$. While $K_{\rm{sat}}$ shows little impact, the peak frequency of the GWs has a strong effective mass dependence due to faster contraction of the PNS for higher values of $m^\star$ owing to a decreased thermal pressure. These more compact PNSs also exhibit more neutrino heating which drives earlier explosions and correlates with the GW amplitude via accretion plumes striking the PNS, exciting the oscillations. We investigate the spatial origin of the GWs and show the agreement between a frequency-radial distribution of the GW emission and a perturbation analysis. We do not rule out overshoot from below via PNS convection as another moderately strong excitation mechanism in our simulations. We also study the combined effect of effective mass and rotation. In all our simulations we find evidence for a power gap near $\sim$1250 Hz, we investigate its origin and report its EOS dependence.

We investigate the luminosity suppression and its effect on the mass-radius relation as well as cooling evolution of highly magnetised white dwarfs. Based on the effect of magnetic field relative to gravitational energy, we suitably modify our treatment of the radiative opacity, magnetostatic equilibrium and degenerate core equation of state to obtain the structural properties of these stars. Although the Chandrasekhar mass limit is retained in the absence of magnetic field and irrespective of the luminosity, strong central fields of about $10^{14}\, {\rm G}$ can yield super-Chandrasekhar white dwarfs with masses up to $1.9\, M_{\odot}$. Smaller white dwarfs tend to remain super-Chandrasekhar for sufficiently strong central fields even when their luminosity is significantly suppressed to $10^{-16}\ L_{\odot}$. Owing to the cooling evolution and simultaneous field decay over $10\ {\rm Gyr}$, the limiting masses of small magnetised white dwarfs can fall to $1.5\ M_{\odot}$ over time. However the majority of these systems still remain practically hidden throughout their cooling evolution because of their high fields and correspondingly low luminosities. Utilising the stellar evolution code $\textit{STARS}$, we obtain close agreement with the analytical mass limit estimates and this suggests that our analytical formalism is physically motivated. Our results argue that super-Chandrasekhar white dwarfs born due to strong field effects may not remain so for long. This explains their apparent scarcity in addition to making them hard to detect because of their suppressed luminosities.

We present a population of 20 radio-luminous supernovae (SNe) with emission reaching $L_{\nu}{\sim}10^{26}-10^{29}\rm{erg s^{-1} Hz^{-1}}$ in the first epoch of the Very Large Array Sky Survey (VLASS) at $2-4$ GHz. Our sample includes one long Gamma-Ray Burst, SN 2017iuk/GRB171205A, and 19 core-collapse SNe detected at $\approx (1-60)$ years after explosion. No thermonuclear explosion shows evidence for bright radio emission, and hydrogen-poor progenitors dominate the sub-sample of core-collapse events with spectroscopic classification at the time of explosion (73%). We interpret these findings into the context of the expected radio emission from the forward shock interaction with the circumstellar medium (CSM). We conclude that these observations require a departure from the single wind-like density profile (i.e., $\rho_{\rm{CSM}}\propto r^{-2}$) that is expected around massive stars and/or a departure from a spherical Newtonian shock. Viable alternatives include the shock interaction with a detached, dense shell of CSM formed by a large effective progenitor mass-loss rate $\dot M \sim (10^{-4}-10^{-1})$ M$_{\odot}$ yr$^{-1}$ (for an assumed wind velocity of $1000\,\rm{km\,s^{-1}}$); emission from an off-axis relativistic jet entering our line of sight; or the emergence of emission from a newly-born pulsar-wind nebula. The relativistic SN\,2012ap that is detected 5.7 and 8.5 years after explosion with $L_{\nu}{\sim}10^{28}$ erg s$^{-1}$ Hz$^{-1}$ might constitute the first detections of an off-axis jet+cocoon system in a massive star. Future multi-wavelength observations will distinguish among these scenarios. Our VLASS source catalogs, which were used to perform the VLASS cross matching, are publicly available at https://doi.org/10.5281/zenodo.4895112.

We present a suite of high-resolution cosmological zoom-in simulations to $z=4$ of a $10^{12}\,{\rm M}_{\odot}$ halo at $z=0$, obtained using seven contemporary astrophysical simulation codes widely used in the numerical galaxy formation community. Physics prescriptions for gas cooling, heating, and star formation, are similar to the ones used in our previous {\it AGORA} disk comparison but now account for the effects of cosmological processes. In this work, we introduce the most careful comparison yet of galaxy formation simulations run by different code groups, together with a series of four calibration steps each of which is designed to reduce the number of tunable simulation parameters adopted in the final run. After all the participating code groups successfully completed the calibration steps, we reach a suite of cosmological simulations with similar mass assembly histories down to $z=4$. With numerical accuracy that resolves the internal structure of a target halo, we find that the codes overall agree well with one another in e.g., gas and stellar properties, but also show differences in e.g., circumgalactic medium properties. We argue that, if adequately tested in accordance with our proposed calibration steps and common parameters, the results of high-resolution cosmological zoom-in simulations can be robust and reproducible. New code groups are invited to join this comparison by generating equivalent models by adopting the common initial conditions, the common easy-to-implement physics package, and the proposed calibration steps. Further analyses of the simulations presented here will be in forthcoming reports from our Collaboration.

Axisymmetric disks of eccentric orbits in near-Keplerian potentials are unstable to an out-of-plane buckling. Recently, Zderic et al. (2020) showed that an idealized disk saturates to a lopsided mode. Here we show that this apsidal clustering also occurs in a primordial scattered disk in the outer solar system which includes the orbit-averaged gravitational influence of the giant planets. We explain the dynamics using Lynden-Bell (1979)'s mechanism for bar formation in galaxies. We also show surface density and line of sight velocity plots at different times during the instability, highlighting the formation of concentric circles and spiral arms in velocity space.

It is still poorly constrained how the densest phase of the interstellar medium varies across galactic environment. A large observing time is required to recover significant emission from dense molecular gas at high spatial resolution, and to cover a large dynamic range of extragalactic disc environments. We present new NOrthern Extended Millimeter Array (NOEMA) observations of a range of high critical density molecular tracers (HCN, HNC, HCO+) and CO isotopologues (13CO, C18O) towards the nearby (11.3 Mpc), strongly barred galaxy NGC 3627. These observations represent the current highest angular resolution (1.85"; 100 pc) map of dense gas tracers across a disc of a nearby spiral galaxy, which we use here to assess the properties of the dense molecular gas, and their variation as a function of galactocentric radius, molecular gas, and star formation. We find that the HCN(1-0)/CO(2-1) integrated intensity ratio does not correlate with the amount of recent star formation. Instead, the HCN(1-0)/CO(2-1) ratio depends on the galactic environment, with differences between the galaxy centre, bar, and bar end regions. The dense gas in the central 600 pc appears to produce stars less efficiently despite containing a higher fraction of dense molecular gas than the bar ends where the star formation is enhanced. In assessing the dynamics of the dense gas, we find the HCN(1-0) and HCO+(1-0) emission lines showing multiple components towards regions in the bar ends that correspond to previously identified features in CO emission. These features are co-spatial with peaks of Halpha emission, which highlights that the complex dynamics of this bar end region could be linked to local enhancements in the star formation.

We analyze VLT/MUSE observations of 37 radio galaxies from the Third Cambridge catalogue (3C) with redshift $<$0.3 searching for nuclear outflows of ionized gas. These observations are part of the MURALES project (a MUse RAdio Loud Emission line Snapshot survey), whose main goal is to explore the feedback process in the most powerful radio-loud AGN. We applied a nonparametric analysis to the [O~III] $\lambda$5007 emission line, whose asymmetries and high-velocity wings reveal signatures of outflows. We find evidence of nuclear outflows in 21 sources, with velocities between $\sim$400 - 1000 km s$^{-1}$, outflowing masses of $\sim 10^5-10^7$ M$_\odot$, and a kinetic energy in the range $\sim 10^{53} - 10^{56}$ erg. In addition, evidence for extended outflows is found in the 2D gas velocity maps of 13 sources of the subclasses of high-excitation (HEG) and broad-line (BLO) radio galaxies, with sizes between 0.4 and 20 kpc. We estimate a mass outflow rate in the range 0.4 - 30 M$_\odot$ yr$^{-1}$ and an energy deposition rate of ${\dot E}_{kin} \sim 10^{42}-10^{45} $ erg s$^{-1}$. Comparing the jet power, the nuclear luminosity of the active galactic nucleus, and the outflow kinetic energy rate, we find that outflows of HEGs and BLOs are likely radiatively powered, while jets likely only play a dominant role in galaxies with low excitation. The low loading factors we measured suggest that these outflows are driven by momentum and not by energy. Based on the gas masses, velocities, and energetics involved, we conclude that the observed ionized outflows have a limited effect on the gas content or the star formation in the host. In order to obtain a complete view of the feedback process, observations exploring the complex multiphase structure of outflows are required.

Two of the most rapidly growing observables in cosmology and astrophysics are gravitational waves (GW) and the neutral hydrogen (HI) distribution. In this work, we investigate the cross-correlation between resolved gravitational wave detections and HI signal from intensity mapping (IM) experiments. By using a tomographic approach with angular power spectra, including all projection effects, we explore possible applications of the combination of the Einstein Telescope and the SKAO intensity mapping surveys. We focus on three main topics: \textit{(i)} statistical inference of the observed redshift distribution of GWs; \textit{(ii)} constraints on dynamical dark energy models as an example of cosmological studies; \textit{(iii)} determination of the nature of the progenitors of merging binary black holes, distinguishing between primordial and astrophysical origin. Our results show that: \textit{(i)} the GW redshift distribution can be calibrated with good accuracy at low redshifts, without any assumptions on cosmology or astrophysics, potentially providing a way to probe astrophysical and cosmological models; \textit{(ii)} the constrains on the dynamical dark energy parameters are competitive with IM-only experiments, in a complementary way and potentially with less systematics; \textit{(iii)} it will be possible to detect a relatively small abundance of primordial black holes within the gravitational waves from resolved mergers. Our results extend towards $\mathrm{GW \times IM}$ the promising field of multi-tracing cosmology and astrophysics, which has the major advantage of allowing scientific investigations in ways that would not be possible by looking at single observables separately.

In a variety of mechanisms generating primordial black holes, each black hole is expected to form along with a surrounding underdense region that roughly compensates the black hole mass. This region will propagate outwards and expand as a shell at the speed of sound in the homogeneous background. Dissipation of the shell due to Silk damping could lead to detectable $\mu$-distortion in the CMB spectrum. While the current bound on the average $\mu$-distortion is $\left| \bar{\mu}\right|\lesssim10^{-4}$, the standard $\Lambda$CDM model predicts $\left| \bar{\mu}\right|\sim10^{-8}$, which could possibly be detected in future missions. It is shown in this work that the non-observation of $\bar{\mu}$ beyond $\Lambda$CDM can place a new upper bound on the density of supermassive primordial black holes within the mass range $10^{6}M_{\odot}\lesssim M\lesssim10^{15}M_{\odot}$. Furthermore, black holes with initial mass $M\gtrsim10^{12}M_{\odot}$ could leave a pointlike distortion with $\mu\gtrsim10^{-8}$ at an angular scale $\sim 1^{\circ}$ in CMB, and its non-observation would impose an even more stringent bound on the population of these stupendously large primordial black holes.

Propagation effects in the interstellar medium and intrinsic profile changes can cause variability in the timing of pulsars, which limits the accuracy of fundamental science done via pulsar timing. One of the best timing pulsars, PSR J1713+0747, has gone through two ``dip'' events in its dispersion measure time series. If these events reflect real changes in electron column density, they should lead to multiple imaging. We show that the events are are well-fit by an underdense corrugated sheet model, and look for associated variability in the pulse profile using principal component analysis. We find that there are transient pulse profile variations, but they vary in concert with the dispersion measure, unlike what is expected from lensing due to a corrugated sheet. The change is consistent in shape across profiles from both the Greenbank and Arecibo radio observatories, and its amplitude appears to be achromatic across the 820 MHz, 1.4 GHz, and 2.3 GHz bands, again unlike expected from interference between lensed images. This result is puzzling. We note that some of the predicted lensing effects would need higher time and frequency resolution data than used in this analysis. Future events appear likely, and storing baseband data or keeping multiple time-frequency resolutions will allow more in-depth study of propagation effects and hence improvements to pulsar timing accuracy.

We study the conditions required for the production of the synchrotron maser emission downstream of a relativistic shock. We show that for weakly magnetized shocks, synchrotron maser emission can be generated at frequencies significantly exceeding the relativistic gyrofrequency. This high-frequency maser emission seems to be the most suitable for interpreting peculiar GHz radio sources. To illustrate this, we consider a magnetar flare model for FRBs. Our analysis shows that the maser emission is radiated away from the central magnetar, which guarantees a short duration of bursts independently of the shock wave radius. If FRBs are produced by the high-frequency maser emission then one can significantly relax the requirements for several key parameters: the magnetic field strength at the production site, luminosity of the flare, and the production site bulk Lorentz factor. To check the feasibility of this model, we study the statistical relation between powerful magnetar flares and the rate of FRBs. The expected ratio is derived by convoluting the redshift-dependent magnetar density with their flare luminosity function above the energy limit determined by the FRB detection threshold. We obtain that only a small fraction, \(\sim10^{-5}\), of powerful magnetar flares trigger FRBs. This ratio agrees surprisingly well with our estimates: we obtained that \(10\%\) of magnetars should be in the evolutionary phase suitable for the production of FRBs, and only \(10^{-4}\) of all flares are expected to be weakly magnetized, which is a necessary condition for the high-frequency maser emission.

We report the discovery of a UHE gamma-ray source, LHAASO J2108+5157, by analyzing the LHAASO-KM2A data of 308.33 live days. Significant excess of gamma-ray induced showers is observed in both energy bands of 25-100 TeV and $\gt$100 TeV with 9.5 sigma and 8.5 sigma, respectively. This source is not significantly favored as an extensive source with the angular extension smaller than the point-spread function of KM2A. The measured energy spectrum from 20 to 200 TeV can be approximately described by a power-law function with an index of -2.83$\pm$ 0.18stat. A harder spectrum is demanded at lower energies considering the flux upper limit set by Fermi-LAT observations. The position of the gamma-ray emission is correlated with a giant molecular cloud, which favors a hadronic origin. No obvious counterparts have been found, deeper multiwavelength observations will help to shed new light on this intriguing UHE source.

We perform for the first time a 3D hydrodynamics simulation of the evolution of the last minutes pre-collapse of the oxygen shell of a fast-rotating massive star. This star has an initial mass of 38 M$_\odot$, a metallicity of $\sim$1/50 Z$_\odot$, an initial rotational velocity of 600 km s$^{-1}$, and experiences chemically homogeneous evolution. It has a silicon- and oxygen-rich (Si/O) convective layer at (4.7-17)$\times 10^{8}$ cm, where oxygen-shell burning takes place. The power spectrum analysis of the turbulent velocity indicates the dominance of the large-scale mode ($\ell \sim 3$), which has also been seen in non-rotating stars that have a wide Si/O layer. Spiral arm structures of density and silicon-enriched material produced by oxygen-shell burning appear in the equatorial plane of the Si/O shell. Non-axisymmetric, large-scale ($m \le 3$) modes are dominant in these structures. The spiral arm structures have not been identified in previous non-rotating 3D pre-supernova models. Governed by such a convection pattern, the angle-averaged specific angular momentum becomes constant in the Si/O convective layer, which is not considered in spherically symmetrical stellar evolution models. Such spiral arms and constant specific angular momentum might affect the ensuing explosion or implosion of the star.

30 Doradus C is a superbubble which emits the brightest nonthermal X- and TeV gamma-rays in the Local Group. In order to explore detailed connection between the high energy radiation and the interstellar medium, we have carried out new CO and HI observations using the Atacama Large Millimeter$/$Submillimeter Array (ALMA), Atacama Submillimeter Telescope Experiment, and the Australia Telescope Compact Array with resolutions of up to 3 pc. The ALMA data of $^{12}$CO($J$ = 1-0) emission revealed 23 molecular clouds with the typical diameters of $\sim$6-12 pc and masses of $\sim$600-10000 $M_{\odot}$. The comparison with the X-rays of $XMM$-$Newton$ at $\sim$3 pc resolution shows that X-rays are enhanced toward these clouds. The CO data were combined with the HI to estimate the total interstellar protons. Comparison of the interstellar proton column density and the X-rays revealed that the X-rays are enhanced with the total proton. These are most likely due to the shock-cloud interaction modeled by the magnetohydrodynamical simulations (Inoue et al. 2012, ApJ, 744, 71). Further, we note a trend that the X-ray photon index varies with distance from the center of the high-mass star cluster, suggesting that the cosmic-ray electrons are accelerated by one or multiple supernovae in the cluster. Based on these results we discuss the role of the interstellar medium in cosmic-ray particle acceleration.

Celestially, Positronium (Ps), has only been observed through gamma-ray emission produced by its annihilation. However, in its triplet state, a Ps atom has a mean lifetime long enough for electronic transitions to occur between quantum states. This produces a recombination spectrum observable in principle at near IR wavelengths, where angular resolution greatly exceeding that of the gamma-ray observations is possible. However, the background in the NIR is dominated by extremely bright atmospheric hydroxyl (OH) emission lines. In this paper we present the design of a diffraction-limited spectroscopic system using novel photonic components - a photonic lantern, OH Fiber Bragg Grating filters, and a photonic TIGER 2-dimensional pseudo-slit - to observe the Ps Balmer alpha line at 1.3122 microns for the first time.

The Davis-Chandrasekhar-Fermi (DCF) method is widely used to indirectly estimate the magnetic field strength from the plane-of-sky field orientation. In this work, we present a set of 3D MHD simulations and synthetic polarization images using radiative transfer of clustered massive star-forming regions. We apply the DCF method on the synthetic polarization maps to investigate its reliability in high-density molecular clumps and dense cores where self-gravity is significant. We investigate the validity of the assumptions of the DCF method step by step and compare the model and estimated field strength to derive the correction factors for the estimated uniform and total (rms) magnetic field strength at clump and core scales. The correction factors in different situations are catalogued. We find the DCF method works well in strong field cases. However, the magnetic field strength in weak field cases could be significantly overestimated by the DCF method when the turbulent magnetic energy is smaller than the turbulent kinetic energy. We investigate the accuracy of the angular dispersion function (ADF, a modified DCF method) method on the effects that may affect the measured angular dispersion and find that the ADF method correctly accounts for the ordered field structure, the beam-smoothing, and the interferometric filtering, but may not be applicable to account for the signal integration along the line of sight in most cases. Our results suggest that the DCF methods should be avoided to be applied below $\sim$0.1 pc scales if the effect of line-of-sight signal integration is not properly addressed.

We present new generation mechanisms of magnetic fields in supernova remnant shocks propagating to partially ionized plasmas in the early universe. Upstream plasmas are dissipated at the collisionless shock, but hydrogen atoms are not dissipated because they do not interact with electromagnetic fields. After the hydrogen atoms are ionized in the shock downstream region, they become cold proton beams that induce the electron return current. The injection of the beam protons can be interpreted as an external force acting on the downstream proton plasma. We show that the effective external force and the electron return current can generate magnetic fields without any seed magnetic fields. The magnetic field strength is estimated to be $B\sim 10^{-14}-10^{-11}~{\rm G}$, where the characteristic lengthscale is the mean free path of charge exchange, $\sim 10^{15}~{\rm cm}$. Since protons are marginally magnetized by the generated magnetic field in the downstream region, the magnetic field could be amplified to larger values and stretched to larger scales by turbulent dynamo and expansion.

Measuring weak lensing cosmic magnification signal is very challenging due to the overwhelming intrinsic clustering in the observed galaxy distribution. In this paper, we modify the Internal Linear Combination (ILC) method to reconstruct the lensing signal with an extra constraint to suppress the intrinsic clustering. To quantify the performance, we construct a realistic galaxy catalogue for the LSST-like photometric survey, covering 20000 $\deg^{2}$ with mean source redshift at $z_{s}\sim 1$. We find that the reconstruction performance depends on the width of the photo-z bin we choose. Due to the correlation between the lensing signal and the source galaxy distribution, the derived signal has smaller systematic bias but larger statistical uncertainty for a narrower photo-z bin. We conclude that the lensing signal reconstruction with the Modified ILC method is unbiased with a statistical uncertainty $<5\%$ for bin width $\Delta z^{P} = 0.2$.

The analysis of stellar spectra depends upon the effective temperature (Teff) and the surface gravity (log g). However, the estimation of log g with high accuracy is challenging. A classical approach is to search for log g that satisfies the ionization balance, i.e., the abundances from neutral and ionized metallic lines being in agreement. We propose a method of using empirical relations between Teff, log g and line-depth ratios, for which we meticulously select pairs of FeI and FeII lines and pairs of CaI and CaII lines. Based on YJ-band (0.97-1.32 micron) high-resolution spectra of 42 FGK stars (dwarfs to supergiants), we selected five FeI-FeII and four CaI-CaII line pairs together with 13 FeI-FeI pairs (for estimating Teff, and derived the empirical relations. Using such relations does not require complex numerical models and tools for estimating chemical abundances. The relations we present allows one to derive Teff and log g with a precision of around 50 K and 0.2 dex, respectively, but the achievable accuracy depends on the accuracy of the calibrators' stellar parameters. It is essential to revise the calibration by observing stars with accurate stellar parameters available, e.g., stars with asteroseismic log g and stars analyzed with complete stellar models taking into account the effects of non-local thermodynamic equilibrium and convection. In addition, the calibrators we used have a limited metallicity range, [Fe/H] between -0.2 and +0.2 dex, and our relations need to be tested and re-calibrated based on a calibrating dataset for a wider range of metallicity.

Long-term gamma-ray variability of a non-blazar Active Galactic Nucleus (AGN) PKS 0521-36 is investigated by using Fermi-LAT pass 8 data covering from 2008 August to 2021 March. The results show that the histogram of the gamma-ray fluxes follows a log-normal distribution. Interestingly, in the analysis of about 5.8-year (from MJD 56317 to 58447) LAT data between two outbursts (occurring during 2012 October and 2019 May respectively), a quasi-periodic oscillation (QPO) with a period of about 1.1 years (about 5 sigma of significance) is found in the Lomb-Scargle Periodogram (LSP), the Weighted Wavelet Z-transform (WWZ) and the REDFIT results. This quasi-periodic signal also appears in the results of Gaussian process modeling the light curve. Therefore, the robustness of the QPO is examined by four different methods. This is the first gamma-ray QPO found in a mildly beamed jet. Our results imply that the gamma-ray outbursts play an important role in the formation of the gamma-ray QPO.

The emergence and development of a power-law tail (PLT) at the high-density end of the observed column-density distribution is thought to be indicative for advanced evolution of star-forming molecular clouds. As shown from many numerical simulations, it corresponds to a morphologically analogous evolution of the mass-density distribution (\rhopdf). The latter may display also a second, shallower PLT at the stage of collapse of the first formed protostellar cores. It is difficult to estimate the parameters of a possible second PLT due to resolution constraints. To address the issue, we extend the method for the extraction of single PLTs from arbitrary density distributions suggested by Veltchev et al.(2019) to detect a second PLT. The technique is elaborated through tests on an analytic \rhopdf{} and applied to a set of hydrodynamical high-resolution simulations of isothermal self-gravitating clouds. In all but one case two PLTs were detected -- the first slope is always steeper and the second one is typically $\partial \ln V /\ln \rho \sim -1$. These results are in a good agreement with numerical and theoretical works and do suggest that the technique extracts correctly double PLTs from smooth PDFs.

The Spectrometer/Telescope for Imaging X-rays (STIX) is the HXR instrument onboard Solar Orbiter designed to observe solar flares over a broad range of flare sizes, between 4-150 keV. We report the first STIX observations of microflares recorded during the instrument commissioning phase in order to investigate the STIX performance at its detection limit. This first result paper focuses on the temporal and spectral evolution of STIX microflares occuring in the AR12765 in June 2020, and compares the STIX measurements with GOES/XRS, SDO/AIA, and Hinode/XRT. For the observed microflares of the GOES A and B class, the STIX peak time at lowest energies is located in the impulsive phase of the flares, well before the GOES peak time. Such a behavior can either be explained by the higher sensitivity of STIX to higher temperatures compared to GOES, or due to the existence of a nonthermal component reaching down to low energies. The interpretation is inconclusive due to limited counting statistics for all but the largest flare in our sample. For this largest flare, the low-energy peak time is clearly due to thermal emission, and the nonthermal component seen at higher energies occurs even earlier. This suggests that the classic thermal explanation might also be favored for the majority of the smaller flares. In combination with EUV and SXR observations, STIX corroborates earlier findings that an isothermal assumption is of limited validity. Future diagnostic efforts should focus on multi-wavelength studies to derive differential emission measure distributions over a wide range of temperatures to accurately describe the energetics of solar flares. Commissioning observations confirm that STIX is working as designed. As a rule of thumb, STIX detects flares as small as the GOES A class. For flares above the GOES B class, detailed spectral and imaging analyses can be performed.

Detecting and understanding rotation in stellar interiors is nowadays one of the unsolved problems in stellar physics. Asteroseismology has been able to provide insights on rotation for the Sun, solar-like stars, and compact objects like white dwarfs. However, this is still very difficult for intermediate-mass stars. These stars are moderate-to-rapid rotators. Rotation splits and shifts the oscillation modes, which makes the oscillation spectrum more complex and harder to interpret. Here we study the oscillation patterns of a sample of benchmark $\delta$~Sct stars belonging to eclipsing binary systems with the objective to find the frequency spacing related to the rotational splitting ($\delta r$). For this task, we combine three techniques: the Fourier transform, the autocorrelation function, and the histogram of frequency differences. The last two showed a similar behaviour. For most of the stars, it was necessary to determine the large separation ($\Delta\nu$) prior to spot $\delta r$. This is the first time we may clearly state that one of the periodicities present in the p~modes oscillation spectra of $\delta$~Sct stars corresponds to the rotational splitting. This is true independently of the stellar rotation rate. These promising results pave the way to find a robust methodology to determine rotational splittings from the oscillation spectra of $\delta$~Sct stars and, thus, understanding the rotational profile of intermediate-mass pulsating stars.

Context. The discovery of an extrasolar planet with an ocean has crucial importance in the search for life beyond Earth. The polarimetric detection of specularly reflected light from a smooth liquid surface is anticipated theoretically, though the polarimetric signature of Earth's ocean has not yet been conclusively detected in disk-integrated planetary light. Aims. We aim to detect and measure the polarimetric signature of the Earth's ocean. Methods. We conducted near-infrared polarimetry for lunar Earthshine and collected data for 32 nights with a variety of ocean fractions in the Earthshine contributing region. Results. A clear positive correlation was revealed between the polarization degree and ocean fraction. We found hourly variations in polarization in accordance with rotational transition of the ocean fraction. The ratios of the variation to the typical polarization degree were as large as ~0.2-1.4. Conclusions. Our observations provide plausible evidence of the polarimetric signature attributed to Earth's ocean. Near-infrared polarimetry may be considered a prospective technique for the search for exoplanetary oceans.

We aim to characterize a diverse selection of dense, potentially star-forming cores, clumps, and clouds within the Milky Way in terms of their dust emission and SF activity. We studied 53 fields that have been observed in the JCMT SCUBA-2 continuum survey SCOPE and have been mapped with Herschel. We estimated dust properties by fitting Herschel observations with modified blackbody functions, studied the relationship between dust temperature and dust opacity spectral index $\beta$, and estimated column densities. We extracted clumps from the SCUBA-2 850 $\mu$m maps with the FellWalker algorithm and examined their masses and sizes. Clumps are associated with young stellar objects found in several catalogs. We estimated the gravitational stability of the clumps with virial analysis. The clumps are categorized as unbound starless, prestellar, or protostellar. We find 529 dense clumps, typically with high column densities from (0.3-4.8)$\times 10^{22}$ cm$^{-2}$, with a mean of (1.5$\pm$0.04)$\times10^{22}$ cm$^{-2}$, low temperatures ($T\sim $10-20 K), and estimated submillimeter $\beta$ =1.7$\pm$0.1. We detect a slight increase in opacity spectral index toward millimeter wavelengths. Masses of the sources range from 0.04 $M_\odot$ to 4259 $M_\odot$. Mass, linear size, and temperature are correlated with distance. Furthermore, the estimated gravitational stability is dependent on distance, and more distant clumps appear more virially bound. Finally, we present a catalog of properties of the clumps.Our sources present a large array of SF regions, from high-latitude, nearby diffuse clouds to large SF complexes near the Galactic center. Analysis of these regions will continue with the addition of molecular line data, which will allow us to study the densest regions of the clumps in more detail.

We present a new, fast, and easy to use tool for modelling light and radial velocity curves of close eclipsing binaries with built-in methods for solving an inverse problem. The main goal of ELISa (Eclipsing binary Learning and Interactive System) is to provide an acceptable compromise between computational speed and precision during the fitting of light curves and radial velocities of eclipsing binaries. The package is entirely written in the Python programming language in a modular fashion, making it easy to install, modify, and run on various operating systems. ELISa implements Roche geometry and the triangulation process to model a surface of the eclipsing binary components, where the surface parameters of each surface element are treated separately. Surface symmetries and approximations based on the similarity between surface geometries were used to reduce the runtime during light curve calculation significantly. ELISa implements the least square trust region reflective algorithm and Markov-chain Monte Carlo optimisation methods to provide the built-in capability to determine parameters of the binary system from photometric observations and radial velocities. The precision and speed of the light curve generator were evaluated using various benchmarks. We conclude that ELISa maintains an acceptable level of accuracy to analyse data from ground-based and space-based observations, and it provides a significant reduction in computational time compared to the current widely used tools for modelling eclipsing binaries.

We present the first on-sky results of the micro-lens ring tip-tilt (MLR-TT) sensor. This sensor utilizes a 3D printed micro-lens ring feeding six multi-mode fibers to sense misaligned light, allowing centroid reconstruction. A tip-tilt mirror allows the beam to be corrected, increasing the amount of light coupled into a centrally positioned single-mode (science) fiber. The sensor was tested with the iLocater acquisition camera at the Large Binocular Telescope in November 2019. The limit on the maximum achieved root mean square reconstruction accuracy was found to be 0.19 $\lambda$/D in both tip and tilt, of which approximately 50% of the power originates at frequencies below 10 Hz. We show the reconstruction accuracy is highly dependent on the estimated Strehl ratio and simulations support the assumption that residual adaptive optics aberrations are the main limit to the reconstruction accuracy. We conclude that this sensor is ideally suited to remove post-adaptive optics non-common path tip tilt residuals. We discuss the next steps for the concept development, including optimizations of the lens and fiber, tuning of the correction algorithm and selection of optimal science cases.

The structure and evolution of protoplanetary disks (PPDs) are largely governed by disk angular momentum transport, mediated by magnetic fields. In the most observable outer disk, PPD gas dynamics is primarily controlled by ambipolar diffusion as the dominant non-ideal magnetohydrodynamic (MHD) effect. In this work, we study the gas dynamics in outer PPDs by conducting a set of global 3D non-ideal MHD simulations with ambipolar diffusion and net poloidal magnetic flux, using the Athena++ MHD code, with resolution comparable to local simulations. Our simulations demonstrate the co-existence of magnetized disk wind and turbulence driven by the magneto-rotational instability (MRI). While MHD wind dominates disk angular momentum transport, the MRI turbulence also contributes significantly. We observe that magnetic flux spontaneously concentrate into axisymmetric flux sheets, leading to radial variations in turbulence levels, stresses, and accretion rates. Annular substructures arise as a natural consequence of magnetic flux concentration. The flux concentration phenomena show diverse properties with different levels of disk magnetization and ambipolar diffusion. The disk generally loses magnetic flux over time, though flux sheets could prevent the leak of magnetic flux in some cases. Our results demonstrate the ubiquity of disk annular substructures in weakly MRI turbulent outer PPDs, and imply a stochastic nature of disk evolution.

We use the problem of dynamical friction within the periodic cube to illustrate the application of perturbation theory in stellar dynamics, testing its predictions against measurements from $N$-body simulation. Our development is based on the explicitly time-dependent Volterra integral equation for the cube's linear response, which avoids the subtleties encountered in analyses based on complex frequency. We obtain an expression for the self-consistent response of the cube to steady stirring by an external perturber. From this we show how to obtain the familiar Chandrasekhar dynamical friction formula and construct an elementary derivation of the Lenard--Balescu equation for the secular quasilinear evolution of an isolated cube composed of $N$ equal-mass stars. We present an alternative expression for the (real-frequency) van Kampen modes of the cube and show explicitly how to decompose any linear perturbation of the cube into a superposition of such modes.

In this thesis, we discuss some of the applications of cosmological perturbation theory in the late universe. We begin by reviewing the tools used to understand the standard model of cosmology theoretically and to compute its observational consequences, including a detailed exposition of cosmological perturbation theory. We then describe the results in this thesis; we present novel analytical solutions for linear-order gravitational waves or tensor perturbations in a flat Friedmann-Robertson-Walker universe containing two perfect fluids -- radiation and pressureless dust -- and allowing for neutrino anisotropic stress. One of the results applies to any sub-horizon gravitational wave in such a universe. Another result applies to gravitational waves of primordial origin (for example, produced during inflation) and works both before and after they cross the horizon. These results improve on analytical approximations previously set out in the literature. Comparison with numerical solutions shows that both these approximations are accurate to within 1% or better, for a wide range of wave-numbers relevant for cosmology. We present a new and independent approach to computing the relativistic galaxy number counts to second order in cosmological perturbation theory. We also derive analytical expressions for the full second-order relativistic observed redshift, for the angular diameter distance and the volume spanned by a survey. We then compare our result with previous works which compute the general distance-redshift relation, finding that our result is in agreement at linear and leading nonlinear order. Lastly, we briefly study a class of almost scale-invariant Gauss-Bonnet modified gravity theory and derive the Einstein-like field equations to first order in cosmological perturbation theory in longitudinal gauge.

Close-in co-orbital planets (in a 1:1 mean motion resonance) can experience strong tidal interactions with the central star. Here, we develop an analytical model adapted to the study of the tidal evolution of those systems. We use a Hamiltonian version of the constant time-lag tidal model, which extends the Hamiltonian formalism developed for the point-mass case. We show that co-orbital systems undergoing tidal dissipation either favour the Lagrange or the anti-Lagrange configurations, depending on the system parameters. However, for all range of parameters and initial conditions, both configurations become unstable, although the timescale for the destruction of the system can be larger than the lifetime of the star. We provide an easy-to-use criterion to determine if an already known close-in exoplanet may have an undetected co-orbital companion.

Global 21cm experiments aim to measure the sky averaged HI absorption signal from cosmic dawn and the epoch of reionisation. However, antenna chromaticity coupling to bright foregrounds can introduce distortions into the observational data of such experiments. We demonstrate a method for guiding the antenna design of a global experiment through data analysis simulations. This is done by performing simulated observations for a range of inserted 21cm signals, then attempting to identify the signals with a data analysis pipeline. We demonstrate this method on five antennae that were considered as potential designs for the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH); a conical log spiral antenna, an inverted conical sinuous antenna and polygonal-, square- and elliptical-bladed dipoles. We find that the log spiral performs significantly better than the other antennae tested, able to correctly and confidently identify every inserted 21cm signal. In second place is the polygonal dipole antenna, which was only unable to detect signals with both very low amplitudes of 0.05K and low centre frequency of 80MHz. The conical sinuous antenna was found to perform least accurately, only able to detect the highest amplitude 21cm signals, and even then with biases. We also demonstrate that, due to the non-trivial nature of chromatic distortion and the processes of correcting for it, these are not the results that could have been expected superficially from the extent of chromatic variation in each antenna.

Exotic electromagnetic energy injection in the early Universe may alter cosmological recombination, and ultimately cosmic microwave background (CMB) anisotropies. Moreover, if energy injection is inhomogeneous, it may induce a spatially-varying ionization fraction, and non-Gaussianity in the CMB. The observability of these signals, however, is contingent upon how far the injected particles propagate and deposit their energy into the primordial plasma, relative to the characteristic scale of energy injection fluctuations. In this study we inspect the spatial properties of energy deposition and perturbed recombination resulting from an inhomogeneous energy injection of sub-10 MeV photons, relevant to accreting primordial black holes (PBHs). We develop a novel Monte-Carlo radiation transport code accounting for all relevant photon interactions in this energy range, and including secondary electron energy deposition efficiency through a new analytic approximation. For a specified injected photon spectrum, the code outputs an injection-to-deposition Green's function depending on time and distance from the injection point. Combining this output with a linearized solution of the perturbed recombination problem, we derive time- and scale-dependent deposition-to-ionization Green's functions. We apply this general framework to accreting PBHs, whose luminosity is strongly spatially modulated by supersonic relative velocities between cold dark matter and baryons. We find that the resulting spatial fluctuations of the free-electron fraction are of the same magnitude as its mean deviation from standard recombination, from which current CMB power spectra constraints are derived. This work suggests that the sensitivity to accreting PBHs might be substantially improved by propagating these inhomogeneities to CMB anisotropy power spectra and non-Gaussian statistics, which we study in subsequent papers.

The High Latitude Survey of the Nancy Grace Roman Space Telescope is expected to measure the positions and shapes of hundreds of millions of galaxies in an area of 2220 deg$^2$. This survey will provide high-quality weak lensing data with unprecedented systematics control. The Roman Space Telescope will survey the sky in near infrared (NIR) bands using Teledyne H4RG HgCdTe photodiode arrays. These NIR arrays exhibit an effect called persistence: charges that are trapped in the photodiodes during earlier exposures are gradually released into later exposures, leading to contamination of the images and potentially to errors in measured galaxy properties such as fluxes and shapes. In this work, we use image simulations that incorporate the persistence effect to study its impact on galaxy shape measurements and weak lensing signals. No significant spatial correlations are found between the galaxy shape changes induced by persistence. On the scales of interest for weak lensing cosmology, the effect of persistence on the weak lensing correlation function is about two orders of magnitude lower than the Roman Space Telescope additive shear error budget, indicating that the persistence effect is expected to be a subdominant contributor to the systematic error budget for weak lensing with the Roman Space Telescope given its current design.

A new concept for the simultaneous detection of primary and secondary scintillation in time projection chambers is proposed. Its core element is a type of very-thick GEM structure supplied with transparent electrodes and machined from a polyethylene naphthalate plate, a natural wavelength-shifter. Such a device has good prospects for scalability and, by virtue of its genuine optical properties, it can improve on the light collection efficiency, energy threshold and resolution of conventional micropattern gas detectors. This, together with the intrinsic radiopurity of its constituting elements, offers advantages for noble gas and liquid based time projection chambers, used for dark matter searches and neutrino experiments. Production, optical and electrical characterization, and first measurements performed with the new device are reported.

A powerful technique to calculate gravitational radiation from binary systems involves a perturbative expansion: if the masses of the two bodies are very different, the "small" body is treated as a point particle of mass $m_p$ moving in the gravitational field generated by the large mass $M$, and one keeps only linear terms in the small mass ratio $m_p/M$. This technique usually yields finite answers, which are often in good agreement with fully nonlinear numerical relativity results, even when extrapolated to nearly comparable mass ratios. Here we study two situations in which the point-particle approximation yields a divergent result: the instantaneous flux emitted by a small body as it orbits the light ring of a black hole, and the total energy absorbed by the horizon when a small body plunges into a black hole. By integrating the Teukolsky (or Zerilli/Regge-Wheeler) equations in the frequency and time domains we show that both of these quantities diverge. We find that these divergences are an artifact of the point-particle idealization, and are able to interpret and regularize this behavior by introducing a finite size for the point particle. These divergences do not play a role in black-hole imaging, e.g. by the Event Horizon Telescope.

Stochastic Inflation is an important framework for understanding the physics of de Sitter space and the phenomenology of inflation. In the leading approximation, this approach results in a Fokker-Planck equation that calculates the probability distribution for a light scalar field as a function of time. Despite its successes, the quantum field theoretic origins and the range of validity for this equation have remained elusive, and establishing a formalism to systematically incorporate higher order effects has been an area of active study. In this paper, we calculate the next-to-next-to-leading order (NNLO) corrections to Stochastic Inflation using Soft de Sitter Effective Theory (SdSET). In this effective description, Stochastic Inflation manifests as the renormalization group evolution of composite operators. The leading impact of non-Gaussian quantum fluctuations appears at NNLO, which is presented here for the first time; we derive the coefficient of this term from a two-loop anomalous dimension calculation within SdSET. We solve the resulting equation to determine the NNLO equilibrium distribution and the low-lying relaxation eigenvalues. In the process, we must match the UV theory onto SdSET at one-loop order, which provides a non-trivial confirmation that the separation into Wilson-coefficient corrections and contributions to initial conditions persists beyond tree level. Furthermore, these results illustrate how the naive factorization of time and momentum integrals in SdSET no longer holds in the presence of logarithmic divergences. It is these effects that ultimately give rise to the renormalization group flow that yields Stochastic Inflation.

In a strong electromagnetic field, gravitational waves are converted into electromagnetic waves of the same frequency, and vice versa. Here we calculate the scattering and conversion cross sections for a planar wave impinging upon a Reissner-Nordstr\"om black hole in vacuum, using the partial-wave expansion and numerical methods. We show that, at long wavelengths, the conversion cross section matches that computed by Feynman-diagram techniques. At short wavelengths, the essential features are captured by a geometric-optics approximation. We demonstrate that the converted flux can exceed the scattered flux at large scattering angles, for highly-charged black holes. In the short-wavelength regime, the conversion effect may be understood in terms of a phase that accumulates along a ray. We compute the scattering angle for which the converted and scattered fluxes are equal, as a function of charge-to-mass ratio. We show that this scattering angle approaches $90$ degrees in the extremal limit.

Weak-scale secluded sector dark matter can reproduce the observed dark matter relic density with thermal freeze-out within that sector. If nature is supersymmetric, three portals to the visible sector - a gauge portal, a Higgs portal, and a gaugino portal - are present. We present gamma ray spectra relevant for indirect detection of dark matter annihilation in such setups. Since symmetries in the secluded sector can stabilize dark matter, $R$-parity is unnecessary, and we investigate the impact of $R$-parity violation on annihilation spectra. We present limits from the Fermi Large Area Telescope observations of dwarf galaxies and projections for Cherenkov Telescope Array observations of the galactic center. Many of our results are also applicable to generic, non-supersymmetric setups.

Global cosmic strings are predicted in many motivated extensions to the Standard Model of particle physics, with close connections to axion dark matter physics. Recent studies suggest that, although subdominant relative to Goldstone emission, gravitational wave (GW) signals from global strings can be detectable with current and planned GW detectors such as LIGO, LISA, DECIGO/BBO, ET/CE and AEDGE/AION, as well as pulsar timing arrays such as PPTA, NANOGrav and SKA. This work is an extensive, updated study on GWs from a global cosmic string network, taking into account of the most recent developments related to the subject. The main analysis is based on the analytical Velocity-dependent One-Scale (VOS) model calibrated with recent simulation results, which provides a generic protocol for such calculations with details given. We also demonstrate how the GW signal can be influenced with variations to the baseline model: this includes considering the uncertainties of model parameters and the potential deviation from the conventional VOS model prediction (i.e. the scaling behavior) as suggested by some of the recent simulation results. Furthermore, we investigated in detail the effect of a non-standard cosmology (e.g. early matter domination or kination) or new particle species on the GW signals from global strings. We demonstrate that the frequency spectrum of GW background from global cosmic strings can be used to probe the cosmic history prior to the Big Bang nucleosynthesis (BBN) (i.e. the primordial dark age) up to a temperature of $T\sim 10^8$ GeV.

The large imbalance in the neutron and proton densities in very neutron rich systems increases the nuclear symmetry energy so that it governs many aspects of neutron stars and their mergers. Extracting the density dependence of the symmetry energy therefore constitutes an important scientific objective. Many analyses have been limited to extracting values for the symmetry energy, $S_0$, and its ``derivative'', $L$, at saturation density $\rho_0 \approx 2.6 \times 10^{14}~\mathrm{g/cm^3}$ $\approx 0.16~\mathrm{nucleons/fm^{3}}$, resulting in constraints that appear contradictory. We show that most experimental observables actually probe the symmetry energy at densities far from $\rho_0$, making the extracted values of $S_0$ or $L$ imprecise. By focusing on the densities these observables actually probe, we obtain a detailed picture of the density dependence of the symmetry energy from $0.25\rho_0$ to $1.5\rho_0$. From this experimentally derived density functional, we extract $L_{01}=53.1\pm6.1 MeV$ at $\rho \approx 0.10~\mathrm{fm^{-3}}$, a neutron skin thickness for $^{208}Pb$ of $R_{np} =$ $0.23\pm0.04$ fm, a symmetry pressure at saturation density of $P_0=3.2\pm1.2 MeV/fm^3$ and suggests a radius for a 1.4 solar mass neutron star of $13.1\pm0.6$ km.

The discovery of interstellar communication signals is complicated by the presence of radio interference. Consequently, interstellar communication signals are hypothesized to have properties that favor discovery in high levels of local planetary radio interference. A hypothesized type of interstellar signal, delta-t delta-f polarized pulse pairs, has properties that are similar to infrequent elements of random noise, while dissimilar from many types of known radio interference. Discovery of delta-t delta-f polarized pulse pairs is aided by the use of interference-filtered receiver systems that are designed to indicate anomalous presence of delta-t delta-f polarized pulse pairs, when pointing a radio telescope to celestial coordinates of a hypothetical transmitter. Observations reported in previous work (ref. arXiv:2105.03727) indicate an anomalous celestial pointing direction having coordinates 5.25 +- 0.15 hours Right Ascension and -7.6 +- 1 degrees Declination. Augmented interference reduction mechanisms used in the current work are described, together with reports of follow-up radio telescope beam transit measurements during 40 days. Conclusions and further work are proposed.

We study how inhomogeneities of the cosmological fluid fields backreact on the homogeneous part of energy density and how they modify the Friedmann equations. In general, backreaction requires to go beyond the pressureless ideal fluid approximation, and this can lead to a reduced growth of cosmological large scale structure. Since observational evidence favours evolution close to the standard growing mode in the linear regime, we focus on two-component fluids in which the non-ideal fluid is gravitationally coupled to cold dark matter and in which a standard growing mode persists. This is realized, e.g. for a baryonic fluid coupled to cold dark matter. We calculate the backreaction for this case and for a wide range of other two-fluid models. Here the effect is either suppressed because the non-ideal matter properties are numerically too small, or because they lead to a too stringent UV cut-off of the integral over the power spectrum that determines backreaction. We discuss then matter field backreaction from a broader perspective and generalize the formalism such that also far-from-equilibrium scenarios relevant to late cosmological times and non-linear scales can be addressed in the future.

Shocks in collisionless plasmas, such as supernovae shocks and shocks driven by coronal mass ejections, are known to be a primary source of energetic particles. Due to their different mass per charge ratio, the interaction of heavy ions with the shock layer differs from that of protons, and injection of these ions into acceleration processes is a challenge. Here we show the first direct observational evidence of magnetic reflection of alpha particles from a high Mach number quasi-perpendicular shock using in-situ spacecraft measurements. The intense magnetic amplification at the shock front associated with nonstationarity modulates the trajectory of alpha particles, some of which travel back upstream as they gyrate in the enhanced magnetic field and experience further acceleration in the upstream region. Our results in particular highlight the important role of high magnetic amplification in seeding heavy ions into the energization processes at nonstationary reforming shocks.

It has been a half-decade since the first direct detection of gravitational waves, which signifies the coming of the era of the gravitational-wave astronomy and gravitational-wave cosmology. The increasing number of the gravitational-wave events has revealed the promising capability of constraining various aspects of cosmology, astronomy, and gravity. Due to the limited space in this review article, we will briefly summarize the recent progress over the past five years, but with a special focus on some of our own works for the Key Project ``Physics associated with the gravitational waves'' supported by the National Natural Science Foundation of China.