Papers for Thursday, Oct 13 2022

We report on results from a one-year soft X-ray observing campaign of the ultraluminous X-ray pulsar NGC 300 ULX-1 by the Neutron star Interior Composition Explorer (NICER) during 2018--2019. Our analysis also made use of data from Swift/XRT and XMM-Newton in order to model and remove contamination from the nearby eclipsing X-ray binary NGC 300 X-1. We constructed and fitted a series of 5-day averaged NICER spectra of NGC 300 ULX-1 in the 0.4--4.0 keV range to evaluate the long-term spectral evolution of the source, and found that an absorbed power-law model provided the best fit overall. Over the course of our observations, the source flux (0.4--4.0 keV; absorbed) dimmed from $2\times10^{-12}$ to below $10^{-13}{\rm\,erg\,s^{-1}\,cm^{-2}}$ and the spectrum softened, with the photon index going from $\Gamma\approx1.6$ to $\Gamma\approx2.6$. We interpret the spectral softening as reprocessed emission from the accretion disk edge coming into view while the pulsar was obscured by the possibly precessing disk. Some spectral fits were significantly improved by the inclusion of a disk blackbody component, and we surmise that this could be due to the pulsar emerging in between obscuration episodes by partial covering absorbers. We posit that we observed a low-flux state of the system (due to line-of-sight absorption) punctuated by the occasional appearance of the pulsar, indicating short-term source variability nested in longer-term accretion disk precession timescales.

Recent work uncovered features in the phase space of the Milky Way's stellar halo which may be attributed to the last major merger. When stellar material from a satellite is accreted onto its host, it phase mixes and appears finely substructured in phase space. For a high-eccentricity merger, this substructure most clearly manifests as numerous wrapping chevrons in $(v_r, r)$ space, corresponding to stripes in $(E, \theta_r)$ space. We introduce the idea of using this substructure as an alternative subhalo detector to cold stellar streams. We simulate an N-body merger akin to the GSE and assess the impact of subhaloes on these chevrons. We examine how their deformation depends on the mass, pericentre, and number of subhaloes. To quantify the impact of perturbers, we utilise the appearance of chevrons in $(E, \theta_r)$ space to introduce a new quantity -- the ironing parameter. We show that: (1) a single flyby of a massive ($\sim 10^{10}$ M$_{\odot}$) subhalo with pericentre comparable to, or within, the shell's apocentre smooths out the substructure, (2) a single flyby of a low mass ($\lesssim 10^8$ M$_{\odot}$) has negligible effect, (3) multiple flybys of subhalos derived from a subhalo mass function between $10^7-10^{10}$ M$_{\odot}$ cause significant damage if deep within the potential, (4) the effects of known perturbers (e.g. Sagittarius) should be detectable and offer constraints on their initial mass. The sensitivity to the populations of subhaloes suggests that we should be able to place an upper limit on the Milky Way's subhalo mass function.

Gas plays an important role in many processes in galaxy formation and evolution, but quantifying the importance of gas has been hindered by the challenge to measure gas masses for large samples of galaxies. Datasets of direct atomic and molecular gas measurements are sufficient to establish simple scaling relations, but often not large enough to quantify three-parameter relations, or second order dependencies. As an alternative approach, we derive here indirect cold gas measurements from optical emission lines using photoionization models for galaxies in the SDSS main galaxy sample and the PHANGS-MUSE survey. We calibrate the gas surface density measurements using xCOLD GASS and PHANGS-ALMA molecular gas measurements to ensure our measurements are reliable. We demonstrate the importance of taking into account the scale-dependence of the relation between optical depth ($\tau_V$) and gas surface density ($\Sigma_{gas}$) and provide a general prescription to estimate $\Sigma_{gas}$ from $\tau_V$, metallicity and the dust-to-metal ratio, at any arbitrary physical resolution. To demonstrate that the indirect cold gas masses are accurate enough to quantify the role of gas in galaxy evolution, we study the mass-metallicity relation (MZR) of SDSS galaxies and show that as a third parameter, gas mass is better than SFR at reducing the scatter of the relation, as predicted by models and simulations.

We characterize Galactic dust filaments by correlating BICEP/Keck and Planck data with polarization templates based on neutral hydrogen (HI) observations. Dust polarization is important for both our understanding of astrophysical processes in the interstellar medium (ISM) and the search for primordial gravitational waves in the cosmic microwave background (CMB). In the diffuse ISM, HI is strongly correlated with the dust and partly organized into filaments that are aligned with the local magnetic field. We analyze the deep BICEP/Keck data at 95, 150, and 220 GHz, over the low-column-density region of sky where BICEP/Keck has set the best limits on inflationary theories. We separate the HI emission into distinct velocity components and demonstrate a detection of dust polarization correlated with the local Galactic HI but not with the HI associated with Magellanic Stream I. We present a robust, multi-frequency detection of polarized dust emission correlated with the filamentary HI morphology template down to 95 GHz. We find the correlation for frequencies between 95-353 GHz to be between ~10-50% over multipoles $20<\ell<200$. We measure the spectral index of the filamentary dust component spectral energy distribution to be $\beta = 1.54 \pm 0.13$. We find no evidence for decorrelation in this region between the filaments and the rest of the dust field or from the inclusion of dust associated with the intermediate velocity HI. Finally, we explore the morphological parameter space in the HI-based filamentary model.

Starburst galaxies are well-motivated astrophysical emitters of high-energy gamma-rays. They are well-known cosmic-ray "reservoirs", thanks to their large magnetic fields which confine high-energy protons for $\sim 10^5$ years. Over such long times, cosmic-ray transport can be significantly affected by scatterings with sub-GeV dark matter. Here we point out that this scattering distorts the cosmic-ray spectrum, and the distortion can be indirectly observed by measuring the gamma-rays produced by cosmic-rays via hadronic collisions. Present gamma-ray data show no sign of such a distortion, leading to stringent bounds on the cross section between protons and dark matter. These are competitive with current bounds, but have large room for improvement with the future gamma-ray measurements in the 0.1--10 TeV range from the Cherenkov Telescope Array, which can strengthen the limits by as much as two orders of magnitude.

Aims. Stars with strong enhancements of r-process elements are rare and tend to be metal-poor, with generally [Fe/H] <-2 dex and found in the halo. In this work we aim to investigate a candidate r-process enriched bulge star with a relatively high metallicity of -0.65 dex, and compare it with a previously published r-rich candidate star in the bulge. Methods. We reconsider the abundance analysis of a high-resolution optical spectrum of the red-giant star 2MASS J18082459-2548444 and determine its europium (Eu) and molybdenum (Mo) abundance, using stellar parameters from five different previous studies. Applying 2MASS photometry, Gaia astrometry and kinematics, we estimate distance, orbits, and population membership of 2MASS J18082459-2548444 and a previously reported r-enriched star 2MASS J18174532-3353235. Results. We find that 2MASS J18082459-2548444 is a relatively metal rich enriched r-process star that is enhanced in Eu and Mo but not substantially enhanced in s-process elements. It has a high probability of membership in the Galactic bulge based on its distance and orbit. We find that both stars show r-process enhancement with elevated [Eu/Fe]-values, even though 2MASS J18174532-3353235 is 1 dex lower in metallicity. Additionally, we find that 2MASS J18174532-3353235 plausibly has a halo or thick disc origin. Conclusions. We conclude that 2MASS J18082459-2548444 represents the first example of a confirmed r-process enhanced star confined to the inner bulge, possibly a relic from a period of enrichment associated with the formation of the bar.

When traveling in the heliosphere, Galactic cosmic rays (GCRs) are subjected to the solar modulation effect, a quasi-periodical change of their intensity caused by the 11-year cycle of solar activity. Here we investigate the association of solar activity and cosmic radiation over five solar cycles, from 1965 to 2020, using a collection of multichannel data from neutron monitors, space missions, and solar observatories. In particular, we focus on the time lag between the monthly sunspot number and the GCR flux variations. We show that the modulation lag is subjected to a 22-year periodical variation, ranging from about 2 to 14 months and following the polarity cycle of the Sun's magnetic field. We also show that the lag is remarkably decreasing with increasing energy of the GCR particles. These results reflect the interplay of basic physics phenomena that cause the GCR modulation effect: the drift motion of charged particles in the interplanetary magnetic field, the latitudinal dependence of the solar wind, the energy dependence of their residence time in the heliosphere. Based on this interpretation, we end up with a global effective formula for the modulation lag and testable predictions for the flux evolution of cosmic particles and antiparticles over the solar cycle.

