Papers for Thursday, Oct 20 2022

We present a data-driven technique to analyze multi-frequency images from upcoming cosmological surveys mapping large sky area. Using full information from the data at the two-point level, our method can simultaneously constrain the large-scale structure (LSS), the spectra and redshift distribution of emitting sources, and the noise in observed data without any prior assumptions beyond the homogeneity and isotropy of cosmological perturbations. In particular, the method does not rely on source detection or photometric or spectroscopic redshift estimates. Here we present the formalism and demonstrate our technique with a mock observation from nine optical and near-infrared photometric bands. Our method can recover the input signal and noise without bias, and quantify the uncertainty on the constraints. Our technique provides a flexible framework to analyze the large-scale structure observation traced by different types of sources, which has potential for wide application to the current or future cosmological datasets such as SPHEREx, Rubin Observatory, Euclid, or the Nancy Grace Roman Space Telescope.

The angular momentum of gas feeding a black hole (BH) is typically misaligned with respect to the BH spin, resulting in a tilted accretion disk. Rotation of the BH drags the surrounding space-time, manifesting as Lense-Thirring torques that lead to disk precession and warping. We study these processes by simulating a thin ($H/r=0.02$), highly tilted ($\mathcal{T}=65^\circ$) accretion disk around a rapidly rotating ($a=0.9375$) BH at extremely high resolutions, which we performed using the general-relativistic magnetohydrodynamic (GRMHD) code H-AMR. The disk becomes significantly warped and continuously tears into two individually precessing sub-disks. We find that mass accretion rates far exceed the standard $\alpha$-viscosity expectations. We identify two novel dissipation mechanisms specific to warped disks that are the main drivers of accretion, distinct from the local turbulent stresses that are usually thought to drive accretion. In particular, we identify extreme scale height oscillations that occur twice an orbit throughout our disk. When the scale height compresses, `nozzle' shocks form, dissipating orbital energy and driving accretion. Separate from this phenomenon, there is also extreme dissipation at the location of the tear. This leads to the formation of low-angular momentum `streamers' that rain down onto the inner sub-disk, shocking it. The addition of low angular momentum gas to the inner sub-disk causes it to rapidly accrete, even when it is transiently aligned with the BH spin and thus unwarped. These mechanisms, if general, significantly modify the standard accretion paradigm. Additionally, they may drive structural changes on much shorter timescales than expected in $\alpha$-disks, potentially explaining some of the extreme variability observed in active galactic nuclei.

Gravitational-wave observations have just started probing the properties of black-hole binary merger populations. The observation of binaries with very massive black holes and small mass ratios motivates the study of dense star clusters as astrophysical environments which can produce such events dynamically. In this paper we present Rapster (for "Rapid cluster evolution"), a new code designed to rapidly model binary black hole population synthesis and the evolution of star clusters based on simple, yet realistic prescriptions. The code can be used to generate large populations of dynamically formed binary black holes. It makes use of existing packages such as SEVN to model the initial black hole mass spectrum, and PRECESSION to model the mass, spin and gravitational recoil of merger remnants. We demonstrate that the event rates and population properties predicted by Rapster are in good agreement with other state-of-the-art codes.

Cosmological inference with large galaxy surveys requires theoretical models that combine precise predictions for large-scale structure with robust and flexible galaxy formation modelling throughout a sufficiently large cosmic volume. Here, we introduce the MillenniumTNG (MTNG) project which combines the hydrodynamical galaxy formation model of IllustrisTNG with the large volume of the Millennium simulation. Our largest hydrodynamic simulation, covering (500 Mpc/h)^3 = (740 Mpc)^3, is complemented by a suite of dark-matter-only simulations with up to 4320^3 dark matter particles (a mass resolution of 1.32 x 10^8 Msun/h) using the fixed-and-paired technique to reduce large-scale cosmic variance. The hydro simulation adds 4320^3 gas cells, achieving a baryonic mass resolution of 2 x 10^7 Msun/h. High time-resolution merger trees and direct lightcone outputs facilitate the construction of a new generation of semi-analytic galaxy formation models that can be calibrated against both the hydro simulation and observation, and then applied to even larger volumes - MTNG includes a flagship simulation with 1.1 trillion dark matter particles and massive neutrinos in a volume of (3000 Mpc)^3. In this introductory analysis we carry out convergence tests on basic measures of non-linear clustering such as the matter power spectrum, the halo mass function and halo clustering, and we compare simulation predictions to those from current cosmological emulators. We also use our simulations to study matter and halo statistics, such as halo bias and clustering at the baryonic acoustic oscillation scale. Finally we measure the impact of baryonic physics on the matter and halo distributions.

Cosmological simulations are an important theoretical pillar for understanding nonlinear structure formation in our Universe and for relating it to observations on large scales. In several papers, we introduce our MillenniumTNG (MTNG) project that provides a comprehensive set of high-resolution, large volume simulations of cosmic structure formation aiming to better understand physical processes on large scales and to help interpreting upcoming large-scale galaxy surveys. We here focus on the full physics box MTNG740 that computes a volume of $(740\,\mathrm{Mpc})^3$ with a baryonic mass resolution of $3.1\times~10^7\,\mathrm{M_\odot}$ using \textsc{arepo} with $80.6$~billion cells and the IllustrisTNG galaxy formation model. We verify that the galaxy properties produced by MTNG740 are consistent with the TNG simulations, including more recent observations. We focus on galaxy clusters and analyse cluster scaling relations and radial profiles. We show that both are broadly consistent with various observational constraints. We demonstrate that the SZ-signal on a deep lightcone is consistent with Planck limits. Finally, we compare MTNG740 clusters with galaxy clusters found in Planck and the SDSS-8 RedMaPPer richness catalogue in observational space, finding very good agreement as well. However, {\it simultaneously} matching cluster masses, richness, and Compton-$y$ requires us to assume that the SZ mass estimates for Planck clusters are underestimated by $0.2$~dex on average. Thanks to its unprecedented volume for a high-resolution hydrodynamical calculation, the MTNG740 simulation offers rich possibilities to study baryons in galaxies, galaxy clusters, and in large scale structure, and in particular their impact on upcoming large cosmological surveys.

Modern redshift surveys are tasked with mapping out the galaxy distribution over enormous distance scales. Existing hydrodynamical simulations, however, do not reach the volumes needed to match upcoming surveys. We present results for the clustering of galaxies using a new, large volume hydrodynamical simulation as part of the MillenniumTNG (MTNG) project. With a computational volume that is $\approx15$ times larger than the next largest such simulation currently available, we show that MTNG is able to accurately reproduce the observed clustering of galaxies as a function of stellar mass. When separated by colour, there are some discrepancies with respect to the observed population, which can be attributed to the quenching of satellite galaxies in our model. We combine MTNG galaxies with those generated using a semi-analytic model to emulate the sample selection of luminous red galaxies (LRGs) and emission line galaxies (ELGs), and show that although the bias of these populations is approximately (but not exactly) constant on scales larger than $\approx10$ Mpc, there is significant scale-dependent bias on smaller scales. The amplitude of this effect varies between the two galaxy types, and also between the semi-analytic model and MTNG. We show that this is related to the distribution of haloes hosting LRGs and ELGs. Using mock SDSS-like catalogues generated on MTNG lightcones, we demonstrate the existence of prominent baryonic acoustic features in the large-scale galaxy clustering. We also demonstrate the presence of realistic redshift space distortions in our mocks, finding excellent agreement with the multipoles of the redshift-space clustering measured in SDSS data.

The early release science results from the $\textit{James Webb Space Telescope (JWST)}$ have yielded an unexpected abundance of high-redshift luminous galaxies that seems to be in tension with current theories of galaxy formation. However, it is currently difficult to draw definitive conclusions form these results as the sources have not yet been spectroscopically confirmed. It is in any case important to establish baseline predictions from current state-of-the-art galaxy formation models that can be compared and contrasted with these new measurements. In this work, we use the new large-volume ($L_\mathrm{box}\sim 740 \, \mathrm{cMpc}$) hydrodynamic simulation of the MilleniumTNG project to make predictions for the high-redshift ($z\gtrsim8$) galaxy population and compare them to recent $\textit{JWST}$ observations. We show that the simulated galaxy population is broadly consistent with observations until $z\sim10$. From $z\approx10-12$, the observations indicate a preference for a galaxy population that is largely dust-free, but is still consistent with the simulations. Beyond $z\gtrsim12$, however, our simulation results underpredict the abundance of luminous galaxies and their star-formation rates by almost an order of magnitude. This indicates either an incomplete understanding of the new $\textit{JWST}$ data or a need for more sophisticated galaxy formation models that account for additional physical processes such as Population~III stars, variable stellar initial mass functions, or even deviations from the standard $\Lambda$CDM model. We emphasise that any new process invoked to explain this tension should only significantly influence the galaxy population beyond $z\gtrsim10$, while leaving the successful galaxy formation predictions of the fiducial model intact below this redshift.

Luminous red galaxies (LRGs) and blue star-forming emission-line galaxies (ELGs) are key tracers of large-scale structure used by cosmological surveys. Theoretical predictions for such data are often done via simplistic models for the galaxy-halo connection. In this work, we use the large, high-fidelity hydrodynamical simulation of the MillenniumTNG project (MTNG) to inform a new phenomenological approach for obtaining an accurate and flexible galaxy-halo model on small scales. Our aim is to study LRGs and ELGs at two distinct epochs, $z = 1$ and $z = 0$, and recover their clustering down to very small scales, $r \sim 0.1 \ {\rm Mpc}/h$, i.e. the one-halo regime, while a companion paper extends this to a two-halo model for larger distances. The occupation statistics of ELGs in MTNG inform us that: (1) the satellite occupations exhibit a slightly super-Poisson distribution, contrary to commonly made assumptions, and (2) that haloes containing at least one ELG satellite are twice as likely to host a central ELG. We propose simple recipes for modeling these effects, each of which calls for the addition of a single free parameter to simpler halo occupation models. To construct a reliable satellite population model, we explore the LRG and ELG satellite radial and velocity distributions and compare them with those of subhalos and particles in the simulation. We find that ELGs are anisotropically distributed within halos, which together with our occupation results provides strong evidence for cooperative galaxy formation (manifesting itself as one-halo galaxy conformity); i.e.~galaxies with similar properties form in close proximity to each other. Our refined galaxy-halo model represents a useful improvement of commonly used analysis tools and thus can be of help to increase the constraining power of large-scale structure surveys.

