Papers for Thursday, Jan 04 2024

We present 3D fully kinetic shearing-box (SB) simulations of pair-plasma magnetorotational turbulence with unprecedented macro-to-microscopic scale separation. We retrieve the expected fluid-model behavior of turbulent magnetic fields and angular-momentum transport, and observe fundamental differences in turbulent fluctuation spectra linked with plasma heating. For the first time, we provide a definitive demonstration of nonthermal particle acceleration in kinetic magnetorotational turbulence agnostically of SB initial conditions, by means of a novel strategy exploiting synchrotron cooling.

We present the mid-infrared (5-12 $\mu$m) phase curve of GJ 367b observed by the Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope (JWST). GJ 367b is a hot (T_eq=1370 K), extremely dense (10.2 +- 1.3 g/cm^3) sub-Earth orbiting a M dwarf on a 0.32 d orbit. We measure an eclipse depth of 79 +- 4 ppm, a nightside planet-to-star flux ratio of 4 +- 8 ppm, and a relative phase amplitude of 0.97 +- 0.10-- all fully consistent with a zero-albedo planet with no heat recirculation. Su ch a scenario is also consistent with the phase offset of 11 +- 5 degrees East to within 2.2$\sigma$. The emission spectrum is likewise consistent with a blackbody with no heat redistribution and a low albedo of A_B ~ 0.1, with the exception of one anomalous wavelength bin that we attribute to unexplained systematics. The emission spectrum puts few constraints on the surface composition, but rules out a CO2 atmosphere >~ 1 bar, an outgassed atmosphere >~10 mbar (under heavily reducing conditions), or an outgassed atmosphere >~0.01 mbar (under heavily oxidizing conditions). The lack of day-night heat recirculation implies that 1 bar atmospheres are ruled out for a wide range of compositions, while 0.1 bar atmospheres are consistent with the data. Taken together with the fact that most of the dayside should be molten, our JWST observations suggest the planet must have lost the vast majority of its initial inventory of volatiles.

Most tidal disruption events (TDEs) are currently found in time-domain optical and soft X-ray surveys, both of which are prone to significant obscuration. The infrared (IR), however, is a powerful probe of dust-enshrouded environments, and hence, we recently performed a systematic search of NEOWISE mid-IR data for nearby, obscured TDEs within roughly 200 Mpc. We identified 18 TDE candidates in galactic nuclei, using difference imaging to uncover nuclear variability amongst significant host galaxy emission. These candidates were selected based on the following IR light curve properties: (1) $L_\mathrm{W2}\gtrsim10^{42}$ erg s$^{-1}$ at peak, (2) fast rise, followed by a slow, monotonic decline, (3) no significant prior variability, and (4) no evidence for AGN activity in WISE colors. The majority of these sources showed no variable optical counterpart, suggesting that optical surveys indeed miss numerous obscured TDEs. Using narrow line ionization levels and variability arguments, we identified 6 sources as possible underlying AGN, yielding a total of 12 TDEs in our gold sample. This gold sample yields a lower limit on the IR-selected TDE rate of $2.0\pm0.3\times10^{-5}$ galaxy$^{-1}$ year$^{-1}$ ($1.3\pm0.2\times10^{-7}$ Mpc$^{-3}$ year$^{-1}$), which is comparable to optical and X-ray TDE rates. The IR-selected TDE host galaxies do not show a green valley overdensity nor a preference for quiescent, Balmer strong galaxies, which are both overrepresented in optical and X-ray TDE samples. This IR-selected sample represents a new population of dusty TDEs that have historically been missed by optical and X-ray surveys and helps alleviate tensions between observed and theoretical TDE rates and the so-called missing energy problem.

Twenty-two extragalactic fast X-ray transients (FXTs) have now been discovered from two decades of Chandra data (analyzing ~259 Ms of data), with 17 associated with distant galaxies (>100 Mpc). Different mechanisms and progenitors have been proposed to explain their properties; nevertheless, after analyzing their timing, spectral parameters, host-galaxy properties, luminosity function, and volumetric rates, their nature remains uncertain. We interpret a sub-sample of nine FXTs that show a plateau or a fast-rise light curve within the framework of a binary neutron star (BNS) merger magnetar model. We fit their light curves and derive magnetar (magnetic field and initial rotational period) and ejecta (ejecta mass and opacity) parameters. This model predicts two zones: an orientation-dependent free zone (where the magnetar spin-down X-ray photons escape freely to the observer) and a trapped zone (where the X-ray photons are initially obscured and only escape freely once the ejecta material becomes optically thin). We argue that six FXTs show properties consistent with the free zone and three FXTs with the trapped zone. This sub-sample of FXTs has a similar distribution of magnetic fields and initial rotation periods to those inferred for short gamma-ray bursts (SGRBs), suggesting a possible association. We compare the predicted ejecta emission fed by the magnetar emission (called merger-nova) to the optical and near-infrared upper limits of two FXTs, XRT 141001 and XRT 210423 where contemporaneous optical observations are available. The non-detections place lower limits on the redshifts of XRT 141001 and XRT 210423 of z>1.5 and >0.1, respectively. If the magnetar remnants lose energy via gravitational waves, it should be possible to detect similar objects with the current advanced LIGO detectors out to a redshift z<0.03, while future GW detectors will be able to detect them out to z=0.5.

