Papers for Friday, Dec 24 2021

Mini-EUSO is a mission of the JEM-EUSO program flying onboard the International Space Station since August 2019. Since the first data acquisition in October 2019, more than 35 sessions have been performed for a total of 52 hours of observations. The detector has been observing Earth at night-time in the UV range and detected a wide variety of transient sources all of which have been modelled through Monte Carlo simulations. Mini-EUSO is also capable of detecting meteors and potentially space debris and we performed simulations for such events to estimate their impact on future missions for cosmic ray science from space. We show here examples of the simulation work done in this framework to analyse the Mini-EUSO data. The expected response of Mini-EUSO with respect to ultra high energy cosmic ray showers has been studied. The efficiency curve of Mini-EUSO as a function of primary energy has been estimated and the energy threshold for Cosmic Rays has been placed to be above 10^{21} eV. We compared the morphology of several transient events detected during the mission with cosmic ray simulations and excluded that they can be due to cosmic ray showers. To validate the energy threshold of the detector, a system of ground based flashers is being used for end-to-end calibration purposes. We therefore implemented a parameterisation of such flashers into the JEM-EUSO simulation framework and studied the response of the detector with respect to such sources.

Astrometric microlensing is a unique tool to measure stellar masses. It allows us to determine the mass of the lensing star with an accuracy of a few per cent. In this paper, we update, extend, and refine our predictions of astrometric-microlensing events based on Gaia's early Data release 3 (eDR3). We selected about 500.000 high-proper-motion stars from Gaia eDR3 with $\mu_{tot}>100\,\mathrm{mas/yr}$ and searched for background sources close to their paths. We applied various selection criteria and cuts in order to exclude spurious sources and co-moving stars. By forecasting the future positions of lens and source we determined epoch of and angular separation at closest approach, and determined an expected positional shift and magnification. Using Gaia~eDR3, we predict 1758 new microlensing events with expected shifts larger than 0.1 mas between the epochs J2010.5 and mid J2066.0. Further we provide more precise information on the angular separation at closest approach for 3084 previously predicted events. This helps to select better targets for observations, especially for events which occur within the next decade. Our search lead to the new prediction of an interesting astrometric-microlensing event by the white dwarf Gaia eDR3-4053455379420641152. In 2025 it will pass by a $G=20.25\,\mathrm{mag}$ star, which will lead to a positional shift of the major image of $\delta\theta_{+}=1.2^{+2.0}_{-0.5}\,\mathrm{mas}$. Since the background source is only $\Delta G=2.45\,\mathrm{mag}$ fainter than the lens, also the shift of the combined center of light will be measurable, especially using a near infrared filter, where the background star is brighter than the lens $\Delta Ks=-1.1\,\mathrm{mag}$

Main-sequence stars with convective envelopes often appear larger and cooler than predicted by standard models of stellar evolution for their measured masses. This is believed to be caused by stellar activity. In a recent study, accurate measurements have been published for the K-type components of the 1.62 day detached eclipsing binary EPIC 219511354, showing the radii and temperatures for both stars to be affected by these discrepancies. This is a rare example of a system in which the age and chemical composition are known, by virtue of being a member of the well-studied open cluster Ruprecht 147 (age $\sim$ 3 Gyr, [Fe/H] = +0.10). Here we report a detailed study of this system with non-standard models incorporating magnetic inhibition of convection. We show that these calculations are able to reproduce the observations largely within their uncertainties, providing robust estimates of the strength of the magnetic fields on both stars: $1600 \pm 130$ G and $1830 \pm 150$ G for the primary and secondary, respectively. Empirical estimates of the magnetic field strengths based on the measured X-ray luminosity of the system are roughly consistent with these predictions, supporting this mechanism as a possible explanation for the radius and temperature discrepancies.

Based on calculations of the Stokes parameters for the Holweger-Mliller model atmosphere, we study sensitivity of the Fe I 525.02 nm line to some kinds of vertical and horizontal magnetic field inhomogeneity. A noticeable asymmetry is shown to appear in the V profile peaks when the vertical gradient is -0.4 mT/km, which is typical of some theoretical flux tube models. The asymmetry is most pronounced in a pure longitudinal magnetic field and at a low macroturbulent velocity. A similar effect is observed for the Q profile in nonlongitudinal fields as well. The Fe I 525.02 nm line is sensitive also to subtelescopic fields of mixed polarity like those observed by Stenflo in IR lines. We argue that the Wilson depression in small-scale flux tubes renders strong-field areas invisible at heliocentric angles greater than 60--65 degrees, since they are screened by surroundings with weaker magnetic fields.

We perform empirical fits to the \emph{Chandra} and \emph{XMM-Newton} spectra of three ultraluminous X-ray sources (ULXs) in the edge-on spiral galaxy NGC 891, monitoring the region over a seventeen year time window. One of these sources has been visible since the early 1990s with \emph{ROSAT} and has been observed multiple times with \emph{Chandra} and \emph{XMM-Newton}. Another has been visible since 2011. We build upon prior analyses of these sources by analyzing all available data at all epochs. Where possible \emph{Chandra} data is used, since its superior spatial resolution allows for more effective isolation of the emission from each individual source, thus providing a better determination of their spectral properties. We also identify a new transient ULX, CXOU J022230.1+421937, which faded from view over the course of a two month period from Nov 2016 to Jan 2017. Modeling of each source at every epoch was conducted using six different models ranging from thermal bremsstrahlung to accretion disk models. Unfortunately, but as is common with many ULXs, no single model yielded a much better fit than the others. The two known sources had unabsorbed luminosities that remained fairly consistent over five or more years. Various possibilities for the new transient ULX are explored.

Dust permeates the interstellar medium, reddening and polarizing background starlight, but dust properties vary with local environment. In order to characterize dust in a highly irradiated diffuse cloud, we measure the reddening and optical polarization towards 27 stars surrounding the mid-latitude $b$=$+$24$^{\circ}$ O9.2IV star $\zeta$ Ophiuchi, using new optical spectroscopy and polarimetry. We incrementally deredden and depolarize with distance, allowing us to distinguish dust components along these sightlines. The data indicate three distinct dust populations: a foreground component characteristic of average Milky Way dust ($R_{\rm{V}}$$\approx$3.1, $d$$\lesssim$180 pc), a highly polarizing mid-distance component in the vicinity of $\zeta$ Oph ($R_{\rm{V}}$$\approx$2.4, 200 pc$<$$d$$<$300 pc), and a non-polarizing distant component ($R_{\rm{V}}$$\approx$3.6, 600 pc$<$$d$$<$2000 pc). Prominent 8 $\mu$m infrared striations spanning the field of view likely have high Polycyclic Aromatic Hydrocarbon content and are illuminated by $\zeta$ Oph. Foreground-subtracted polarizations roughly align with these striations, which, we argue, lie immediately behind $\zeta$ Oph and constitute the highly-polarizing mid-distance dust. This component polarizes very efficiently ($P_{\rm{V}}$$>$9.1$E(B-V)$), implying a high degree of grain alignment and suggesting that the bulk of the polarization occurs in a small fraction of the volume. The large $R_{\rm{V}}$ in the distant component reveals that dust above the Galactic Plane ($z$$>$250 pc) may contain a greater fraction of large grains than the Milky Way average.