We present a detailed analysis of the ionised gas distribution and kinematics in the inner ~ 200 pc of NGC 4546, host of a Low Luminosity Active Galactic Nucleus (LLAGN). Using GMOS-IFU observations, with a spectral coverage of 4736-6806 A and an angular resolution of 0.7 arcsec, we confirm that the nuclear emission is consistent with photoionisation by an AGN, while the gas in the circumnuclear region may be ionised by hot low-mass evolved stars. The gas kinematics in the central region of NGC 4546 presents three components: (i) a disc with major axis oriented along a position angle of 43deg +/- 3deg, counter rotating relative to the stellar disc; (ii) non-circular motions, evidenced by residual velocities of up to 60 km/s, likely associated to a previous capture of a dwarf satellite by NGC 4546; and (iii) nuclear outflows in ionised gas, identified as a broad component (sigma ~ 320 km/s) in the line profiles, with a mass outflow rate of dMout = 0.3 +/- 0.1 Msun/yr and a total mass of Mout = (9.2 +/- 0.8)E3 Msun in ionised gas, corresponding to less than 3 per cent of the total mass of ionised gas in the inner 200 pc of NGC 4546. The kinetic efficiency of the outflow is roughly 0.1 per cent, which is smaller than the outflow coupling efficiencies predicted by theoretical studies to AGN feedback become efficient in suppressing star formation in the host galaxy.

Approximate methods to populate dark matter halos with galaxies are of great utility to large galaxy surveys. However, the limitations of simple halo occupation models (HODs) preclude a full use of small-scale galaxy clustering data and call for more sophisticated models. We study two galaxy populations, luminous red galaxies (LRGs) and star-forming emission-line galaxies (ELGs), at two epochs, $z=1$ and $z=0$, in the large volume, high-resolution hydrodynamical simulation of the MillenniumTNG project. In a partner study we concentrated on the small-scale, one-halo regime down to $r\sim 0.1{\rm Mpc}/h$, while here we focus on modeling galaxy assembly bias in the two-halo regime, $r\gtrsim 1{\rm Mpc}/h$. Interestingly, the ELG signal exhibits scale dependence out to relatively large scales ($r\sim 20{\rm Mpc}/h$), implying that the linear bias approximation for this tracer is invalid on these scales, contrary to common assumptions. The 10-15\% discrepancy present in the standard halo model prescription is only reconciled when we augment our halo occupation model with a dependence on extrinsic halo properties ("shear" being the best-performing one) rather than intrinsic ones (e.g., concentration, peak mass). We argue that this fact constitutes evidence for two-halo galaxy conformity. Including tertiary assembly bias (i.e. a property beyond mass and "shear") is not an essential requirement for reconciling the galaxy assembly bias signal of LRGs, but the combination of external and internal properties is beneficial for recovering ELG the clustering. We find that centrals in low-mass haloes dominate the assembly bias signal of both populations. Finally, we explore the predictions of our model for higher-order statistics such as nearest-neighbor counts. The latter supplies additional information about galaxy assembly bias and can be used to break degeneracies between halo model parameters.

We introduce a novel technique for constraining cosmological parameters and galaxy assembly bias using non-linear redshift-space clustering of galaxies. We scale cosmological N-body simulations and insert galaxies with the SubHalo Abundance Matching extended (SHAMe) empirical model to generate over 175,000 clustering measurements spanning all relevant cosmological and SHAMe parameter values. We then build an emulator capable of reproducing the projected galaxy correlation function at the monopole, quadrupole and hexadecapole level for separations between $0.1\,h^{-1}{\rm Mpc}$ and $25\,h^{-1}{\rm Mpc}$. We test this approach by using the emulator and Monte Carlo Markov Chain (MCMC) inference to jointly estimate cosmology and assembly bias parameters both for the MTNG740 hydrodynamic simulation and for a semi-analytical galaxy formation model (SAM) built on the MTNG740-DM dark matter-only simulation, obtaining unbiased results for all cosmological parameters. For instance, for MTNG740 and a galaxy number density of $n\sim 0.01 h^{3}{\rm Mpc}^{-3}$, we obtain $\sigma_{8}=0.799^{+0.039}_{-0.044}$ ($\sigma_{8,{\rm MTNG}} =$ 0.8159), and $\Omega_\mathrm{M}h^2= 0.138^{+ 0.025}_{- 0.018}$ ($\Omega_{\mathrm{M}} h^2_{\rm MTNG} =$ 0.142). For fixed Hubble parameter ($h$), the constraint becomes $\Omega_\mathrm{M}h^2= 0.137^{+ 0.011}_{- 0.012}$. Our method performs similarly well for the SAM and for other tested sample densities. We almost always recover the true amount of galaxy assembly bias within one sigma. The best constraints are obtained when scales smaller than $2\,h^{-1}{\rm Mpc}$ are included, as well as when at least the projected correlation function and the monopole are incorporated. These methods offer a powerful way to constrain cosmological parameters using galaxy surveys.

Using data from the Galactic Arecibo L-band Feed Array Continuum Transit Survey (GALFACTS), we report the discovery of two previously unidentified, very compressed, thin, and straight polarized filaments approximately centred at Galactic coordinates, $(l,b)=(182.5^\circ,-4.0^\circ)$, which we call G182.5--4.0. Using data from the Isaac Newton Telescope Galactic Plane Survey (IGAPS), we also find straight, long, and extremely thin H$\alpha$ filaments coincident with the radio emission. These filaments are positioned in projection at the edge of the Orion-Eridanus superbubble and we find evidence indicating that the filaments align with the coherent magnetic field of the outer Galaxy. We find a lower limit on the total radio flux at 1.4~GHz to be $0.7\pm0.3$~Jy with an average linearly polarized fraction of $40\substack{+30 \\ -20}\%$. We consider various scenarios that could explain the origin of these filaments, including a shell-type supernova remnant (SNR), a bow shock nebula associated with a pulsar, or relic fragments from one or more supernova explosions in the adjacent superbubble, with a hybrid scenario being most likely. This may represent an example of a new class of objects that is neither an SNR nor a bow shock. The highly compressed nature of these filaments and their alignment with Galactic plane suggests a scenario where this object formed in a magnetic field that was compressed by the expanding Orion-Eridanus superbubble, suggesting that the object is related to this superbubble and implying a distance of $\sim$400~pc.

Evidence for almost spatial flatness of the Universe has been provided from several observational probes, including the Cosmic Microwave Background (CMB) and Baryon Acoustic Oscillations (BAO) from galaxy clustering data. However, other than inflation, and in this case only in the limit of infinite time, there is no strong a priori motivation for a spatially flat Universe. Using the renormalization group (RG) technique in curved spacetime, we present in this work a theoretical motivation for spatial flatness. Starting from a general spacetime, the first step of the RG, coarse-graining, gives a Friedmann-Lema\^itre-Robertson-Walker (FLRW) metric with a set of parameters. Then, we study the rescaling properties of the curvature parameter, and find that zero spatial curvature of the FLRW metric is singled out as the unique scale-free, non-singular background for cosmological perturbations.

The NISP (Near Infrared Spectrometer and Photometer) is one of the two Euclid instruments. It operates in the near-IR spectral region (950-2020nm) as a photometer and spectrometer. The instrument is composed of: a cold (135 K) optomechanical subsystem consisting of a Silicon carbide structure, an optical assembly, a filter wheel mechanism, a grism wheel mechanism, a calibration unit, and a thermal control system, a detection system based on a mosaic of 16 H2RG with their front-end readout electronic, and a warm electronic system (290 K) composed of a data processing / detector control unit and of an instrument control unit that interfaces with the spacecraft via a 1553 bus for command and control and via Spacewire links for science data. This paper presents: the final architecture of the flight model instrument and subsystems, and the performance and the ground calibration measurement done at NISP level and at Euclid Payload Module level at operational cold temperature.

With the most trans-iron elements detected of any star outside the Solar System, HD 222925 represents the most complete chemical inventory among r-process-enhanced, metal-poor stars. While the abundance pattern of the heaviest elements identified in HD 222925 agrees with the scaled Solar r-process residuals, as is characteristic of its r-process-enhanced classification, the newly measured lighter r-process elements display marked differences from their Solar counterparts. In this work, we explore which single astrophysical site (if any) produced the entire range of elements ($34\leq Z\leq 90$) present in HD 222925. We find that the abundance pattern of lighter r-process elements newly identified in HD 222925 presents a challenge for our existing nucleosynthesis models to reproduce. The most likely astrophysical explanation for the elemental pattern of HD 222925 is that its light r-process elements were created in rapidly expanding ejecta (e.g., from shocked, dynamical ejecta of compact object merger binaries). However, we find that the light r-process-element pattern can also be successfully reproduced by employing different nuclear mass models, indicating a need for a fresh investigation of nuclear input data for elements with $46\lesssim Z\lesssim 52$ by experimental methods. Either way, the new elemental abundance pattern of HD 222925 -- particularly the abundances obtained from space-based, ultraviolet (UV) data -- call for a deeper understanding of both astrophysical r-process sites and nuclear data. Similar UV observations of additional r-process-enhanced stars will be required to determine whether the elemental abundance pattern of HD 222925 is indeed a canonical template for the r-process at low metallicity.

