Papers for Wednesday, Aug 03 2022

Using one of the largest volumes of the hydrodynamical cosmological simulation suit Magneticum, we study the evolution of protoclusters identified at redshift = 4, with properties similar to SPT2349-56. We identify 42 protoclusters in the simulation, as massive and equally rich in substructures as observed, confirming that these structures are already virialized. The dynamics of the internally fast rotating member galaxies within these protoclusters resembles observations, merging rapidly to form the cores of the BCGs of the assembling clusters. Half of the gas reservoir of these structures is in a hot phase, with the metal-enrichment at a very early stage. These systems show a good agreement with the observed amount of cold star-forming gas, largely enriched to solar values. We predict that some of the member galaxies are already quenched at z = 4, rendering them undetectable through measurements of their gas reservoir. Tracing the evolution of protoclusters reveals that none of the typical mass indicators at high redshift are good tracers to predict the present-day mass of the system. We find that none of the simulated protoclusters with properties as SPT2349-56 at z = 4.3, are among the top ten most massive clusters at redshift z = 0, with some barely reaching masses of M = 2 x 10^14Msun. Although the average star-formation and mass-growth rates in the simulated galaxies match observations at high redshift reasonably well, the simulation fails to reproduce the extremely high total star-formation rates within observed protoclusters, indicating that the sub-grid models are lacking the ability to reproduce higher star-formation efficiency (or lower depletion timescales).

We apply the clustering algorithm HDBSCAN on the Gaia early third data release astrometry combined with the Gaia second data release radial velocity measurements of almost 5.5 million stars to identify the local stellar kinematic substructures in the solar neighborhood. Understanding these structures helps build a more complete picture of the formation of the Milky Way, as well as an empirical phase space distribution of dark matter that would inform detection experiments. The main goal of this study is to provide a list of the most stable clusters, by taking into account the measurement uncertainties and studying the stability of the clustering results. We apply the clustering algorithm in two spaces, in velocity space in order to study recently accreted structures, and in action-angle space to find phase-mixed structures. We find 23 (6) robust clusters in velocity space (action-angle space) that are consistently not associated with noise. They are attributed to the known structures: the Gaia Sausage-Enceladus, the Helmi Stream, and globular cluster NGC 3201 are found in both spaces, while NGC 104 and the thick disk (Sequoia) are identified in velocity space (action-angle space). We discuss the kinematic properties of these structures and study whether many of the small clusters belong to a similar larger cluster based on their chemical abundances. Although we do not identify any new structures, we find that the HDBSCAN member selection of already known structures is unstable to input kinematics of the stars when resampled within their uncertainties. We therefore present the most stable subset of local kinematic structures, which are consistently identified by the clustering algorithm, and emphasize the need to take into account error propagation during both the manual and automated identification of stellar structures, both for existing ones as well as future discoveries. (abridged)

We implement a new observational method for mapping the aliphatic hydrocarbon content in the solid phase in our Galaxy, based on spectrophotometric imaging of the 3.4 $\mu$m absorption feature from interstellar dust. We previously demonstrated this method in a field including the Galactic Centre cluster. We applied the method to a new field in the Galactic centre where the 3.4 $\mu$m absorption feature has not been previously measured and we extended the measurements to a field in the Galactic plane to sample the diffuse local interstellar medium, where the 3.4 $\mu$m absorption feature has been previously measured. We have analysed 3.4 $\mu$m optical depth and aliphatic hydrocarbon column density maps for these fields. Optical depths are found to be reasonably uniform in each field, without large source-to-source variations. There is, however, a weak trend towards increasing optical depth in a direction towards $b=0^{\circ}$ in the Galactic centre. The mean value of column densities and abundances for aliphatic hydrocarbon were found to be about several $\rm \times 10^{18} \, cm^{-2}$ and several tens $\times 10^{-6}$, respectively for the new sightlines in the Galactic plane. We conclude that at least 10-20% of the carbon in the Galactic plane lies in aliphatic form.

The magnetorotational instability (MRI) is an important process in sufficiently ionized accretion disks, as it can create turbulence that acts as an effective viscosity, mediating angular momentum transport. Due to its local nature, it is often analyzed in the shearing box approximation with Eulerian methods, which otherwise would suffer from large advection errors in global disk simulations. In this work, we report on an extensive study that applies the quasi-Lagrangian, moving-mesh code ${\rm \small AREPO}$, combined with the Dedner cleaning scheme to control deviations from $\nabla \cdot B = 0$, to the problem of magnetized flows in shearing boxes. We find that we can resolve the analytical linear growth rate of the MRI with mean background magnetic field well. In the zero net flux case, there is a threshold value for the strength of the divergence cleaning above which the turbulence eventually dies out, and in contrast to previous Eulerian simulations, the strength of the MRI does not decrease with increasing resolution. If we increase the vertical aspect ratio of our box we find the mean-field dynamo described in Shi et al. (2016), as well as an active shear current effect that can sustain MRI turbulence for at least 200 orbits. In stratified simulations, we obtain an active $\alpha \omega$ dynamo and the characteristic butterfly diagram, again for at least 200 orbits. Our results compare well with previous results obtained with static grid codes such as ${\rm\small ATHENA}$. We thus conclude that ${\rm \small AREPO}$ represents a particularly attractive alternative for global disk simulations, where the method benefits from its quasi-Lagrangian nature, as well as for shearing box simulations with large density variations, where ${\rm \small AREPO}$'s continuously adaptive resolution is advantageous.

In this paper we present COMET, a Gaussian process emulator of the galaxy power spectrum multipoles in redshift-space. The model predictions are based on one-loop perturbation theory and we consider two alternative descriptions of redshift-space distortions: one that performs a full expansion of the real- to redshift-space mapping, as in recent effective field theory models, and another that preserves the non-perturbative impact of small-scale velocities by means of an effective damping function. The outputs of COMET can be obtained at arbitrary redshifts (up to $z \sim 3$), for arbitrary fiducial background cosmologies, and for a large parameter space that covers the shape parameters $\omega_c$, $\omega_b$, and $n_s$, as well as the evolution parameters $h$, $A_s$, $\Omega_K$, $w_0$, and $w_a$. This flexibility does not impair COMET's accuracy, since we exploit an exact degeneracy between the evolution parameters that allows us to train the emulator on a significantly reduced parameter space. While the predictions are sped up by at least two orders of magnitude, validation tests reveal an accuracy of $0.1\,\%$ for the monopole and quadrupole ($0.3\,\%$ for the hexadecapole), or alternatively, better than $0.25\,\sigma$ for all three multipoles in comparison to statistical uncertainties expected for the Euclid survey with a tenfold increase in volume. We show that these differences translate into shifts in mean posterior values that are at most of the same size, meaning that COMET can be used with the same confidence as the exact underlying models. COMET is a publicly available Python package that also provides the tree-level bispectrum multipoles in redshift-space and Gaussian covariance matrices.

Many gamma-ray bursts (GRBs) have been observed from radio wavelengths, and a few at very-high energies (VHEs, > 100GeV). The HAWC gamma-ray observatory is well suited to study transient phenomena at VHEs due to its large field of view and duty cycle. These features allow for searches of VHE emission and can probe different model assumptions of duration and spectra. In this paper, we use data collected by HAWC between December 2014 and May 2020 to search for emission in the energy range from 80 to 800 GeV coming from a sample 47 short GRBs that triggered the Fermi, Swift and Konus satellites during this period. This analysis is optimized to search for delayed and extended VHE emission within the first 20 s of each burst. We find no evidence of VHE emission, either simultaneous or delayed, with respect to the prompt emission. Upper limits (90% confidence level) derived on the GRB fluence are used to constrain the synchrotron self-Compton forward-shock model. Constraints for the interstellar density as low as $10^{-2}$ cm$^{-3}$ are obtained when assuming z=0.3 for bursts with the highest keV-fluences such as GRB 170206A and GRB 181222841. Such a low density makes observing VHE emission mainly from the fast cooling regime challenging.

Constraining the Galactic foregrounds with multi-frequency Cosmic Microwave Background (CMB) observations is an essential step towards ultimately reaching the sensitivity to measure primordial gravitational waves (PGWs), the sign of inflation after the Big-Bang that would be imprinted on the CMB. The BICEP Array telescope is a set of multi-frequency cameras designed to constrain the energy scale of inflation through CMB B-mode searches while also controlling the polarized galactic foregrounds. The lowest frequency BICEP Array receiver (BA1) has been observing from the South Pole since 2020 and provides 30 GHz and 40 GHz data to characterize the Galactic synchrotron in our CMB maps. In this paper, we present the design of the BA1 detectors and the full optical characterization of the camera including the on-sky performance at the South Pole. The paper also introduces the design challenges during the first observing season including the effect of out-of-band photons on detectors performance. It also describes the tests done to diagnose that effect and the new upgrade to minimize these photons, as well as installing more dichroic detectors during the 2022 deployment season to improve the BA1 sensitivity. We finally report background noise measurements of the detectors with the goal of having photon noise dominated detectors in both optical channels. BA1 achieves an improvement in mapping speed compared to the previous deployment season.