Recent work has shown that aside from the classical view of collisions by increasingly massive planetesimals, the accretion of mm- to m-sized 'pebbles' can also reproduce the mass-orbit distribution of the terrestrial planets. Here, we perform N-body simulations to study the effects of pebble accretion onto growing planetesimals of different diameters located in the inner Solar System. The simulations are run to occur during the lifetime of the gas disc while also simultaneously taking Jupiter's growth into account. We find that pebble accretion can increase the mass in the solid disc by at least a few times its initial mass with reasonable assumptions that pebbles fragment to smaller-sized grains at the snow line and that gas-disc-induced orbital migration effects are in force. Such a large contribution in mass by pebbles would seem to imply that the isotopic composition of the inner Solar System should be similar to the pebble source (i.e. outer Solar System). This implication appears to violate the observed nucleosynthetic isotopic dichotomy of the sampled Solar System. Thus, pebble accretion played little or no role in terrestrial planet formation.

The collimation of relativistic jets launched from the vicinity of supermassive black holes (SMBHs) at the centers of active galactic nuclei (AGN) is one of the key questions to understand the nature of AGN jets. However, little is known about the detailed jet structure for AGN like quasars since very high angular resolutions are required to resolve these objects. We present very long baseline interferometry (VLBI) observations of the archetypical quasar 3C 273 at 86 GHz, performed with the Global Millimeter VLBI Array, for the first time including the Atacama Large Millimeter/submillimeter Array. Our observations achieve a beam-formed angular resolution down to $\sim$60 ${\rm \mu}$as, resolving the innermost part of the jet ever on scales of $\sim 10^5$ Schwarzschild radii. Our observations, including close-in-time High Sensitivity Array observations of 3C 273 at 15, 22, and 43 GHz, suggest that the inner jet collimates parabolically, while the outer jet expands conically, similar to jets from other nearby low luminosity AGN. We discovered the jet collimation break at $\sim8\times10^{6}$ Schwarzschild radii, providing the first compelling evidence for structural transition in a quasar jet. The location of the collimation break for 3C 273 is farther downstream the sphere of gravitational influence (SGI) from the central SMBH. With the results for other AGN jets, our results show that the end of the collimation zone in AGN jets is governed not only by the SGI of the SMBH but also by the more diverse properties of the central nuclei.

The way supermassive black holes (SMBH) in galactic centers accumulate their mass is not completely determined. At large scales, it is governed by galactic encounters, mass inflows connected to spirals arms and bars, or due to expanding shells from supernova (SN) explosions in the central parts of galaxies. The investigation of the latter process requires an extensive set of gas dynamical simulations to explore the muti-dimensional parameter space needed to frame the phenomenon. The aims of this paper are to extend our investigation of the importance of supernovae for inducing accretion onto a SMBH and carry out a comparison between the fully hydrodynamic code Flash and the much less computationally intensive code Ring, which uses the thin shell approximation. We simulate 3D expanding shells in a gravitational potential similar to that of the Galactic Center with a variety of homogeneous and turbulent environments. In homogeneous media, we find convincing agreement between Flash and Ring in the shapes of shells and their equivalent radii throughout their whole evolution until they become subsonic. In highly inhomogeneous, turbulent media, there is also a good agreement of shapes and sizes of shells, and of the times of their first contact with the central 1 pc sphere, where we assume that they join the accretion flow. The comparison supports the proposition that a SN occurring at a galactocentric distance of 5 pc typically drives 1 - 3 $M_\odot$ into the central 1 pc around the galactic center.

We present stellar mass fractions and composite luminosity functions (LFs) for a sample of \Ncl\ clusters from the Massive and Distant Clusters of WISE Survey (MaDCoWS) at a redshift range of $0.951 \leq z \leq 1.43$. Using SED fitting of optical and deep mid-infrared photometry, we establish the membership of objects along the lines-of-sight to these clusters and calculate the stellar masses of member galaxies. We find stellar mass fractions for these clusters largely consistent with previous works, including appearing to display a negative correlation with total cluster mass. We measure a composite $3.6~\mathrm{\mu m}$ LF down to $m^*+2.5$ for all 12 clusters. Fitting a Schechter function to the LF, we find a characteristic $3.6~\mathrm{\mu m}$ magnitude of $m^*=19.83\pm0.12$ and faint-end slope of $\alpha=-0.81\pm0.10$ for the full sample at a mean redshift of $\bar{z} = 1.18$. We also divide the clusters into high- and low-redshift bins at $\bar{z}=1.29$ and $\bar{z}=1.06$ respectively and measure a composite LF for each bin. We see a small, but statistically significant evolution in $m^*$ and $\alpha$ -- consistent with passive evolution -- when we study the joint fit to the two parameters, which is probing the evolution of faint cluster galaxies at $z\sim1$. This highlights the importance of deep IR data in studying the evolution of cluster galaxy populations at high-redshift.

We present the blind Westerbork Coma Survey probing the HI content of the Coma galaxy cluster with the Westerbork Synthesis Radio Telescope. The survey covers the inner $\sim$ 1 Mpc around the cluster centre, extending out to 1.5 Mpc towards the south-western NGC 4839 group. The survey probes the atomic gas in the entire Coma volume down to a sensitivity of $\sim$ 10$^{19}$ cm$^{-2}$ and 10$^8$ M$_{\odot}$. Combining automated source finding with source extraction at optical redshifts and visual verification, we obtained 40 HI detections of which 24 are new. Over half of the sample displays perturbed HI morphologies indicative of an ongoing interaction with the cluster environment. With the use of ancillary UV and mid-IR, data we measured their stellar masses and star formation rates and compared the HI properties to a set of field galaxies spanning a similar stellar mass and star formation rate range. We find that $\sim$ 75 \% of HI-selected Coma galaxies have simultaneously enhanced star formation rates (by $\sim$ 0.2 dex) and are HI deficient (by $\sim$ 0.5 dex) compared to field galaxies of the same stellar mass. According to our toy model, the simultaneous HI deficiency and enhanced star formation activity can be attributed to either HI stripping of already highly star forming galaxies on a very short timescale, while their H$_2$ content remains largely unaffected, or to HI stripping coupled to a temporary boost of the HI-to-H$_2$ conversion, causing a brief starburst phase triggered by ram pressure before eventually quenching the galaxy.

The Wavelength-shifting Optical Module (WOM) is a novel photosensor concept for the instrumentation of large detector volumes with single-photon sensitivity. The key objective is to improve the signal-to-noise ratio which is achieved by decoupling the photosensitive area of a sensor from the cathode area of its photomultiplier tube (PMT). The WOM consists of a transparent tube with two PMTs attached to its ends. The tube is coated with wavelength-shifting paint absorbing ultra-violet photons with nearly $100\,\%$ efficiency. Depending on the environment, e.g. air (ice), up to $73\,\%$ $(41\,\%)$ of the subsequently emitted optical photons can be captured by total internal reflection and propagate towards the PMTs where they are recorded. The optical properties of the paint, the geometry of the tube and the coupling of the tube to the PMTs have been optimized for maximal sensitivity based on theoretical derivations and experimental evaluations. Prototypes were built to demonstrate the technique and to develop a reproducible construction process. Important measurable characteristics of the WOM are the wavelength dependent effective area, the transit time spread of detected photons and the signal-to-noise ratio. The WOM outperforms bare PMTs especially with respect to the low signal-to-noise ratio with an increase of a factor up to 8.9 in air (5.2 in ice). Since the gain in sensitivity is mostly in the UV-regime, the WOM is an ideal sensor for Cherenkov and scintillation detectors.