Liller 1 and Terzan 5 are two massive systems in the Milky-Way bulge hosting populations characterized by significantly different ages ($\Delta t>7-8$ Gyr) and metallicities ($\Delta$[Fe/H]$\sim1$ dex). Their origin is still strongly debated in the literature and all formation scenarios proposed so far require some level of fine-tuning. The detailed star formation histories (SFHs) of these systems may represent an important piece of information to assess their origin. Here we present the first attempt to perform such an analysis for Liller 1. The first key result we find is that Liller 1 has been forming stars over its entire lifetime. More specifically, three broad SF episodes are clearly detected: 1) a dominant one, occurred some 12-13 Gyr ago with a tail extending for up to $\sim3$ Gyr, 2) an intermediate burst, between 6 and 9 Gyr ago, 3) and a recent one, occurred between 1 and 3 Gyr ago. The old population contributes to about $70\%$ of the total stellar mass and the remaining fraction is almost equally split between the intermediate and young populations. If we take these results at a face value, they would suggest that this system unlikely formed through the merger between an old globular cluster and a Giant Molecular Cloud, as recently proposed. On the contrary, our findings provide further support to the idea that Liller 1 is the surviving relic of a massive primordial structure that contributed to the Galactic bulge formation, similarly to the giant clumps observed in star-forming high-redshift galaxies.

We present a comprehensive study of 29 short gamma-ray bursts (SGRBs) observed $\approx 0.8-60$ days post-burst using $Chandra$ and $XMM-Newton$. We provide the inferred distributions of SGRB jet opening angles and true event rates to compare against neutron star merger rates. We perform uniform analysis and modeling of their afterglows, obtaining 10 opening angle measurements and 19 lower limits. We report on two new opening angle measurements (SGRBs 050724A and 200411A) and eight updated values, obtaining a median value of $\langle \theta_{\rm j} \rangle \approx 6.1^{\circ}$ [-3.2$^{\circ}$,+9.3$^{\circ}$] (68\% confidence on the full distribution) from jet measurements alone. For the remaining events, we infer $\theta_{\rm j}\gtrsim 0.5-26^{\circ}$. We uncover a population of SGRBs with wider jets of $\theta_{\rm j} \gtrsim 10^{\circ}$ (including two measurements of $\theta_{\rm j} \gtrsim 15^{\circ}$), representing $\sim 28\%$ of our sample. Coupled with multi-wavelength afterglow information, we derive a total true energy of $\langle E_{\rm true, tot} \rangle \approx 10^{49}-10^{50}$\,erg which is consistent with MHD jet launching mechanisms. Furthermore, we determine a range for the beaming-corrected event rate of $\mathfrak{R}_{\rm true} \approx360-1800$ Gpc$^{-3}$ yr$^{-1}$, set by the inclusion of a population of wide jets on the low end, and the jet measurements alone on the high end. From a comparison with the latest merger rates, our results are consistent with the majority of SGRBs originating from binary neutron star mergers. However, our inferred rates are well above the latest neutron star-black hole merger rates, consistent with at most a small fraction of SGRBs originating from such mergers.

Starobinsky inflation is an attractive, fundamental model to explain the Planck measurements, and its higher-order extension may allow us to probe quantum gravity effects. We show that future CMB data combined with the 21cm intensity map from SKA will meaningfully probe such an extended Starobinsky model. A combined analysis will provide a precise measurement and intriguing insight into inflationary dynamics, even accounting for correlations with astrophysical parameters.

Galactic cosmic rays (GCRs) inside the heliosphere are affected by magnetic turbulence and Solar wind disturbances which result in the so-called solar modulation effect. To investigate this phenomenon, we have performed a data-driven analysis of the temporal dependence of the GCR flux over the solar cycle. With a global statistical inference of GCR data collected in space by AMS-02, PAMELA, and CRIS on monthly basis, we have determined the dependence of the GCR diffusion parameters upon time and rigidity. In this conference, we present our results for GCR protons and nuclei, we discuss their interpretation in terms of basic processes of particle transport and their relations with the dynamics of the heliospheric plasma.

The observed variability of the cosmic-ray intensity in the interplanetary space is driven by the evolution of the Sun's magnetic activity over its 11-year quasiperiodical cycle. Investigating the relationship between solar activity indices and cosmic-ray intensity measurements is then essential for understanding the fundamental processes of particle transport in the heliosphere. Here we have performed a global characterization the solar modulation of cosmic rays over the solar activity cycle and for different energies of the cosmic particles. We present our cross-correlation studies using data from space experiments, neutron monitors and solar observatories collected over several solar cycles.

In this work we present updated forecasts on parameterised modifications of gravity that can capture deviations of the behaviour of cosmological density perturbations beyond $\Lambda$CDM. For these forecasts we adopt the SKA Observatory (SKAO) as a benchmark for future cosmological surveys at radio frequencies, combining a continuum survey for weak lensing and angular galaxy clustering with an HI galaxy survey for spectroscopic galaxy clustering that can detect baryon acoustic oscillations and redshift space distortions. Moreover, we also add 21cm HI intensity mapping, which provides invaluable information at higher redshifts, and can complement tomographic resolution, thus allowing us to probe redshift-dependent deviations of modified gravity models. For some of these cases, we combine the probes with other optical surveys, such as the Dark Energy Spectroscopic Instrument (DESI) and the Vera C. Rubin Observatory (VRO). We show that such synergies are powerful tools to remove systematic effects and degeneracies in the non-linear and small-scale modelling of the observables. Overall, we find that the combination of all SKAO radio probes will have the ability to constrain the present value of the functions parameterising deviations from $\Lambda$CDM ($\mu$ and $\Sigma$) with a precision of $2.7\%$ and $1.8\%$ respectively, competitive with the constraints expected from optical surveys and with constraints we have on gravitational interactions in the standard model. Exploring the radio-optical synergies, we find that the combination of VRO with SKAO can yield extremely tight constraints on $\mu$ and $\Sigma$ ($0.9\%$ and $0.7\%$ respectively), which are further improved when the cross-correlation between intensity mapping and DESI galaxies is included.

The stochastic gravitational wave background (SGWB) is a rich resource of cosmological information, encoded both in its source statistics and its anisotropies induced by propagation effects. We provide a theoretical description of it, without employing the Boltzmann equation which is exclusively valid in the geometric optics limit. Our formalism is based on the so-called classical matter approximation and it is able to capture wave-optics effects, such as interference and diffraction. We show that the interaction between the gravitational waves and the cosmic structures along the line-of-sight produce an observable scalar and vector polarization. We build the two point correlation function describing the statistics of the tensor modes of the SGWB, and introduce the gravitational Stokes parameters. In the case of an unpolarized, Gaussian, statistically homogeneous and isotropic SGWB, we show that the interaction with matter does not generate a net difference between left- and right- helicity tensor modes, as expected, but it may produce Q- and U- polarization modes.

We draw attention to two unanswered questions in radio pulsars' theory. The first one concerns the inclination angles between magnetic and rotation axes. We show that the very existence of interpulse pulsars indicates that all possible angles are realized. The second one is the question of the break of the death line on the PPdot-diagram. We show that this break can be easily explained by a very mechanism of secondary plasma generation.

Precision cosmology with gravitational wave (GW) sources requires understanding the interplay between GW source population and cosmological parameters governing the dynamics of the Universe. With the fast increase of GW detections, for exploring many aspects of cosmology and fundamental physics it is necessary to develop a tool which can simulate GW mock samples for several population and cosmological models with and without a galaxy catalog. We have developed a new code called GWSim, allowing to make GW mock events from a large range of configurations, varying the cosmology, the merger rate, and the GW source parameters (mass and spin distributions in particular), for a given network of GW detectors. We restrict the cosmology to spatially flat universes, including models with varying dark energy equation of state. GWSim provides each mock event with a position in the sky and a redshift; these values can be those of random host galaxies coming from an isotropic and homogeneous simulated Universe or a user-supplied galaxy catalog. We use realistic detector configurations of the LIGO and Virgo network of detectors to show the performance of this code for the latest observation runs and the upcoming observation run.

The fate of stars in the zero-age main-sequence (ZAMS) range $\approx 8-12$ Msun is unclear. They could evolve to form white dwarfs or explode as electron-capture supernovae (SNe) or iron core-collapse SNe (CCSNe). Even though the initial mass function indicates that this mass range should account for over 40% of all CCSNe progenitors, few have been observationally confirmed, likely owing to the faintness and rapid evolution of these transients. In this paper, we present a sample of nine Ca-rich/O-poor Type IIb SNe detected by the Zwicky Transient Facility with progenitors likely in this mass range. We perform a holistic analysis of the spectroscopic and photometric properties of the sample. These sources have a flux ratio of [Ca II] $\lambda \lambda$7291, 7324 to [O I] $\lambda \lambda$6300, 6364 of $\gtrsim$ 2 in their nebular spectra. Comparing the measured [O I] luminosity ($\lesssim 10^{39} \mathrm{erg\ s^{-1}}$) and derived oxygen mass ($\lesssim 0.1$ Msun) with theoretical models, we infer that the progenitor ZAMS mass for these explosions is less than 12 Msun. These correspond to He-stars with core masses less than around 3 Msun. We find that the ejecta properties (Mej $\lesssim 1$ Msun) are also consistent with those expected for such low mass He-stars. The low ejecta mass of these sources indicates a class of strongly-stripped SNe that is a transition between the regular stripped-envelope SNe and ultra-stripped SNe. The progenitor could be stripped by a main sequence companion and result in the formation of a neutron star $-$ main sequence binary. Such binaries have been suggested to be progenitors of neutron star $-$ white dwarf systems that could merge within a Hubble time, and be detectable with LISA.