We use our cluster population model, cBHBd, to explore the mass distribution of merging black hole binaries formed dynamically in globular clusters. We include in our models the effect of mass growth through hierarchical mergers and compare the resulting distributions to those inferred from the third gravitational wave transient catalogue. We find that none of our models can reproduce the peak at $m_1\simeq 10M_\odot$ in the primary black hole mass distribution that is inferred from the data. This disfavours a scenario where most of the detected sources are formed in globular clusters. On the other hand, a globular cluster origin can account for the inferred secondary peak at $m_1\simeq 35M_\odot$, which requires that the most massive clusters form with half-mass densities $\rho_{\rm h,0} \gtrsim 10^4 M_\odot \rm pc^{-3}$. Finally, we find that the lack of a high mass cut--off in the inferred mass distribution can be also explained by the repopulation of an initial mass gap through hierarchical mergers. Matching the inferred merger rate above $\simeq 50M_\odot$ requires both initial cluster densities $\rho_{\rm h,0} \gtrsim 10^4 M_\odot \rm pc^{-3}$, and that black holes form with nearly zero spin. A hierarchical merger scenario makes specific predictions for the appearance and position of multiple peaks in the black hole mass distribution, which can be tested against future data.

Single electron Sensitive Read Out (SiSeRO) is a novel on-chip charge detector output stage for charge-coupled device (CCD) image sensors. Developed at MIT Lincoln Laboratory, this technology uses a p-MOSFET transistor with a depleted internal gate beneath the transistor channel. The transistor source-drain current is modulated by the transfer of charge into the internal gate. At Stanford, we have developed a readout module based on the drain current of the on-chip transistor to characterize the device. Characterization was performed for a number of prototype sensors with different device architectures, e.g. location of the internal gate, MOSFET polysilicon gate structure, and location of the trough in the internal gate with respect to the source and drain of the MOSFET (the trough is introduced to confine the charge in the internal gate). Using a buried-channel SiSeRO, we have achieved a charge/current conversion gain of >700 pA per electron, an equivalent noise charge (ENC) of around 6 electrons root mean square (RMS), and a full width half maximum (FWHM) of approximately 140 eV at 5.9 keV at a readout speed of 625 Kpixel/s. In this paper, we discuss the SiSeRO working principle, the readout module developed at Stanford, and the characterization test results of the SiSeRO prototypes. We also discuss the potential to implement Repetitive Non-Destructive Readout (RNDR) with these devices and the preliminary results which can in principle yield sub-electron ENC performance. Additional measurements and detailed device simulations will be essential to mature the SiSeRO technology. However, this new device class presents an exciting technology for next generation astronomical X-ray telescopes requiring fast, low-noise, radiation hard megapixel imagers with moderate spectroscopic resolution.

New experiments that target the B-mode polarization signals in the Cosmic Microwave Background require more sensitivity, more detectors, and thus larger-aperture millimeter-wavelength telescopes, than previous experiments. These larger apertures require ever larger vacuum windows to house cryogenic optics. Scaling up conventional vacuum windows, such as those made of High Density Polyethylene (HDPE), require a corresponding increase in the thickness of the window material to handle the extra force from the atmospheric pressure. Thicker windows cause more transmission loss at ambient temperatures, increasing optical loading and decreasing sensitivity. We have developed the use of woven High Modulus Polyethylene (HMPE), a material 100 times stronger than HDPE, to manufacture stronger, thinner windows using a pressurized hot lamination process. We discuss the development of a specialty autoclave for generating thin laminate vacuum windows and the optical and mechanical characterization of full scale science grade windows, with the goal of developing a new window suitable for BICEP Array cryostats and for future CMB applications.

We here study the transfer process of material from one hemisphere to the other (deposition of airfall material) on an active comet nucleus, specifically 67P/Churyumov-Gerasimenko. Our goals are to: 1) quantify the thickness of the airfall debris layers and how it depends on the location of the target area, 2) determine the amount of $\mathrm{H_2O}$ and $\mathrm{CO_2}$ ice that are lost from icy dust assemblages of different sizes during transfer through the coma, and 3) estimate the relative amount of vapor loss in airfall material after deposition in order to understand what locations are expected to be more active than others on the following perihelion approach. We use various numerical simulations, that include orbit dynamics, thermophysics of the nucleus and of individual coma aggregates, coma gas kinetics and hydrodynamics, as well as dust dynamics due to gas drag, to address these questions. We find that the thickness of accumulated airfall material varies substantially with location, and typically is of the order $0.1$-$1\,\mathrm{m}$. The airfall material preserves substantial amounts of water ice even in relatively small (cm-sized) coma aggregates after a rather long ($12\,\mathrm{h}$) residence in the coma. However, $\mathrm{CO_2}$ is lost within a couple of hours even in relatively large (dm-sized) aggregates, and is not expected to be an important component in airfall deposits. We introduce reachability and survivability indices to measure the relative capacity of different regions to simultaneously collect airfall and to preserve its water ice until the next perihelion passage, thereby grading their potential of contributing to comet activity during the next perihelion passage.

In this work, we analyze a sample of $\sim$4000 massive ($M_*\geq 10^{11} M_\odot$ at $z=0$) galaxies in TNG300, the $(300 \mathrm{Mpc})^3$ box of the IllustrisTNG simulation suite. We characterize the shape and kinematics of these galaxies with a focus on the kinematic misalignment ($\Psi_\mathrm{int}$) between the angular momentum (AM) and morphological major axis. We find that the traditional purely shape- or kinematics-based classifications are insufficient to characterize the diversity of our sample and define a new set of classes based on the rates of change of the galaxies' morphological and kinematic axes. We show that these classes are mostly stable over time and correspond to six distinct populations of galaxies: the rapid AM reorienters (58% of our sample), unsettled galaxies (10%), spinning disks (10%), twirling cigars (16%), misaligned slow reorienters (3%), and regular prolate rotators (galaxies that display major axis rotation; 2%). We demonstrate that the most-recent significant (mass-ratio $\mu>1/10$) mergers of these galaxies are the primary cause for their present-day properties and find that these mergers are best characterized at the point of the satellite's final infall -- that is, much closer to the final coalescence than has been previously thought. We show that regular prolate rotators evolve from spinning disk progenitors that experience a radial merger along their internal AM direction. Finally, we argue that these regular prolate rotators are distinct from the similarly-sized population of rapid AM reorienters with large $\Psi_\mathrm{int}$, implying that a large $\Psi_\mathrm{int}$ is not a sufficient condition for major axis rotation.

We study the impact of stellar cooling due to light axion emission on the formation and evolution of black hole binaries, via stable mass transfer and the common envelope scenario. We find that in the presence of light axion emission, no binary black hole mergers are formed with black holes in the lower mass gap ($M_{\rm BH} < 4 {\rm M}_\odot $) via the common envelope formation channel. In some systems, this happens because axions prevent Roche lobe overflow. In others, they prevent the common envelope from being ejected. Our results apply to axions with couplings $ g_{a \gamma} \gtrsim 10^{10}\, \rm GeV^{-1}$ (to photons) or $\alpha_{ae} \gtrsim 10^{-26} $ (to electrons) and masses $ m_a \ll 10 \, \rm keV$. Light, weakly coupled particles may therefore apparently produce a mass gap $2 {\rm M}_\odot < M_{\rm BH} < 4 {\rm M}_\odot $ in the LIGO/Virgo/KAGRA data, when no mass gap is present in the stellar remnant population.

We present the first measurements of HI galaxy scaling relations from a blind survey at $z>0.15$. We perform spectral stacking of 9023 spectra of star-forming galaxies undetected in HI at $0.23<z<0.49$, extracted from MIGHTEE-HI Early Science datacubes, acquired with the MeerKAT radio telescope. We stack galaxies in bins of galaxy properties ($M_*$, SFR, and sSFR, with ${\rm sSFR}\equiv M_*/{\rm SFR}$), obtaining $\gtrsim 5\sigma$ detections in most cases, the strongest HI-stacking detections to date in this redshift range. With these detections, we are able to measure scaling relations in the probed redshift interval, finding evidence for a moderate evolution from the median redshift of our sample $z_{\rm med}\sim 0.37$ to $z\sim 0$. In particular, low-$M_*$ galaxies ($\log_{10}(M_*/{\rm M_\odot})\sim 9$) experience a strong HI depletion ($\sim 0.5$ dex in $\log_{10}(M_{\rm HI}/{\rm M}_\odot)$), while massive galaxies ($\log_{10}(M_*/{\rm M_\odot})\sim 11$) keep their HI mass nearly unchanged. When looking at the star formation activity, highly star-forming galaxies evolve significantly in $M_{\rm HI}$ ($f_{\rm HI}$, where $f_{\rm HI}\equiv M_{\rm}/M_*$) at fixed SFR (sSFR), while at the lowest probed SFR (sSFR) the scaling relations show no evolution. These findings suggest a scenario in which low-$M_*$ galaxies have experienced a strong HI depletion during the last $\sim4$ Gyr, while massive galaxies have undergone a significant HI replenishment through some accretion mechanism, possibly minor mergers. Interestingly, our results are in good agreement with the predictions of the SIMBA simulation. We conclude that this work sets novel important observational constraints on galaxy scaling relations.

A science goal of many future X-ray observatories is mapping the cosmic web through deep exposures of faint diffuse sources. Such observations require low background and the best possible knowledge of the remaining unrejected background. The dominant contribution to the background above 1-2 keV is from Galactic Cosmic Ray protons. Their flux and spectrum are modulated by the solar cycle but also by solar activity on shorter timescales. Understanding this variability may prove crucial to reducing background uncertainty for ESA's Athena X-ray Observatory and other missions with large collecting area. We examine of the variability of the particle background as measured by ACIS on the Chandra X-ray Observatory and compare that variability to that measured by the Alpha Magnetic Spectrometer (AMS), a precision particle detector on the ISS. We show that cosmic ray proton variability measured by AMS is well matched to the ACIS background and can be used to estimate proton energies responsible for the background. We discuss how this can inform future missions.