Regular large-scale polarimetric observations in Crimean astrophysical observatory began in the early 1960s. In 2002 - 2017 the single-channel aperture photometer-polarimeter with a quarter-wave plate at the 2.6-m Shajn mirror telescope (SMT) was used. We accumulated a large homogeneous data set of polarimetric observations of different types of objects that are to be published separately. Correct polarimetric data processing requires high polarization standards and zero-polarization stars. We aim to improve the data reduction and calibration process to obtain further results with highest possible accuracy. High time resolution broad-band (WR, R, V, B, U) polarization observations are made of 98 known standard stars (527 time series with total duration about 184 hours). We determined values of linear and circular polarization for 98 nearby Northern bright stars. This catalogue is not compilative, but obtained using the same instrument and technique during large time interval. It will be used for our future research and it may be used by other authors. We implemented the least squares approach for determination of the Stokes parameters. It allowed us to obtain results with the accuracy better then we obtained using previously used methods. We report suspicious or variable stars that are not suitable as standards for high precision polarimetry.

The total event number of fast radio bursts (FRBs) is accumulating rapidly with the improvement of existing radio telescopes and the completion of new facilities. Especially, the Five-hundred-meter Aperture Spherical radio Telescope (FAST) Collaboration has just reported more than one thousand bursts in a short observing period of 47 days \citep{LiD2021}. The interesting bimodal distribution in their work motivates us to revisit the definition of FRBs. In this work, we ascribe the bimodal distribution to two physical kinds of radio bursts, which may have different radiation mechanisms. We propose to use brightness temperature to separate two subtypes. For FRB 20121102A, the critical brightness temperature is $T_{\rm B,cri}\simeq10^{33}\,\rm K$. Bursts with $T_{\rm B}\geq T_{\rm B,cri}$ are denoted as "classical" FRBs, and further we find a tight pulse width-fluence relation ($T\propto\mathcal{F_\nu}^{0.306}$) for them. On the contrary, the other bursts are considered as "atypical" bursts that may originate from a different physical process. We suggest that for each FRB event, a similar dividing line should exist but $T_{\rm B,cri}$ is not necessarily the same. Its exact value depends on FRB radiation mechanism and properties of the source.

The surface temperature distributions of central compact objects (CCOs) are powerful probes of their crustal magnetic field strengths and geometries. Here we model the surface temperature distribution of RX J0822$-$4300, the CCO in the Puppis A supernova remnant (SNR), using $471$ ks of XMM-Newton data. We compute the energy-dependent pulse profiles in sixteen energy bands, fully including the general relativistic effects of gravitational redshift and light bending, to accurately model the two heated surface regions of different temperatures and areas, in addition to constraining the viewing geometry. This results in precise measurements of the two temperatures: $kT_{\rm warm} = (1+z) \times 0.222_{-0.019}^{+0.018}$ keV and $kT_{\rm hot} = (1+z) \times 0.411\pm0.011$ keV. For the first time, we are able to measure a deviation from a pure antipodal hot-spot geometry, with a minimum value of $1.\!^{\circ}1 \pm 0.\!^{\circ}2$, and an expectation value of $9.\!^{\circ}35 \pm 0.\!^{\circ}17$ among the most probable geometries. The discovery of this asymmetry, along with the factor of $\approx2$ temperature difference between the two emitting regions, may indicate that RX J0822$-$4300 was born with a strong, tangled crustal magnetic field.

Combining data from the ACS Virgo Cluster Survey (ACSVCS) and the Next Generation Virgo cluster Survey (NGVS), we extend previous studies of color gradients of the globular cluster (GC) systems of the two most massive galaxies in the Virgo cluster, M87 and M49, to radii of $\sim 15~R_e$ ($\sim 200$ kpc for M87 and $\sim 250$ kpc for M49). We find significant negative color gradients, i.e., becoming bluer with increasing distance, out to these large radii. The gradients are driven mainly by the outwards decrease of the ratio of red to blue GC numbers. The color gradients are also detected out to $\sim 15~R_e$ in the red and blue sub-populations of GCs taken separately. In addition, we find a negative color gradient when we consider the satellite low-mass elliptical galaxies as a system, i.e., the satellite galaxies closer to the center of the host galaxy usually have redder color indices, both for their stars and GCs. According to the "two phase" formation scenario of massive early-type galaxies, the host galaxy accretes stars and GCs from low-mass satellite galaxies in the second phase. So the accreted GC system naturally inherits the negative color gradient present in the satellite population. This can explain why the color gradient of the GC system can still be observed at large radii after multiple minor mergers.

A scheme of Bayesian rotation inversion, which allows us to compute the probability of a model of a stellar rotational profile, is developed. The validation of the scheme with simple rotational profiles and the corresponding sets of artificially generated rotational shifts has been successfully carried out, and we can correctly distinguish the (right) rotational model, prepared beforehand for generating the artificial rotational shifts, with the other (wrong) rotational model. The Bayesian scheme is applied to a gamma Dor-delta Sct type hybrid star, KIC 11145123, leading to a result that the convective core of the star might be rotating much faster (~ 10 times faster) than the other regions of the star. The result is consistent with that previously suggested by Hatta et al. (2019) based on a 3-zone modeling, further strengthening their argument from a Bayesian point of view.

A photometric and spectroscopic investigation is performed on five W Ursae Majoris eclipsing binaries (EWs) J015818.6+260247 (hereinafter as J0158b), J073248.4+405538 (hereinafter as J0732), J101330.8+494846 (hereinafter as J1013), J132439.8+130747 (hereinafter as J1324) and J152450.7+245943 (hereinafter as J1524). The photometric data are collected with the help of the 1.3\,m Devasthal Fast Optical Telescope (DFOT), the 1.04\,m Sampurnanand Telescope (ST) and the TESS space mission. The low-resolution spectra of the 4\,m Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) are used for spectroscopic analysis. The orbital period change of these systems is determined using our and previously available photometric data from different surveys. The orbital period of J1013 and J1524 is changing with a rate of $-2.552(\pm0.249)\times 10^{-7}$ days $yr^{-1}$ and $-6.792(\pm0.952)\times 10^{-8}$ days $yr^{-1}$, respectively, while others do not show any orbital period change. The orbital period change of J1013 and J1524 corresponds to a mass transfer rate of $2.199\times10^{-7} M_{\odot}\,yr^{-1}$ and $6.151\times10^{-8}M_{\odot}\,yr^{-1}$ from the primary to the secondary component in these systems. It is likely that angular momentum loss via magnetic braking may also be responsible for the observed orbital period change in the case of J1524. All systems have a mass-ratio lower than 0.5, except J0158b with a mass-ratio of 0.71. All the systems are shallow type contact binaries. The J0158b and J1524 are A-subtype while others are W-subtype. The $H_{\alpha}$ emission line region is compared with template spectra prepared using two inactive stars with the help of STARMOD program. The J0158, J1324 and J1524 systems show excess emission in the residual spectra after subtraction of the template.

The dynamics of irregularly-shaped particles subjected to the combined effect of gas drag and radiative forces and torques in a cometary environment is investigated. The equations of motion are integrated over distances from the nucleus surface up to distances where the gas drag is negligible. The aerodynamic forces and torques are computed assuming a spherically symmetric expanding gas. The calculations are limited to particle sizes in the geometric optics limit, which is the range of validity of our radiative torque calculations. The dynamical behaviour of irregular particles is quite different to those exhibited by non-spherical but symmetric particles such as spheroids. An application of the dynamical model to comet 67P/Churyumov-Gerasimenko, the target of the Rosetta mission, is made. We found that, for particle sizes larger than about 10 micrometer, the radiative torques are negligible in comparison with the gas-driven torques up to a distance of about 100 km from the nucleus surface. The rotation frequencies of the particles depend on their size, shape, and the heliocentric distance, while the terminal velocities, being also dependent on size and heliocentric distance, show only a very weak dependence on particle shape. The ratio of the sum of the particles projected areas in the sun-to-comet direction to that of the sum of the particles projected areas in any direction perpendicular to it is nearly unity, indicating that the interpretation of the observed u-shaped scattering phase function by Rosetta/OSIRIS on comet 67P coma cannot be linked to mechanical alignment of the particles.