The Pierre Auger Observatory has been detecting ultra-high energy cosmic rays (UHECRs) for more than fifteen years. An essential feature of the Observatory is its hybrid design: cosmic rays above $100~$PeV are detected through the observation of the associated air showers with different and complementary techniques, from surface detector arrays and fluorescence telescopes to radio antennas. The analyses of the multi-detector data have enabled high-statistics and high-precision studies of the energy spectrum, mass composition and distribution of arrival directions of UHECRs. The resulting picture is summarized in this contribution. While no discrete source of UHECRs has been identified so far, the extragalactic origin of the particles has been recently determined from the arrival directions above 8~EeV, and the ring is closing around nearby astrophysical sites at higher energies. Also, the established upper limits on fluxes of UHE neutrinos and photons have implications on dark matter and cosmological aspects that are also presented in this contribution.

We construct a family of simple single-field inflation models consistent with Planck / BICEP Keck bounds which have a parametrically small tensor amplitude and no running of the scalar spectral index. The construction consists of a constant-roll hilltop inflaton potential with the end of inflation left as a free parameter induced by higher-order operators which become dominant late in inflation. This construction directly demonstrates that there is no lower bound on the tensor/scalar ratio for simple single-field inflation models.

We present the results of a {\it James Webb Space Telescope} NIRCam investigation into the young massive star cluster (YMC) population in the luminous infrared galaxy VV 114. We identify 374 compact YMC candidates with a $S/N \geq 3$, 5, and 5 at F150W, F200W, and F356W respectively. A direct comparison with our {\it HST} cluster catalog reveals that $\sim 20\%$ of these sources are undetected at optical wavelengths. Based on {\it yggdrasil} stellar population models, we identify 16 YMCs in our {\it JWST} catalog with F150W-F200W and F200W-F356W colors suggesting they are all very young ($1-4$ Myr), dusty ($A_{V} \geq 5$), and massive ($10^{5.8} < M_{\odot} < 10^{6.1}$). The discovery of these `hidden' sources, many of which are found in the `overlap' region between the two nuclei, triples the number of $t < 4 $ Myr clusters detected in VV 114. Extending the cluster age distribution ($dN/d\tau \propto \tau^{\gamma}$) to the youngest ages, we find a slope of $\gamma = -1.27 \pm 0.32$ for $10^{6} < \tau (\mathrm{yr}) < 10^{7}$, which is consistent with the previously determined value from $10^{7} < \tau (\mathrm{yr}) < 10^{8.5}$, and confirms that VV 114 has a steep age gradient for all massive star clusters across the entire range of cluster ages observed. Finally, the consistency between our {\it JWST}- and {\it HST}-derived age distribution slopes indicates that the balance between cluster formation and destruction has not been significantly altered in VV 114 over the last 0.5 Gyr.

We study the stability of Hall MHD with strong magnetic fields in which Landau quantization of electrons is important. We find that the strong-field Hall modes can be destabilized by the dependence of the differential magnetic susceptibility on magnetic field strength. This instability is studied using linear perturbation theory, and is found to have typical growth time of order $\lesssim 10^3$ yrs, with the growth time decreasing as a function of wavelength of the perturbation. The instability is self-limiting, turning off following a period of local field growth by a few percent of the initial value. Finite temperature is also shown to limit the instability, with sufficiently high temperatures eliminating it altogether. Alfv\'{e}n waves can show similar unstable behaviour on shorter timescales. We find that Ohmic heating due to the large fields developed via the instability and magnetic domain formation is not large enough to account for observed magnetar surface temperatures. However, Ohmic heating is enhanced by the oscillatory differential magnetic susceptibility of Landau-quantized electrons, which could be important to magneto-thermal simulations of neutron star crusts.

We conduct hydrodynamical cosmological zoom simulations of fourteen voids to study the ability of haloes to accrete gas at different locations throughout the voids at z = 0. Measuring the relative velocity of haloes with respect to their ambient gas, we find that a tenth of the haloes are expected to be unable to accrete external gas due to its fast flow passed them (so called 'fast flow haloes'). These are typically located near void walls. We determine that these haloes have recently crossed the void wall and are still moving away from it. their motion counter to that of ambient gas falling towards the void wall results in fast flows that make external gas accretion very challenging, and often cause partial gas loss via the resultant ram pressures. Using an analytical approach, we model the impact of such ram pressures on the gas inside haloes of different masses. A halo's external gas accretion is typically cut off, with partial stripping of halo gas. For masses below a few times 10$^{9}$ M$_{\odot}$, their halo gas is heavily truncated but not completely stripped. We identify numerous examples of haloes with a clear jelly-fish like gas morphology, indicating their surrounding gas is being swept away, cutting them off from further external accretion. These results highlight how, even in the relatively low densities of void walls, a fraction of galaxies can interact with large-scale flows in a manner that has consequences for their gas content and ability to accrete gas.

Understanding the chemistry of the interstellar medium (ISM) is fundamental for the comprehension of the Galactic and stellar evolution. X-rays provide an excellent way to study the dust chemical composition and crystallinity along different sight-lines in the Galaxy. In this work we study the dust grain chemistry in the diffuse regions of the interstellar medium in the soft X-ray band (<1 keV). We use newly calculated X-ray dust extinction cross sections, obtained from laboratory data, in order to investigate the oxygen K and iron L shell absorption. We explore the XMM-Newton and Chandra spectra of 5 low-mass X-ray binaries located in the Galactic plane, and we model the gas and dust features of oxygen and iron simultaneously. The dust samples used for this study include silicates with different Mg:Fe ratios, sulfides, iron oxides and metallic iron. Most dust samples are in both amorphous and crystalline lattice configuration. The extinction cross sections have been computed using Mie scattering approximation and assuming a power law dust size distribution. We find that the Mg-rich amorphous pyroxene (Mg0.75Fe0.25SiO3) represents the largest fraction of dust towards most of the X-ray sources, about 70% on average. Additionally, we find that ~15% of the dust column density in our lines of sight is in Fe metallic. We do not find strong evidence for ferromagnetic compounds, such as Fe3O4 or iron sulfides (FeS, FeS2). Our study confirms that the iron is heavily depleted from the gas phase into solids; more than 90% of iron is in dust. The depletion of neutral oxygen is mild, between 10-20% depending on the line of sight.

We present the technical design, construction and testing of the Colibri telescope array at Elginfield Observatory near London, Ontario, Canada. Three 50-cm telescopes are arranged in a triangular array and are separated by 110-160 metres. During operation, they will monitor field stars at the intersections of the ecliptic and galactic plane for serendipitous stellar occultations (SSOs) by trans-Neptunian objects (TNOs). At a frame rate of 40 frames per second (fps), Fresnel diffraction in the occultation light curve can be resolved and, with coincident detections, be used to estimate basic properties of the occulting object. Using off-the-shelf components, the Colibri system streams imagery to disk at a rate of 1.5 GB/s for next-day processing by a custom occultation detection pipeline. The imaging system has been tested and is found to perform well, given the moderate site conditions. Limiting magnitudes at 40 fps are found to be about 12.1 (temporal SNR=5, visible light Gaia G band) with time-series standard deviations ranging from about 0.035 mag to >0.2 mag. SNR is observed to decrease linearly with magnitude for stars fainter than about G = 9.5 mag. Brighter than this limit, SNR is constant, suggesting that atmospheric scintillation is the dominant noise source. Astrometric solutions show errors typically less than approximately 0.3 pixels (0.8 arc seconds) without a need for high-order corrections.

Scalar field dark matter offers an interesting alternative to the traditional WIMP dark matter picture. Astrophysical and cosmological simulations are useful to constraining the mass of the dark matter particle in this model. This is particularly true at low mass where the wavelike nature of the dark matter particle manifests on astrophysical scales. These simulations typical use a classical field approximation. In this work, we look at extending these simulations to include quantum corrections. We look into both the ways in which large corrections impact the predictions of scalar field dark matter, and the timescales on which these corrections grow large. Corrections tend to lessen density fluctuations and increase the effect of "quantum pressure". During collapse, these corrections grow exponentially, quantum corrections would become important in about ~30 dynamical times. This implies that the predictions of classical field simulations may differ from those with quantum corrections for systems with short dynamical times.