Scientists have long speculated about the potential habitability of Venus, not at the 700K surface, but in the cloud layers located at 48-60 km altitudes, where temperatures match those found on Earth's surface. However, the prevailing belief has been that Venus' clouds cannot support life due to the cloud chemical composition of concentrated sulfuric acid-a highly aggressive solvent. In this work, we study 20 biogenic amino acids at the range of Venus' cloud sulfuric acid concentrations (81% and 98% w/w, the rest water) and temperatures. We find 19 of the biogenic amino acids we tested are either unreactive (13 in 98% w/w and 12 in 81% w/w) or chemically modified in the side chain only, after four weeks. Our major finding, therefore, is that the amino acid backbone remains intact in concentrated sulfuric acid. These findings significantly broaden the range of biologically relevant molecules that could be components of a biochemistry based on a concentrated sulfuric acid solvent.

In ultraviolet (UV) astronomical observations, photons from the sources are very few compared to the visible or infrared (IR) wavelength ranges. Detectors operating in the UV usually employ a photon-counting mode of operation. These detectors usually have an image intensifier sensitive to UV photons and a readout mechanism that employs photon counting. The development of readouts for these detectors is resource-intensive and expensive. In this paper, we describe the development of a low-cost UV photon-counting detector processing unit that employs a Raspberry Pi with its in built readout to perform the photon-counting operation. Our system can operate in both 3x3 and 5x5 window modes at 30 frames per sec (fps), where 5x5 window mode also enables the provision of detection of double events. The system can be built quickly from readily available custom-off-the-shelf (COTS) components and is thus used in inexpensive CubeSats or small satellite missions. This low-cost solution promises to broaden access to UV observations, advancing research possibilities in space-based astronomy.

We present a comparative study of the molecular gas in two galaxies from the LEGUS sample: barred spiral NGC 1313 and flocculent spiral NGC 7793. These two galaxies have similar masses, metallicities, and star formation rates, but NGC 1313 is forming significantly more massive star clusters than NGC 7793, especially young massive clusters (<10 Myr, >10^4 Msol). Using ALMA CO(2-1) observations of the two galaxies with the same sensitivities and resolutions of 13 pc, we directly compare the molecular gas in these two similar galaxies to determine the physical conditions responsible for their large disparity in cluster formation. By fitting size-linewidth relations for the clouds in each galaxy, we find that NGC 1313 has a higher intercept than NGC 7793, implying that its clouds have higher kinetic energies at a given size scale. NGC 1313 also has more clouds near virial equilibrium than NGC 7793, which may be connected to its higher rate of massive cluster formation. However, these virially bound clouds do not show a stronger correlation with young clusters than that of the general cloud population. We find surprisingly small differences between the distributions of molecular cloud populations in the two galaxies, though the largest of those differences are that NGC 1313 has higher surface densities and lower free-fall times.

We compare the molecular cloud properties in sub-galactic regions of two galaxies, barred spiral NGC 1313, which is forming many massive clusters, and flocculent spiral NGC 7793, which is forming significantly fewer massive clusters despite having a similar star formation rate to NGC 1313. We find that there are larger variations in cloud properties between different regions within each galaxy than there are between the galaxies on a global scale, especially for NGC 1313. There are higher masses, linewidths, pressures, and virial parameters in the arms of NGC 1313 and center of NGC 7793 than in the interarm and outer regions of the galaxies. The massive cluster formation of NGC 1313 may be driven by its greater variation in environments, allowing more clouds with the necessary conditions to arise, although no one parameter seems primarily responsible for the difference in star formation. Meanwhile NGC 7793 has clouds that are as massive and have as much kinetic energy as clouds in the arms of NGC 1313, but have densities and pressures more similar to the interarm regions and so are less inclined to collapse and form stars. The cloud properties in NGC 1313 and NGC 7793 suggest that spiral arms, bars, interarm regions, and flocculent spirals each represent distinct environments with regard to molecular cloud populations. We see surprisingly little difference in surface densities between the regions, suggesting that the differences in surface densities frequently seen between arm and interarm regions of lower-resolution studies are indicative of the sparsity of molecular clouds, rather than differences in their true surface density.