This paper studies the long-term migration of disturbed regolith materials on the surface of Solar System small bodies from the viewpoint of nonlinear dynamics. We propose an approximation model for secular mass movement, which combines the complex topography and irregular gravitational field. Choosing asteroid 101955 Bennu as a representative, the global change of the dynamical environment is examined, which presents a division of the creeping-sliding-shedding regions for a spun-up asteroid. In the creeping region, the dynamical equation of disturbed regolith grains is established based on the assumption of "trigger-slide" motion mode. The equilibrium points, local manifolds and large-scale trajectories of the system are calculated to clarify the dynamical characteristics of long-term regolith movement. Generally, we find for a low spin rate, the surface regolith grains flow toward the middle latitudes from the polar/equatorial regions, which is dominated by the gradient of the geopotential. While spun up to a high rate, regolith grains tend to migrate toward the equator, which happens in parallel with a topological shift of the local equilibria at low latitudes. From a long-term perspective, we find the equilibrium points dominate the global trends of regolith movements. Using the methodology developed in this paper, we give a prospect or retrospect to the secular motion of regolith materials during the spin-up process, and the results reveal a significant regulatory role of the equilibria. Through a detailed look at the dynamical scheme under different spin rates, we achieve a macro forecast of the global trends of regolith motion during the spin-up process, which explains the global geologic evolution driven by the long-term movements of regolith materials.

We examine the effectiveness of identifying distinct evolutionary histories in IllustrisTNG-100 galaxies using unsupervised machine learning with Gaussian Mixture Models on galaxy observables. We focus on how clustering compressed metallicity histories and star formation histories produces subpopulations of galaxies with distinct evolutionary properties (for both halo mass assembly and ex-situ mass fraction). By contrast, clustering with photometric colours fails to resolve such histories. We identify several populations of interest that reflect a variety of evolutionary scenarios supported by the literature. Notably, we identify a population of galaxies inhabiting the upper-red sequence, $M_{*}>10^{10} M_{\odot}$ that has a significantly higher ex-situ merger mass fraction present at fixed masses and a star formation history that has yet to fully quench, in contrast to an overlapping satellite-dominated population along the red sequence, which is distinctly quiescent. Extending the clustering to study four clusters instead of three further divides quiescent galaxies, while star forming ones are mostly constrained to a single cluster, demonstrating a variety of supported pathways to quenching. In addition to these populations, we identify a handful of populations from our other clusters that are readily applicable to observational surveys, allowing for possible extensions of this work in an observational context, and to corroborate results within the IllustrisTNG ecosystem.

In the time domain, the radio sky in particular along the Galactic plane direction may vary significantly because of various energetic activities associated with stars, stellar and supermassive black holes. Using multi-epoch Very Large Array surveys of the Galactic plane at 5.0 GHz, Becker et al. (2010) presented a catalogue of 39 variable radio sources in the flux density range 1-70 mJy. To probe their radio structures and spectra, we observed 17 sources with the very-long-baseline interferometric (VLBI) imaging technique and collected additional multi-frequency data from the literature. We detected all of the sources at 5 GHz with the Westerbork Synthesis Radio Telescope, but only G23.6644-0.0372 with the European VLBI Network (EVN). Together with its decadal variability and multi-frequency radio spectrum, we interpret it as an extragalactic peaked-spectrum source with a size of <~10 pc. The remaining sources were resolved out by the long baselines of the EVN because of either strong scatter broadening at the Galactic latitude <1 deg or intrinsically very extended structures on centi-arcsec scales. According to their spectral and structural properties, we find that the sample has a diverse nature. We notice two young H II regions and spot a radio star and a candidate planetary nebula. The rest of the sources are very likely associated with radio active galactic nuclei (AGN). Two of them also displays arcsec-scale faint jet activity. The sample study indicates that AGN are commonplace even among variable radio sources in the Galactic plane.

The standard cosmological model assumes a homogeneous and isotropic universe as the background spacetime on large scales called the cosmological principle. However, some observations suggest the possibility of an inhomogeneous and anisotropic universe at large scales. In this paper, we investigate a model of the universe with random inhomogeneities and anisotropies on very large scales, motivated by the supercurvature dark energy model in [Phys. Rev. D 99 103512 (2019)]. In this model, the authors introduced a scalar field with $\mathcal{O}(1)$ inhomogeneities on a scale sufficiently larger than the current horizon scale (superhorizon scale), and the potential energy of the scalar field explains the accelerating expansion, with slight deviations from the cosmological principle. We aim at clarifying the theoretical prediction on the large-scale structure (LSS) of the matter component in this model. Based on the work on the superhorizon scale fluctuations (superhorizon mode) presented in [arXiv:2111.14174], we derive the equations that the perturbative components to the LSS obey as a generalization of the cosmological perturbations theory, which is solved to find the influence of the dark energy inhomogeneities on the formation of the LSS. Finally, we show that the model can be consistent with observations by comparing the $\sigma_8$ predicted by the numerical solution of the model with the $\sigma_8$ indicated by observations such as Planck and SDSS.

The Laser Interferometer Space Antenna (LISA) aims to observe gravitational waves in the mHz regime over its 10-year mission time. LISA will operate laser interferometers between three spacecrafts. Each spacecraft will utilize independent clocks which determine the sampling times of onboard phasemeters to extract the interferometric phases and, ultimately, gravitational wave signals. To suppress limiting laser frequency noise, signals sampled by each phasemeter need to be combined in post-processing to synthesize virtual equal-arm interferometers. The synthesis in turn requires a synchronization of the independent clocks. This article reports on the experimental verification of a clock synchronization scheme down to LISA performance levels using a hexagonal optical bench. The development of the scheme includes data processing that is expected to be applicable to the real LISA data with minor modifications. Additionally, some noise coupling mechanisms are discussed.

We present a reassessment of the radial abundance gradients of He, C, N, O, Ne, S, Cl, and Ar in the Milky Way using the deep optical spectra of 42 HII regions presented in Arellano-C\'ordova et al. (2020, 2021) and M\'endez-Delgado et al. (2020) exploring the impact of: (1) new distance determinations based on Gaia EDR3 parallaxes and (2) the use of Peimbert's temperature fluctuations paradigm ($t ^ 2> 0$) for deriving ionic abundances. We find that distances based on Gaia EDR3 data are more consistent with kinematic ones based on Galactic rotation curves calibrated with radio parallaxes, which give less dispersion and uncertainties than those calibrated with spectrophotometric stellar distances. The distances based on the Gaia parallaxes --DR2 or EDR3-- eliminate the internal flattening observed in previous determinations of the Galactic gradients at smaller distances than $\sim 7$ kpc. Abundances and gradients determined assuming $ t ^ 2> 0 $ -- not only for O but also for the rest of elements -- are not affected by the abundance discrepancy problem and give elemental abundances much consistent with the solar ones for most elements. We find that our radial abundance gradient of He is consistent with the most accurate estimates of the primordial He abundance. We do not find evidence of azimuthal variations in the chemical abundances of our sample. Moreover, the small dispersion in the O gradient -- indicator of metallicity in photoionized regions -- indicate that the gas of the HII regions is well mixed in the sampled areas of the Galaxy.