The 2020 Decadal Survey on Astronomy and Astrophysics endorsed space-based high contrast imaging for the detection and characterization of habitable exoplanets as a key priority for the upcoming decade. To advance the maturity of starlight suppression techniques in a space-like environment, we are developing the Space Coronagraph Optical Bench (SCoOB) at the University of Arizona, a new thermal vacuum (TVAC) testbed based on the Coronagraphic Debris Exoplanet Exploring Payload (CDEEP), a SmallSat mission concept for high contrast imaging of circumstellar disks in scattered light. When completed, the testbed will combine a vector vortex coronagraph (VVC) with a Kilo-C microelectromechanical systems (MEMS) deformable mirror from Boston Micromachines Corp (BMC) and a self-coherent camera (SCC) with a goal of raw contrast surpassing $10^{-8}$ at visible wavelengths. In this proceedings, we report on our wavefront sensing and control efforts on this testbed in air, including the as-built performance of the optical system and the implementation of algorithms for focal-plane wavefront control and digging dark holes (regions of high contrast in the focal plane) using electric field conjugation (EFC) and related algorithms.

The development of spaceborne coronagraphic technology is of paramount importance to the detection of habitable exoplanets in visible light. In space, coronagraphs are able to bypass the limitations imposed by the atmosphere to reach deeper contrasts and detect faint companions close to their host star. To effectively test this technology in a flight-like environment, a high-contrast imaging testbed must be designed for operation in a thermal vacuum (TVAC) chamber. A TVAC-compatible high-contrast imaging testbed is undergoing development at the University of Arizona inspired by a previous mission concept: The Coronagraphic Debris and Exoplanet Exploring Payload (CDEEP). The testbed currently operates at visible wavelengths and features a Boston Micromachines Kilo-C DM for wavefront control. Both a vector vortex coronagraph and a knife-edge Lyot coronagraph operating mode are under test. The optics will be mounted to a 1 x 2 meter pneumatically isolated optical bench designed to operate at 10^-8 torr and achieve raw contrasts of 10^-8 or better. The validation of our optical surface quality, alignment procedure, and first light results are presented. We also report on the status of the testbed's integration in the vaccum chamber.

In asteroseismology, the surface effect is a disparity between the observed and the modelled oscillation frequencies. It originates from improper modelling of the surface layers in stars with solar-like oscillations. Correcting the surface effect usually requires using functions with free parameters, which are conventionally fitted to the observed frequencies. On the basis that the correction should vary smoothly across the H--R diagram, we parameterize it as a simple function of three stellar surface properties: surface gravity, effective temperature, and metallicity. We determine this function by fitting stars ranging from main-sequence dwarfs to red-giant-branch stars. The absolute amount of the surface correction increases with surface gravity, but the ratio between it and $\nu_{\rm max}$ decreases. Applying the prescription has an advantage of eliminating unrealistic surface correction, which improves parameter estimations with stellar modelling. Using two open clusters, we found that adopting the prescription can help reduce the scatter of the model-derived ages for each star in the same cluster. For an application, we provide a new revision for the $\Delta\nu$ scaling relation, using our prescription to account for the surface effect in models. The values of the correction factor, $f_{\Delta\nu}$, are up to 2\% smaller than those determined without the surface effect considered, suggesting decreases of up to 4\% in asteroseismic scaling radii and up to 8\% in asteroseismic scaling masses. This revision brings the asteroseismic properties into agreement with those determined from eclipsing binaries. Finally, the new correction factor and the stellar models with the corrected frequencies are made publicly available.

We investigated a wealth of X-ray and gamma-ray spectral energy distribution (SED) and multi-band light curve (LC) data of the gamma-ray binary HESS J0632+057 using a phenomenological intrabinary shock (IBS) model. Our baseline model assumes that the IBS is formed by colliding winds from a putative pulsar and its Be companion, and particles accelerated in the IBS emit broadband radiation via synchrotron (SY) and inverse-Compton upscattering (ICS) processes. Adopting the latest orbital solution and system geometry (Tokayer et al. 2021), we reproduced the global X-ray and TeV LC features, two broad bumps at $\phi \sim 0.3$ and $\sim0.7$, with the SY and ICS model components. We found these TeV LC peaks originate from ICS emission caused by the enhanced seed photon density near periastron and superior conjunction or Doppler-beamed emission of bulk-accelerated particles in the IBS at inferior conjunction. While our IBS model successfully explained most of the observed SED and LC data, we found that phase-resolved SED data in the TeV band require an additional component associated with ICS emission from pre-shock particles (produced by the pulsar wind). This finding indicates a possibility of delineating the IBS emission components and determining the bulk Lorentz factors of the pulsar wind at certain orbital phases.

Ultra-fast outflows (UFOs) with mildly relativistic velocities are measured using the X-ray spectra of radio-quiet and -loud active galactic nuclei (AGNs). In general, UFOs are believed to be generated from the accretion disk around a black hole (BH). A line-force driving model is suggested to be the mechanism to drive UFOs from the accretion disk. In this paper, we use the non-hydrodynamic approach to examine the influences of radiation-drag effects on the line-force-driven winds generated from the accretion disk. We find that the radiation-drag effects can significantly weaken the line-force-driven winds. Compared with the case without the radiation-drag effects, when the radiation-drag effects are considered, the maximum speed of winds is reduced by $\sim$60\%--70\%, the mass outflow rate is reduced by $\sim$50\%--80\%, and the kinetic power is reduced by about an order of magnitude. The radiation-drag effects narrow the area where the winds are generated.

One critical remaining issue to unclear the initial mass function of the first (Population III) stars is the final fate of secondary protostars formed in the accretion disk, specifically whether they merge or survive. We focus on the magnetic effects on the first star formation under the cosmological magnetic field. We perform a suite of ideal magnetohydrodynamic simulations until 1000 years after the first protostar formation. Instead of the sink particle technique, we employ a stiff equation of state approach to represent the magnetic field structure connecting to protostars. Ten years after the first protostar formation in the cloud initialized with $B_0 = 10^{-20}$ G at $n_0 = 10^4\,{\rm cm^{-3}}$, the magnetic field strength around protostars amplifies from pico- to kilo-gauss, which is the same strength as the present-day star. The magnetic field rapidly winds up since the gas in the vicinity of the protostar ($\leq\!10$ au) has undergone several tens orbital rotations in the first decade after protostar formation. As the mass accretion progresses, the vital magnetic field region extends outward, and the magnetic braking eliminates fragmentation of the disk that would form in the unmagnetized model. On the other hand, assuming a gas cloud with small angular momentum, this amplification might not work because the rotation would be slower. However, disk fragmentation would not occur in that case. We conclude that the exponential amplification of the cosmological magnetic field strength, about $10^{-18}$ G, eliminates disk fragmentation around the Population III protostars.

Instrumenting a gigaton of ice at the geographic South Pole, the IceCube Neutrino Observatory has been at the forefront of groundbreaking scientific discoveries over the past decade. These include the observation of a flux of TeV-PeV astrophysical neutrinos, detection of the first astrophysical neutrino on the Glashow resonance and evidence of the blazar TXS 0506+056 as the first known astronomical source of high-energy neutrinos. Several questions, however, remain, pertaining to the precise origins of astrophysical neutrinos, their production mechanisms at the source and in Earth's atmosphere and in the context of physics beyond the Standard Model. This proceeding highlights some of our latest results, from new constraints on neutrino interactions and oscillations to the latest measurements of the astrophysical neutrino flux and searches for their origins to future prospects with IceCube-Gen2.

Current, state-of-the-art CCDs are close to being able to deliver all key performance figures for future strategic X-ray missions except for the required frame rates. Our Stanford group is seeking to close this technology gap through a multi-pronged approach of microelectronics, signal processing and novel detector devices, developed in collaboration with the Massachusetts Institute of Technology (MIT) and MIT Lincoln Laboratory (MIT-LL). Here we report results from our (integrated) readout electronics development, digital signal processing and novel SiSeRO (Single electron Sensitive Read Out) device characterization.

An X-ray flash, expected in a very early phase of a nova outburst, was at last detected with the {\it SRG}/eROSITA in the classical nova YZ Reticuli 2020. The observed flash timescale, luminosity, and blackbody temperature substantially constrain the nova model. We present light curve models of the X-ray flash for various white dwarf (WD) masses and mass accretion rates. We have found the WD mass in YZ Ret to be as massive as $M_{\rm WD}\sim 1.3 ~M_\odot$ with mass accretion rates of $\dot M_{\rm acc}\sim 5 \times 10^{-10}- 5\times 10^{-9} ~M_\odot$ yr$^{-1}$ including the case that the mass accretion rate is changing between them, to be consistent with the {\it SRG}/eROSITA observation. The X-ray observation confirms the luminosity to be close to the Eddington limit at the X-ray flash. The occurrence of optically thick winds, with the photospheric radius exceeding $\sim 0.1~R_\odot$, terminated the X-ray flash of YZ Ret by strong absorption. This sets a constrain on the starting time of wind mass loss. A slight contamination of the core material into the hydrogen rich envelope seems to be preferred to explain a very short duration of the X-ray flash.