We report the discovery of a two-dimensional Galaxy Manifold within the multi-dimensional luminosity space of local galaxies. The multi-dimensional luminosity space is constructed using 11 bands that span from far ultraviolet to near-infrared for redshift < 0.1 galaxies observed with GALEX, SDSS, and UKIDSS. The two latent parameters are sufficient to express 93.2% of the variance in the galaxy sample, suggesting that this Galaxy Manifold is one of the most efficient representations of galaxies. The transformation between the observed luminosities and the manifold parameters as an analytic mapping is provided. The manifold representation provides accurate (85%) morphological classifications with a simple linear boundary, and galaxy properties can be estimated with minimal scatter (0.12 dex and 0.04 dex for star formation rate and stellar mass, respectively) by calibrating with the two-dimensional manifold location. Under the assumption that the manifold expresses the possible parameter space of galaxies, the evolution on the manifold is considered. We find that constant and exponentially decreasing star formation histories form almost orthogonal modes of evolution on the manifold. Through these simple models, we understand that the two modes are closely related to gas content, which suggests the close relationship of the manifold to gas accretion. Without assuming a star formation history, a gas-regulated model reproduces an exponentially declining star formation history with a timescale of $\sim$1.2 Gyrs on the manifold. Lastly, the found manifold suggests a paradigm where galaxies are characterized by their mass/scale and specific SFR, which agrees with previous studies of dimensionality reduction.

I discuss the feasibility of a conceptual space-based neutrino detector that utilizes the Ice Giants as Targets for Galactic Neutrinos. The purpose of this research stems from the concept of wanting to find a new method of observing the Galactic Core (GC) of the Milky Way and the Supermassive black hole, Sag A*. Observations of the GC have been made in every accessible wavelength except for the regions of space that are too dense for photons to probe. In these regions, we may instead use neutrinos. Neutrinos from the Active Galactic Nucleus are emitted at extreme energies, 10 GeV to EeV scales, but have an extremely low flux measured here at Earth. Neutrino telescopes such as the IceCube Observatory have only been able to measure a handful of neutrinos that might correlate to the GC. But using Gravitational lensing, our sun can be used as a lens which increases the 'light' collection power for neutrinos by a factor of $10^{13}$, with the trade-off that the minimum focal point is located at 22 AU. This means that Uranus and Neptune are suitable natural targets for these neutrinos to interact with and observe the effects of a spacecraft in orbit. Initial studies use GEANT4, a particle physics simulation toolbox developed by CERN, to facilitate the propagation of energetic particles passing through the atmosphere of Neptune. Various aspects are studied to understand the behaviors of these particle interactions. For each of these aspects, we modify several variables such as particle type, energy, interaction depth, and orbital distance from the surface. I also discuss the versatility of this neutrino detector which has the possibility of mapping out the inner structure of the Ice Giants, in-depth studies of the neutrinos coming from the GC, and possibilities to use this method for other cosmic neutrino sources.

In this paper, we report the multiwavelength observations of the partial filament eruption associated with a C1.2 class flare in NOAA active region 11236 on 13 June 2011. The event occurred at the eastern limb in the field of view (FOV) of Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO) spacecraft and was close to the disk center in the FOV of Extreme-UltraViolet Imager (EUVI) on board the behind Solar Terrestrial Relations Observatory (STEREO) spacecraft. During eruption, the filament splits into two parts: the major part and runaway part. The major part flows along closed loops and experiences bifurcation at the loop top. Some of the materials move forward and reach the remote footpoint, while others return back to the original footpoint. The runaway part flows along open field lines, which is evidenced by a flare-related type III radio burst. The runaway part also undergoes bifurcation. The upper branch of escapes the corona and evolves into a jet-like narrow coronal mass ejection (CME) at a speed of 324 km s-1, while the lower branch falls back to the solar surface. A schematic cartoon is proposed to explain the event and provides a new mechanism of partial filament eruptions

We present and analyze observations of polarized dust emission at 850 $\mu$m towards the central 1 pc $\times$ 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the BISTRO (B-fields in Star-forming Region Observations) survey. The orientations of the magnetic field follow the spiral structure of Mon R2, which are well-described by an axisymmetric magnetic field model. We estimate the turbulent component of the magnetic field using the angle difference between our observations and the best-fit model of the underlying large-scale mean magnetic field. This estimate is used to calculate the magnetic field strength using the Davis-Chandrasekhar-Fermi method, for which we also obtain the distribution of volume density and velocity dispersion using a column density map derived from $Herschel$ data and the C$^{18}$O ($J$ = 3-2) data taken with HARP on the JCMT, respectively. We make maps of magnetic field strengths and mass-to-flux ratios, finding that magnetic field strengths vary from 0.02 to 3.64 mG with a mean value of 1.0 $\pm$ 0.06 mG, and the mean critical mass-to-flux ratio is 0.47 $\pm$ 0.02. Additionally, the mean Alfv\'en Mach number is 0.35 $\pm$ 0.01. This suggests that in Mon R2, magnetic fields provide resistance against large-scale gravitational collapse, and magnetic pressure exceeds turbulent pressure. We also investigate the properties of each filament in Mon R2. Most of the filaments are aligned along the magnetic field direction and are magnetically sub-critical.

Obscuring winds driven away from active supermassive black holes are rarely seen due to their transient nature. They have been observed with multi-wavelength observations in a few Seyfert 1 galaxies and one broad absorption line radio-quiet quasar so far. An X-ray obscuration event in MR 2251-178 was caught in late 2020, which triggered multi-wavelength (NIR to X-ray) observations targeting this radio-quiet quasar. In the X-ray band, the obscurer leads to a flux drop in the soft X-ray band from late 2020 to early 2021. X-ray obscuration events might have a quasi-period of two decades considering earlier events in 1980 and 1996. In the UV band, a forest of weak blueshifted absorption features emerged in the blue wing of Ly$\alpha$ $\lambda1216$ in late 2020. Our XMM-Newton, NuSTAR, and HST/COS observations are obtained simultaneously, hence, the transient X-ray obscuration event is expected to account for the UV outflow, although they are not necessarily caused by the same part of the wind. Both blueshifted and redshifted absorption features were found for He {\sc i} $\lambda10830$, but no previous NIR spectra are available for comparison. The X-ray observational features of MR 2251-178 shared similarities with some other type 1 AGNs with obscuring wind. However, observational features in the UV to NIR bands are distinctly different from those seen in other AGN with obscuring winds. A general understanding of the observational variety and the nature of obscuring wind is still lacking.

It has been investigated the possibility of the various atmospheres over water oceans. We have considered the H$_2$ atmosphere and He atmosphere concerning to N$_2$ atmosphere over oceans. One of the main subjects in astrobiology is to estimate the habitable zone. If there is an ocean on the planet with an atmosphere, there is an upper limit to the outgoing infrared radiation called the Komabayashi-Ingersoll limit (KI-limit). This limit depends on the components of the atmospheres. We have investigated this dependence under the simple model, using the one-dimensional gray radiative-convective equilibrium model adopted by Nakajima et al. (1992). The outgoing infrared radiation ($F_{IRout}$) with the surface temperature ($T_s$) has shown some peculiar behavior. The examples for H$_2$, He, and N$_2$ background gas for H$_2$O vapour are investigated. There is another limit called the Simpson-Nakajima limit (SN-limit) mainly composed of vapour. This steam limit does not depend on the background atmosphere components. Under super-Earth case ($g=2\times$9.8 m/s$^2$), several cases are also calculated. The KI-limit dependence on the initial pressure is presented. The various emission rates by Koll & Cronin (2019) are investigated.

Very fast novae are novae which evolve exceptionally quickly (on timescales of only days). Due to their rapid evolution, very fast novae are challenging to detect and study, especially at early times. Here we report the discovery, which was made as part of our Transient UV Objects project, of a probable very fast nova in the nearby spiral galaxy NGC 300. We detected the rise to the peak (which is rarely observed for very fast novae) in the near-ultraviolet (NUV), with the first detection just ~2 hours after the eruption started. The peak and early stages of the decay were also observed in UV and optical bands. The source rapidly decayed by 2 NUV magnitudes within 3.5 days, making it one of the fastest novae known. In addition, a likely quiescent counterpart was found in archival near-infrared Spitzer and VIRCAM images but not in any deep optical and UV observations, indicating a very red spectral shape in quiescence. The outburst and quiescence properties suggest that the system is likely a symbiotic binary. We discuss this new transient in the context of very fast novae in general and specifically as a promising supernova Type Ia progenitor candidate, due to its very high inferred WD mass (~1.35 Ms; determined by comparing this source to other very fast novae).