The science behind galaxy interaction and mergers has a fundamental role and gives us an insight into galaxy formation and its evolution. Fluctuating angular momentum is responsible for extraordinary events like polar rings, tidal tails, and ripples. To study different phenomena related to galaxy interactions, various parameters like the mass ratio of the interacting galaxy, orbital parameters, mass distribution, morphologies are required. Convolutional Neural Networks (CNN) are widely used to classify image data. Thus, we used CNN as our approach to the problem. In this work, we will be using data from state-of-the-art magneto-hydrodynamic simulations of galaxy mergers from the GalMer database at different dynamical parameters using image snapshots of merging pairs of galaxies and feeding them to our Deep Learning model (ResNet). The dynamical parameters we are aiming for; would be spin, relative inclination ($i$), viewing angle ($\theta$), and azimuthal angle ($\phi$). We aim to download bulk data using the web scraping method. The first approach is to create different combinations of these parameters to form 60 classes. Feeding the data into the model, we achieved 93.63% accuracy. As we received good results in minute classification, we moved to our second approach, regression. Here the model can predict the continuous and exact values of the dynamical parameters. We have achieved a 99.86% R-squared value and the mean squared error of 0.0833 on testing data. In the end, we used data from Sloan Digital Sky Survey to test our trained model on some real images.

Thirty years after the discovery of the first very-high-energy {\gamma}-ray source by the Whipple telescope, the field experienced a revolution mainly driven by the third generation of imaging atmospheric Cherenkov telescopes (IACTs). The combined use of large mirrors and the invention of the imaging technique at the Whipple telescope, stereoscopic observations, developed by the HEGRA array and the fine-grained camera, pioneered by the CAT telescope, led to a jump by a factor of more than ten in sensitivity. The advent of advanced analysis techniques led to a vast improvement in background rejection, as well as in angular and energy resolutions. Recent instruments already have to deal with a very large amount of data (petabytes), containing a large number of sources often very extended (at least within the Galactic plane) and overlapping each other, and the situation will become even more dramatic with future instruments. The first large catalogues of sources have emerged during the last decade, which required numerous, dedicated observations and developments, but also made the first population studies possible. This paper is an attempt to summarize the evolution of the field towards the building up of the source catalogues, to describe the first population studies already made possible, and to give some perspectives in the context of the upcoming, new generation of instruments.

We demonstrate two new approaches that have been developed to aid the production of future hard X-ray catalogs, and specifically to reduce the reliance on human intervention during the detection of faint excesses in maps that also contain systematic noise. A convolutional neural network has been trained on data from the INTEGRAL/ISGRI telescope to create a source detection tool that is more sensitive than previous methods, whilst taking less time to apply to the data and reducing the human subjectivity involved in the process. This new tool also enables searches on smaller observation timescales than was previously possible. We show that a method based on Bayesian reasoning is better able to combine the detections from multiple observations than previous methods. When applied to data from the first 1000 INTEGRAL revolutions these improved techniques detect 25 sources (about 5% of the total sources) which were previously undetected in the stacked images used to derive the published catalog made using the same dataset.

The composition of giant planets reflects their formation history. Planetesimal accretion during the phase of planetary migration could lead to the delivery of heavy elements into giant planets. In our previous paper (Shibata et al. 2020) we showed that planetesimal accretion during planetary migration occurs in a rather narrow region of the protoplanetary disk, which we refer as "the sweet spot for accretion". The goal of this paper is to reveal the nature of the sweet spot and investigate the role of the sweet spot in determining the composition of gas giant planets. We analytically derive the required conditions for the sweet spot. Then, we compare the derived equations with the numerical simulations. We find that the conditions required for the sweet spot can be expressed by the ratio of the gas damping timescale of the planetesimal orbits and the planetary migration timescale. If the planetary migration timescale depends on the surface density of disk gas inversely, the location of the sweet spot does not change with the disk evolution. The mass of planetesimals accreted by the planet depends on the amount of planetesimals that are shepherded by mean motion resonances. Our analysis suggests that tens Earth-mass of planetesimals can be shepherded into the sweet spot without planetesimal collisions. However, as more planetesimals are trapped into mean motion resonances, collisional cascade can lead to fragmentation of planetesimals. This could affect the location of the sweet spot and the population of small objects in planetary systems. We conclude that the composition of gas giant planets depends on whether the planets crossed the sweet spot during their formation. Constraining the metallicity of cold giant planets, that are expected to be outer than the sweet spot, would reveal key information for understanding the origin of heavy elements in giant planets.

In this work, we study the effects of nuclear symmetry energy slope on neutron star dense matter equation of state and its impact on neutron star observables (mass-radius, tidal response). We construct the equation of state within the framework of covariant density functional theory implementing coupling schemes of non-linear and density-dependent models with viability of heavier non-nucleonic degrees of freedom. The slope of symmetry energy parameter ($L_{\text{sym}}$) is adjusted following density-dependence of isovector meson coupling to baryons. We find that smaller values of $L_{\text{sym}}$ at saturation favour early appearance of $\Delta$-resonances in comparison to hyperons leading to latter's threshold at higher matter densities. We also investigate the dependence of $L_{\text{sym}}$ on tidal deformability and compactness parameter of a $1.4~M_\odot$ neutron star for different equation of states and observe similar converging behaviour for larger $L_{\text{sym}}$ values.

Simulations of idealised star-forming filaments of finite length typically show core growth which is dominated by two cores forming at its respective end. The end cores form due to a strong increasing acceleration at the filament ends which leads to a sweep-up of material during the filament collapse along its axis. As this growth mode is typically faster than any other core formation mode in a filament, the end cores usually dominate in mass and density compared to other cores forming inside a filament. However, observations of star-forming filaments often do not show this prevalence of cores at the filament edges. We use numerical simulations of accreting filaments forming in a finite converging flow to explore a possible mechanism which leads to a suppression of the end cores. While such a setup still leads to end cores that soon begin to move inwards, the continued accumulation of material outside of these makes a key difference: their positions now lie within the larger filamentary structure and not at its edges. This softens their inward gravitational acceleration as they are embedded by new material further out. As a result, these two cores do not grow as fast as expected for the edge effect and thus do not dominate over other core formation modes in the filament.

Magnetic fields play an important role in the evolution of interstellar medium and star formation. As the only direct probe of interstellar field strength, credible Zeeman measurements remain sparse due to the lack of suitable Zeeman probes, particularly for cold, molecular gas. Here we report the detection of a magnetic field of $+$3.8 $\pm$ 0.3 $\mu$G through the HI narrow self-absorption (HINSA) toward L1544, a well-studied prototypical prestellar core in an early transition between starless and protostellar phases characterized by high central number density and low central temperature. A combined analysis of the Zeeman measurements of quasar HI absorption, HI emission, OH emission, and HINSA reveals a coherent magnetic field from the atomic cold neutral medium (CNM) to the molecular envelope. The molecular envelope traced by HINSA is found to be magnetically supercritical, with a field strength comparable to that of the surrounding diffuse, magnetically subcritical CNM despite a large increase in density. The reduction of the magnetic flux relative to the mass, necessary for star formation, thus seems to have already happened during the transition from the diffuse CNM to the molecular gas traced by HINSA, earlier than envisioned in the classical picture where magnetically supercritical cores capable of collapsing into stars form out of magnetically subcritical envelopes.