The European Space Agency's Euclid mission is one of the upcoming generation of large-scale cosmology surveys, which will map the large-scale structure in the Universe with unprecedented precision. The development and validation of the SGS pipeline requires state-of-the-art simulations with a high level of complexity and accuracy that include subtle instrumental features not accounted for previously as well as faster algorithms for the large-scale production of the expected Euclid data products. In this paper, we present the Euclid SGS simulation framework as applied in a large-scale end-to-end simulation exercise named Science Challenge 8. Our simulation pipeline enables the swift production of detailed image simulations for the construction and validation of the Euclid mission during its qualification phase and will serve as a reference throughout operations. Our end-to-end simulation framework starts with the production of a large cosmological N-body & mock galaxy catalogue simulation. We perform a selection of galaxies down to I_E=26 and 28 mag, respectively, for a Euclid Wide Survey spanning 165 deg^2 and a 1 deg^2 Euclid Deep Survey. We build realistic stellar density catalogues containing Milky Way-like stars down to H<26. Using the latest instrumental models for both the Euclid instruments and spacecraft as well as Euclid-like observing sequences, we emulate with high fidelity Euclid satellite imaging throughout the mission's lifetime. We present the SC8 data set consisting of overlapping visible and near-infrared Euclid Wide Survey and Euclid Deep Survey imaging and low-resolution spectroscopy along with ground-based. This extensive data set enables end-to-end testing of the entire ground segment data reduction and science analysis pipeline as well as the Euclid mission infrastructure, paving the way to future scientific and technical developments and enhancements.

We present a comprehensive analysis of the Hubble Space Telescope observations of the atmosphere of WASP-121 b, a ultra-hot Jupiter. After reducing the transit, eclipse, and phase-curve observations with a uniform methodology and addressing the biases from instrument systematics, sophisticated atmospheric retrievals are used to extract robust constraints on the thermal structure, chemistry, and cloud properties of the atmosphere. Our analysis shows that the observations are consistent with a strong thermal inversion beginning at ~0.1 bar on the dayside, solar to subsolar metallicity Z (i.e., -0.77 < log(Z) < 0.05), and super-solar C/O ratio (i.e., 0.59 < C/O < 0.87). More importantly, utilizing the high signal-to-noise ratio and repeated observations of the planet, we identify the following unambiguous time-varying signals in the data: i) a shift of the putative hotspot offset between the two phase-curves and ii) varying spectral signatures in the transits and eclipses. By simulating the global dynamics of WASP-121 b atmosphere at high-resolution, we show that the identified signals are consistent with quasi-periodic weather patterns, hence atmospheric variability, with signatures at the level probed by the observations (~5% to ~10%) that change on a timescale of ~5 planet days; in the simulations, the weather patterns arise from the formation and movement of storms and fronts, causing hot (as well as cold) patches of atmosphere to deform, separate, and mix in time.

In order to study exoplanets, a comprehensive characterization of the fundamental properties of the host stars, such as angular diameter, temperature, luminosity, and age, is essential, as the formation and evolution of exoplanets are directly influenced by the host stars at various points in time. In this paper, we present interferometric observations taken of directly imaged planet host 51 Eridani at the CHARA Array. We measure the limb-darkened angular diameter of HD 29391 to be $\theta_{\rm LD} = 0.450\pm 0.004 \rm ~mas$ and combining with the Gaia zero-point corrected parallax, we get a stellar radius of $1.45 \pm 0.01 \rm~R_{\odot}$. We use the PARSEC isochrones to estimate an age of $23.2^{+1.7}_{-1.6} \rm ~Myr$ and a mass of $1.550 \pm 0.005 \rm ~M_{\odot}$. The age and mass agree well with values in the literature, determined through a variety of methods ranging from dynamical age trace-backs to lithium depletion boundary methods. We derive a mass of $4.1^{+0.5}_{-0.4} \rm ~M_{Jup}$ for 51 Eri b using the Sonora Bobcat models, which further supports the possibility of 51 Eri b forming under either the hot-start formation model or the warm-start formation model.