Disc winds and planet formation are considered to be two of the most important mechanisms that drive the evolution and dispersal of protoplanetary discs and in turn define the environment in which planets form and evolve. While both have been studied extensively in the past, we combine them into one model by performing three-dimensional radiation-hydrodynamic simulations of giant planet hosting discs that are undergoing X-ray photo-evaporation, with the goal to analyse the interactions between both mechanisms. In order to study the effect on observational diagnostics, we produce synthetic observations of commonly used wind-tracing forbidden emission lines with detailed radiative transfer and photo-ionisation calculations. We find that a sufficiently massive giant planet carves a gap in the gas disc that is deep enough to affect the structure and kinematics of the pressure-driven photo-evaporative wind significantly. This effect can be strong enough to be visible in the synthetic high-resolution observations of some of our wind diagnostic lines, such as the [OI] 6300 \r{A} or [SII] 6730 \r{A} lines. When the disc is observed at inclinations around 40{\deg} and higher, the spectral line profiles may exhibit a peak in the redshifted part of the spectrum, which cannot easily be explained by simple wind models alone. Moreover, massive planets can induce asymmetric substructures within the disc and the photo-evaporative wind, giving rise to temporal variations of the line profiles that can be strong enough to be observable on timescales of less than a quarter of the planet's orbital period.

AGN feedback during galaxy merger has been the most favoured model to explain black hole-galaxy co-evolution. However, how the AGN-driven jet/wind/radiation is coupled with the gas of the merging galaxies, which leads to positive feedback, momentarily enhanced star formation, and subsequently negative feedback, a decline in star formation, is poorly understood. Only a few cases are known where the jet and companion galaxy interaction leads to minor off-axis distortions in the jets and enhanced star formation in the gas-rich minor companions. Here, we briefly report one extraordinary case, RAD12, discovered by RAD@home citizen science collaboratory, where for the first time a radio jet-driven bubble ~137 kpc is showing a symmetric reflection after hitting the incoming galaxy which is not a gas-rich minor but a gas-poor early-type galaxy in a major merger. Surprisingly, neither positive feedback nor any radio lobe on the counter jet side, if any, is detected. It is puzzling if RAD12 is a genuine one-sided jet or a case of radio lobe trapped, compressed and re-accelerated by shocks during the merger. This is the first imaging study of RAD12 presenting follow-up with the GMRT, archival MeerKAT radio data and CFHT optical data.

Decaying Cold Dark Matter (DCDM) is a model that is currently under investigation regarding primarily the $S_8$ tension between cosmic microwave background (CMB) and certain large-scale structure measurements. The decay into one massive and one (or more) massless daughter particle(s) leads to a suppression of the power spectrum in the late universe that depends on the relative mass splitting $\epsilon=(1-m^2/M^2)/2$ between the mother and massive daughter particle as well as the lifetime $\tau$. In this work we investigate the impact of the BOSS DR14 one-dimensional Lyman-$\alpha$ forest flux power spectrum on the DCDM model using a conservative effective model approach to account for astrophysical uncertainties. Since the suppression of the power spectrum due to decay builds up at low redshift, we find that regions in parameter space that address the $S_8$ tension can be well compatible with the Lyman-$\alpha$ forest. Nevertheless, for values of the degeneracy parameter $\epsilon\sim 0.1-0.5\%$, for which the power suppression occurs within the scales probed by BOSS Lyman-$\alpha$ data, we find improved constraints compared to previous CMB and galaxy clustering analyses, obtaining $\tau\gtrsim 18$ Gyrs for small mass splitting. Furthermore, our analysis of the BOSS Lyman-$\alpha$ flux power spectrum allows for values $\tau\sim 10^2$ Gyrs, $\epsilon\sim 1\%$, that have been found to be preferred by a combination of Planck and galaxy clustering data with a KiDS prior on $S_8$, and we even find a marginal preference within this regime.

The spectral index $n_s$ of scalar perturbation is the significant initial condition set by inflation theory for our observable Universe. According to Planck results, current constraint is $n_s = 0.965\pm 0.004$, while an exact scale-invaiant Harrison-Zeldovich spectrum, i.e. $n_s=1$, has been ruled out at $8.4\sigma$ significance level. However, it is well-known that the standard $\Lambda$CDM model is suffering from the Hubble tension, which is at $\sim 5\sigma$ significance level. This inconsistency likely indicates that the comoving sound horizon at last scattering surface is actually lower than expected, which so seems to be calling for the return of $n_s=1$. Here, in light of recent observations we find strong evidence for a $n_s=1$ Universe. And we show that if so, it would be confirmed conclusively by CMB-S4 experiment.

A particular open problem in cosmology is whether dark matter on small scales is clumpy, forming gravitationally-bound halos distributed within the Galaxy. The practical difficulties inherent in testing this hypothesis stem from the fact that, on astrophysical scales, dark matter is solely observable via its gravitational interaction with other objects. This thesis presents a gravitational-lensing-based solution for the mapping and characterisation of low-mass, dark matter halos via their signature in millisecond pulsar observations. This involves: first, determining the time delay and magnification surfaces generated in the frame of reference of the halo; second, obtaining the corresponding pulsar signature in the reference frame of the observer; and last, generalising the method to multiple halos at varying distances. We discuss whether the delay is observationally detectable for both single and multiple lenses. The key dependency of the time delay is the density profile adopted for the halo. I utilise a variety of proposed halo mass profiles -- elliptical, Schwarzschild, horizontal-disc lenses and the Navarro-Frenk-White (NFW) density profile -- which are applicable over a broad range of halo masses. I demonstrate the use of Hankel transforms to increase the efficiency of the relativistic time delay calculation. The observational signatures of such halos are best identified using millisecond pulsars due to their high rotational frequencies and period stability. My method does not require major adjustments when searching for signs of lensing, thus it is unnecessary to implement specialist data reduction pipelines. Thus we can leverage data from both existing and future surveys easily. This method is readily extensible to nearby globular clusters and galaxies, pending improvements in pulsar detection at such distances.

This paper is dedicated to the new implicit unstructured coronal code COCONUT, which aims at providing fast and accurate inputs for space weather forecast as an alternative to empirical models. We use all 20 available magnetic maps of the solar photosphere covering the date of the 2nd of July 2019 which corresponds to a solar eclipse on Earth. We use the same standard pre-processing on all maps, then perform coronal MHD simulations with the same numerical and physical parameters. In the end, we quantify the performance for each map using three indicators from remote-sensing observations: white-light total solar eclipse images for the streamers' edges, EUV synoptic maps for coronal holes and white-light coronagraph images for the heliospheric current sheet. We discuss the performance for space weather forecasts and we show that the choice of the input magnetic map has a strong impact. We find performances between 24% to 85% for the streamers' edges, 24% to 88% for the coronal hole boundaries and a mean deviation between 4 to 12 degrees for the heliospheric current sheet position. We find that the HMI runs are globally performing better on all indicators, with the GONG-ADAPT being the second-best choice. HMI runs perform better for the streamers' edges, GONG-ADAPT for polar coronal holes, HMI synchronic for equatorial coronal holes and for the streamer belt. We especially showcase the importance of the filling of the poles. This demonstrates that the solar poles have to be taken into account even for ecliptic plane previsions.

We investigate the cosmological evolution of interstellar dust with a semi-analytical galaxy formation model ($\nu^2$GC), focusing on the evolution of grain size distribution. The model predicts the statistical properties of dust mass and grain size distribution in galaxies across cosmic history. We confirm that the model reproduces the relation between dust-to-gas ratio and metallicity in the local Universe, and that the grain size distributions of the Milky Way (MW)-like sample become similar to the so-called MRN distribution that reproduces the observed MW extinction curve. Our model, however, tends to overpredict the dust mass function at the massive end at redshift $z\lesssim 0.8$ while it reproduces the abundance of dusty galaxies at higher redshifts. We also examine the correlation between grain size distribution and galaxy properties (metallicity, specific star formation rate, gas fraction, and stellar mass), and observe a clear trend of large-grain-dominated, small-grain-dominated, and MRN-like grain size distributions from unevolved to evolved stages. As a consequence, the extinction curve shapes are flat, steep, and intermediate (MW-like) from the unevolved to evolved phases. At a fixed metallicity, the grain size distribution tends to have larger fractions of small grains at lower redshift; accordingly, the extinction curve tends to be steeper at lower redshift. We also predict that supersolar-metallicity objects at high redshift have flat extinction curves with weak 2175 \AA bump strength.