Binaries consisting of a hot subdwarf star and an accreting white dwarf (WD) are sources of gravitational wave radiation at low frequencies and possible progenitors of type Ia supernovae if the WD mass is large enough. Here, we report the discovery of the third binary known of this kind: it consists of a hot subdwarf O (sdO) star and a WD with an orbital period of 3.495 hours and an orbital shrinkage of 0.1 s in 6 yr. The sdO star overfills its Roche lobe and likely transfers mass to the WD via an accretion disk. From spectroscopy, we obtain an effective temperature of $T_{\mathrm{eff}}=54\,240\pm1\,840$ K and a surface gravity of $\log{g}=4.841\pm0.108$ for the sdO star. From the light curve analysis, we obtain a sdO mass of $M_{\mathrm{sdO}}=0.55$ ${\mathrm{M_{\odot}}}$ and a mass ratio of $q=M_{\mathrm{WD}}/M_{\mathrm{sdO}}=0.738\pm0.001$. Also, we estimate that the disk has a radius of $\sim 0.41R_\odot$ and a thickness of $\sim 0.18R_\odot$. The origin of this binary is probably a common envelope ejection channel, where the progenitor of the sdO star is either an RGB star or, more likely, an early AGB star; the sdO star will subsequently evolve into a WD and merge with its WD companion, likely resulting in an R CrB star. The outstanding feature in the spectrum of this object is strong Ca H&K lines, which are blueshifted by $\sim$200 km/s and likely originate from the recently ejected common envelope, and we estimated that the remnant CE material in the binary system has a density $\sim 6\times 10^{-10} {\rm g/cm^3}$.

Elemental abundances of extrinsic carbon stars provide insight into the poorly understood origin and evolution of elements in the early Galaxy. In this work, we present the results of a detailed spectroscopic analysis of four potential carbon star candidates from the Hamburg/ESO Survey (HES) HE~0457$-$1805, HE~0920$-$0506, HE~1241$-$0337, and HE~1327$-$2116. This analysis is based on the high-resolution spectra obtained with Mercator/HERMES (R$\sim$86,000) and SUBARU/HDS (R$\sim$50,000). Although the abundances of a few elements, such as, Fe, C, and O are available from medium-resolution spectra, we present the first ever detailed high-resolution spectroscopic analysis for these objects. The object HE~0457$-$1805 and HE~1241$-$0337 are found to be CEMP-s stars, HE~0920$-$0506 a CH star, and HE~1327$-$2116 a CEMP-r/s star. The object HE~0457$-$1805 is a confirmed binary, whereas the binary status of the other objects is unknown. The locations of program stars on the absolute carbon abundance A(C) vs [Fe/H] diagram point at their binary nature. We have examined various elemental abundance ratios of the program stars and confirmed the low-mass nature of their former AGB companions. We have shown that the i-process models could successfully reproduce the observed abundance pattern in HE~1327$-$2116. The parametric model based analysis performed for HE~0457$-$1805, HE~0920$-$0506, and HE~1241$-$0337 based on the FRUITY models confirmed that the surface chemical composition of these three objects are influenced by pollution from low-mass AGB companions.

Dense gas is important for galaxy evolution and star formation. Optically-thin dense-gas tracers, such as isotopologues of HCN, HCO+, etc., are very helpful to diagnose excitation conditions of dense molecular gas. However, previous studies of optically-thin dense-gas tracers were mostly focusing on average properties of galaxies as a whole, due to limited sensitivity and angular resolution. M82, a nearby prototype starburst galaxy, offers a unique case for spatially-resolved studies with single-dish telescopes. With the IRAM 30-m telescope, we observed the J = 1 - 0 transition of H13CN, HC15N, H13CO+, HN13C, H15NC, and SiO J = 2 - 1, HC3N J= 10 - 9, H2CO J = 2 - 1 toward five positions along the major axis of M82. The intensity ratios of I(HCN)/I(H13CN) and I(HCO+)/I(H13CO+) show a significant spatial variation along the major axis, with lower values in the central region than those on the disk, indicating higher optical depths in the central region. The optical depths of HCO+ lines are found to be systematically higher than those of HCN lines at all positions. Futhermore, we find that the 14N/15N ratios have an increasing gradient from the center to the outer disk.

We investigate the star-forming main sequence of the host galaxies of a large, well-defined sample of 453 redshift $\sim$0.3 quasars with previously available star formation rates by deriving stellar masses from modeling their broad-band ($grizy$) spectral energy distribution. We perform two-dimensional, simultaneous, multi-filter decomposition of Pan-STARRS1 3$\pi$ Steradian Survey images to disentangle the active galactic nucleus (AGN) from its host galaxy, by explicitly considering, for the first time, the wavelength variation of galaxy structures. We quantify the S\'ersic profiles and sizes of the host galaxies from mock AGNs generated from both real and idealized galaxies. Detailed morphological classifications of the calibration galaxy sample with Hubble Space Telescope images enable us to estimate crude morphological types of the quasars. Although the majority ($\sim$60%) of the quasars are hosted by bulge-dominated, early-type galaxies, a substantial fraction ($\sim$40%) reside in disk-dominated, late-type galaxies, suggesting that at least in these systems major mergers have not played a significant role in regulating their AGN activity, in agreement with recent simulations and observations of nearby quasars. The vast majority ($\sim$90%) of the quasars have star formation rates that place them on or above the galaxy star-forming main sequence, with more rapidly accreting AGNs displaced further above the main sequence. Quasar host galaxies generally follow the stellar mass-size relation defined by inactive galaxies, both for late-type and early-type systems, but roughly 1/3 of the population has smaller sizes at a given stellar mass, reminiscent of compact star-forming galaxies at higher redshift.

Recent laboratory efforts and telescopic observations of Europa have shown the relevance of a yellow colouration of sodium chloride (NaCl) caused by crystal defects generated by irradiation. We further investigate this process by irradiating (with energetic electrons) different types of analogues where NaCl is associated in different ways to water ice. We produce two types of icy analogues: compact slabs and granular particles where we investigate two particle sizes (5 and 70 $\mu$m). We perform electron irradiation at cryogenic temperatures (100 K) and under high vacuum (10-7 mbar) conditions, with energies of 1 and 5 keV. We observe the formation of two different types of colour centres. The so-called F-centres (460 nm) were formed in every sample, but the intensity of the absorption band within the compact slabs surpassed any other icy analogues and was comparable to the intensity of the absorption band within pure NaCl grains. M-centres (720 nm) have not been detected at the surface of Europa so far, and were close to the detection limit during our irradiation of compact slabs. The slabs could be good analogues for Europa_s surface as they produce mainly F-centres. Other notable differences have been observed between compact slabs and granular samples, such as the presence of an absorption band at 580 nm attributed to colloids of Na, exclusively within granular samples. Such absorptions have not been reported in previous studies.

Time-domain astronomy is entering a new era as wide-field surveys with higher cadences allow for more discoveries than ever before. The field has seen an increased use of machine learning and deep learning for automated classification of transients into established taxonomies. Training such classifiers requires a large enough and representative training set, which is not guaranteed for new future surveys such as the Vera Rubin Observatory, especially at the beginning of operations. We present the use of Gaussian processes to create a uniform representation of supernova light curves from multiple surveys, obtained through the Open Supernova Catalog for supervised classification with convolutional neural networks. We also investigate the use of transfer learning to classify light curves from the Photometric LSST Astronomical Time Series Classification Challenge (PLAsTiCC) dataset. Using convolutional neural networks to classify the Gaussian process generated representation of supernova light curves from multiple surveys, we achieve an AUC score of 0.859 for classification into Type Ia, Ibc, and II. We find that transfer learning improves the classification accuracy for the most under-represented classes by up to 18% when classifying PLAsTiCC light curves, and is able to achieve an AUC score of 0.945 when including photometric redshifts for classification into six classes (Ia, Iax, Ia-91bg, Ibc, II, SLSN-I). We also investigate the usefulness of transfer learning when there is a limited labelled training set to see how this approach can be used for training classifiers in future surveys at the beginning of operations.

Laser Guide Star Adaptive Optics (LGS-AO) is becoming routine in several astronomical observatories. The use of powerful lasers generates sensible Raman emissions on the uplink laser beam path, plus secondary Rayleigh scattering from atmospheric molecules and Mie scattering from aerosols. This paper reports the results of a campaign done with the 10.4m Gran Telescopio CANARIAS (GTC); this campaign was undertaken to assess the spectral and photometric contamination coming from a 589 nm laser uplink beam scattering and Raman emission induced on the GTC spectro-imager OSIRIS by laser launched about 1 km off-axis. The photometric contamination is due to primary and secondary scattering of the uplink photons, as well by the Raman inelastic scattering. We have propagated the laser beam creating a mesospheric LGS, then pointed and focused the GTC telescope toward the uplink laser beam, at different heights and up to the LGS, taking into account the observing geometry. In our observations, the Raman emissions for O2 and N2 vibrational lines are visible at 20 km, weakening with altitude and becoming undetectable above 30 km. The scattering of the focused uplink beam is detectable at less than +/-0.2 arcmin from the center of the beam, while for the focused LGS the scattering is narrower, being detectable at less than +/-0.1 arcmin around the plume. Recommendations for Laser Traffic Control Systems (LTCS) are given accordingly.

Faraday rotation studies of distant radio sources can constrain the evolution and the origin of cosmic magnetism. We use data from the LOFAR Two Metre Sky Survey: Data Release 2 (LoTSS DR2) to study the dependence of the Faraday rotation measure (RM) on redshift. By focusing on radio sources that are close in terms of their projection on the sky, but physically unrelated (random pairs), we measure the RM difference, $\Delta$RM, between the two sources. Thus, we isolate the extragalactic contribution to $\Delta$RM from other contributions. We present a statistical analysis of the resulting sample of random pairs and find a median absolute RM difference |$\Delta$RM| $ = (1.79 \pm 0.09)$ rad/m$^{2}$ , with |$\Delta$RM| uncorrelated both with respect to the redshift difference of the pair and the redshift of the nearer source, and a median excess of random pairs over physical pairs of $(1.65 \pm 0.10)$ rad/m$^{2}$. We seek to reproduce this result with Monte Carlo simulations assuming a non vanishing seed cosmological magnetic field and a redshift evolution of the comoving magnetic field strength that varies as $1/(1 + z)^{\gamma}$. We find the best fitting results $B_0 \equiv B_{\rm comoving}(z = 0) \lesssim (2.0 \pm 0.2)$ nG and $\gamma \lesssim 4.5 \pm 0.2$ that we conservatively quote as upper limits due to an unmodelled but non vanishing contribution of local environments to the RM difference. A comparison with cosmological simulations shows our results to be incompatible with primordial magnetogenesis scenarios with uniform seed fields of order nG.