We report the observations of two self-lensing pulses from KIC 12254688 in Transiting Exoplanet Survey Satellite (TESS) light curves. This system, containing a F2V star and white-dwarf companion, was amongst the first self-lensing binary systems discovered by the Kepler Space Telescope over the past decade. Each observed pulse occurs when the white dwarf transits in front of its companion star, gravitationally lensing the star's surface, thus making it appear brighter to a distant observer. These two pulses are the very first self-lensing events discovered in TESS observations. We describe the methods by which the data were acquired and detrended, as well as the best-fit binary parameters deduced from our self-lensing+radial velocity model. We highlight the difficulties of finding new self-lensing systems with TESS, and we discuss the types of self-lensing systems that TESS may be more likely to discover in the future.

PSR J1012+5307 is a millisecond pulsar with an extremely low-mass (ELM) white dwarf (WD) companion in an orbit of 14.5 hours. Magnetic braking (MB) plays an important role in influencing the orbital evolution of binary systems with a low-mass ($\lt 1-2~M_{\odot}$) donor star. At present, there exist several different MB descriptions. In this paper, we investigate the formation of PSR J1012+5307 as a probe to test the plausible MB model. Employing a detailed stellar evolution model by the MESA code, we find that the Convection And Rotation Boosted MB and the 'Intermediate' MB models can reproduce the WD mass, WD radius, WD surface gravity, neutron-star mass, and orbital period observed in PSR J1012+5307. However, our simulated WD has higher effective temperature than the observation. Other three MB mechanisms including the standard MB model are too weak to account for the observed orbital period in a Hubble time. A long cooling timescale caused by H-shell flashes of the WD may alleviate the discrepancy between the simulated effective temperature and the observed value.

Linear RMS-flux relation has been well established in different spectral states of all accreting systems. In this work, we study the evolution of the frequency-dependent RMS-flux relation of MAXI J1820+070 during the initial decaying phase of the 2018 outburst with Insight-HXMT over a broad energy range 1-150 keV. As the flux decreases, we first observe a linear RMS-flux relation at frequencies from 2 mHz to 10 Hz, while such a relation breaks at varying times for different energies, leading to a substantial reduction in the slope. Moreover, we find that the low-frequency variability exhibits the highest sensitivity to the break, which occurs prior to the hard-to-hard state transition time determined through time-averaged spectroscopy, and the time deviation increases with energy. The overall evolution of the RMS-flux slope and intercept suggests the presence of a two-component Comptonization system. One component is radially extended, explaining the strong disk-corona coupling before the break, while the other component extends vertically, contributing to the reduction of the disk-corona coupling after the break. A further vertical expansion of the latter component is required to accommodate the dynamic evolution observed in the RMS-flux slope. In conclusion, we suggest that the RMS-flux slope in 1-150 keV band can be employed as an indicator of the disk-corona coupling and the hard-to-hard state transition in MAXI J1820+070 could be partially driven by the changes in the corona geometry.

We compare, with data from the quasars, the Hubble parameter measurements, and the Pantheon+ type Ia supernova, three different relations between X-ray luminosity ($L_X$) and ultraviolet luminosity ($L_{UV}$) of quasars. These three relations consist of the standard and two redshift-evolutionary $L_X$-$L_{UV}$ relations which are constructed respectively by considering a redshift dependent correction to the luminosities of quasars and using the statistical tool called copula. By employing the PAge approximation for a cosmological-model-independent description of the cosmic background evolution and dividing the quasar data into the low-redshift and high-redshift parts, we find that the constraints on the PAge parameters from the low-redshift and high-redshift data, which are obtained with the redshift-evolutionary relations, are consistent with each other, while they are not when the standard relation is considered. If the data are used to constrain the coefficients of the relations and the PAge parameters simultaneously, then the observations support the redshift-evolutionary relations at more than $3\sigma$. The Akaike and Bayes information criteria indicate that there is strong evidence against the standard relation and mild evidence against the redshift-evolutionary relation constructed by considering a redshift dependent correction to the luminosities of quasars. This suggests that the redshift-evolutionary $L_X$-$L_{UV}$ relation of quasars constructed from copula is favored by the observations.

The illumination phase functions for Starlink Mini satellites are determined for times of twilight and darkness. Those functions are then evaluated to give apparent magnitudes over a grid of points across the sky and over a range of solar angles below the horizon. Sky maps and a table of satellite magnitude distributions are presented. The largest areas of sky with satellites brighter than magnitudes 6 and 7 both occur during twilight. Brightness surges, known as flares, are also characterized.