The warm-hot intergalactic medium (warm-hot IGM, or WHIM) pervades the filaments of the Cosmic Web and harbours half of the Universe's baryons. The WHIM's thermodynamic properties are notoriously hard to measure. Here we estimate a galaxy group - WHIM boundary temperature using a new method. In particular, we use a radio image of the giant radio galaxy (giant RG, or GRG) created by NGC 6185, a massive nearby spiral. We analyse this extraordinary object with a Bayesian 3D lobe model and deduce an equipartition pressure $P_\mathrm{eq} = 6 \cdot 10^{-16}\ \mathrm{Pa}$ -- among the lowest found in RGs yet. Using an X-ray-based statistical conversion for Fanaroff-Riley II RGs, we find a true lobe pressure $P = 1.5\substack{+1.7\\-0.4}\cdot 10^{-15}\ \mathrm{Pa}$. Cosmic Web reconstructions, group catalogues, and MHD simulations furthermore imply an $\mathrm{Mpc}$-scale IGM density $1 + \delta_\mathrm{IGM} = 40\substack{+30\\-10}$. The buoyantly rising lobes are crushed by the IGM at their inner side, where an approximate balance between IGM and lobe pressure occurs: $P_\mathrm{IGM} \approx P$. The ideal gas law then suggests an IGM temperature $T_\mathrm{IGM} = 11\substack{+12\\-5} \cdot 10^6\ \mathrm{K}$, or $k_\mathrm{B}T_\mathrm{IGM} = 0.9\substack{+1.0\\-0.4}\ \mathrm{keV}$, at the virial radius -- consistent with X-ray-derived temperatures of similarly massive groups. Interestingly, the method is not performing at its limit: in principle, estimates $T_\mathrm{IGM} \sim 4 \cdot 10^6\ \mathrm{K}$ are already possible -- rivalling the lowest X-ray measurements available. The technique's future scope extends from galaxy group outskirts to the WHIM. In conclusion, we demonstrate that observations of GRGs in Cosmic Web filaments are finally sensitive enough to probe the thermodynamics of galaxy groups and beyond.

3D simulations of high mass young stellar object (HMYSO) growth show that their circumstellar discs fragment onto multiple self-gravitating objects. Accretion of these by HMYSO may explain episodic accretion bursts discovered recently. We post-process results of a previous 3D simulation of a HMYSO disc with a 1D code that resolves the disc and object dynamics down to the stellar surface. We find that burst-like deposition of material into the inner disc seen in 3D simulations by itself does not always signify powerful accretion bursts. Only high density post-collapse clumps crossing the inner computational boundary may result in observable bursts. The rich physics of the inner disc has a significant impact on the expected accretion bursts: (1) In the standard turbulent viscosity discs, migrating objects can stall at a migration trap at the distance of a few au from the star. However, in discs powered by magnetised winds, the objects are able to cross the trap and produce bursts akin to those observed so far. (2) Migrating objects may interact with and modify the thermal (hydrogen ionisation) instability of the inner disc, which can be responsible for longer duration and lower luminosity bursts in HMYSOs. (3) If the central star is bloated to a fraction of an au by a previous episode of high accretion rate, or if the migrating object is particularly dense, a merger rather than a disc-mediated accretion burst results; (4) Object disruption bursts may be super-Eddington, leading to episodic feedback on HMYSO surroundings via powerful outflows.

Primordial magnetic fields (PMFs) are possible candidates for explaining the observed magnetic fields in galaxy clusters. Two competing scenarios of primordial magnetogenesis have been discussed in the literature: inflationary and phase-transitional. We study the amplification of both large- and small-scale correlated magnetic fields corresponding to inflation- and phase-transition-generated PMFs in a massive galaxy cluster. We employ high-resolution magnetohydrodynamic cosmological zoom-in simulations to resolve turbulent motions in the intracluster medium. We find that the turbulent amplification is more efficient for the large-scale, inflationary models while the phase-transition-generated seed fields show moderate growth. The differences between the models are imprinted on the spectral characteristics of the field (such as amplitude and shape of the magnetic power spectrum) and therewith on the final correlation length. We find one order of magnitude difference in the final strengths between the inflation- and phase-transition-generated magnetic field and a factor of $1.5$ difference in their final coherence scales. Thus, the final configuration of the magnetic field retains information on the PMF generation scenarios. Our findings have implications for future extragalactic Faraday rotation surveys with the possibility of distinguishing between different magnetogenesis scenarios.

Many core collapse supernovae (SNe) with hydrogen-poor and low-mass ejecta, such as ultra-stripped SNe and type Ibn SNe, are observed to interact with dense circumstellar material (CSM). These events likely arise from the core-collapse of helium stars which have been heavily stripped by a binary companion and ejected significant mass during the last weeks to years of their lives. In helium star models run to days before core-collapse, we identify a range of helium core masses $\approx 2.5 -3 M_{\odot}$ whose envelopes expand substantially due to helium shell burning while the core undergoes neon and oxygen burning. When modeled in binary systems, the rapid expansion of these helium stars induces extremely high rates of late-stage mass transfer ($\dot{M} \gtrsim 10^{-2} M_\odot/{\rm yr}$) beginning weeks to decades before core-collapse. We consider two scenarios for producing CSM in these systems: either mass transfer remains stable and mass loss is driven from the system in the vicinity of the accreting companion, or mass transfer becomes unstable and causes a common envelope event (CEE) through which the helium envelope is unbound. The ensuing CSM properties are consistent with the CSM masses ($\sim 10^{-2}-1 M_\odot$) and radii ($\sim 10^{13}-10^{16} {\rm cm}$) inferred for ultra-stripped SNe and several type Ibn SNe. Furthermore, systems that undergo a CEE could produce short-period NS binaries that merge in less than 100 Myr.

The first Pop III stars formed out of primordial, metal free gas, in minihalos at z>20, and kickstarted the cosmic processes of reionizaton and enrichment. While these stars are likely more massive than their enriched counterparts, the current unknowns of their astrophysics include; when the first Pop III stars ignited, how massive they were, and when and how the era of the first stars ended. Investigating these questions requires an exploration of a multi-dimensional parameter space, including the slope of the Pop III stellar initial mass function (IMF) and the strength of the non-ionizing UV background. In this work, we present a novel model which treats both the slope and maximum mass of Pop III stars as truly free parameters while including the physics of the fragmentation of primordial gas. Our results also hint at a non-universal Pop III IMF which is dependent on the efficiency of primordial gas fragmentation. Our relatively simple model reproduces the results from hydrodynamic simulations, but with a computational efficiency which allows us to investigate the observable differences between a wide range of potential Pop III IMFs. In addition, the evolution of the number density of Pop III stars may provide insight into the evolution of the H2 dissociating background. While the slope of the Pop III IMF does not significantly affect the predicted number density of the first stars, more top heavy IMFs produce Pop III star clusters which are 2-3 magnitudes brighter than their more bottom heavy counterparts. While the Pop III star clusters are too dim for direct detection by JWST, we find they are within the reach of gravitational lensing.