The AMS-02 collaboration recently released cosmic-ray F/Si data of unprecedented accuracy. CR F is predominantly produced by fragmentation of heavier progenitors, while Si is mostly accelerated at source. This ratio is thus maximally sensitive to CR propagation. We study the compatibility of the transport parameters derived from the F/Si ratio with those obtained from the lighter Li/C, Be/C, and B/C ratios. We also inspect the CR source abundance of F, one of the few elements with a high first ionisation potential but only moderately volatile, and a potentially key element to study the acceleration mechanism of CRs. We use the 1D diffusion model implemented in the USINE code and perform $\chi^2$ analyses accounting for several systematic effects (energy correlations in data, nuclear cross sections and solar modulation uncertainties). We also take advantage of the EXFOR nuclear database to update the F production cross sections for its most important progenitors (identified to be 56Fe, 32S, 28Si, 27Al, 24Mg, 22Ne, and 20Ne). The transport parameters obtained from AMS-02 F/Si data are compatible with those obtained from AMS-02 (Li,Be,B)/C data. The combined fit of all these ratios leads to a $\chi^2/dof\approx1.1$, with $\lesssim 10\%$ adjustments of the B and F production cross sections (the latter are based on very few nuclear data points, and would strongly benefit from new measurements). The F/Si ratio is compatible with a pure secondary origin of F, with a best-fit relative source abundance 19F/28Si$_{CRS}\sim 10^{-3}$ and an upper limit of $\sim 5\times 10^{-3}$. Unfortunately, this limit is not sufficient to test global acceleration models of CR nuclei, for which values at the level of $\sim 10^{-4}$ are required. Such levels could be attained with F/Si data of a few percent accuracy at a few tens of TV, possibly within reach of the next generation of CR experiments.

We study the effect of two Modified Gravity (MG) theories, $f(R)$ and nDGP, on three probes of large-scale structure, the real space power spectrum estimator $Q_0$, bispectrum and voids, and validate fast approximate COLA simulations against full $N$-body simulations for the prediction of these probes. We find that using the first three even multipoles of the redshift space power spectrum to estimate $Q_0$ is enough to reproduce the MG boost factors of the real space power spectrum for both halo and galaxy catalogues. By analysing the bispectrum and reduced bispectrum of Dark Matter (DM), we show that the strong MG signal present in the DM bispectrum is mainly due to the enhanced power spectrum. We warn about adopting screening approximations in simulations as this neglects non-linear contributions that can source a significant component of the MG bispectrum signal at the DM level, but we argue that this is not a problem for the bispectrum of galaxies in redshift space where the signal is dominated by the non-linear galaxy bias. Finally, we perform void-finding on our galaxy mock catalogues by the ZOBOV watershed algorithm. To apply a linear model for Redshift-Space Distortion (RSD) in the void-galaxy cross-correlation function, we first examine the effects of MG on the void profiles entering into the RSD model. We find relevant MG signals in the integrated-density, velocity dispersion and radial velocity profiles in the nDGP theory. Fitting the RSD model for the linear growth rate, we recover the linear theory prediction in an nDGP model, which is larger than the $\Lambda$CDM prediction at the $3 \sigma$ level. In $f(R)$ theory we cannot naively compare the results of the fit with the linear theory prediction as this is scale-dependent, but we obtain results that are consistent with the $\Lambda$CDM prediction.

Recent analyses have found intriguing correlations between the colour ($c$) of type Ia supernovae (SNe Ia) and the size of their mass-step, the relationship between host galaxy stellar mass and Hubble residual. These analyses suggest that the underlying cause of this relationship is dust. Using a sample of 675 photometrically-classified SNe Ia from the Dark Energy Survey 5-year sample, we study the differences in Hubble residual for a variety of host and local properties for subsamples split by their colour ($c$). We find a $3\sigma$ difference for the size of the mass-step when comparing blue ($c < 0$) and red ($c > 0$) SNe. We observe the lowest r.m.s. scatter ($\sim 0.14$) in Hubble residual for blue SNe in low mass or blue environments, suggesting that these objects provide the most homogeneous sample for cosmological analyses. By fitting for $c$-dependent relationships between Hubble residuals and $M_\mathrm{stellar}$, approximating existing dust models, we remove the mass-step from the data but find significant remaining steps in rest-frame $U-R$, indicating that current dust modelling based on $M_\mathrm{stellar}$ may not fully explain the remaining dispersion in SN luminosity. The most dispersion is removed by instead accounting for a $c$-dependent relationship between Hubble residuals and global $U-R$, resulting in $\leq 1\sigma$ remaining steps in other environmental properties, suggesting that $U-R$ provides different information about the environment of SNe Ia to $M_\mathrm{stellar}$. This $c$-dependent $U-R$ relation implies that $U-R$ may be more closely linked to dust, motivating the future inclusion of galaxy $U-R$ colour in the correction for SN distance biases.

We present the results of a broadband (0.5-78 keV) X-ray spectral study of the persistent Galactic black hole X-ray binary GRS 1758-258 observed simultaneously by Swift and NuSTAR. Fitting with an absorbed power-law model revealed a broad Fe line and reflection hump in the spectrum. We used different flavours of the relativistic reflection model for the spectral analysis. All models indicate the spin of the black hole in GRS 1758-258 is >0.92. The source was in the low hard state during the observation, with the hot electron temperature of the corona estimated to be kT$_e$ ~ 140 keV. The black hole is found to be accreting at ~1.5 % of the Eddington limit during the observation, assuming the black hole mass of 10 $M_{\odot}$ and distance of 8 kpc.

X-ray studies of jellyfish galaxies opened a window in the physics of the interplay between intracluster medium (ICM) and interstellar medium (ISM). In this paper, we present the study of an archival \textit{Chandra} observation of the GASP jellyfish galaxy JO194. We observe X-ray emission extending from the stellar disk to the unwinding spiral arms with an average temperature of $kT=0.79\pm0.03$ keV. To investigate the origin of the X-ray emission, we compare the observed X-ray luminosities with those expected from the star formation rates (SFR) obtained from H$\alpha$ emission. We estimate an X-ray luminosity excess of a factor $\sim2-4$ with respect to the SF, therefore we conclude that SF is not the main responsible for the extended X-ray emission of JO194. The metallicity in the spiral arms ($Z=0.24^{+0.19}_{-0.12} Z_{\odot}$) is consistent with that of the ICM around JO194 ($Z=0.35\pm0.07$), thus we suggest that the ICM radiative cooling dominates the X-ray emission of the arms. We speculate that the X-ray plasma results from the ISM-ICM interplay, although the nature of this interplay is still mostly unknown. Finally, we observe that the X-ray properties of JO194 are consistent with those of two other GASP galaxies with different stellar mass, phase-space conditions in their hosting clusters, and local ICM conditions. We suggest that the conditions required to induce extended X-ray emission in jellyfish galaxies are established at the beginning of the stripping, and they can persist on long time scales so that galaxies in different clusters and evolutionary stages can present similar extended X-ray emission.

The processes regulating protoplanetary disk evolution are constrained by studying how mass accretion rates scale with stellar and disk properties. The spread in these relations can be used as a constraint to the models of disk evolution, but only if the impact of accretion variability is correctly accounted for. While the effect of variability might be substantial in the embedded phases of star formation, it is often considered limited at later stages. Here we report on the observed large variation in the accretion rate for one target, XX Cha, and we discuss the impact on population studies of classical T Tauri stars. The mass accretion rate determined by fitting the UV-to-near-infrared spectrum in recent X-Shooter observations is compared with the one measured with the same instrument 11 years before. XX Cha displays an accretion variability of almost 2 dex between 2010 and 2021. Although the timescales on which this variability happens are uncertain, XX Cha displays an extreme accretion variability for a classical T Tauri star. If such behavior is common among classical T Tauri stars, possibly on longer timescales than previously probed, it could be relevant for discussing the disk evolution models constrained by the observed spread in accretion rates. Finally, we remark that previous studies of accretion variability based on spectral lines may have underestimated the variability of some targets.

We present the results from a coordinated XMM-Newton $+$ NuSTAR observation of the Seyfert 1 Galaxy ESO 511$-$G030. With this joint monitoring programme, we conduct a detailed variability and spectral analysis. The source remained in a low flux and very stable state throughout the observation period, although there are slight fluctuations of flux over long timescales. The broadband (0.3-78~keV) spectrum shows the presence of a power-law continuum with a soft excess below 2~keV, a relatively narrow iron K$\alpha$ emission ($\sim$6.4~keV), and an obvious cutoff at high energies. We find that the soft excess can be modeled by two different possible scenarios: a warm ($kT_{\rm e} \sim$ 0.19~keV) and optically thick ($\tau - 18\sim25$) Comptonizing corona or a relativistic reflection from a high-density ($\log [n_{\rm e}/{\rm cm}^{-3}]=17.1 \sim 18.5$) inner disc. All models require a low temperature ($kT_{\rm e} \sim$ 13~keV) for the hot corona.