There is growing evidence that high-mass star formation (HMSF) is a multiscale, dynamical process in molecular clouds, where filaments transport gas material between larger and smaller scales. We analyze here multiscale gas dynamics in an HMSF filamentary cloud, G034.43+00.24 (G34), using APEX observations of the C18O (2-1), HCO+/H13CO+ (3-2), and HCN/H13CN (3-2) lines. We find large-scale, filament-aligned velocity gradients from C18O emission, which drive filamentary gas inflows onto dense clumps in the middle ridge of G34. The nature of these inflows is gravity driven. We also find clump-scale gas infall in the middle ridge of the MM2, MM4, and MM5 clumps from other lines. Their gas infall rates could depend on large-scale filamentary gas inflows since the infall/inflow rates on these two scales are comparable. We confirm that the multiscale, dynamical HMSF scenario is at work in G34. It could be driven by gravity up to the filament scale, beyond which turbulence originating from several sources, including gravity, could be in effect in G34.

For long wavelength gravitational wave (GW), it is easy to diffract when it is lensed by celestial objects. Traditional diffractive integral formula has ignored large angle diffraction, which is adopted in most of cases. However, in some special cases (e. g. a GW source lensed by its companion in a binary system, where the lens is very close to the source), large angle diffraction could be important. Our previous works have proposed a new general diffractive integral formula which has including large angle diffraction case. In this paper, we have investigated how much difference between this general diffractive formula and traditional diffractive integral formula could be under these special cases with different parameters. We find that the module of amplification factor for general diffractive formula could become smaller than that of traditional diffractive integral basically with a factor $r_F\simeq0.674$ when the distance between lens and sources is $D_{\rm LS}=1$ AU and lens mass $M_{\rm L}=1M_\odot$. Their difference is so significant that it is detectable. Furthermore, we find that the proportionality factor $r_F$ is gradually increasing from 0.5 to 1 with increasing $D_{\rm LS}$ and it is decreasing with increasing $M_{\rm L}$. As long as $D_{\rm LS}\lesssim3$ AU (with $M_{\rm L}=1M_\odot$) or $M_{\rm L}\gtrsim0.1M_\odot$ (with $D_{\rm LS}=1$ AU ), the difference between new and traditional formulas is enough significant to be detectable. It is promising to test this new general diffractive formula by next-generation GW detectors in the future GW detection.

An important aspect of solar energetic particle (SEP) events is their source populations. Elemental abundance enhancements of impulsive SEP events, originating in presumed coronal reconnection episodes, can be fitted to steep power laws of A/Q, where A and Q are the atomic mass and ionic charge. Since thermal electron energies are enhanced and nonthermal electron distributions arise in the reconnection process, we might expect that ionic charge states Q would be increased through ionization interactions with those electron populations during the acceleration process. The temperature estimated from the SEPs corresponds to the charge state during the acceleration process, while the actual charge state measured in situ may be modified as the SEPs pass through the corona. We examine whether the temperature estimation from the A/Q would differ with various kappa values in a kappa function representing high-energy tail deviating from a Maxwellian velocity distribution. We find that the differences in the A/Q between a Maxwellian and an extreme kappa distribution are about 10-30. We fit power-law enhancement of element abundances as a function of their A/Q with various kappa values. Then, we find that the derived source region temperature is not significantly affected by whether or not the electron velocity distribution deviates from a Maxwellian, i.e., thermal, distribution. Assuming that electrons are heated in the acceleration region, the agreement of the SEP charge state during acceleration with typical active region temperatures suggests that SEPs are accelerated and leave the acceleration region in a shorter time than the ionization time scale.

I identify a point-symmetric morphology of the supernova remnant (SNR) G352.7-0.1 and propose that the outer axially-symmetric structure is the remnant of a common envelope evolution (CEE) of the progenitor system, while the inner structure is the ejecta of a thermonuclear explosion triggered by the merger of a white dwarf (WD) and the core of an asymptotic giant branch (AGB) star. The main radio structure of SNR G352.7-0.1 forms an outer (large) ellipse. The bright X-ray emitting gas forms a smaller ellipse with a symmetry axis inclined to the symmetry axis of the large radio ellipse. The high abundance of iron and the energy of its X-ray lines suggest a type Ia supernova (SN Ia). The massive swept-up gas suggests a relatively massive progenitor system. I propose a scenario with progenitors of initial masses of M1=5-7Mo and M2=4-5 Mo. At a later phase, the WD remnant of the primary star and the AGB secondary star experience a CEE that ejects the circumstellar material that swept up more ISM to form the large elliptical radio structure. An explosion during the merger of the WD with the core of the AGB star triggered a super-Chandrasekhar thermonuclear explosion that formed the inner structure that is bright in X-ray. A tertiary star in the system caused the misalignment of the two symmetry axes. This study adds to the rich variety of evolutionary routes within the different scenarios of normal and peculiar SNe Ia.