In many black hole systems, the accretion disk is expected to be misaligned with respect to the black hole spin axis. If the scale height of the disk is much smaller than the misalignment angle, the spin of the black hole can tear the disk into multiple, independently precessing `sub-disks'. This is most likely to happen during outbursts in black hole X-Ray binaries (BHXRBs) and in active galactic nuclei (AGN) accreting above a few percent of the Eddington limit, because the disk becomes razor-thin. Disk tearing has the potential to explain variability phenomena including quasi-periodic oscillations (QPOs) in BHXRBs and changing-look phenomena in AGN. Here, we present the first radiative two-temperature GRMHD simulation of a strongly tilted ($65^{\circ}$) accretion disk around a $M_{BH}=10M_{\odot}$ black hole, which tears and precesses. This leads to luminosity swings between a few percent and $50 \%$ of the Eddington limit on sub-viscous timescales. Surprisingly, even where the disk is radiation pressure dominated, the accretion disk is thermally stable over $t \gtrsim 14,000 r_g/c$. This suggests warps play an important role in stabilizing the disk against thermal collapse. The disk forms two nozzle shocks perpendicular to the line of nodes where the scale height of the disk decreases $10$-fold and the electron temperature reaches $T_e \sim 10^8-10^9 K$. In addition, optically thin gas crossing the tear between the inner and outer disk gets heated to $T_e \sim 10^8 K$. This suggests that warped disks may emit a Comptonized spectrum that deviates substantially from idealized models.

We present a new $\texttt{python}$ package SARABANDE for measuring 3 & 4 Point Correlation Functions (3/4 PCFs) in $\mathcal{O}(N_{\rm g} \log N_{\rm g})$ time using Fast Fourier Transforms (FFTs), with $N_{\rm g}$ the number of grid points used for the FFT. SARABANDE can measure both projected and full 3 and 4 PCFs on gridded 2D and 3D datasets. The general technique is to generate suitable angular basis functions on an underlying grid, radially bin these to create kernels, and convolve these kernels with the original gridded data to obtain expansion coefficients about every point simultaneously. These coefficients are then combined to give us the 3/4 PCF as expanded in our basis. We apply SARABANDE to simulations of the Interstellar Medium (ISM) to show the results and scaling of calculating both the full and projected 3/4 PCFs.

The astorb database at Lowell Observatory is an actively curated catalog of all known asteroids in the Solar System. astorb has heritage dating back to the 1970's and has been publicly accessible since the 1990's. Beginning in 2015 work began to modernize the underlying database infrastructure, operational software, and associated web applications. That effort has involved the expansion of astorb to incorporate new data such as physical properties (e.g. albedo, colors, spectral types) from a variety of sources. The data in astorb are used to support a number of research tools hosted at https://asteroid.lowell.edu. Here we present a full description of the software tools, computational foundation, and data products upon which the astorb ecosystem has been built.

Radio galaxies are luminous structures created by the jets of supermassive black holes, and consist of atomic nuclei, relativistic electrons, and magnetic fields. In exceptional cases, radio galaxies attain cosmological, megaparsec extents - and thus turn into giants. Giants embody the most extreme known mechanism through which galaxies can impact the Cosmic Web around them. The triggers of giant growth remain a mystery. Excitingly, new sensitive low-frequency sky surveys hold promise to change this situation. In this work, we perform a precision measurement of the distribution of giant growth's central dynamical quantity: total length. We first construct a statistical geometric framework for radio galaxies that is both rigorous and practical. We then search the LOFAR Two-metre Sky Survey DR2 for giants, discovering 2050 previously unknown specimina: more than have been found in all preceding literature combined. Spectacular discoveries include the longest giant hosted by an elliptical galaxy, the longest giant hosted by a spiral galaxy, and 13 giants with an angular length larger than that of the full Moon. By combining theory and observations - and carefully forward modelling selection effects - we infer that giant radio galaxy lengths are well described by a Pareto distribution with tail index $-3.5 \pm 0.5$. This finding is a new observational constraint for models and simulations of radio galaxy growth. In addition, for the first time, we determine the comoving number density of giants, $5 \pm 2\ (100\ \mathrm{Mpc})^{-3}$, and the volume-filling fraction of giant radio galaxy lobes in clusters and filaments, $5\substack{+8\\-2}\cdot 10^{-6}$. We conclude that giants are truly rare - not only from an observational perspective, but also from a cosmological one. At any moment in time, most clusters and filaments - the building blocks of the modern Cosmic Web - do not harbour giants.

We present ALMA band 6 images of the 12CO, 13CO and C18O J=2-1 line emissions for the circumstellar disc around HD 169142, at ~8 au spatial resolution. We resolve a central gas depleted cavity, along with two independent near-symmetric ring-like structures in line emission: a well-defined inner gas ring [~25 au] and a second relatively fainter and diffuse outer gas ring [~65 au]. We identify a localised super-Keplerian feature or vertical flow with a magnitude of ~75m/s in the 12CO map. This feature has the shape of an arc that spans azimuthally across a PA range of ~60 to 45 degrees and radially in between the B1[26au] and B2[59au] dust rings. Through reconstruction of the gas surface density profile, we find that the magnitude of the background perturbations by the pressure support and self-gravity terms are not significant enough to account for the kinematic excess. If of planetary origin, the relative depletion in the gas-density profile would suggest a 1 Mj planet. In contrast, the central cavity displays relatively smooth kinematics, suggesting either a low mass companion and/or a binary orbit with a minimal vertical velocity component.

Modeling the 21cm global signal from the Cosmic Dawn is challenging due to the many poorly constrained physical processes that come into play. We address this problem using the semi-analytical code "Cosmic Archaeology Tool" (CAT). CAT follows the evolution of dark matter halos tracking their merger history and provides an ab initio description of their baryonic evolution, starting from the formation of the first (Pop III) stars and black holes (BHs) in mini-halos at z > 20. The model is anchored to observations of galaxies and AGN at z < 6 and predicts a reionization history consistent with constraints. In this work we compute the evolution of the mean global 21cm signal between $4\leq z \leq 40$ based on the rate of formation and emission properties of stars and accreting black holes. We obtain an absorption profile with a maximum depth $\delta {\rm T_b} = -95$ mK at $z \sim 26.5$ (54 MHz). This feature is quickly suppressed turning into an emission signal at $z = 20$ due to the contribution of accreting BHs that efficiently heat the IGM at $z < 27$. The high-$z$ absorption feature is caused by the early coupling between the spin and kinetic temperature of the IGM induced by Pop III star formation episodes in mini-halos. Once we account for an additional radio background from early BHs, we are able to reproduce the timing and the depth of the EDGES signal only if we consider a smaller X-ray background from accreting BHs, but not the shape.

We investigate the physical mechanism that decides the saturation level of turbulence in collapsing gas clouds. We perform a suite of high-resolution numerical simulations following the collapse of turbulent gas clouds with various effective polytropic exponents $\gamma_{\rm eff}$, initial Mach numbers $\mathcal{M}_0$, and initial turbulent seeds. Equating the energy injection rate by gravitational contraction and the dissipation rate of turbulence, we obtain an analytic expression of the saturation level of turbulence, and compare it with the numerical results. Consequently, the numerical results are well described by the analytic model, given that the turbulent driving scale in collapsing gas clouds is one-third of Jeans length of collapsing core. These results indicate that the strength of turbulence at the first core formation in the early universe/present-day star-formation process can be estimated solely by $\gamma_{\rm eff}$.

We present Gemini-NIFS, VLT-SINFONI and Keck-OSIRIS observations of near-infrared [Fe II] emission associated with the well-studied jets from three active T Tauri stars; RW Aur A, RY Tau and DG Tau taken from 2012-2021. We primarily covered the redshifted jet from RW Aur A, and the blueshifted jets from RY Tau and DG Tau, to investigate long-term time variabilities potentially related to the activities of mass accretion and/or the stellar magnetic fields. All of these jets consist of several moving knots with tangential velocities of 70-240 km s-1, ejected from the star with different velocities and at irregular time intervals. Via comparison with literature, we identify significant differences in tangential velocities for the DG Tau jet between 1985-2008 and 2008-2021. The sizes of the individual knots appear to increase with time, and in turn, their peak brightnesses in the 1.644-micron emission decreased up to a factor of ~30 during the epochs of our observations. A variety of the decay timescales measured in the [Fe II] 1.644 micron emission can be attributed to different pre-shock conditions if the moving knots are unresolved shocks. However, our data do not exclude the possibility that these knots are due to non-uniform density/temperature distributions with another heating mechanism, or in some cases due to stationary shocks without proper motions. Spatially resolved observations of these knots with significantly higher angular resolutions are necessary to better understand their physical nature.

Probing the solar corona is crucial to study the coronal heating and solar wind acceleration. However, the transient and inhomogeneous solar wind flows carry large-amplitude inherent Alfven waves and turbulence, which make detection more difficult. We report the oscillation and propagation of the solar wind at 2.6 solar radii (Rs) by observation of China Tianwen and ESA Mars Express with radio telescopes. The observations were carried out on Oct.9 2021, when one coronal mass ejection (CME) passed across the ray paths of the telescope beams. We obtain the frequency fluctuations (FF) of the spacecraft signals from each individual telescope. Firstly, we visually identify the drift of the frequency spikes at a high spatial resolution of thousands of kilometers along the projected baselines. They are used as traces to estimate the solar wind velocity. Then we perform the cross-correlation analysis on the time series of FF from different telescopes. The velocity variations of solar wind structure along radial and tangential directions during the CME passage are obtained. The oscillation of tangential velocity confirms the detection of streamer wave. Moreover, at the tail of the CME, we detect the propagation of an accelerating fast field-aligned density structure indicating the presence of magnetohydrodynamic waves. This study confirm that the ground station-pairs are able to form particular spatial projection baselines with high resolution and sensitivity to study the detailed propagation of the nascent dynamic solar wind structure.

We present MMT spectroscopic observations of two massive galaxy cluster candidates at redshift $z\sim0.07$ that show extended and diffuse X-ray emission in the ROSAT All Sky Survey (RASS) images. The targets were selected from a previous catalog of 303 newly-identified cluster candidates with the similar properties using the intra-cluster medium emission. Using the new MMT Hectospec data and SDSS archival spectra, we identify a number of member galaxies for the two targets and confirm that they are galaxy clusters at $z=0.079$ and 0.067, respectively. The size of the two clusters, calculated from the distribution of the member galaxies, is roughly 2 Mpc in radius. We estimate cluster masses using three methods based on their galaxy number overdensities, galaxy velocity dispersions, and X-ray emission. The overdensity-based masses are $(6\sim8) \rm \times10^{14}\ M_\odot$, comparable to the masses of large clusters at low redshift. The masses derived from velocity dispersions are significantly lower, likely due to their diffuse and low concentration features. Our result suggests the existence of a population of large clusters with very diffuse X-ray emission that have been missed by most previous searches using the RASS images. If most of the 303 candidates in the previous catalog are confirmed to be real clusters, this may help to reduce the discrepancy of cosmological results between the CMB and galaxy cluster measurements.

The spin-down limit of the continuous gravitational wave strain (assuming pulsars as triaxial rotators) depends upon the value of the intrinsic spin frequency derivative of the pulsar, among other parameters. In order to get more accurate intrinsic spin frequency derivative values, dynamical effects contributing to the measured spin frequency derivative values must be estimated via more realistic approaches. In this work, we calculate improved values of the spin-down limit of the gravitational wave strain (assuming pulsars as triaxial rotators) for a set of 237 pulsars for which a targeted search for continuous gravitational waves was carried out by the LVK Collaboration, recently. We use `GalDynPsr', a python-based public package, to calculate more realistic values of intrinsic spin frequency derivatives, and consequently, we get, more realistic values of the spin-down limit. For 136 pulsars, we obtain a higher value of the spin-down limit as compared to recently estimated values by the LVK Collaboration.

Magnetic reconnection is one of the major particle acceleration processes in space and astrophysical plasmas. Low-energy supra-thermal particles emitted by magnetic reconnection are a source of ionisation for circumstellar discs, influencing their chemical, thermal and dynamical evolution. The aim of this work is to propose a first investigation to evaluate how energetic particles can propagate in the circumstellar disc of a T Tauri star and how they affect the ionisation rate of the disc plasma. To that end, we have collected experimental and theoretical cross sections for the production of H$^+$, H$_2^+$ and He$^+$ by electrons and protons. Starting from theoretical injection spectra of protons and electrons emitted during magnetic reconnection events, we have calculated the propagated spectra in the circumstellar disc considering the relevant energy loss processes. We have considered fluxes of energetic particles with different spectral indices and different disc magnetic configurations, generated at different positions from the star considering the physical properties of the flares as deduced from the observations obtained by the Chandra Orion Ultra Deep point source catalogue. We have then computed the ionisation rates for a disc whose structure has been calculated with the radiation thermo-chemical code {\tt ProDiMo}. We find that energetic particles are potentially a very strong source of local ionisation with ionisation rates exceeding by several orders of magnitude the contribution due to X-rays, stellar energetic particles and radioactivity in the inner disc.

Since its launch, the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) satellite has discovered hundreds of X-ray sources, many of which lack proper classification. This mission also led to the discovery of new categories of high mass X-ray binaries (HMXB). We use the spectra of the X-Shooter instrument at the Very Large Telescope (VLT) of the European Southern Observatory (ESO) to better understand the nature of 3 accreting binaries (IGR J10101-5654, IGR J11435-6109 and IGR J12489-6243) discovered by INTEGRAL. We mainly focused on the lines and continuum from the X-Shooter spectra. We used atlases to constrain the nature of the sources and also complemented the spectra with measurements taken by Spitzer and the Wide-field Infrared Survey Explorer (WISE) in infrared, and parallaxes from Gaia for the distances. We determined the nature of each binary system: a BeHMXB system with a companion star of spectral type B0.5Ve with peculiar carbon emission for IGR J10101-5654 and IGR J11435-6109, and a CV system with an evolved K star (K0IV-K2IV) for IGR J12489-6243. We also estimated some geometrical parameters of the decretion disk and neutron star's orbit in the case of IGR J11435-6109.