The recent increase in computing power and the vast amount of coming data have allowed the incursion of machine learning methods as analysis tools in observational cosmology. The main aim of this work is to introduce the artificial neural networks and some of its applications to cosmology, starting by the basic mathematical theory and followed by its computational implementation. We build a neural network from scratch with any desired architecture, to generate models from two datasets based on the Friedmann equation: one describing the function with their respective uncertainties and other from different cosmological parameters. We also trained a multilayer perceptron from solutions of the differential equations that describe the Universe components, to provide a model that is able to accurately approximate the real evolution and to reduce the computational time by 53\%. Finally, a neural network was trained to perform a classification task for stellar objects: stars, quasars and galaxies. The performance of the resulting network was tested on a set of 200 objects, where it obtained an accuracy rate of 97.65\%. These applications are just a small sample of the potential that Deep Learning can have in physics, in particular to cosmology.

Betelgeuse is an important variable star with many observations in the AAVSO database, but there is an annual gap of about four months where Betelgeuse is close to the sun and not observable at night. This gap could be filled with daylight observations. The star is bright enough to be imaged with small telescopes during the day, so photometry is possible when the sun is up. We present V band photometry of Alpha Ori taken with an amateur telescope equipped with an interline-transfer CCD camera and neutral density filter. These data compare favorably with contemporaneous nighttime photometry. The method used is a variation on ensemble photometry (using other bright daytime stars), and involves large stacks of very short exposures. The ensemble method provided V magnitudes of Betelgeuse with calculated errors of 0.020 +-0.008 mag from February to April 2021. From May to July, at the closest distances to the sun, the photometry of Betelgeuse could be continued with mean errors of 0.040 +-0.013 mag.

In the curvaton scenario, the curvature perturbation is generated after inflation at the curvaton decay, which may have a prominent non-Gaussian effect. For a model with a non-trivial kinetic term, an enhanced curvature perturbation on a small scale can be realized, which can lead to copious production of primordial black holes (PBHs) and induce secondary gravitational waves (GWs). We find that under the assumption that thus formed PBHs would not overclose the universe, the non-Gaussianity of the curvature perturbation can be well approximated by the local quadratic form. When the curvaton energy fraction is small at the moment of curvaton decay, both the PBH abundance and the induced GW spectrum depends only on the power spectrum of the curvature perturbation. If asteroid-mass PBHs are the cold dark matter of the universe, we have $\Omega_\text{IGW}\gtrsim10^{-11}$ at $\sim10^{-2}~\text{Hz}$, which is detectable in the planned space GW detectors such as LISA.

We report the discovery of a nearby large, diffuse galaxy that shows star formation, using Ultra Violet Imaging Telescope (UVIT) far-UV observations, archival optical data from Multi-Unit Spectroscopic Explorer (MUSE) and Dark Energy Camera Legacy Survey (DECaLS), and InfraRed Survey Facility (IRSF) near-infrared observations. The galaxy was not detected earlier due to its superposition with the background galaxy, NGC 6902A. They were together mistakenly classified as an interacting system. NGC 6902A is at a redshift of 0.05554, but MUSE observations indicate that the interacting tail is a separate star-forming, foreground galaxy at a redshift of 0.00980. We refer to the new galaxy as UVIT J202258.73-441623.8 (UVIT J2022). The near-infrared observations show that UVIT J2022 has a stellar mass of 8.7$\times$10$^{8}$M$_{\odot}$. Its inner disk (R$<$4 kpc) shows UV and H$\alpha$ emission from ongoing massive star formation. The rest of the disk is extremely low luminosity, has a low stellar surface density, and extends out to a radius of R$\sim$9 kpc. The velocity and metallicity distribution maps and the star formation history indicate that UVIT J2022 has undergone three bursts of star formation. The latest episode is ongoing, which is supported by the presence of widespread H$\alpha$ and UV emission in its inner disk. The galaxy also shows patchy spiral arms in far-UV, and there is a metallicity enhancement along a bar-like feature. UVIT J2022 is thus a unique example of triggered star formation in a diffuse galaxy, resulting in the growth of its inner stellar disk. Our study raises the intriguing possibility that (i) there could be similar diffuse galaxies that have been mistakenly interpreted as interacting galaxies due to their superposition, and (ii) UV or H$\alpha$ could be a way to detect such diffuse galaxies in our local universe.

Matter power spectra emulators, such as the Euclid Emulator and CosmicEmu, are trained on simulations to correct the non-linear part of the power spectrum. Map-based analyses retrieve additional non-Gaussian information from the density field, whether through human-designed statistics such as peak counts, or via machine learning methods such as convolutional neural networks. The simulations required for these methods are very resource-intensive, both in terms of computing time and storage. Map-level density field emulators, based on deep generative models, have recently been proposed to address these challenges. In this work, we present a novel mass map emulator of the KiDS-1000 survey footprint, which generates noise-free spherical maps in a fraction of a second. It takes a set of cosmological parameters $(\Omega_M, \sigma_8)$ as input and produces a consistent set of 5 maps, corresponding to the KiDS-1000 tomographic redshift bins. To construct the emulator, we use a conditional generative adversarial network architecture and the spherical convolutional neural network $\texttt{DeepSphere}$, and train it on N-body-simulated mass maps. We compare its performance using an array of quantitative comparison metrics: angular power spectra $C_\ell$, pixel/peaks distributions, $C_\ell$ correlation matrices, and Structural Similarity Index. Overall, the agreement on these summary statistics is $<10\%$ for the cosmologies at the centre of the simulation grid, and degrades slightly on grid edges. Finally, we perform a mock cosmological parameter estimation using the emulator and the original simulation set. We find a good agreement in these constraints, for both likelihood and likelihood-free approaches. The emulator is available at https://tfhub.dev/cosmo-group-ethz/models/kids-cgan/1.

The past decade has brought impressive advances in astrophysics of cosmic rays (CRs) and multi-wavelength astronomy -- thanks to the new instrumentation launched into space and built on the ground. Understanding the astrophysical backgrounds to better precision than the observed data is vital in moving to this new territory. The state-of-the-art CR propagation code called GALPROP is designed to address exactly this challenge. Having 25 years of development behind it, the GALPROP framework has become a de-facto standard in astrophysics of CRs, diffuse photon emissions (radio- to gamma-rays), and searches for new physics. GALPROP uses information from astronomy, particle physics, and nuclear physics to predict CRs and their associated emissions self-consistently, providing a unifying modelling framework. The range of its physical validity covers 18 orders of magnitude in energy, from sub-keV to PeV energies for particles and from micro-eV to PeV energies for photons. The framework and the datasets are public and are extensively used by many experimental collaborations, and by thousands of individual researchers worldwide for interpretation of their data and for making predictions. This paper details the latest release of the GALPROP framework, further developments of its initially auxiliary datasets that grew into independent studies of the Galactic structure -- distributions of gas, dust, radiation and magnetic fields -- as well as the extension of its modelling capabilities. Example applications included with the distribution illustrating usage of the new features are also described.