Star formation takes place in filamentary molecular clouds which arise by physical processes that take place in the cold, neutral medium (CNM). We address the necessary conditions for this diffuse ($n \approx 30$ cm$^{-3}$), cold (T $\approx$ 60 K), magnetized gas undergoing shock waves and supersonic turbulence, to produce filamentary structures capable of fragmenting into cluster forming regions. Using RAMSES and a magnetized CNM environment as our initial conditions, we simulate a 0.5 kpc turbulent box to model a uniform gas with magnetic field strength of 7 $\mu G$, varying the 3D velocity dispersion via decaying turbulence. We use a surface density of $320 M_{\odot} pc^{-2}$, representative of the inner 4.0 kpc CMZ of the Milky Way and typical luminous galaxies. Filamentary molecular clouds are formed dynamically via shocks within a narrow range of velocity dispersions in the CNM of 5 - 10 km/s with a preferred value at 8 km/s. Cluster sink particles appear in filaments which exceed their critical line mass, occurring optimally for velocity dispersions of 8 km/s. Tracking the evolution of magnetic fields, we find that they lead to double the dense star forming gas than in purely hydro runs. Perpendicular orientations between magnetic field and filaments can increase the accretion rates onto filaments and hence their line masses. Because magnetic fields help support gas, MHD runs result in average temperatures an order of magnitude higher than unmagnetized counterparts. Finally, we find magnetic fields delay the onset of cluster formation by $\propto 0.4$ Myr.

Photodissociated gas bears the signature of the dynamical evolution of the ambient interstellar medium impacted by the mechanical and radiative feedback from an expanding HII region. Here we present an analysis of the kinematics of the young Trifid nebula, based on velocity-resolved observations of the far-infrared fine-structure lines of [C II] at 158 micron and [O I] at 63 micron. The distribution of the photodissociated regions (PDRs) surrounding the nebula is consistent with a shell-like structure created by the HII region expanding with a velocity of 5 km/s. Comparison of ratios of [C II] and [O I] 63 micron intensities for identical velocity components with PDR models indicate a density of 1e4 /cm^3. The red- and blue-shifted PDR shells with a combined mass of 516 Msun have a kinetic energy of ~1e47 erg. This is consistent with the thermal energy of the HII region as well as with the energy deposited by the stellar wind luminosity from HD 169442A, an O7 V star, over the 0.5 Myr lifetime of the star. The observed momentum of the PDR shell is lower than what theoretical calculations predict for the radial momentum due to the shell being swept up by an expanding HII region, which suggests that significant mass loss has occurred in M20 due to the dispersal of the surrounding gas by the advancing ionization front.

We use the VERITAS imaging air Cherenkov Telescope (IACT) array to obtain the first measured angular diameter of $\beta$ UMa at visual wavelengths using stellar intensity interferometry (SII) and independently constrain the limb-darkened angular diameter. The age of the Ursa Major moving group has been assessed from the ages of its members, including nuclear member Merak ($\beta$ UMa), an A1-type subgiant, by comparing effective temperature and luminosity constraints to model stellar evolution tracks. Previous interferometric limb-darkened angular-diameter measurements of $\beta$ UMa in the near-infrared (CHARA Array, $1.149 \pm 0.014$ mas) and mid-infrared (Keck Nuller, $1.08 \pm 0.07$ mas), together with the measured parallax and bolometric flux, have constrained the effective temperature. This paper presents current VERITAS-SII observation and analysis procedures to derive squared visibilities from correlation functions. We fit the resulting squared visibilities to find a limb-darkened angular diameter of $1.07 \pm 0.04 {\rm (stat)} \pm 0.05$ (sys) mas, using synthetic visibilities from a stellar atmosphere model that provides a good match to the spectrum of $\beta$ UMa in the optical wave band. The VERITAS-SII limb-darkened angular diameter yields an effective temperature of $9700\pm200\pm 200$ K, consistent with ultraviolet spectrophotometry, and an age of $390\pm 29 \pm 32 $ Myr, using MESA Isochrones and Stellar Tracks (MIST). This age is consistent with $408 \pm 6$ Myr from the CHARA Array angular diameter.