We present the discovery of the simultaneous flux variation of a massive young stellar object (MYSO) G036.70+00.09 (G036.70) both in the maser emission and mid-infrared (MIR;$\lambda=3$--$5$~$\mu$m) bands. Utilizing the ALLWISE and NEOWISE archival databases covering a long time span of approximately 10 years with a cadence of 6 months, we confirmed that G036.70 indicates a stochastic year-long MIR variability with no signs of the WISE band color change of W1 (3.4~$\mu$m) $-$W2 (4.6~$\mu$m). Cross-matching the MIR data set with the high-cadence 6.7~GHz class II methanol maser flux using a Hitachi 32-m radio telescope that discovered its periodicity in the methanol maser of 53.0--53.2 days, we also determine the flux correlations between the two bands at two different timescales, year-long and day-long, both of which have never been reported in MYSOs except when they are in a state of the accretion burst phase. The results of our study support the scenario that a class II methanol maser is pumped up by infrared emission from accreting disks of MYSOs. We also discuss the possible origins of MIR and maser variability. To explain the two observed phenomena, a stochastic year-long MIR variability with no signs of significant color change and maser-MIR variability correlation, change in mass accretion rate and line-of-sight extinction because of nonaxisymmetric dust density distribution in a rotating accretion disk are possible origins. Observations through spectroscopic monitoring of accretion-related emission lines are essential for determining the origin of the observed variability in G036.70.

Roughly 10% of quasars are "radio-loud", producing copious radio emission in large jets. The origin of the low-level radio emission seen from the remaining 90% of quasars is unclear. Observing a sample of eight radio-quiet quasars with the Very Long Baseline Array, we discovered that their radio properties depend strongly on their Eddington ratio (r_Edd=L_AGN/L_Edd). At lower Eddington ratios (r_Edd < 0.3), the total radio emission of the AGN predominately originates from an extremely compact region, possibly as small as the accretion disk. At higher Eddington ratios (r_Edd > 0.3), the relative contribution of this compact region decreases significantly, and though the total radio power remains about the same, the emission now originates from regions >100 pc large. The change in the physical origin of the radio-emitting plasma region with r_Edd is unexpected, as the properties of radio-loud quasars show no dependence with Eddington ratio. Our results suggest that at lower Eddington ratios the magnetised plasma is likely confined by the accretion disk corona, and only at higher Eddington ratios escapes to larger scales. Stellar-mass black holes show a similar dependence of their radio properties on the accretion rate, supporting the paradigm which unifies the accretion onto black holes across the mass range.

The shape of the faint-end of the high-z galaxy luminosity function (LF) informs early star formation and reionization physics during the Cosmic Dawn and Epoch of Reionization. Until recently, based on the strong gravitational lensing cluster deep surveys, the Hubble Frontier Fields (HFF) has found a potential turnover in the ultraviolet (UV) LF at z$\sim$6. In this paper, we analyze the contribution of extremely faint galaxies with the magnitude larger than the turnover magnitude in LF to cosmic reionization. We apply the measurement from HFF to our suppressed star formation efficiency model, including three free parameters: halo mass threshold $M_t$, curvature parameter $\beta$ and a UV conversion factor $l_{\rm UV}$. According to our fit of 68\% confidence level, the high-redshift star formation in haloes smaller than $ M_t=1.82^{+2.86}_{-1.08}\times10^{10} \rm M_{\odot}$ is found to be dampened. The turnover magnitude $\rm \gtrsim -13.99-2.45$, correspondingly the halo mass $\lesssim(4.57+20.03)\times10^{9} \rm M_{\odot}$. We find that the absorption trough in the global 21-cm signal is sensitive to our SFE model parameters. Together with ($\beta$, $l_{\rm UV}$) = ($2.17^{+2.42}_{-1.72}$, $9.33^{+0.43}_{-0.42} \rm ~erg~yr ~s^{-1}M_{\odot}^{-1})$, the trough locates at $\sim$ $134^{+10}_{-17}$ $\rm MHz$ with an amplitude of $\sim$ $-237^{-6}_{+7}$ $\rm mK$, compared to (106\rm MHz, -212\rm mK) in the absence of turnover. Besides, we find that the star formation of faint galaxies has also an impact on the 21-cm power spectra. The best fitting peak power decreases by $\sim4\%$ and shifts towards smaller scales from $0.88 h \rm Mpc^{-1}$ to $0.91 h \rm Mpc^{-1}$. According to our calculation, such impact is distinguishable with the forthcoming Square Kilometre Array.

Gamma-ray bursts are flashes of light from distant exploding stars. Cube satellites that monitor photons across different energy bands are used to detect these bursts. There is a need for computationally efficient algorithms, able to run using the limited computational resource onboard a cube satellite, that can detect when gamma-ray bursts occur. Current algorithms are based on monitoring photon counts across a grid of different sizes of time window. We propose a new algorithm, which extends the recently developed FOCuS algorithm for online change detection to Poisson data. Our algorithm is mathematically equivalent to searching over all possible window sizes, but at half the computational cost of the current grid-based methods. We demonstrate the additional power of our approach using simulations and data drawn from the Fermi gamma-ray burst catalogue.

Topological defects are a fossil relic of early Universe phase transitions, with cosmic strings being the best motivated example. While in most cases one studies Nambu-Goto or Abelian-Higgs strings, one also expects that cosmologically realistic strings should have additional degrees of freedom in their worldsheets, one specific example being superstrings from Type IIB superstring theory. Here we continue the scientific exploitation of our recently developed multi-GPU field theory cosmic strings code to study the evolution of U(1)$\times$U(1) multitension networks, which are a numerically convenient proxy: these contain two lowest-tension strings networks able to interact and form bound states, providing a convenient first approximation to the behaviour expected from cosmic superstrings. (...) We rely on the largest field theory simulations of this model so far, specifically $4096^3$, $\Delta x = 0.5$ boxes. We present robust evidence of scaling for the lightest strings, measured through a complete and self-consistent set of correlation length and velocity diagnostics. We also find a linearly growing average length of the bound state segments, consistent with a scaling behaviour. (In previously reported lower-resolution simulations, such behaviour had only been identified with carefully engineered initial conditions, rich in those segments.) Finally, while we see no evidence of a large population of bound states forming at early stages of the network evolution, we do present tentative evidence for an asymptotic constant value of the fraction of bound states, with this value being different in the radiation and the matter eras. Our work demonstrates that our GPU-accelerated field theory code can by successfully extended beyond the simple Abelian-Higgs approximation, and enables future detailed studies of realistic string networks and of their observational signatures.

The ESA/KU Leuven CubeSpec mission is specifically designed to provide low-cost space-based high-resolution optical spectroscopy. Here we highlight the science requirements and capabilities of CubeSpec. The primary science goal is to perform pulsation mode identification from spectroscopic line profile variability and empower asteroseismology of massive stars.

Given an increasing number of gamma-ray bursts accompanied by potential kilonovae there is a growing importance to advance modelling of kilonova afterglows. This electromagnetic signature might play an important role in the multi-messenger picture of binary neutron star mergers and offer a possibility to infer the ejecta properties and, consequently, the binary properties. In this work, we investigate how the presence of two electron populations that follow a Maxwellian (thermal) and a power-law (non-thermal) distributions affect kilonova afterglow light curves. The modelling is done using the semi-analytic afterglow model, $\texttt{PyBlastAfterglow}$. We consider a set of kilonova ejecta profiles from ab-initio numerical relativity binary neutron star merger simulations, targeted to GW170817. We find that the emission from thermal electrons dominates at early times. If the interstellar medium density is sufficiently high (${\simeq}0.1\,$cm$^{-3}$) it adds an early time peak to the light curve. As ejecta decelerates the spectral and temporal indexes change in a characteristic way that, if observed, can be used to reconstruct the ejecta velocity distribution. For the low interstellar medium density, inferred for GRB170817A, the emission from the thermal electron population is dimmer than that from the non-thermal. We also assess how kilonova afterglow light curves change if the interstellar medium has been altered by laterally expanding gamma-ray burst ejecta. For the latter we consider properties informed by observations of GRB170817A. We find that the main effect is the emission suppression at early time ${\lesssim}10^{3}\,$days, and at its maximum it reaches ${\sim}40\%$. The subsequent rebrightening, when these ejecta break through and shocks form, is very mild (${\lesssim}10\%$), smeared over time and may not be observable.

Whether it be due to rapid rotation or binary interactions, deviations from spherical symmetry are common in massive stars. These deviations from spherical symmetry are known to cause non-uniform distributions of various parameters across the surface including temperature, which can drive internal mixing processes within the envelopes of these massive stars. Despite how common these 3D distortions are, they are often neglected in spectroscopic analyses. We present a new spectral analysis code called SPAMMS (Spectroscopic PAtch Model for Massive Stars) specifically designed to analyze non-spherical systems. We discuss how the code works and discuss its assumptions. Furthermore, we demonstrate how SPAMMS can be applied to a variety of different types of systems and we show how it can model 3D effects in a way that current analysis techniques are not able to.

Large inner dust gaps in transition disks are frequently posited as evidence of giant planets sculpting gas and dust in the disk, or the opening of a gap by photoevaporative winds. Although the former hypothesis is strongly supported by the observations of planets and deep depletions in gas within the gap some disks, many T Tauri stars hosting transition disks accrete at rates typical for an undepleted disk, raising the question of how gap opening occurs in these objects. We thus present an analysis of the structure of the transition disk around the T Tauri star DM Tau, which is strongly accreting ($\sim 10^{-8.3}~\mathrm{M}_\odot~ \mathrm{yr}^{-1}$) and turbulent ($\alpha=0.078 \pm 0.02$). Using the DALI thermochemical code, we fit disk models to simultaneously reproduce the accretion rate, high level of turbulence, the gas traced by ALMA band 6 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O J=2--1 lines, and the observed dust emission from the mm continuum and spectral energy distribution. We find a shallow depletion in gas surface density of $\sim 10$ relative to the outer disk and a gas rich inner disk is consistent with the observations. The planet mass of $<1$ M$_\mathrm{Jup}$ implied by the gap depth is in tension with predictions for dust trapping in a highly viscous disk, which requires a more massive planet of of $\sim10$M$_\mathrm{Jup}$. Photoevaporative models including a dead zone can qualitatively reproduce some features of the DM Tau disk, but still struggle to explain the high accretion rates and the observed mm continuum flux.