Nielsen-Olesen vortices in the Abelian-Higgs (AH) model are the simplest realisations of cosmic strings in a gauge field theory. Large-scale numerical solutions show that the dominant decay channel of a network of AH strings produced from random initial conditions is classical field radiation. However, they also show that with special initial conditions, loops of string can be created for which classical field radiation is suppressed, and which behave like Nambu-Goto (NG) strings with a dominant decay channel into gravitational radiation. This indicates that cosmic strings are generically sources of both high-energy particles and gravitational waves. Here we adopt a simple parametrisation of the AH string network allowing for both particle and gravitational wave production, which sets the basis for a "multi-messenger" investigation of this model. We find that, in order to explain the NANOGrav detection of a possible gravitational wave background, while satisfying the constraint on NG-like loop production from simulations and bounds from the cosmic microwave background, the tension of the AH string in Planck units $G\mu$ and the fraction of the NG-like loops $f_{\rm NG}$ should satisfy $G\mu f_{\rm NG}^{2.6} \gtrsim 3.2\times 10^{-13}$ at 95$\%$ confidence. On the other hand, for such string tensions, constraints from the diffuse gamma-ray background (DGRB) indicate that more than 97$\%$ of the total network energy should be converted to dark matter (DM) or dark radiation. We also consider joint constraints on the annihilation cross-section, the mass, and the relic abundance of DM produced by decays of strings. For example, for a DM mass of 500 GeV, the observed relic abundance can be explained by decaying AH strings that also account for the NANOGrav signal.

Flare-associated quasi-periodic pulsations (QPPs) in radio and X-ray wavelengths, particularly those related to nonthermal electrons, contain important information about the energy release and transport processes during flares. However, the paucity of spatially resolved observations of such QPPs with a fast time cadence has been an obstacle for us to further understand their physical nature. Here, we report observations of such a QPP event occurred during the impulsive phase of a C1.8-class eruptive solar flare using radio imaging spectroscopy data from the Karl G. Jansky Very Large Array (VLA) and complementary X-ray imaging and spectroscopy data. The radio QPPs, observed by the VLA in the 1--2 GHz with a sub-second cadence, are shown as three spatially distinct sources with different physical characteristics. Two radio sources are located near the conjugate footpoints of the erupting magnetic flux rope with opposite senses of polarization. One of the sources displays a QPP behavior with a ~5-s period. The third radio source, located at the top of the post-flare arcade, coincides with the location of an X-ray source and shares a similar period of ~25--45 s. We show that the two oppositely polarized radio sources are likely due to coherent electron cyclotron maser (ECM) emission. On the other hand, the looptop QPP source, observed in both radio and X-rays, is consistent with incoherent gyrosynchrotron and bremsstrahlung emission, respectively. We conclude that the concurrent, but spatially distinct QPP sources must involve multiple mechanisms which operate in different magnetic loop systems and at different periods.

Measuring the magnetic field in cosmic filaments reveals how the Universe is magnetised and the process that magnetised it. Using the Rotation Measures (RM) at 144-MHz from the LoTSS DR2 data, we analyse the rms of the RM extragalactic component as a function of redshift to investigate the evolution with redshift of the magnetic field in filaments. From previous results, we find that the extragalactic term of the RM rms at 144-MHz is dominated by the contribution from filaments (more than 90 percent). Including an error term to account for the minor contribution local to the sources, we fit the data with a model of the physical filament magnetic field, evolving as $B_f = B_{f,0}\,(1+z)^\alpha$ and with a density drawn from cosmological simulations of five magnetogenesis scenarios. We find that the best-fit slope is in the range $\alpha = [-0.2, 0.1]$ with uncertainty of $\sigma_\alpha = 0.4$--0.5, which is consistent with no evolution. The comoving field decreases with redshift with a slope of $\gamma = \alpha - 2 = [-2.2, -1.9]$. The mean field strength at $z=0$ is in the range $B_{f,0}=39$--84~nG. For a typical filament gas overdensity of $\delta_g=10$ the filament field strength at $z=0$ is in the range $B_{f,0}^{10}=8$--26~nG. A primordial stochastic magnetic field model with initial comoving field of $B_{\rm Mpc} = 0.04$--0.11~nG is favoured. The primordial uniform field model is rejected.

Modeling the internal chemistry of molecular clouds is critical to accurately simulating their evolution. To reduce computational expense, 3D simulations generally restrict their chemical modeling to species with strong heating and cooling effects. We address this by post-processing tracer particles in the SILCC-Zoom molecular cloud simulations. Using a chemical network of 39 species and 299 reactions (including freeze-out of CO and H$_2$O), and a novel iterative algorithm to reconstruct a filled density grid from sparse tracer particle data, we produce time-dependent density distributions for various species. We focus upon the evolution of HCO$^+$, which is a critical formation reactant of CO but is not typically modeled on-the-fly. We analyse the evolution of the tracer particles to assess the regime in which HCO$^+$ production preferentially takes place. We find that the HCO$^+$ content of the cold molecular gas forms in situ around $n_\textrm{HCO$^+$}\simeq10^3$-$10^4$ cm$^{-3}$, over a time-scale of approximately 1 Myr, rather than being distributed to this density regime via turbulent mixing from deeper in the cloud. We further show that the dominant HCO$^+$ formation pathway is dependent on the visual extinction, with the reaction H$_3^+$ + CO contributing 90% of the total HCO$^+$ production flux above $A_\textrm{V,3D}=3$. Using our novel grid reconstruction algorithm, we produce the very first maps of the HCO$^+$ column density, $N$(HCO$^+$), and show that it reaches values as high as $10^{15}$ cm$^{-2}$. We find that 50% of the HCO$^+$ mass is located in an $A_\textrm{V}$-range of $\sim$10-30, and in a density range of $10^{3.5}$-$10^{4.5}$ cm$^{-3}$. Finally, we compare our $N$(HCO$^+$) maps to recent observations of W49A and find good agreement.

Solar active regions (ARs) play a fundamental role in driving many of the geo-effective eruptions which propagate into the Solar System. However, we are still unable to consistently predict where and when ARs will occur across the solar disk by identifying pre-emergence signatures in observables such as the Doppler velocity (without using Helioseismic methods). Here we aim to determine the earliest time at which pre-emergence signatures, specifically the Horizontal Divergent Flow (HDF), can be confidently detected using data from the Solar Dynamics Observatory's Helioseismic and Magnetic Imager (SDO/HMI). Initially, we follow previous studies using the thresholding method, which searches for significant increases in the number of pixels that display a specific line-of-sight velocity. We expand this method to more velocity windows and conduct a basic parameter study investigating the effect of cadence on the inferred results. Our findings agree with previous studies with $37.5$% of ARs displaying a HDF, with average lead times between the HDF and flux emergence of $58$ minutes. We present a new potential signature of flux emergence which manifests as cadence-independent transient disruptions to the amplitudes of multiple velocity windows and recover potential pre-emergence signatures for 10 of the 16 ARs studied, with lead times of 60-156 minutes. Several effects can influence both the estimated times of both HDF and flux emergence suggesting that one may need to combine Doppler and magnetic field data to get a reliable indicator of continued flux emergence.

The quasiparticle propagation away from the track of a highly ionizing particle in superfluid helium at low temperatures has previously been shown to exhibit anisotropy. We discuss the mechanism responsible for this behavior and show that it occurs for nuclear scattering by dark matter for recoil energies down to a few keV, and perhaps lower. This makes it possible to extend WIMP searches with interaction cross sections that reach into the neutrino floor in a meaningful energy range.

CPD-54 810 is a double-lined detached eclipsing binary containing two mid-F type dwarfs on an eccentric 26-day orbit. We perform a combined analysis of the extensive photometry obtained by the TESS space mission along with previously published observations to obtain a full orbital and physical solution for the system. We measure the following model-independent masses and radii: M1 = 1.3094+/-0.0051 Msun, M2 = 1.0896+/-0.0034 Msun, R1 = 1.9288+/-0.0030 Rsun, and R2 = 1.1815+/-0.0037 Rsun. We employ a Bayesian approach to obtain the bolometric flux for both stars from observed magnitudes, colours, and flux ratios. These bolometric fluxes combined with the stars' angular diameters (from R1, R2 and the parallax from Gaia EDR3) lead directly to the stars' effective temperatures: Teff,1 = 6462+/-43 K, and Teff,2 = 6331+/-43 K, with an additional systematic error of 0.8% (13 K) from the uncertainty in the zero-point of the flux scale. Our results are robust against the choice of model spectra and other details of the analysis. CPD-54 810 is an ideal benchmark system that can be used to test stellar parameters measured by large spectroscopic surveys or derived from asteroseismology, and calibrate stellar models by providing robust constraints on the measured parameters. The methods presented here can be applied to many other detached eclipsing binary systems to build a catalogue of well-measured benchmark stars.