WISE J224607.6-052634.9 (W2246-0526) is a hot dust-obscured galaxy at $z$ = 4.601, and the most luminous obscured quasar known to date. W2246-0526 harbors a heavily obscured supermassive black hole that is most likely accreting above the Eddington limit. We present observations with the Atacama Large Millimeter/submillimeter Array (ALMA) in seven bands, including band 10, of the brightest far-infrared (FIR) fine-structure emission lines of this galaxy: [OI]$_{63\mu m}$, [OIII]$_{88\mu m}$, [NII]$_{122\mu m}$, [OI]$_{145\mu m}$, [CII]$_{158\mu m}$, [NII]$_{205\mu m}$, [CI]$_{370\mu m}$, and [CI]$_{609\mu m}$. A comparison of the data to a large grid of Cloudy radiative transfer models reveals that a high hydrogen density ($n_{H}\sim3\times10^3$ cm$^{-3}$) and extinction ($A_{V}\sim300$ mag), together with extreme ionization ($log(U)=-0.5$) and a high X-ray to UV ratio ($\alpha_{ox}\geq-0.8$) are required to reproduce the observed nuclear line ratios. The values of $\alpha_{ox}$ and $U$ are among the largest found in the literature and imply the existence of an X-ray-dominated region (XDR). In fact, this component explains the a priori very surprising non-detection of the [OIII]$_{88\mu m}$ emission line, which is actually suppressed, instead of boosted, in XDR environments. Interestingly, the best-fitted model implies higher X-ray emission and lower CO content than what is detected observationally, suggesting the presence of a molecular gas component that should be further obscuring the X-ray emission over larger spatial scales than the central region that is being modeled. These results highlight the need for multiline infrared observations to characterize the multiphase gas in high redshift quasars and, in particular, W2246-0526 serves as an extreme benchmark for comparisons of interstellar medium conditions with other quasar populations at cosmic noon and beyond.

Hot dust in the proximity of AGNs strongly emits in the Near Infrared producing a red excess that, in type 2 sources, can be modeled to measure its temperature. In the era of high spatial-resolution multi-wavelength data, mapping the hot dust around Supermassive Black Holes is important for the efforts to achieve a complete picture of the dust role and distribution around these compact objects. In this work we propose a methodology to detect the hot dust emission in the proximity of Type 2 AGNs and measure its temperature using K-band spectra ($\lambda_c$ = 2.2\,$\mu$m). To achieve this, we have developed NIRDust, a Python package for modeling K-band spectra, estimate the dust temperature and characterize the involved uncertainties. We tested synthetic and real spectra in order to check the performance and suitability of the physical model over different types of data. Our tests on synthetic spectra demonstrated that the obtained results are influenced by the signal-to-noise ratio (S/N) of the input spectra. However, we accurately characterized the uncertainties, which remained below $\sim$150 K for an average S/N per pixel exceeding 20. Applying NIRDust to NGC 5128 (Centaurus A), observed with the Gemini South Telescope, we estimated a dust temperature of 662 and 667 K from Flamingos-2 spectra and 697 and 607 K from GNIRS spectra using two different approaches.

Relativistic jets accompany the collapse of massive stars, the merger of compact objects, or the accretion of gas in active galactic nuclei. They carry information about the central engine and generate electromagnetic radiation. No self-consistent simulations have been able to follow these jets from their birth at the black hole scale to the Newtonian dissipation phase, making the inference of central engine property through astronomical observations undetermined. We present the general relativistic moving-mesh framework to achieve the continuity of jet simulations throughout space-time. We implement the general relativistic extension for the moving-mesh relativistic hydrodynamic code-JET, and develop a tetrad formulation to utilize the HLLC Riemann solver in the general relativistic moving mesh code. The new framework is able to trace the radial movement of the relativistic jets from the central region where strong gravity holds all the way to distances of jet dissipation.