With JWST, new opportunities to study the formation and evolution of galaxies in the early Universe are now emerging. Spitzer constraints on rest-optical properties of $z \gtrsim 7$ galaxies demonstrated the power of using stellar masses and star formation histories (SFHs) of galaxies to indirectly infer the star formation history of the Universe. However, only the brightest individual objects at $z \gtrsim 8$ could be detected with Spitzer, making it difficult to robustly constrain past activity at $z \gtrsim 10$. Here, we leverage the greatly improved rest-optical sensitivity of JWST at $z \gtrsim 8$ to constrain the ages and SFHs of eleven UV-bright ($M_\text{UV} \lesssim -19.5$) galaxies selected to lie at $z \sim 8.5 - 11$, then investigate implications for star formation activity at $z \gtrsim 15$. We infer the properties of individual objects in our sample with two spectral energy distribution modelling codes, then infer a distribution of ages for bright $z \sim 8.5 - 11$ galaxies. We find a median age of $\sim 30$ Myr, younger than that inferred at $z \sim 7$ with a similar analysis, which is consistent with an evolution towards larger specific star formation rates at early times. The age distribution suggests that only $\sim 9$ percent of bright $z \sim 8.5 - 11$ galaxies would be similarly luminous at $z \gtrsim 15$, implying that the number density of bright galaxies declines by approximately an order of magnitude between $z \sim 8.5 - 11$ and $z \sim 15$. This evolution is challenging to reconcile with some early JWST results suggesting that the abundance of bright galaxies does not significantly decrease towards very early times, but we suggest this tension may be eased if young stellar populations form on top of older stellar components, or if bright galaxies at $z \sim 15$ are observed during a burst of star formation.

Despite having long been the standard for quantifying relief on Earth and beyond, elevation has its limitations. The zero-elevation datum is defined by arbitrary and inconsistent conventions, especially on planets without a sea level, hence the lack of a universally standardized way to quantify relief. Furthermore, when quantifying relief on such planets, the elevation of a point is rather meaningless on its own, deriving most of its value when compared to the elevation of other points. In light of these considerations, this paper introduces a universally consistent framework for quantifying relief that does not require a datum altogether, and is instead based on physically meaningful concepts. Designed to be mathematically elegant and free of arbitrary parameters, the so-called datumless measures are divided into the datumless point measures and the datumless surface measures. As opposed to elevation, which describes the vertical position of a point relative to a datum, the datumless point measures directly describe the vertical position of a point relative to local terrain, making them useful for comparing the relief of features such as mountains across different planets. In the meantime, the datumless surface measures quantify various aspects of relief within a region, as opposed to that of a single point. This is done through datumless formulations of surface area and surface mean value, which can be directly applied to the fractal-like planetary surface without projecting it onto a reference ellipsoid. Altogether, this paper lays the groundwork for a datumless framework that enables future topographic tasks to transcend the limitations of elevation.

High-contrast imaging presents us with the opportunity to study circumstellar disks and the planets still embedded within them -- providing key insights into the formation and evolution of planetary systems. However, the post-processing techniques that are often needed to suppress stellar halo light typically result in significant and variable loss of circumstellar light, even when using relatively conservative approaches like reference star differential imaging (RDI). We introduce ``constrained reference star differential imaging" (constrained RDI), a new class of RDI point spread function (PSF)-subtraction techniques for systems with circumstellar disks. Constrained RDI utilizes either high-resolution polarized intensity (PI) images or disk models to severely limit or even eliminate the signal loss due to oversubtraction that is common to RDI. We demonstrate the ability of constrained RDI utilizing polarimetric data to yield an oversubtraction-free detection of the AB Aurigae protoplanetary disk in total intensity. PI-constrained RDI allows us to decisively recover the spectral signature of the confirmed, recently-discovered protoplanet, AB Aurigae b (Currie et al. 2022). We further demonstrate that constrained RDI can be a powerful analysis tool for soon-to-be-acquired James Webb Space Telescope coronagraphic imaging of disks. In both cases, constrained RDI provides analysis-ready products that enable more detailed studies of disks and more robust verification of embedded exoplanets.

The ultraviolet (UV) continuum slope ($\beta$ where f$_\lambda\propto \lambda^\beta$) of galaxies is sensitive to a variety of properties, from the metallicity and age of the stellar population to the attenuation from dust through the galaxy. Considerable attention has focused on identifying reionization-era galaxies with very blue UV slopes ($\beta<-3$). Not only do such systems provide a signpost of low metallicity stars, but they also identify galaxies that likely have ionizing photons leaking from their HII regions as such blue UV slopes can only be seen if the reddening effect of nebular continuum has been diminished. In this paper we present a search for reionization-era galaxies with very blue UV colors in recent JWST/NIRCam imaging of the EGS field. We characterize UV slopes for a large sample of $z\simeq 7-11$ galaxies, finding a median value of $\beta =-2.1$. Three of the lower luminosity (M$_{\rm{UV}}\simeq -19.5$) and lower stellar mass (5-6$\times10^7$M$_\odot$) systems exhibit both extremely blue UV slopes ($\beta=-3.1$ to $-3.2$) and rest-optical photometry indicating weak nebular line emission. Each system is very compact (r$_e<$260 pc) with very high star formation rate surface densities. We model the SEDs with a suite of BEAGLE models with varying levels of ionizing photon escape. The SEDs cannot be reproduced with our fiducial (f$_{\rm{esc,HII}}$=0) or alpha enhanced (Z$_*<Z_{\rm{ISM}}$) models. The combined blue UV slopes and weak nebular emission are best-fit by models with significant ionizing photon escape from HII regions (f$_{\rm{esc,HII}}$=0.6-0.8) and extremely low metallicity massive stars (Z$_*$=0.01-0.06 Z$_\odot$). The discovery of these galaxies highlights the potential for JWST to identify large numbers of candidate Lyman Continuum leaking galaxies in the reionization era and suggests low metallicity stellar populations may be very

We conduct a comprehensive study on dropout galaxy candidates at $z\sim 9-17$ using the first 90 arcmin$^2$ JWST/NIRCam images taken by the early release observations (ERO) and early release science (ERS) programs. With the JWST simulation images, we find that a number of foreground interlopers are selected with a weak photo-$z$ determination ($\Delta \chi^2>4$). We thus carefully apply a secure photo-$z$ selection criterion ($\Delta \chi^2>9$) and conventional color criteria with confirmations of the ERO NIRSpec spectroscopic redshifts, and obtain a total of 25 dropout galaxies at $z\sim 9-17$, including two candidates at $z_\mathrm{phot}=16.45_{-0.32}^{+0.09}$ and $16.66_{-0.34}^{+1.86}$. We perform thorough comparisons of dropout galaxies found in our work with recent JWST studies, and conclude that our galaxy sample is reliable enough for statistical analyses. We derive the UV luminosity functions at $z\sim 9-17$, and confirm that our UV luminosity functions at $z\sim 9$ and $12$ agree with those determined by previous HST and JWST studies. The cosmic star-formation rate density decreases from $z\sim 9$ to $12$, and perhaps to $17$, but the densities at $z\sim12-17$ are higher than the constant star formation efficiency model. Interestingly, there are six bright galaxy candidates at $z\sim 11-17$ with $M_{\rm UV}<-19.5$ whose stellar masses are very high, $10^{8-9} M_\odot$. Because a majority ($\sim 70\%$) of these galaxies shows no signatures of AGNs in their morphologies, the high cosmic star-formation rate densities and the existence of these stellar massive galaxies are explained by no suppression of star-formation by the UV background radiation at the pre-reionization epoch or an efficient UV radiation production by Population III-like star formation.

We present a self-consistent cross-calibration of the three main molecular gas mass tracers in galaxies, the $\rm ^{12}CO$(1-0), [CI]($^3P_1$-$^3P_0$) lines, and the submm dust continuum emission, using a sample of 407 galaxies, ranging from local disks to submillimetre-selected galaxies (SMGs) up to $z \approx 6$. A Bayesian method is used to produce galaxy-scale universal calibrations of these molecular gas indicators, that hold over 3-4 orders of magnitude in infrared luminosity, $L_{\rm IR}$. Regarding the dust continuum, we use a mass-weighted dust temperature, $T_{\rm mw}$, determined using new empirical relations between temperature and luminosity. We find the average $L/M_{\rm mol}$ gas mass conversion factors to be $\alpha_{850}= 6.9\times10^{12}\,\rm W\,Hz^{-1}\,M_{\odot}^{-1}$, $\alpha_{\rm CO} = \rm 4\,M_{\odot} (K\,km\,s^{-1}\,pc^2)^{-1}$ and $\alpha_{\rm CI} = \rm 17.0 \,M_{\odot} (K\,km\,s^{-1}\,pc^2)^{-1}$, based on the assumption that the mean dust properties of the sample ($\kappa_H$ = gas-to-dust ratio/dust emissivity) will be similar to those of local metal rich galaxies and the MW. The tracer with the least intrinsic scatter is [CI](1-0), while CO(1-0) has the highest. The conversion factors show a weak but significant correlation with $L_{\rm IR}$. Assuming dust properties typical of metal-rich galaxies, we infer a neutral carbon abundance $X_{\rm CI} = [C^0/\rm mol]=1.6\times 10^{-5}$, similar to that in the MW. We find no evidence for bimodality of $\alpha_{\rm CO}$ between main-sequence (MS) galaxies and those with extreme star-formation intensity, i.e. ULIRGs and SMGs. The means of the three conversion factors are found to be similar between MS galaxies and ULIRGs/SMGs, to within 10-20%. We show that for metal-rich galaxies, near-universal average values for $\alpha_{\rm CO}$, $X_{\rm CI}$ and $\kappa_H$ are adequate for global molecular gas estimates.