It is reported that the Large High Altitude Air Shower Observatory (LHAASO) observed very high energy photons from GRB 20221009A, with the highest energy reaching 18 TeV. We find that observation of such high energy photons is a quite nontrivial fact since extragalactic background light could absorb these photons severely and the flux is too weak to be observed. Therefore we discuss the potential new mechanism for us to observe these photons, and suggest that Lorentz invariance violation induced threshold anomaly of the process (\gamma\gamma\to e^-e^+) provides a candidate to explain this phenomenon.

The LIGO-Virgo gravitational wave detectors have confidently observed 4 events involving neutron stars: two binary neutron star (BNS) mergers (GW170817 and GW190425), and two neutron star-black hole mergers (GW200105 and GW200115). However, our theoretical understanding of the remnant properties of such systems is incomplete due to the complexities related to the modeling of matter effects and the very high computational cost of corresponding numerical relativity simulations. An important such property is the recoil velocity, which is imparted onto the remnant due to the anisotropic emission of gravitational radiation and the dynamical ejection of matter in the kilonova. In this work, we combine gravitational radiation as well as dynamical ejecta distributions, computed by the Computational Relativity numerical simulations, to get accurate estimates for BNS remnant recoil velocities. We find that recoils due to ejection of matter dominate those caused by gravitational wave emission. Knowledge of BNS remnant recoil velocities is important in determining if the remnant is retained by its environment for future hierarchical mergers which, in turn, can form binaries with black holes in the so-called lower mass gap of 3 to 5 solar masses.

We present a new open-source data-reduction pipeline to reconstruct spectral data cubes from raw SPHERE integral-field spectrograph (IFS) data. The pipeline is written in Python and based on the pipeline that was developed for the CHARIS IFS. It introduces several improvements to SPHERE data analysis that ultimately produce significant improvements in postprocessing sensitivity. We first used new data to measure SPHERE lenslet point spread functions (PSFs) at the four laser calibration wavelengths. These lenslet PSFs enabled us to forward-model SPHERE data, to extract spectra using a least-squares fit, and to remove spectral crosstalk using the measured lenslet PSFs. Our approach also reduces the number of required interpolations, both spectral and spatial, and can preserve the original hexagonal lenslet geometry in the SPHERE IFS. In the case of least-squares extraction, no interpolation of the data is performed. We demonstrate this new pipeline on the directly imaged exoplanet 51 Eri b and on observations of the hot white dwarf companion to HD 2133. The extracted spectrum of HD 2133B matches theoretical models, demonstrating spectrophotometric calibration that is good to a few percent. Postprocessing on two 51 Eri b data sets demonstrates a median improvement in sensitivity of 80% and 30% for the 2015 and 2017 data, respectively, compared to the use of cubes reconstructed by the SPHERE Data Center. The largest improvements are seen for poorer observing conditions. The new SPHERE pipeline takes less than three minutes to produce a data cube on a modern laptop, making it practical to reprocess all SPHERE IFS data.

We investigate the synergy of upcoming galaxy surveys and gravitational wave (GW) experiments in constraining late-time cosmology, examining the cross-correlations between the weak lensing of gravitational waves (GW-WL) and the galaxy fields. Without focusing on any specific GW detector configuration, we benchmark the requirements for the high-precision measurement of cosmological parameters by considering several scenarios, varying the number of detected GW events and the uncertainty on the inference of the source luminosity distance and redshift. We focus on $\Lambda$CDM and scalar-tensor cosmologies, using the Effective Field Theory formalism as a unifying language. We find that, in some of the explored setups, GW-WL contributes to the galaxy signal by doubling the accuracy on non-$\Lambda$CDM parameters, allowing in the most favourable scenarios to reach even percent and sub-percent level bounds. Though the most extreme cases presented here are likely beyond the observational capabilities of currently planned individual GW detectors, we show nonetheless that - provided that enough statistics of events can be accumulated - GW-WL offers the potential to become a cosmological probe complementary to LSS surveys, particularly for those parameters that cannot be constrained by other GW probes such as standard sirens.

This work reports the detection of a multi peaked colour pattern in the integrated colours distribution of globular clusters associated to the giant elliptical galaxy NGC 4486, using Next Generation Virgo Survey data. This feature is imprinted on the well known bimodal colour distribution of these clusters. Remarkably, the pattern is similar to that found in previous works based on photometry from the HST Advanced Camera Virgo Survey, in less massive Virgo galaxies. This characteristic can be traced up to to 45 arcmin (217 Kpc) in galactocentric radius. This suggests that globular cluster formation in Virgo has been regulated, at least partially, by a collective process composed by several discrete events, working on spatial scales comparable to the size of the galaxy cluster. Furthermore, the presence of a similar colour pattern in NGC 5128, at the outskirsts of the Virgo Super-cluster, poses an intriguing question about the spatial scale of the phenomenon. The nature of the process that leads to the colour pattern is unknown. However, energetic events connected with galaxy or sub-galaxy cluster mergers and SMBH activity, in the early Universe, appear as possible candidates to explain an eventual enhancement/quenching of the globular clusters formation, reflected in the modulation of their integrated colours. Such events, presumably, may also have had an impact on the whole star formation history in Virgo galaxies.

Future Gaia and Legacy Survey of Space and Time data releases, together with wide area spectroscopic surveys, will deliver large samples of resolved binary stars with phase space coordinates, albeit with low-cadence. Given an eccentricity law $f(\epsilon)$, we derive properties of (i) the velocity distribution $v/\sqrt{G M/r}$ normalised by the value for a circular orbit at the measured separation $r$; (ii) the astrometric acceleration distributions $a/\left(G M/r^2\right)$ again normalised to the circular orbit value. Our formulation yields analytic predictions for the full statistical distribution for some commonly used eccentricity laws, if the timescale of data-sampling is comparable to or exceeds the binary period. In particular, the velocity distribution for the linear eccentricity law is surprisingly simple. With Bayesian analysis, we suggest a method to infer the eccentricity distribution based on the measured velocity distribution.

Massive colliding-wind binaries that host a Wolf-Rayet (WR) star present a potentially important source of dust and chemical enrichment in the interstellar medium (ISM). However, the chemical composition and survival of dust formed from such systems is not well understood. The carbon-rich WR (WC) binary WR~140 presents an ideal astrophysical laboratory for investigating these questions given its well-defined orbital period and predictable dust-formation episodes every 7.93 years around periastron passage. We present observations from our Early Release Science program (ERS1349) with the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) Medium-Resolution Spectrometer (MRS) and Imager that reveal the spectral and spatial signatures of nested circumstellar dust shells around WR~140. MIRI MRS spectroscopy of the second dust shell and Imager detections of over 17 shells formed throughout the past $\gtrsim130$ years confirm the survival of carbonaceous dust grains from WR~140 that are likely carriers of "unidentified infrared" (UIR)-band features at 6.4 and 7.7 $\mu$m. The observations indicate that dust-forming WC binaries can enrich the ISM with organic compounds and carbonaceous dust.

Finite temperature effects in the Standard Model tend to restore the electroweak symmetry in the early universe, but new fields coupled to the higgs field may as well reverse this tendency, leading to the so-called electroweak symmetry non-restoration (EW SNR) scenario. Previous works on EW SNR often assume that the reversal is due to the thermal fluctuations of new fields with negative quartic couplings to the higgs, and they tend to find that a large number of new fields are required. We observe that EW SNR can be minimally realized if the field(s) coupled to the higgs field develop(s) a stable condensate. We show that one complex scalar field with a sufficiently large global-charge asymmetry can develop a condensate as an outcome of thermalization and keep the electroweak symmetry broken up to temperatures well above the electroweak scale. In addition to providing a minimal benchmark model, our work hints on a class of models involving scalar condensates that yield electroweak symmetry non-restoration in the early universe.

The detection of gravitational waves has opened a new era for astronomy, allowing for the combined use of gravitational wave and electromagnetic emissions to directly probe the physics of compact objects, still poorly understood. So far, the theoretical modelling of these sources has mainly relied on standard numerical techniques as grid-based methods or smoothed particle hydrodynamics, with only a few recent attempts at using new techniques as moving-mesh schemes. Here, we introduce a general relativistic extension to the mesh-less hydrodynamic schemes in the code GIZMO, which benefits from the use of Riemann solvers and at the same time perfectly conserves angular momentum thanks to a generalised leap-frog integration scheme. We benchmark our implementation against many standard tests for relativistic hydrodynamics, either in one or three dimensions, and also test the ability to preserve the equilibrium solution of a Tolman-Oppenheimer-Volkoff compact star. In all the presented tests, the code performs extremely well, at a level at least comparable to other numerical techniques.