Mergers of binary compact objects, accompanied with electromagnetic (EM) counterparts, offer excellent opportunities to explore varied cosmological models, since gravitational waves (GW) and EM counterparts always carry the information of luminosity distance and redshift, respectively. $f(T)$ gravity, which alters the background evolution and provides a friction term in the propagation of GW, can be tested by comparing the modified GW luminosity distance with the EM luminosity distance. Considering the third-generation gravitational-wave detectors, Einstein Telescope and two Cosmic Explorers, we simulate a series of GW events of binary neutron stars (BNS) and neutron-star-black-hole (NSBH) binary with EM counterparts. These simulations can be used to constrain $f(T)$ gravity (specially the Power-law model $f(T)=T+\alpha(-T)^\beta$ in this work) and other cosmological parameters, such as $\beta$ and Hubble constant. In addition, combining simulations with current observations of type Ia supernovae and baryon acoustic oscillations, we obtain tighter limitations for $f(T)$ gravity. We find that the estimated precision significantly improved when all three data sets are combined ($\Delta \beta \sim 0.03$), compared to analyzing the current observations alone ($\Delta \beta \sim 0.3$). Simultaneously, the uncertainty of the Hubble constant can be reduced to approximately $1\%$.

In this article we explore the holographic approach to neutron stars in the realm of Quantum Field Theory (QFT). We delve into the structures of neutron stars, emphasizing the application of the AdS/CFT duality in modeling them. We discuss both "bottom-up" and "top-down" holographic models, comparing their predictions with astrophysical observations. Finally, we demonstrate the potential broader applications of the holography method in areas like superconductivity, highlighting the methodological significance of string theory and QFT in astrophysics.

Rotating neutron stars that support long-lived, non-axisymmetric deformations known as mountains have long been considered potential sources of gravitational radiation. However, the amplitude from such a source is very weak and current gravitational-wave interferometers have yet to witness such a signal. The lack of detections has provided upper limits on the size of the involved deformations, which are continually being constrained. With expected improvements in detector sensitivities and analysis techniques, there is good reason to anticipate an observation in the future. This review concerns the current state of the theory of neutron-star mountains. These exotic objects host the extreme regimes of modern physics, which are related to how they sustain mountains. We summarise various mechanisms that may give rise to asymmetries, including crustal strains built up during the evolutionary history of the neutron star, the magnetic field distorting the star's shape and accretion episodes gradually constructing a mountain. Moving beyond the simple rotating model, we also discuss how precession affects the dynamics and modifies the gravitational-wave signal. We describe the prospects for detection and the challenges moving forward.

Possible forms of black hole images, viewed by a distant observer, are calculated basing on general relativity and equations of motion in the Kerr-Newman metric. Black hole image is a gravitationally lensed image of the black hole event horizon. It may be viewed as a black spot on the celestial sphere, projected inside the position of classical black hole shadow. In the nearest future it would be possible to verify modified gravity theories by observations of astrophysical black hole with Space Observatory Millimetron.

We resolve a long standing question regarding the suitable effective diffusion coefficient of the spherically-symmetric transport equation, which is valid at long times. To that end, we generalize a transport solution in three dimensions for homogeneous media, to include general collisional properties, including birth-death events and linearly anisotropic scattering. This is done by introducing an exact scaling law relating the Green function of the pure-scattering case with the general collision case, which is verified using deterministic and Monte-Carlo simulations. Importantly, the effective diffusion coefficient is identified by inspecting the transport solution at long times.

The geodesic method has played a crucial role in understanding the circular orbits generated by compact objects, culminating in the definition of the photon sphere, which was later generalized to a photon surface in arbitrary spacetimes. This new formulation extends the concept of the photon sphere in a broader sense, including dynamical spacetimes, as shown by the Vaidya solution. The photon surface essentially defines the null geodesics, which are originally tangent to the temporal surface, and keeps them confined to this surface. However, this formalism does not cover all classes of particles, and to overcome this limitation, a more comprehensive approach, denoted as the "massive particle surface", has been proposed that also accounts for charged massive particles. Indeed, the photon surface concept is recovered when the charge and mass of the particles are zero. In this work, we use these three formalisms to check the consistency of the results for the values of the radius of the photon sphere ($r_{ps}$) and the radius of the "innermost stable circular orbit" (ISCO) ($r_{\rm ISCO}$) for some gravitational models. In our results, the first model is described by conformal gravity, with the peculiarity that $g_{00}\neq-g_{11}^{-1}$. The second model, i.e. the Culetu solution, is developed by coupling General Relativity with nonlinear electrodynamics, which requires the consideration of the effective metric ($g_{\rm eff}^{\mu\nu}$) for geodesic approaches. Furthermore, we have also analysed the expressions for $r_{ps}$ and $r_{\rm ISCO}$ in a general static and spherically symmetric metric. Under these circumstances, we have found a discrepancy of $r_{ps}$ and $r_{\rm ISCO}$ obtained by the massive particle surface formalism as compared to the geodesic and photon surface formalisms.