Optically quiescent quasars (OQQs) represent a recently systematised class of infrared-luminous active galactic nuclei (AGN) which have galaxy-like optical continua. They may represent an interesting, brief phase in the AGN life cycle, e.g. either cocooned within high-covering-factor media or indicative of recent triggering, though their nature remains unclear. Here, we present the first targeted simultaneous X-ray observations of an OQQ, our previously identified prototype, SDSS J075139.06+402811.2 at $z$=0.587. The source is significantly detected over 0.5-16 keV with XMM-Newton and NuSTAR, unambiguously confirming the presence of current accretion activity. Spectral modelling yields an intrinsic luminosity $L_{\rm 2-10 keV}$ $\approx$4.4$ \times$10$^{43}$ erg s$^{-1}$, well within the AGN regime, but underluminous relative to its infrared power. It is lightly obscured, with log $N_{\rm H}$ [cm$^{-2}$] $\approx$22.

With just a month of data, JWST is already transforming our view of the Universe, revealing and resolving starlight in unprecedented populations of galaxies. Although ``HST-dark" galaxies have previously been detected at long wavelengths, these observations generally suffer from a lack of spatial resolution which limits our ability to characterize their sizes and morphologies. Here we report on a first view of starlight from a subset of the HST-dark population that are bright with JWST/NIRCam (4.4$\mu$m<24.5mag) and very faint or even invisible with HST ($<$1.6$\mu$m). In this Letter we focus on a dramatic and unanticipated population of physically extended galaxies ($\gtrsim$0.17''). These 12 galaxies have photometric redshifts $2<z<6$, high stellar masses $M_{\star}\gtrsim 10^{10}~M_{\odot}$, and significant dust-attenuated star formation. Surprisingly, the galaxies have elongated projected axis ratios at 4.4$\mu$m, suggesting that the population is disk-dominated or prolate. Most of the galaxies appear red at all radii, suggesting significant dust attenuation throughout. We refer to these red, disky, HST-dark galaxies as Ultra-red Flattened Objects (UFOs). With $r_e$(F444W)$\sim1-2$~kpc, the galaxies are similar in size to compact massive galaxies at $z\sim2$ and the cores of massive galaxies and S0s at $z\sim0$. The stellar masses, sizes, and morphologies of the sample suggest that some could be progenitors of lenticular or fast-rotating galaxies in the local Universe. The existence of this population suggests that our previous censuses of the universe may have missed massive, dusty edge-on disks, in addition to dust-obscured starbursts.

A new sufficient condition for photons emitted near future null infinity to reach future null infinity is derived by studying null geodesics in the Bondi coordinates in asymptotically flat spacetimes. In our previous works [arXiv:2106.03150, arXiv:2110.10917], such a condition was established for photons emitted in outward or tangential directions to constant radial surfaces. This paper improves our previous result by including photons emitted in inward directions. In four dimensions, imposing the same assumptions on the metric functions as [arXiv:2106.03150, arXiv:2110.10917], we prove that photons reach future null infinity if their initial values of $|dr/du|$ are smaller than a certain quantity, where $r$ and $u$ are the radial and retarded time coordinates, respectively. This quantity is determined by the asymptotic properties of the metric and is connected to the conjectured maximal luminosity. In higher dimensions, photons emitted with $dr/du>-(1-1/\sqrt{3})\approx -0.423$ are shown to reach future null infinity without the assumptions on the metric functions.

Based on entropy considerations and the arrow of time Penrose argued that the universe must have started in a special initial singularity with vanishing Weyl curvature. This is often interpreted to be at odds with inflation. Here we argue just the opposite, that Penrose's persuasions are in fact consistent with inflation. Using the example of power law inflation, we show that inflation begins with a past null singularity, where Weyl tensor vanishes when the metric is initially exactly conformally flat. This initial state precisely obeys Penrose's conditions. The initial null singularity breaks $T$-reversal spontaneously and picks the arrow of time. It can be regulated and interpreted as a creation of a universe from nothing, initially fitting in a bubble of Planckian size when it materializes. Penrose's initial conditions are favored by the initial $O(4)$ symmetry of the bubble, selected by extremality of the regulated Euclidean action. The predicted observables are marginally in tension with the data, but they can fit if small corrections to power law inflation kick in during the last 60 efolds.

We investigate a new possible solution to the Hubble constant tension. we propose a simple resolution to the problem assuming that a first-order phase transition related to H0 transition occurred in the early Universe. The early evolution of the Universe is a result of hybrid inflation that has lasted for a specific period until symmetry breaking takes place. Fitting our model to measurements from Planck and SH0ES data provides a key explanation of discrepancies of H0 measurements. The quantum fluctuations calculated in this model have significant results on the reheating parameters Nre and Tre. Therefore, new constraints must be taken into consideration to fit these parameters to recent results.

In this paper we discuss a detection method for the Cosmic Neutrino Background using bremsstrahlung from a neutrino scattering process which has no kinematic threshold, does not rely on a resonance and would in principle allow to measure the velocity distribution of the relic neutrinos. As a concrete example we calculate the rate for solar neutrinos scattering from a relic neutrino emitting a photon. We also provide the energy and angular distributions of the emitted photons.

There is growing interest in using current and future gravitational-wave interferometers to search for anisotropies in the gravitational-wave background. One guaranteed anisotropic signal is the kinematic dipole induced by our peculiar motion with respect to the cosmic rest frame, as measured in other full-sky observables such as the cosmic microwave background. Our prior knowledge of the amplitude and direction of this dipole is not explicitly accounted for in existing searches by LIGO/Virgo/KAGRA, but could provide crucial information to help disentangle the sources which contribute to the gravitational-wave background. Here we develop a targeted search pipeline which uses this prior knowledge to enable unbiased and minimum-variance inference of the dipole magnitude. Our search generalises existing methods to allow for a time-dependent signal model, which captures the annual modulation of the dipole due to the Earth's orbit. We validate our pipeline on mock data, demonstrating that neglecting this time dependence can bias the inferred dipole by as much as $\sim10\%$. We then run our analysis on the full LIGO/Virgo O1+O2+O3 dataset, obtaining upper limits on the dipole amplitude that are consistent with existing anisotropic search results.

Electro-thermal coupling in semiconductor bolometers is known to create nonlinearities in transient detector response, particularly when such detectors are biased outside of their ideal regions (i.e. past the turnover point in their IV curves). This effect is further compounded in the case where a stray capacitance in the bias circuit is present, for example in long cryogenic cabling. We present a physical model of the influence of such electro-thermal coupling and stray capacitance in a composite NTD germanium bolometer, in which previous experimental data at high $V_{\rm bias}$ resulted in oscillations of the impulse response of the detector to irradiation by alpha particles. The model reproduces the transient oscillations seen in the experimental data, depending both on electro-thermal coupling and stray capacitance. This is intended as an experimental and simulated example of such oscillations, demonstrated for the specific case of this bolometric detector.

We investigate transient radiation processes in the non-linear superradiance (SR) regime of the Doppler broadened Maxwell-Bloch equations when the velocity distribution is of total bandwidth greatly exceeding that of the transient process itself. We demonstrate the formation of global polarisation phase correlation and the quenching of temporal structure if a smooth distribution is inverted above the critical threshold required to enter the non-linear SR regime. We propose candidate stochastic velocity distributions capable of sustaining finite temporal structure in the non-linear emission process. We develop a novel algorithm for simulating the Doppler broadened Maxwell-Bloch equations which is $O(n)$ complex in the number of velocity channels $n$ whenever the emerging polarisation correlation is of moderate bandwidth, and we apply it to a stochastic velocity distribution in order to demonstrate sustained delay and duration of peak intensity in the widely Doppler broadened limit. We discuss the transverse inversion process and recognise an autoregulation mechanism on the number of molecules cooperatively participating in SR emission. This mechanism has the effect of limiting the temporal duration of the intensity pulse to a lower bound proportional to the length of the sample, which we confirm through simulation.

We investigate how GWs pass through the spacetime of a Schwarzschild black hole using time-domain numerical simulations. Our work is based on the perturbed 3+1 Einstein's equations up to the linear order. We show explicitly that our perturbation equations are covariant under infinitesimal coordinate transformations. Then we solve a symmetric second-order hyperbolic wave equation with a spatially varying wave speed. As the wave speed in our wave equation vanishes at the horizon, our formalism can naturally avoid boundary conditions at the horizon. Our formalism also does not contain coordinate singularities and, therefore, does not need regularity conditions. Then, based on our code, we simulate both finite and continuous initially plane-fronted wave trains passing through the Schwarzschild black hole. We find that for the finite wave train, the wave zone of GWs is wildly twisted by the black hole. While for the continuous wave train, unlike geometric optics, GWs can not be sheltered by the back hole. A strong beam and an interference pattern appear behind the black hole along the optical axis. Moreover, we find that the back-scattering due to the interaction between GWs and the background curvature is strongly dependent on the direction of the propagation of the trailing wavefront relative to the black hole.