In this paper, we extend our previous study \cite{DelPopolo2021} on the Lemaitre-Tolman (LT) model showing how the prediction of the model changes when the equation of state parameter ($w$) of dark energy is modified. In the previous study, it was considered that dark energy was merely constituted by the cosmological constant. In this paper, as in the previous study, we also took into account the effect of angular momentum and dynamical friction ($J\eta$ LT model) that modifies the evolution of a perturbation, initially moving with the Hubble flow. As a first step, solving the equation of motion, we calculated the relationship between mass, $M$, and the turn-around radius, $R_0$. If one knows the value of the turn-around radius $R_0$, it is possible to obtain the mass of the studied objects. As a second step, we build up, as in the previous paper, a relationship between the velocity, $v$, and radius, $R$. The relation was fitted to data of groups and clusters. Since the relationship $v-R$ depends on the Hubble constant and the mass of the object, we obtained optimized values of the two parameters of the objects studied. The mass decreases of a factor of maximum 25\% comparing the $J\eta$ LT results (for which $w=-1$) and the case $w=-1/3$, while the Hubble constant increases going from $w=-1$ to $w=-1/3$. Finally, the obtained values of the mass, $M$, and $R_0$ of the studied objects can put constraints to the dark energy equation of state parameter, $w$.

Upcoming large galaxy surveys will subject the standard cosmological model, $\Lambda$CDM, to new precision tests. These can be tightened considerably if theoretical models of galaxy formation are available that can predict galaxy clustering and galaxy-galaxy lensing on the full range of measurable scales throughout volumes as large as those of the surveys and with sufficient flexibility that uncertain aspects of the underlying astrophysics can be marginalised over. This, in particular, requires mock galaxy catalogues in large cosmological volumes that can be directly compared to observation, and can be optimised empirically by Monte Carlo Markov Chains or other similar schemes to eliminate or estimate astrophysical parameters related to galaxy formation when constraining cosmology. Semi-analytic galaxy formation methods implemented on top of cosmological dark matter simulations offer a computationally efficient approach to construct physically based and flexibly parametrised galaxy formation models, and as such they are more potent than still faster, but purely empirical models. Here we introduce an updated methodology for the semi-analytic L-GALAXIES code, allowing it to be applied to simulations of the new MillenniumTNG project, producing galaxies directly on fully continuous past lightcones, potentially over the full sky, out to high redshift, and for all galaxies more massive than $\sim 10^8\,{\rm M}_\odot$. We investigate the numerical convergence of the resulting predictions, and study the projected galaxy clustering signals of different samples. The new methodology can be viewed as an important step towards more faithful forward-modelling of observational data, helping to reduce systematic distortions in the comparison of theory to observations.

We modelled dust grain-size distributions for carbonaceous and silicates dust, as well as for free-flying iron nanoparticles in the environment of a $\gamma$-ray burst (GRB) afterglow, GRB 180325A. This GRB, at $z=2.2486$, has an unambiguous detection of the 2175 \r{A} extinction feature with $R_V=4.58$ and $A_V=1.58$. In addition to silicates, polycyclic aromatic hydrocarbons (PAH), and graphite, we used iron nanoparticles grain-size distributions for the first time to model the observed extinction curve of GRB 180325A. We fit the observed extinction for four model permutations, using 232 sets of silicates, graphite, carbon abundance in hydrocarbon molecules ($b_C$), and fraction of iron abundance in free-flying nanoparticles ($b_{\text{Fe}}$). These four different permutations were chosen to test iron nanoparticles significance and carbon abundance in hydrocarbons. Our results indicate that iron nanoparticles contribution is insignificant and there is a degeneracy of carbon abundances, with the range $(0.0 \leq b_C \leq 0.7)\times10^{-5}$ providing the best-fit to the observed extinction curve of GRB 180325A. We therefore favour the simplest model of silicates and polycyclic aromatic hydrocarbons. The silicates are dominant and contribute to the entire wavelength range of the GRB extinction curve while graphite contributes towards both the 2175 \r{A} bump and the UV extinction. The afterglow peak luminosity ($1.5\times10^{51}$ ergs/s) indicates dust destruction may have taken place. We conclude that further investigations into other potential contributors of extinction are warranted, particularly for steep UV extinction.

In order to extract maximal information about cosmology from the large-scale structure of the Universe, one needs to use every bit of signals that can be observed. Beyond the spatial distributions of astronomical objects, the spatial correlations of tensor fields, such as galaxy spins and shapes, are ones of promising sources that we can access in the era of large surveys in near future. The perturbation theory is a powerful tool to analytically describe the behaviors and evolutions of correlation statistics for a given cosmology. In this paper, we formulate a nonlinear perturbation theory of tensor fields in general, based on the existing formulation of integrated perturbation theory for the scalar-valued bias, generalizing it to include the tensor-valued bias. To take advantage of rotational symmetry, the formalism is constructed on the basis of irreducible decomposition of tensors, identifying physical variables which are invariant under the rotation of the coordinates system. Describing fundamental formulations and calculation techniques, this paper is expected to serve as a useful reference for future applications of perturbation theory to tensor fields in general.

We study the origin of the anomalous deep absorption in a spectrum of the SDSS J110511.15+530806.5 distant quasar (z=1.929) obtained by the Sloan Digital Sky Survey in Data Release 7 of the optical catalog. We aim to estimate the velocity of absorbing material, and we show that this material considerably affects our measurements of the black hole (BH) mass in massive quasars with the use of common virial mass estimators. The spectral shape of the quasar was modeled assuming that the accretion disk emission is influenced by a hot corona, warm skin, and absorbing material located close to the nucleus. The whole analysis was undertaken with XSPEC models and tools. The overall spectral shape was represented with the AGNSED model, while the deep absorption is well described by two Gaussians. The observed spectrum and the fitting procedure allowed us to estimate the BH mass in the quasar as $3.52 \pm 0.01 \times 10^9 M_{\odot}$, the nonzero BH spin is $a_* = 0.32 \pm 0.04$, and the accretion rate is $\dot m=0.274 \pm 0.001$. The velocities of the detected absorbers lie in the range of 6330-108135 km/s. When we consider that absorption is caused by the CIV ion, one absorber is folding toward the nucleus with a velocity of 73887 km/s. We derived a BI index of about 20300 km/s and a mass outflow rate up to 38.5 % of the source accretion rate. The high absorption observed in SDSS J110511.15+530806.5 is evidence of fast winds that place the source in the group of objects on the border with UFO (ultra-fast outflows), strong broad absorption line (BAL), and fast failed radiatively accelerated dusty outflow (FRADO). This absorption affects the BH mass measurement by two orders of magnitude as compared to virial mass estimation.

At a distance of 50 kpc, Supernova 1987A is an ideal target to study how a young supernova (SN) evolves in time. Its equatorial ring, filled with material expelled from the progenitor star about 20,000 years ago, has been engulfed with SN blast waves. Shocks heat dust grains in the ring, emitting their energy at mid-infrared (IR) wavelengths We present ground-based 10--18$\mu$m monitoring of the ring of SN 1987A from day 6067 to 12814 at a resolution of 0.5", together with SOFIA photometry at 10-30 $\mu$m. The IR images in the 2000's (day 6067-7242) showed that the shocks first began brightening the east side of the ring. Later, our mid-IR images from 2017 to 2022 (day 10952-12714) show that dust emission is now fading in the east, while it has brightened on the west side of the ring. Because dust grains are heated in the shocked plasma, which can emit X-rays, the IR and X-ray brightness ratio represent shock diagnostics. Until 2007 the IR to X-ray brightness ratio remained constant over time, and during this time shocks seemed to be largely influencing the east side of the ring. However, since then, the IR to X-ray ratio has been declining, due to increased X-ray brightness. Whether the declining IR brightness is because of dust grains being destroyed or being cooled in the post-shock regions will require more detailed modelling.

KT Eridani was a very fast nova in 2009 peaking at V=5.42 mag. We marshal large data sets of photometry to finally work out the nature of KT Eri. From the TESS light curve, as confirmed with our radial velocity curve, we find an orbital period of 2.61595 days. With our 272 spectral energy distributions from simultaneous BVRIJHK measures, the companion star has a temperature of 6200$\pm$500 K. Our century-long average in quiescence has V=14.5. With the Gaia distance (5110$^{+920}_{-430}$ parsecs), the absolute magnitude is +0.7$\pm$0.3. We converted this absolute magnitude (corrected to the disc light alone) to accretion rates, with a full integration of the alpha-disc model. This accretion rate is very high at 3.5x10$^{-7}$ solar masses per year. Our search and analysis of archival photographs shows that no eruption occurred from 1928--1954 or after 1969. With our analysis of the optical light curve, the X-ray light curve, and the radial velocity curve, we derive a white dwarf mass of 1.25$\pm$0.03 solar masses. With the high white dwarf mass and very-high accretion rate, KT Eri must require a short time to accumulate the required mass to trigger the next nova event. Our detailed calculations give a recurrence time-scale of 12 years with a total range of 5--50 years. When combined with the archival constraints, we conclude that the recurrence time-scale must be between 40--50 years. So, KT Eri is certainly a recurrent nova, with the prior eruption remaining undiscovered in a solar gap of coverage from 1959 to 1969.

Here, we present a systematic study of 59 stripped-envelope supernovae (SESNe) (including Type IIb, Ib, Ic, and transitional events) to map a possible reason for the so-called mass-discrepancy problem. In this scenario, we assume the tension between the estimated ejected masses from early- and late-time light curves (LC) is due to approximations generally used in analytical models. First, we examine the assumption that the R-band light curve is indeed a good approximation of the bolometric light curve. Next, we test the generally used assumption that rise-time to maximum brightness is equal to the effective diffusion time-scale that can be used to derive the ejecta mass from the early LC. In addition, we analyze the effect of gamma-ray and positron-leakage, which play an important role in forming the shape of the tails of SESNe, and also can be crucial to gaining the ejecta masses from the late-time LC data. Finally, we consider the effect of the different definitions of velocity that are needed for the ejecta mass calculations.