The 21~cm transition from neutral Hydrogen promises to be the best observational probe of the Epoch of Reionisation. The main difficulty in measuring the 21 cm signal is the presence of bright foregrounds that require very accurate interferometric calibration. Closure quantities may circumvent the calibration requirements but may be, however, affected by direction dependent effects, particularly antenna primary beam responses. This work investigates the impact of antenna primary beams affected by mutual coupling on the closure phase and its power spectrum. Our simulations show that primary beams affected by mutual coupling lead to a leakage of foreground power into the EoR window, which can be up to $\sim4$ orders magnitude higher than the case where no mutual coupling is considered. This leakage is, however, essentially confined at $k < 0.3$~$h$~Mpc$^{-1}$ for triads that include 29~m baselines. The leakage magnitude is more pronounced when bright foregrounds appear in the antenna sidelobes, as expected. Finally, we find that triads that include mutual coupling beams different from each other have power spectra similar to triads that include the same type of mutual coupling beam, indicating that beam-to-beam variation within triads (or visibility pairs) is not the major source of foreground leakage in the EoR window.

We present a novel framework for the equation of state of dense and hot Quantum Chromodynamics (QCD), which focuses on the region of the phase diagram relevant for neutron star mergers and core-collapse supernovae. The model combines predictions from the gauge/gravity duality with input from lattice field theory, QCD perturbation theory, chiral effective theory and statistical modeling. It is therefore, by construction, in good agreement with theoretical constraints both at low and high densities and temperatures. The main ingredients of our setup are the non-perturbative V-QCD model based on the gauge/gravity duality, a van der Waals model for nucleon liquid, and the DD2 version of the Hempel-Schaffner-Bielich statistical model of nuclear matter. By consistently combining these models, we also obtain a description for the nuclear to quark matter phase transition and its critical endpoint. The parameter dependence of the model is represented by three (soft, intermediate and stiff) variants of the equation of state, all of which agree with observational constraints from neutron stars and their mergers. We discuss resulting constraints for the equation of state, predictions for neutron stars and the location of the critical point.

We derive some of the central equations governing quantum fluctuations in gravitational waves, making use of general relativity as a sensible effective quantum theory at large distances. We begin with a review of classical gravitational waves in general relativity, including the energy in each mode. We then form the quantum ground state and coherent state, before then obtaining an explicit class of squeezed states. Since existing gravitational wave detections arise from merging black holes, and since the quantum nature of black holes remains puzzling, one can be open-minded to the possibility that the wave is in an interesting quantum mechanical state, such as a highly squeezed state. We compute the time and space two-point correlation functions for the quantized metric perturbations. We then constrain its amplitude with LIGO observations. Using existing LIGO data, we place a bound on the (exponential) squeezing parameter of the quantum gravitational wave state of $\zeta<41$.

The increasing observational pressure on the standard cosmological model motivates analyses going beyond the paradigm of the collision-less cold dark matter (DM). Since the only clear evidence for the existence of DM is based on gravitational interactions, it seems particularly fitting to study them in this sector where extensions to the standard model can be naturally introduced. A promising avenue can be obtained using modifications of the space-time metric coupled to DM and induced by the presence of a new ultra-light scalar field $\phi$. The $\phi$ field can contribute to the DM density and can couple all the matter species universally, including additional heavy DM particles. We present a simple two-component DM model employing derivative conformal interactions between the two DM species. This can simultaneously: 1) guarantee the necessary symmetries to stabilize the dark species, 2) predict a subdominant thermal relic density of the heavy DM component, and 3) alleviate small-scale structure tensions of the cold DM scenario due to the possible co-interactions in the dark sector. The scenario is highly predictive with future observational prospects ranging from the Large Hadron Collider (LHC) to gravitational-wave searches, and can be generalized to more rich and realistic dark sector models.

We produce gravitational waveforms for nonspinning compact binaries undergoing a quasicircular inspiral. Our approach is based on a two-timescale expansion of the Einstein equations in second-order self-force theory, which allows first-principles waveform production in milliseconds. Although the approach is designed for extreme mass ratios, our waveforms agree remarkably well with those from full numerical relativity, even for comparable-mass systems. Our results will be invaluable in accurately modelling extreme-mass-ratio inspirals for the LISA mission and intermediate-mass-ratio systems currently being observed by the LIGO-Virgo-KAGRA Collaboration.

Heat always flows from hotter to a colder temperature until thermal equilibrium be finally restored in agreement with the usual (zeroth, first and second) laws of thermodynamics. However, Tolman and Ehrenfest demonstrated that the relation between inertia and weight uniting all forms of energy in the framework of general relativity implies that the standard equilibrium condition is violated in order to maintain the validity of the first and second law of thermodynamics. Here we demonstrate that the thermal equilibrium condition for a static self-gravitating fluid, besides being violated, is also heavily dependent on the underlying spacetime geometry (whether Riemannian or non-Riemannian). As a particular example, a new equilibrium condition is deduced for a large class of Weyl and f(R) type gravity theories. Such results suggest that experiments based on the foundations of the heat theory (thermal sector) may also be used for confronting gravity theories and prospect the intrinsic geometric nature of the spacetime structure.

We study a warm inflationary model for different expansions assuming an anisotropic universe described by Bianchi I metric. The universe is filled with a scalar field or inflaton, radiation, and bulk viscous pressure. We carry out the inflationary analysis for different solutions of such universe in two different cases $\Gamma$ and $\xi$ as constant and variable parameters, respectively. We compare the obtained results with the recent observations, in order to find the observational constraints on the parameters space of the models. Moreover, we attempt to present a better judgment among the considered models by calculation of the non-linear parameter $f_{NL}$ describing the non-Gaussianity property of the models. Additionally, we investigate the warm inflationary models with viscous pressure from the Weak Gravity Conjecture approach, considering the swampland criteria.

In the presence of a Weyl scaling invariant cosmological action, black holes no longer have an event horizon and an apparent horizon. So they have no trapped surfaces, and may be "leaky", emitting a "black hole wind" which could lead to star formation in the neighborhood of the hole. In this paper we formulate and analyze a one-dimensional model for star formation resulting from a postulated outgoing black hole particle flux, incident on a distant spherical surface modeled as a set of planar disks surrounding the hole. Using the Toomre analysis of the Jeans instability of a disk, we compute conditions for a disk collapse instability and estimate the collapse time. We suggest a mechanism for giving the disk angular momentum around the central black hole. This gives a possible explanation for the puzzling observation of young stars forming in the vicinity of the black hole Sagittarius A* central to the Milky Way galaxy.

We numerically study the primordial black hole (PBH) formation by an isocurvature perturbation of a massless scalar field on super Hubble scales in the radiation-dominated universe. As a first step we perform simulations of spherically symmetric configurations. For the initial condition, we employ the spatial gradient expansion and provide the general form of the growing mode solutions valid up through the second order in this expansion. The initial scalar field profile is assumed to be Gaussian with a characteristic comoving wavenumber $k$; $\sim\exp(-k^2R^2)$, where $R$ is the radial coordinate. We find that a PBH is formed for a sufficiently large amplitude of the scalar field profile. Nevertheless, we find that the late time behavior of the gravitational collapse is dominated by the dynamics of the fluid but not by the scalar field, which is analogous to the PBH formation from an adiabatic perturbation in the radiation-dominated universe.