We combine amortized neural posterior estimation with importance sampling for fast and accurate gravitational-wave inference. We first generate a rapid proposal for the Bayesian posterior using neural networks, and then attach importance weights based on the underlying likelihood and prior. This provides (1) a corrected posterior free from network inaccuracies, (2) a performance diagnostic (the sample efficiency) for assessing the proposal and identifying failure cases, and (3) an unbiased estimate of the Bayesian evidence. By establishing this independent verification and correction mechanism we address some of the most frequent criticisms against deep learning for scientific inference. We carry out a large study analyzing 42 binary black hole mergers observed by LIGO and Virgo with the SEOBNRv4PHM and IMRPhenomXPHM waveform models. This shows a median sample efficiency of $\approx 10\%$ (two orders-of-magnitude better than standard samplers) as well as a ten-fold reduction in the statistical uncertainty in the log evidence. Given these advantages, we expect a significant impact on gravitational-wave inference, and for this approach to serve as a paradigm for harnessing deep learning methods in scientific applications.

The flux of cosmic rays in the heliosphere is subjected to variations that are related to the Sun's magnetic activity. To study this effect, updated time series of multichannel observations are needed. Here we present a web application that collects real-time data on solar activity proxies, interplanetary plasma parameters, and charged cosmic-ray data. The data are automatically retrieved on daily basis from several space missions or observatories. With this application, the data can be visualized and download into a common format. Along with observational data, the application aims to provide real-time calculations for the solar modulation of cosmic rays in the heliosphere.

The physics of de Sitter space is essential to our understanding of our cosmological past, present, and future. It forms the foundation for the statistical predictions of inflation in terms of quantum vacuum fluctuations that are being tested with cosmic surveys. In addition, the current expansion of the universe is dominated by an apparently constant vacuum energy and we again find our universe described by a de Sitter epoch. Despite the success of our predictions for cosmological observables, conceptual questions of the nature of de Sitter abound and are exacerbated by technical challenges in quantum field theory and perturbative quantum gravity in curved backgrounds. In recent years, significant process has been made using effective field theory techniques to tame these breakdowns of perturbation theory. We will discuss how to understand the long-wavelength fluctuations produced by accelerating cosmological backgrounds and how to resolve both the UV and IR obstacles that arise. Divergences at long wavelengths are resummed by renormalization group (RG) flow in the EFT. For light scalar fields, the RG flow manifests itself as the stochastic inflation formalism. In single-field inflation, long-wavelength metric fluctuations are conserved outside the horizon to all-loop order, which can be understood easily in EFT terms from power counting and symmetries.

Dual-phase xenon time projection chamber (TPC) detectors have demonstrated superior search sensitivities to dark matter over a wide range of particle masses. To extend their sensitivity to include low-mass dark matter interactions, it is critical to characterize both the light and charge responses of liquid xenon to sub-keV nuclear recoils. In this work, we report a new nuclear recoil calibration in the LUX detector $\textit{in situ}$ using neutron events from a pulsed Adelphi Deuterium-Deuterium neutron generator. We demonstrate direct measurements of light and charge yields down to 0.45 keV (1.4 scintillation photons) and 0.27 keV (1.3 ionization electrons), respectively, approaching the physical limit of liquid xenon detectors. We discuss the implication of these new measurements on the physics reach of dual-phase xenon TPCs for nuclear-recoil-based low-mass dark matter detection.

Quantum treatment of physical reference frame leads to the Ricci flow of quantum spacetime, which is a quite rigid framework to quantum and renormalization effect of gravity. The theory has a low characteristic energy scale described by a unique constant: the critical density of the universe. At low energy long distance (cosmic or galactic) scale, the theory modifies Einstein's gravity which naturally gives rise to a cosmological constant as a counter term of the Ricci flow at leading order and an effective scale dependent Einstein-Hilbert action. In the weak and static gravity limit, the framework naturally contains a MOdified Newtonian Dynamics (MOND)-like theory first proposed by Milgrom. When local curvature is large, Newtonian gravity is recovered. When local curvature is lower or comparable with the asymptotic background curvature corresponding to the characteristic energy scale, the baryonic Tully-Fisher relation can be obtained. For intermediate general curvature, the interpolating Lagrangian function gives a similar curve to the observed radial acceleration relation of galaxies. The critical acceleration constant $a_{0}$ introduced in MOND is related to the low characteristic energy scale of the theory. The cosmological constant gives a universal leading order contribution to $a_{0}$ and the flow effect gives the next order scale dependent contribution beyond MOND, which equivalently induces the "cold dark matter" to the theory. $a_{0}$ is consistent with galaxian data when the "dark matter" is about 5 times the baryonic matter. The Ricci flow of quantum spacetime is proposed as a possible underlying theory of MOND and related acceleration discrepancy phenomenon at long distance scale.

The leading order Green's function solutions of the general relativistic thin disc equations are computed, using a pseudo-Newtonian potential and asymptotic Laplace mode matching techniques. This solution, valid for a vanishing ISCO stress, is constructed by ensuring that it reproduces the leading order asymptotic behaviour of the near-ISCO, Newtonian, and global WKB limits. Despite the simplifications used in constructing this solution, it is typically accurate, for all values of the Kerr spin parameter $a$ and at all radii, to less than a percent of the full numerically calculated solutions of the general relativistic disc equations. These solutions will be of use in studying time-dependent accretion discs surrounding Kerr black holes.

Likelihood-to-evidence ratio estimation is usually cast as either a binary (NRE-A) or a multiclass (NRE-B) classification task. In contrast to the binary classification framework, the current formulation of the multiclass version has an intrinsic and unknown bias term, making otherwise informative diagnostics unreliable. We propose a multiclass framework free from the bias inherent to NRE-B at optimum, leaving us in the position to run diagnostics that practitioners depend on. It also recovers NRE-A in one corner case and NRE-B in the limiting case. For fair comparison, we benchmark the behavior of all algorithms in both familiar and novel training regimes: when jointly drawn data is unlimited, when data is fixed but prior draws are unlimited, and in the commonplace fixed data and parameters setting. Our investigations reveal that the highest performing models are distant from the competitors (NRE-A, NRE-B) in hyperparameter space. We make a recommendation for hyperparameters distinct from the previous models. We suggest a bound on the mutual information as a performance metric for simulation-based inference methods, without the need for posterior samples, and provide experimental results.

In this paper, we investigate the thermalization of Hawking radiation from primordial black holes (PBHs) in the early Universe, taking into account the interference effect on thermalization of high energy particles, known as Landau-Pomeranchuk-Migdal (LPM) effect. Small PBHs with masses $ \lesssim 10^9 \, \mathrm{g} $ completely evaporate before the big bang nucleosynthesis (BBN). The Hawking radiation emitted from these PBHs heats up the ambient plasma with temperature lower than the Hawking temperature, which results in a non-trivial temperature profile around the PBHs, namely a hot spot surrounding a PBH with a broken power-law tail. We find that the hot spot has a core with a radius much larger than the black hole horizon and its highest temperature is independent of the initial mass of the PBH such as $2 \times 10^{9} \, {\rm GeV} \times (\alpha/0.1)^{19/3}$, where $\alpha$ generically represents the fine-structure constants. We also briefly discuss the implications of the existence of the hot spot for phenomenology.

I discuss how one can apply the covariant formalism developed by Vilkovisky and DeWitt to obtain frame invariant fifth force calculations for scalar-tensor theories. Fifth forces are severely constrained by astrophysical measurements. It was shown previously that for scale-invariant Higgs-dilaton gravity, in a particular choice of Jordan frame, the dilaton fifth force is dramatically suppressed, evading the observational constraints. Using a geometric approach I extend this result to all frames, and show that the usual dichotomy of "Jordan frame" versus "Einstein frame" is better understood as a continuum of frames: submanifold slices of a more general field space.

CP violation and the violation of baryon-minus-lepton number B-L do not necessarily have to occur simultaneously in order to accomplish successful leptogenesis. Instead, it suffices if new CP-violating interactions at high energies result in primordial charge asymmetries, which are then reprocessed into a nonvanishing B-L asymmetry by right-handed neutrinos (RHNs) at lower energies. In this paper, we study this novel mechanism known as "wash-in leptogenesis", utilizing axion inflation as the source of high-scale CP violation. We specifically consider axion inflation coupled to the Standard Model hypercharge sector, which results in the dual production of hypermagnetic helicity and fermionic charge asymmetries. Although the survival of these charges is endangered by sphaleron processes, magnetic diffusion, and the chiral plasma instability, we find a large range of viable scenarios. We consistently account for RHN flavor effects and coherence among the Standard Model lepton flavors across a wide range of RHN masses. We find a lower bound of 10^(5...9) GeV on the mass of the lightest RHN involved in wash-in leptogenesis, depending on the onset of turbulence in the chiral plasma and the Hubble scale of inflation. Our model is representative of a broader class of new leptogenesis scenarios and suggests interesting observational signatures with regard to intergalactic magnetic fields, primordial black holes, and gravitational waves.