To date, over 500 short-period rocky planets with equilibrium temperatures above 1500 K have been discovered. Such planets are expected to support magma oceans, providing a direct interface between the interior and atmosphere. This provides a unique opportunity to gain insight into their interior compositions through atmospheric observations. A key process in doing such work is the vapor outgassing from the lava surface. LavAtmos is an open-source code that calculates the equilibrium chemical composition of vapor above a melt for a given composition and temperature. Results show that the produced output is in good agreement with the partial pressures obtained from experimental laboratory data as well as with other similar codes from literature. LavAtmos allows for the modeling of vaporisation of a wide range of different mantle compositions of hot-rocky exoplanets. In combination with atmospheric chemistry codes, this enables the characterization of interior compositions through atmospheric signatures.

We investigate the dynamics of clumps that coexisted with/in advection-dominated accretion flows by considering thermal conductivity. Thermal conduction can be one of the effective factors in the energy transportation of ADAFs; hence it may indirectly affect the dynamics of clumps by means of a contact force between them and their host medium. We first study the ensemble of clumps by assuming them as collision-less particles and secondly we find the orbital motion of these clouds as individuals. For both parts, clumps are subject to the gravity of the central object and a drag force. The strong coupling between clumps and ADAF leads to equality between the average treatment of the clumps and the dynamics of their background. By employing the collision-less Boltzmann equation we calculate the velocity dispersion of the clumps which turns out approximately one order of magnitude higher than the ADAF. In fact, involving drag force in such a system causes the angular momentum of the clumps can be transported outwards by the ADAF, and hence the clouds eventually will be captured at the tidal radius. The results show that the presence of thermal conduction increases the root of the averaged radial velocity square and this in turn speeds up the process of capturing the clouds through the tidal force. In the end, we focus on a typical individual cloud, the spiral orbits appear thanks to only the toroidal component of friction force. The parametric study again proves that the operation of thermal conduction helps for decreasing the lifetime of clumps.

In the near future, wide field surveys will discover 1000's of new strongly lensed quasars, and these will be monitored with unprecedented cadence by the Legacy Survey of Space and Time (LSST). Many of these will undergo caustic-crossing microlensing events over the 10-year LSST survey, in which a sharp caustic feature from a stellar body in the lensing galaxy crosses the inner accretion disk. Caustic-crossing events offer the unique opportunity to probe the vicinity of the central supermassive black hole for 100s of quasars with multi-platform follow-up triggered by LSST monitoring. To prepare for these observations, we have developed detailed simulations of caustic-crossing light curves. These employ a realistic analytic model of the inner accretion disk that reveals the strong surface brightness asymmetries introduced when fully accounting for both special- and general-relativistic effects. We demonstrate that an inflection in the caustic-crossing light curve due to the innermost stable circular orbit (ISCO) can be detected in reasonable follow-up observations and can be analyzed to constrain ISCO size. We also demonstrate that a convolutional neural network can be trained to predict ISCO size more reliably than traditional approaches and can also recover source orientation with high accuracy.

The photometric analysis and spectroscopic study of the long period low mass ratio deep contact binary KN Per were executed. The light curves of BV(RI)$_c$-band were from the Ningbo Bureau of Education and Xinjiang Observatory Telescope (NEXT) at the Xingming Observatory. Through the analysis of Wilson-Devinney (W-D) program, KN Per was found as an A-type low mass ratio deep contact binary (q=0.236, f=53.4\%). A cool spot applied on the primary component was introduced to explain the unequal maxima of the light curve. Based on the O-C analysis, we found that the rate of the increasing orbital period is $\dot{P}$ = 5.12 $\pm$ (0.30) $\times$ 10$^{-7}$ d/yr, meaning the mass transfer from the secondary component to the primary one. By analyzing the spectroscopic data, we detected chromospheric activity emission line indicators, which is corresponding to the light-curve analysis. 71 long period (P $>$ 0.5 days) contact binaries including our target were collected. The evolutionary states of all collected stars were investigated by the illustrations of mass-radius, mass-luminosity, and log M$_{T}$ - log J$_{o}$. The relations of some physical parameters were also determined. With the instability parameters of KN Per, we determined that it is a stable contact binary system at present.

Radiation pressure-driven outflows from luminous accreting supermassive black holes are an important part of active galactic nucleus (AGN) feedback. The effective Eddington limit, based on absorption of radiation by dust, not electron scattering, is revealed in the plane of AGN absorption column density $N_{\mathrm{H}}$ as a function of Eddington fraction $\lambda_{\mathrm{Edd}} = L_{\mathrm{bol}}/L_{\mathrm{Edd}}$, where a lack of objects is seen in the region where the effective limit is exceeded. Here, we conduct radiation simulation using the CLOUDY code to deduce the radiative force applied onto dusty gas at the nucleus and compare to the gravitational force to reveal the outflow region and its boundary with long-lived absorption clouds. We also investigate how the outflow condition is affected by various AGN and dust properties and distribution. As expected, the dust abundance has the largest effect on the $N_{\mathrm{H}} - \lambda_{\mathrm{Edd}}$ diagram since the higher the abundance, the more effective the radiative feedback, while the impact of the inner radius of the dusty gas shell, the shell width and the AGN spectral shape are relatively negligible. The presence of other central masses, such as a nuclear star cluster, can also make the feedback less effective. The AGN spectral energy distribution depends on the mass of the black hole and its spin. Though the effects of the AGN SED on the diagram are relatively small, the fraction of ionizing ultraviolet (UV) photons from the blackbody accretion disc is affected more by black hole mass than spin, and can influence the efficiency of radiation pressure.

We report a spectroscopic study on red supergiant stars (RSGs) in the irregular dwarf galaxy IC 1613 in the Local Group. We derive the effective temperatures ($T_\mathrm{eff}$) and metallicities of 14 RSGs by synthetic spectral fitting to the spectra observed with the MMIRS instrument on the MMT telescope for a wavelength range from 1.16 $\mu$m to 1.23 $\mu$m. A weak bimodal distribution of the RSG metallicity centered on the [Fe/H]=$-0.65$ is found, which is slightly lower than or comparable to that of the Small Magellanic Cloud (SMC). There is no evidence for spatial segregation between the metal rich ([Fe/H]$>-0.65$) and poor ([Fe/H]$<-0.65$) RSGs throughout the galaxy. The mean effective temperature of our RSG sample in IC 1613 is higher by about 250 K than that of the SMC. However, no correlation between $T_\mathrm{eff}$ and metallicity within our RSG sample is found. We calibrate the convective mixing length ($\alpha_{\mathrm{MLT}}$) by comparing stellar evolutionary tracks with the RSG positions on the HR diagram, finding that models with $\alpha_{\mathrm{MLT}}=2.2-2.4 H_P$ can best reproduce the effective temperatures of the RSGs in IC 1613 for both Schwarzschild and Ledoux convection criteria. This result supports our previous study that a metallicity dependent mixing length is needed to explain the RSG temperatures observed in the Local Group, but we find that this dependency becomes relatively weak for RSGs having a metallicity equal to or less than the SMC metallicity.

Direct imaging of exoplanets is a challenging task due to the small angular distance and high contrast relative to their host star, and the presence of quasi-static noise. We propose a new statistical method for direct imaging of exoplanets based on a likelihood ratio detection map, which assumes that the noise after the background subtraction step obeys a Laplacian distribution. We compare the method with two detection approaches based on signal-to-noise ratio (SNR) map after performing the background subtraction by the widely used Annular Principal Component Analysis (AnnPCA). The experimental results on the Beta Pictoris data set show the method outperforms SNR maps in terms of achieving the highest true positive rate (TPR) at zero false positive rate (FPR).

Space systems miniaturization has been increasingly popular for the past decades, with over 1600 CubeSats and 300 sub-CubeSat sized spacecraft estimated to have been launched since 1998. This trend towards decreasing size enables the execution of unprecedented missions in terms of quantity, cost and development time, allowing for massively distributed satellite networks, and rapid prototyping of space equipment. Pocket-sized spacecraft can be designed in-house in less than a year and can reach weights of less than 10g, reducing the considerable effort typically associated with orbital flight. However, while Systems Engineering methodologies have been proposed for missions down to CubeSat size, there is still a gap regarding design approaches for picosatellites and smaller spacecraft, which can exploit their potential for iterative and accelerated development. In this paper, we propose a Systems Engineering methodology that abstains from the classic waterfall-like approach in favor of agile practices, focusing on available capabilities, delivery of features and design "sprints". Our method, originating from the software engineering disciplines, allows quick adaptation to imposed constraints, changes to requirements and unexpected events (e.g. chip shortages or delays), by making the design flexible to well-defined modifications. Two femtosatellite missions, currently under development and due to be launched in 2023, are used as case studies for our approach, showing how miniature spacecraft can be designed, developed and qualified from scratch in 6 months or less. We claim that the proposed method can simultaneously increase confidence in the design and decrease turnaround time for extremely small satellites, allowing unprecedented missions to take shape without the overhead traditionally associated with sending cutting-edge hardware to space.

In this paper, we modeled the dynamics and radiation physics of the rarity event GRB 221009A afterglow in detail. Based on the analysis results of the {\tt ASGARD} package we developed, the afterglow data of GRB 221009A strongly favors the origin of a relativistic jet propagating in a stellar-wind-dominated environment. Therefore, GRB 221009A is a typical lesson for the very high energy (VHE) afterglow of long gamma-ray burst in a wind environment. We also show the broadband spectral energy distribution (SED) analysis results of GRB 221009A, and find that the synchrotron self-Compton (SSC) radiation component of GRB 221009A is very bright in the $0.1-10$ TeV band. The integrated SED of $0-2000$ s after {\em Fermi}/GBM trigger shows that the bright SSC component can be easily observed in TeV band, above the detection sensitivity of LHASSO, MAGIC and CTA. We predict that the SSC radiation peak flux of GRB 221009A in the first 2000 second integral SED is $\sim 10^{-7}~\rm erg~cm^{-2}~s^{-1}$, with a peak energy close to 300 GeV. Furthermore, we find that the inclusion of GeV observations could break the degeneracy between model parameters, highlighting the significance of high-energy observations in determining accurate parameters for GRB afterglows.