We present an analysis that reconciles the gravitational-wave signal GW190521 observed by the Advanced LIGO and Advanced Virgo detectors with the electromagnetic flare ZTF19abanrhr observed by the Zwicky Transient Facility. We analyze GW190521 under a mass-ratio prior uniform in $Q \in [1,4]$ and using the state-of-the-art waveform model for black-hole mergers \texttt{NRSur7dq4}. We find a $90\%$ credible region for the black-hole masses extending far outside what originally reported by \cite{GW190521D}, where our maximum likelihood masses reside. We find a $15\%$ probability that both black holes avoid the pair-instability supernova gap. We infer a three-dimensional sky-location highly consistent with ZTF19abanrhr, obtaining an odds-ratio ${\cal{O}}_{C/R}=72:1$ that strongly favors the hypothesis of a true coincidence over a random one. Combining this event with the neutron-star merger GW170817, we estimate a Hubble constant H$_0=72.1^{+10.6}_{-6.4}\mathrm{km\,s^{-1}\,Mpc^{-1}}$ at the $68\%$ credible level.

The Newtonian and general relativistic equations for radial motion in the field of a central massive object are identical, and have the property that particles infalling from rest at infinity, and black hole "wind" particles with relativistic velocity leaking out of a Schwarzschild-like black hole nominal horizon, both have the same magnitude of velocity at any radius from the hole. Hence when equally massive infalling and wind particles collide at any radius, they yield collision products with zero center of mass velocity, which can then nucleate star formation at the collision radius. We suggest that this gives a general mechanism by which a central black hole can catalyze galaxy formation.

The solutions of the Einstein-Maxwell system of equations, useful for the modeling of compact stars have a long and rich history. It took just a year since Einstein's general theory of relativity was published for Karl Schwarzschild to obtain the first exact solution of Einstein's field equations. The number of viable exact solutions has been growing since then. Different models have been constructed for a variety of applications. In this article, we have discussed different ways of generating a static spherically symmetric anisotropic fluid model. The purpose of the present article is to present a simple classification scheme for static and spherically symmetric anisotropic fluid solutions. The known solutions are reviewed and compartmentalized according to this scheme so that we can illustrate general ideas about these solutions without being exhaustive.

The extraordinary neutrino flux produced in extreme astrophysical environments like the early universe, core-collapse supernovae and neutron star mergers may produce coherent quantum neutrino oscillations on macroscopic length scales. The Hamiltonian describing this evolution can be mapped into quantum spin models with all-to-all couplings arising from neutrino-neutrino forward scattering. To date many studies of these oscillations have been performed in a mean-field limit where the neutrinos time evolve in a product state. In this paper we examine a simple two-beam model evolving from an initial product state and compare the mean-field and many-body evolution. The symmetries in this model allow us to solve the real-time evolution for the quantum many-body system for hundreds or thousands of spins, far beyond what would be possible in a more general case with an exponential number ($2^N$) of quantum states. We compare mean-field and many-body solutions for different initial product states and ratios of one- and two-body couplings, and find that in all cases in the limit of infinite spins the mean-field (product state) and many-body solutions coincide for simple observables. This agreement can be understood as a consequence of the fact that the typical initial condition represents a very local but dense distribution about a mean energy in the spectrum of the Hamiltonian. We explore quantum information measures like entanglement entropy and purity of the many-body solutions, finding intriguing relationships between the quantum information measures and the dynamical behavior of simple physical observables.

We investigate spectral properties of turbulence in the solar wind that is a weakly collisional astrophysical plasma, accessible to in-situ observations. Using the Helios search coil magnetometer measurements in the fast solar wind, in the inner heliosphere, we focus on properties of the turbulent magnetic fluctuations at scales smaller than the ion characteristic scales, the so-called kinetic plasma turbulence. At such small scales, we show that the magnetic power spectra between 0.3 and 0.9 AU from the Sun have a generic shape $\sim f^{-8/3}\exp{(-f/f_d)}$ where the dissipation frequency $f_d$ is correlated with the Doppler shifted frequency $f_{\rho e}$ of the electron Larmor radius. This behavior is statistically significant: all the observed kinetic spectra are well described by this model, with $f_d = f_{\rho e}/1.8$. Our results indicate that the electron gyroradius plays the role of the dissipation scale and marks the end of the electromagnetic cascade in the solar wind.

In unimodular gravity the stress-energy tensor of matter is not necessarily conserved, and so the theory offers a natural framework for interacting dark energy models, where dark energy has a constant equation of state $w=-1$. We derive the equations of motion for linear cosmological perturbations in interacting dark energy models of this class, focusing on the scalar sector. Then, we consider a specific model where the energy-momentum transfer potential is proportional to the energy density of cold dark matter; this transfer potential has the effect of inducing an effective equation of state $w_{\rm eff}\neq0$ for cold dark matter. We analyze in detail the evolution of perturbations during radiation domination on super-Hubble scales, finding that the well-known large-scale instability that affects a large class of interacting dark energy models is absent in this model. To avoid a gradient instability, energy must flow from dark matter to dark energy. Finally, we show that interacting dark energy models with $w=-1$ are equivalent to a class of generalized dark matter models.

Monitoring of high energy cosmic ray neutrons is of particular interest for cosmic ray water Cherenkov detectors as intense bundles of delayed neutrons have been found to arrive after the initial passage of a high energy air shower. In this paper we explore the possibility of building large-area high-energy neutron monitors using gadolinium-loaded Water Cherenkov Detectors (WCDs). GEANT4 simulations of photon production in WCDs are used to estimate the maximum detection efficiency for a hypothetical system. Requiring a series of neutron induced gamma ray flashes distributed over an extended period of time (up to 20{\mu}s) was shown to be an effective way to discriminate high energy neutron interactions from other backgrounds. Results suggest that neutron detection efficiencies of 4-15% may be possible using a gadolinium-loaded detection system above 200 MeV. The magnitude of gadolinium loading was also shown to significantly modify the timing response of the simulated detector.

Axions, hypothetical particles theorized to solve the strong CP-problem, are presently being considered as strong candidates as cold dark matter constituents. The signal power of resonant-based axion detectors, known as haloscopes, is directly proportional to their quality factor $Q$. In this paper, the impact of the use of superconductors in the performances of the haloscopes is studied by evaluating the obtainable $Q$. In particular, the surface resistance $R_s$ of NbTi, Nb$_3$Sn, YBa$_2$Cu$_3$O$_{7-\delta}$ and FeSe$_{0.5}$Te$_{0.5}$ is computed in the frequency, magnetic field and temperature ranges of interest, starting from the measured vortex motion complex resistivity and screening lengths of these materials. From $R_s$ the quality factor $Q$ of a cylindrical haloscope with copper conical bases and superconductive lateral wall, operating with the TM$_{010}$ mode, is evaluated and used to perform a comparison of the performances of the different materials. Both YBa$_2$Cu$_3$O$_{7-\delta}$ and FeSe$_{0.5}$Te$_{0.5}$ are shown to improve the measurement sensitivity by almost an order of magnitude with respect to a whole Cu cavity, while NbTi is shown to be suitable only at lower frequencies (<10 GHz). Nb$_3$Sn can give an intermediate improvement in the whole spectrum of interest.

Although capillary and permeability are the two most important physical properties controlling fluid distribution and flow in nature, the interconnectivity function between them was a pressing challenge. Because knowing permeability leads to determining capillary pressure. Geodynamics (e.g., subsurface water, CO2 sequestration) and organs (e.g., plants, blood vessels) depend on capillary pressure and permeability. The first determines how far the fluid can reach, while the second determines how fast the fluid can flow in porous media. They are also vital to designing synthetic materials and micro-objects like membranes and micro-robotics. Here, we reveal the capillary and permeability intertwined behavior function. And demonstrate the unique physical connectors: pore throat size and network, linking capillary pressure and permeability. Our discovery quantifies the inverse relationship between capillary pressure and permeability for the first time, which we analytically derived and experimentally proved.