The multi-messenger detection of GW170817 showed that binary neutron star (BNS) mergers are progenitors of (at least some) short gamma-ray bursts (GRBs), and that short GRB jets (and their afterglows) can have structures (and observational properties) more complex than predicted by the standard top-hat jet scenario. Indeed, the emission from the structured jet launched in GW170817 peaked in the radio band (cm wavelengths) at about 100 d since merger - a timescale much longer than the typical time span of radio follow-up observations of short GRBs. Moreover, radio searches for a potential late-time radio flare from the fast tail of the neutron-rich debris that powered the kilonova associated with GW170817 (AT2017gfo) have extended to even longer timescales (years after the merger). In light of this, here we present the results of an observational campaign targeting a sample of seven, years-old GRBs in the Swift/BAT sample with no redshift measurements and no promptly-identified X-ray counterpart. Our goal is to assess whether this sample of short GRBs could harbor nearby BNS mergers, searching for the late-time radio emission expected from their ejecta. We found one radio candidate counterpart for one of the GRBs in our sample, GRB111126A, though an origin related to emission from star formation or from an AGN in its host galaxy cannot be excluded without further observations.

When Type Ia supernovae are used to infer cosmological parameters, their luminosities are compared to those from a homogeneous cosmology. In this note we propose a test to examine to what degree SN Ia have been observed on lines of sight where the average matter density is \textit{not} representative of the homogeneous background. We use a modification to the Tripp estimator to explicitly de-lens SN Ia magnitudes, and show that it reduces scatter. We apply our test to the Pantheon SN Ia compilation, and find two redshift bins which indicate a moderate bias to over-density at $\sim 2\sigma$. Using our revised estimator, the effect on cosmological parameters from Pantheon in $\Lambda$CDM is however small with a change in mean value from $\Omega_{\rm m} = 0.317 \pm 0.027$ (baseline) to $\Omega_{\rm m} = 0.312 \pm 0.025$ (de-lensed). For the Flat $w$CDM case it is $\Omega_{\rm m} = 0.332 \pm 0.049$ and $w = -1.16 \pm 0.16$ (baseline) versus $\Omega_{\rm m} = 0.316 \pm 0.048$ and $w = -1.12 \pm 0.15$ (de-lensed). The change is larger if maximum likelihood values are used. We note that the effect of lensing on cosmological parameters may be larger for future high-z surveys.

We are using archived data from HST of transiting exoplanet L~98-59~b to place constraints on its potentially hot atmosphere. We analyze the data from five transit visits and extract the final combined transmission spectrum using Iraclis. Then we use the inverse atmospheric retrieval code TauREx to analyze the combined transmission spectrum. There is a weak absorption feature near 1.40~$\mu m$ and 1.55~$\mu m$ in the transmission spectrum, which can be modeled by a cloudy atmosphere with abundant HCN. However, the unrealistically high abundance of HCN derived cannot be explained by any equilibrium chemical model with reasonable assumptions. Thus, the likeliest scenario is that L~98-59~b has a flat, featureless transmission spectrum in the WFC3/G141 bandpass due to a thin atmosphere with high mean molecular weight, an atmosphere with an opaque aerosol layer, or no atmosphere, and it is very unlikely for L~98-59~b to have a clear hydrogen-dominated primary atmosphere. Due to the narrow wavelength coverage and low spectral resolution of HST/WFC3 G141 grism observation, we cannot tell these different scenarios apart. Our simulation shows future higher precision measurements over wider wavelengths from the James Webb Space Telescope (JWST) can be used to better characterize the planetary atmosphere of L~98-59~b.

The dynamics of galaxies in an expanding universe is often determined for gravitational and dark matter in an Einstein-de Sitter universe, or alternatively by modifying the gravitational long-range attractions in the Newtonian dynamics (MOND). Here the time evolution of galaxies is determined by simulations of systems with pure gravitational forces by classical Molecular Dynamic simulations. A time reversible algorithm for formation and aging of gravitational systems by self-assembly of baryonic objects, recently derived (Eur. Phys. J. Plus 2022, 137:99), is extended to include the Hubble expansion of the space. The algorithm is stable for billions of time steps without any adjustments. The algorithm is used to simulate simple models of the Milky Way with the Hubble expansion of the universe, and the galaxies are simulated for times which corresponds to more than 25 Gyr. The rotating galaxies lose bound objects from time to time, but they are still stable at the end of the simulations. The simulations indicate that the explanation for the dynamics of galaxies may be that the universe is very young in cosmological times. Although the models of the Milky Way are rather stable at 13-14 Gyr, which corresponds to the cosmological time of the universe, the Hubble expansion will sooner or later release the objects in the galaxies. But the simulations indicate that this will first happen in a far away future.

We cross-matched the 4XMM-DR10 catalog with the HyperLEDA database and obtained the new sample of galaxies that contain X-ray sources. Excluding duplicate observations and false matches, we present a total of 7759 galaxies with X-ray sources. In the current work, we present general properties of the sample: namely the distribution in equatorial coordinates, radial velocity distribution, morphological type, and X-ray fluxes. The sample includes morphological classification for 5241 galaxies with X-ray emission, almost half of which, 42\%, are elliptical (E, E-S0). Most galaxies in the sample have nuclear X-ray emission (6313 or 81\%), and the remaining 1443 (19\%) present X-ray emission from the host galaxy. This sample can be used for future deep studies of multi wavelengths properties of the galaxies with X-ray emission.

Knowledge of the interior density distribution of an asteroid can reveal its composition and constrain its evolutionary history. However, most asteroid observational techniques are not sensitive to interior properties. We investigate the interior constraints accessible through monitoring variations in angular velocity during a close encounter. We derive the equations of motion for a rigid asteroid's orientation and angular velocity to arbitrary order and use them to generate synthetic angular velocity data for a representative asteroid on a close Earth encounter. We develop a toolkit AIME (Asteroid Interior Mapping from Encounters) which reconstructs asteroid density distribution from these data, and we perform injection-retrieval tests on these synthetic data to assess AIME's accuracy and precision. We also perform a sensitivity analysis to asteroid parameters (e.g., asteroid shape and orbital elements), observational set-up (e.g., measurement precision and cadence), and the mapping models used. We find that high precision in rotational period estimates (<~0.27 seconds) are necessary for each cadence, and that low perigees (<~18 Earth radii) are necessary to resolve large-scale density non-uniformities to uncertainties ~0.1% of the local density under some models.

Our aim is to understand the effect of environment to galaxy quenching in various local and global environments. We focus on galaxies with very old stellar populations (VO galaxies), typically found in the centers of clusters and groups, and search for such galaxies in the lowest global density environments, watersheds between superclusters. We use the Sloan Digital Sky Survey MAIN galaxy sample to calculate the luminosity-density field and get global density field, to determine groups and filaments, and to obtain data on galaxy properties. We divide groups into low- and high-luminosity groups based on the highest luminosity of groups in the watershed region, $L_{gr} = 15 \times10^{12} h^{-2} L_{sun}$. Our study shows that the global density is most strongly related to the richness of galaxy groups. Its influence on the overall star formation quenching in galaxies is less strong. Correlations between the morphological properties of galaxies and the global density field are the weakest. The watershed regions are populated mostly by single galaxies (70% of all galaxies there), and by low-luminosity groups. Still, approximately one-third of all galaxies in the watershed regions are VO galaxies. They have lower stellar masses, smaller stellar velocity dispersions, and stellar populations that are up to 2Gyr younger than those of VO galaxies in other global environments. In higher density global environments, the morphological properties of galaxies are very similar. Differences in galaxy properties are the largest between satellites and brightest group galaxies. Our results suggest that galaxy evolution is determined by the birthplace of galaxies in the cosmic web, and mainly by internal processes which lead to the present-day properties of galaxies. This may explain the similarity of (VO) galaxies in extremely different environments.

The dark matter annihilation cross section can be amplified by orders of magnitude if the annihilation occurs into a narrow resonance, or if the dark-matter particles experience a long-range force before annihilation (Sommerfeld effect). We show that when both enhancements are present they factorize completely, that is, all long-distance non-factorizable effects cancel at leading order in the small-velocity and narrow-width expansion. We then investigate the viability of ``super-resonant'' annihilation from the coaction of both mechanisms in Standard Model Higgs portal and simplified MSSM-inspired dark-matter scenarios.

The reconstruction of an inflationary universe in the context of one $f(\phi)T$ gravity, in which $T$ corresponds to the trace of energy momentum tensor is studied. To realize this reconstruction during the inflationary epoch, we consider as attractor the scalar spectral index $n_s$ in terms of of the number of $e$-folds $N$, in the framework of the slow-roll approximation. By assuming a specific function $f(\phi)$ together with the simplest attractor $n_s(N)$, we find different expressions for the reconstructed effective potential $V(\phi)$. Additionally, we analyze the era of reheating occurs after of the reconstruction obtained during the inflationary epoch. In this scenario we determine the duration and temperature during the reheating epoch, in terms of the equation of state parameter and the observational parameters. In this context, the different parameters associated to the reconstructed model are restricted during the scenarios of inflation and reheating by considering the recent astronomical observations.

In this work we consider the effects of Lorentz Invariance Violation over the observed flux of very high-energy neutrinos. For that, we study the neutrino propagation in a Modified Dispersion Relation scenario with a superluminal velocity. This makes the neutrino unstable and causes a cut-off in the flux of detected neutrinos. Using simple models, one can approximate the location of the cut-off as function of the parameters of new physics and the closest source.

We study the electrical and thermal conductivity of iron at high pressures using time-dependent density functional theory. In doing so, we particularly consider the impact of a Hubbard correction (+\textit{U}) specifically for regions where strong electron correlations are present. Using the TDDFT+U methodology, we examine the anisotropy in the thermal conductivity of HCP iron, which may provide insights into the transport properties at conditions relevant to the core-mantle boundary and the interior of the Earth.

We review the distortions of spectra of relic neutrinos due to the interactions with electrons, positrons, and neutrinos in the early universe. We solve integro-differential kinetic equations for the neutrino density matrix, including vacuum three-flavor neutrino oscillations, oscillations in electron and positron background, a collision term and finite temperature corrections to electron mass and electromagnetic plasma up to the next-to-leading order $\mathcal{O}(e^3)$. After that, we estimate the effects of the spectral distortions in neutrino decoupling on the number density and energy density of the Cosmic Neutrino Background (C$\nu$B) in the current universe, and discuss the implications of these effects on the capture rates in direct detection of the C$\nu$B on tritium, with emphasis on the PTOLEMY-type experiment. In addition, we find a precise value of the effective number of neutrinos, $N_{\rm eff}=3.044$. However, QED corrections to weak interaction rates at order $\mathcal{O}(e^2 G_F^2)$ and forward scattering of neutrinos via their self-interactions have not been precisely taken into account in the whole literature so far. Recent studies suggest that these neglections might induce uncertainties of $\pm(10^{-3} - 10^{-4})$ in $N_{\rm eff}$.

We extend the Cecotti-Kallosh model of Starobinsky inflation in supergravity by adding a holomorphic function to the superpotential in order to generate a large peak in the power spectrum of scalar (curvature) perturbations. In our approach, the singular non-canonical kinetic terms are largely responsible for inflation (as an attractor solution), whereas the superpotential is engineered to generate a production of primordial black holes. We study the cases with (i) a linear holomorphic function, (ii) a quadratic holomorphic function, and (iii) an exponential holomorphic function, as regards the dependence of inflation and primordial black holes production upon parameters of those functions and initial conditions, as well as verify viability of inflation with our superpotentials. We find that an efficient production of primordial black holes consistent with CMB measurements is only possible in the second (ii) case. We calculate the masses of the produced primordial black holes and find that they are below the Hawking (black hole) evaporation limit, so that they cannot be part of the current dark matter in our Universe.

Solar Energetic Particles (SEPs) are charged particles accelerated within the solar atmosphere or the interplanetary space by explosive phenomena such as solar flares or Coronal Mass Ejections (CMEs). Once injected into the interplanetary space, they can propagate towards Earth, causing space weather related phenomena. For their analysis, interplanetary in-situ measurements of charged particles are key. The recently expanded spacecraft fleet in the heliosphere not only provides much-needed additional vantage points, but also increases the variety of missions and instruments for which data loading and processing tools are needed. This manuscript introduces a series of Python functions that will enable the scientific community to download, load, and visualize charged particle measurements of the current space missions that are especially relevant to particle research as time series or dynamic spectra. In addition, further analytical functionality is provided that allows the determination of SEP onset times as well as their inferred injection times. The full workflow, which is intended to be run within Jupyter Notebooks and can also be approachable for Python laymen, will be presented with scientific examples. All functions are written in Python, with the source code publicly available at GitHub under a permissive license. Where appropriate, available Python libraries are used, and their application is described.

The natural cycles of the surface-to-atmosphere fluxes of carbon dioxide (CO$_2$) and other important greenhouse gases are changing in response to human influences. These changes need to be quantified to understand climate change and its impacts, but this is difficult to do because natural fluxes occur over large spatial and temporal scales. To infer trends in fluxes and identify phase shifts and amplitude changes in flux seasonal cycles, we construct a flux-inversion system that uses a novel spatially varying time-series decomposition of the fluxes, while also accommodating physical constraints on the fluxes. We incorporate these features into the Wollongong Methodology for Bayesian Assimilation of Trace-gases (WOMBAT, Zammit-Mangion et al., Geosci. Model Dev., 15, 2022), a hierarchical flux-inversion framework that yields posterior distributions for all unknowns in the underlying model. We apply the new method, which we call WOMBAT v2.0, to a mix of satellite observations of CO$_2$ mole fraction from the Orbiting Carbon Observatory-2 (OCO-2) satellite and direct measurements of CO$_2$ mole fraction from a variety of sources. We estimate the changes to CO$_2$ fluxes that occurred from January 2015 to December 2020, and compare our posterior estimates to those from an alternative method based on a bottom-up understanding of the physical processes involved. We find substantial trends in the fluxes, including that tropical ecosystems trended from being a net source to a net sink of CO$_2$ over the study period. We also find that the amplitude of the global seasonal cycle of ecosystem CO$_2$ fluxes increased over the study period by 0.11 PgC/month (an increase of 8%), and that the seasonal cycle of ecosystem CO$_2$ fluxes in the northern temperate and northern boreal regions shifted earlier in the year by 0.4-0.7 and 0.4-0.9 days, respectively (2.5th to 97.5th posterior percentiles).

Two-field mimetic gravity was recently realized by looking at the singular limit of the conformal transformation between the auxiliary metric and the physical metric with two scalar fields involved. In this paper, we reanalyze the singular conformal limit and find a more general solution for the conformal factor A, which greatly broadens the form of two-field mimetic constraint and thus extends the two-field mimetic gravity. We find the general setup still mimics the roles of dark matter at the cosmological background level. Moreover, we extend the action by introducing extra possible term for phenomenological interests. Surprisingly, some special cases are found to be equivalent to general relativity, k-essence theory and Galileon theory. Finally, we further extend the theory by allowing the expression of mimetic constraint to be arbitrary without imposed condition, and show that the dark matter-like behavior is unaffected even in this extension.

I show that a gravitational slingshot using a stellar-mass black hole (BH) orbiting SgrA* could launch robotic spacecraft toward M31 at $0.1\,c$, a speed that is ultimately limited by the tensile strength of steel and the BH mass, here conservatively estimated as $m_{\rm bh}=5\,M_\odot$. The BH encounter must be accurate to $\lesssim 1\,$km, despite the fact that the BH is dark. Navigation guided by gravitational microlensing can easily achieve this. Deceleration into M31 would rely on a similar engagement (but in reverse) with an orbiting BH near the M31 center. Similarly for a return trip, if necessary. Colonization of M31 planets on 50 Myr timescales is therefore feasible provided that reconstruction of humans, trans-humans, or androids from digital data becomes feasible in the next few Myr. The implications for Fermi's Paradox (FP) are discussed. FP is restated in a more challenging form. The possibility of intergalactic colonization on timescales much shorter than the age of Earth significantly tightens FP. It can thereby impact our approach to astrobiology on few-decade timescales. I suggest using a network of tight white-dwarf-binary "hubs" as the backbone of a $0.002\,c$ intra-Galactic transport system, which would enable complete exploration of the Milky Way (hence full measurement of all non-zero terms in the Drake equation) on 10 Myr timescales. Such a survey would reveal the reality and/or severity of an "intergalactic imperative".

An expanding civilization could rapidly spread through the galaxy, so the absence of extraterrestrial settlement in the solar system implies that such expansionist civilizations do not exist. This argument, often referred to as the Fermi paradox, typically assumes that expansion would proceed uniformly through the galaxy, but not all stellar types may be equally useful for a long-lived civilization. We suggest that low-mass stars, and K-dwarf stars in particular, would be ideal migration locations for civilizations that originate in a G-dwarf system. We use a modified form of the Drake Equation to show that expansion across all low-mass stars could be accomplished in 2 Gyr, which includes waiting time between expansion waves to allow for a close approach of a suitable destination star. This would require interstellar travel capabilities of no more than ~0.3 ly to settle all M-dwarfs and ~2 ly to settle all K-dwarfs. Even more rapid expansion could occur within 2 Myr, with travel requirements of ~10 ly to settle all M-dwarfs and ~50 ly to settle all K-dwarfs. The search for technosignatures in exoplanetary systems can help to place constraints on the presence of such a "low-mass Galactic Club" in the galaxy today.

We present a model of Cosmological Electroweak Symmetry Breaking (CEWSB), where a Higgs-like field and a cosmological background of weak boson gauge fields interact with gravity to realize the epoch of cosmic inflation, which is then followed by a Higgs resonance preheating. As a result, the scale of electroweak symmetry breaking is linked with the end of inflation. The theory is equipped with a shift symmetry that can protect the Higgs mass, and it has close semblance to natural inflation and its variants. As the Higgs field's amplitude decays at the end of inflation, its mass emerges. The model has a built in Higgs self-resonance preheating mechanism which leads to the possible emergence of the cosmic microwave background (CMB) due to resonant Higgs, quark and lepton production after inflation. We provide a pathway to implement a similar mechanism with the realistic Higgs-doublet of the standard electroweak theory and discuss phenomenological considerations.

We examine thermal warm dark matter (WDM) models that are being probed by current constraints, and the relationship between the particle dark matter spin and commensurate thermal history. Two primary classes are examined: spin-1/2 particles (e.g., thermalized sterile neutrinos, axinos) and partially to fully thermally populated spin-3/2 particles (e.g., gravitinos). We present new transfer function fits for thermal WDM candidates in particle mass regimes beyond the range of previous work, and at the scales of current constraints and beyond. Importantly, we find that the standard, predominantly used, spin-1/2, thermal WDM particle produces a "cooler" transfer function than that determined in previous work, so that current work may slightly overstate their constraints on particle mass. We also analyze the entropy requirements for these WDM models to successfully produce observed dark matter densities. We explore the early Universe physics of gravitinos as either partially thermalized or fully thermalized species, which considerably changes the particle dark matter candidates' thermalization history and effects on structure formation.

Accelerating supermassive black holes, connected to cosmic strings, could contribute to structure formation and get captured by galaxies if their velocities are small. This would mean that the acceleration of these black holes is small too. Such a slow acceleration has no significant effect on the shadow of such supermassive black holes. We also show that, for slowly accelerating black holes, the angular position of images in the gravitational lensing effects do not change significantly. We propose a method to observe the acceleration of these black holes through the gravitational lensing. The method is based on the observation that differential time delays associated with the images are substantially different with respect to the case of non-accelerating black holes.