Papers for Thursday, Nov 18 2021

Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the frontier of precision astrometry from evidence of Earth-mass habitable worlds around the nearest stars, to distant Milky Way objects, and out to the Local Group of galaxies. As we enter the era of the James Webb Space Telescope and the new ground-based, adaptive-optics-enabled giant telescopes, by obtaining these high precision measurements on key objects that Gaia could not reach, a mission that focuses on high precision astrometry science can consolidate our theoretical understanding of the local Universe, enable extrapolation of physical processes to remote redshifts, and derive a much more consistent picture of cosmological evolution and the likely fate of our cosmos. Already several missions have been proposed to address the science case of faint objects in motion using high precision astrometry missions: NEAT proposed for the ESA M3 opportunity, micro-NEAT for the S1 opportunity, and Theia for the M4 and M5 opportunities. Additional new mission configurations adapted with technological innovations could be envisioned to pursue accurate measurements of these extremely small motions. The goal of this White Paper is to address the fundamental science questions that are at stake when we focus on the motions of faint sky objects and to briefly review instrumentation and mission profiles.

Early-type galaxies (ETGs) possess total radial density profiles that are largely described by singular isothermal spheres, which can lead to non-Gaussian line-of-sight velocity dispersion (LOSVD) under anisotropic stellar orbits. However, recent observations of local ETGs in the MASSIVE Survey reveal outer kinematic structures at $1.5 R_{eff}$ that are inconsistent with fixed isothermal density profiles; the authors proposed varying density profiles as an explanation. We aim to verify this conjecture and understand the influence of stellar assembly on the formation of these non-canonical kinematic features through mock ETGs with well-studied density profiles in IllustrisTNG. We create 2D maps mimicking integral-field-unit (IFU) observations to extract projected stellar kinematic features for 207 ETGs in the TNG100-1 box with stellar mass $M_{\ast} > 10^{11} M_{\odot}$. The mock observations reproduce the key outer ($1.5R_{eff}$) kinematic structures in the MASSIVE ETGs, including the puzzling positive correlation between velocity dispersion profile outer slope $\gamma_{outer}$ and the kurtosis $h_4$'s gradient. We find that $h_4$ is uncorrelated with stellar orbital anisotropy beyond $R_{eff}$; instead we find that the variations in $\gamma_{outer}$ and outer $h_4$ (which dominates the $h_4$ gradient) are both driven by variations of the density profile at the outskirts across different ETGs. These findings corroborate the proposed conjecture and rule out velocity anisotropy as the dominant driver of non-Gaussian outer kinematic structure in ETGs. We also find that the outer kurtosis and anisotropy correlate with different stellar assembly components, with the former related to minor mergers or flyby interactions while the latter is mainly driven by major mergers, suggesting distinct stellar assembly origins that decorrelates these two properties.

Secular oscillations in multi-planet systems can drive chaotic evolution of an inner body through non-linear resonant perturbations. This "secular chaos" readily pushes the particle to an extreme eccentricity, triggering tidal interactions or collision with the central star. We present a numerical study of secular chaos in systems with two planets and a test particle using the ring-averaging method, with emphasis on the relationship between the planets' properties and the time-scale of chaotic diffusion. We find that secular chaos can excite extreme eccentricities on time-scales spanning several orders of magnitude in a given system. We apply our results to the evolution of planetary systems around white dwarfs (WDs), specifically the tidal disruption and high-eccentricity migration of planetesimals and planets. We find that secular chaos driven by large ($\gtrsim 10 M_{\oplus}$), distant ($\gtrsim 10 \, {\rm au}$) planets can sustain metal accretion onto the WD from tidal disruption of planetesimals over Gyr time-scales. The accretion rate is consistent with that inferred for WDs with atmospheric pollution along the cooling sequence if the total mass of planetesimals in the chaotic zone at the beginning of the WD stage is similar to that of the Solar System's main asteroid belt. Based on the occurrence of long-period exoplanets and exo-asteroid belts, we conclude that secular chaos can be a significant (perhaps dominant) channel for polluting WDs. Secular chaos can also produce short-period planets and planetesimals around WDs in concert with various circularization mechanisms. We discuss prospects for detecting exoplanets driving secular chaos using direct imaging and microlensing.

We present LEGWORK (LISA Evolution and Gravitational Wave Orbit Kit), an open-source Python package for making predictions about stellar-origin gravitational wave sources and their detectability in LISA or other space-based gravitational wave detectors. LEGWORK can be used to evolve the orbits of sources due to gravitational wave emission, calculate gravitational wave strains (using post-Newtonian approximations), compute signal-to-noise ratios and visualise the results. It can be applied to a variety of potential sources, including binaries consisting of white dwarfs, neutron stars and black holes. Although we focus on double compact objects, in principle LEGWORK can be used for any system with a user-specified orbital evolution, such as those affected by a third object or gas drag. We optimised the package to make it efficient for use in population studies which can contain tens-of-millions of sources. This paper describes the package and presents several potential use cases. We explain in detail the derivations of the expressions behind the package as well as identify and clarify some discrepancies currently present in the literature. We hope that LEGWORK will enable and accelerate future studies triggered by the rapidly growing interest in gravitational wave sources.

This paper aims to quantify how the lowest halo mass that can be detected with galaxy-galaxy strong gravitational lensing depends on the quality of the observations and the characteristics of the observed lens systems. Using simulated data, we measure the lowest detectable NFW mass at each location of the lens plane, in the form of detailed \emph{sensitivity maps}. In summary, we find that: (i) the lowest detectable mass $M_{\rm low}$ decreases linearly as the signal-to-noise ratio (SNR) increases and the sensitive area is larger when we decrease the noise; (ii) a moderate increase in angular resolution (0.07" vs 0.09") and pixel scale (0.01" vs 0.04") improves the sensitivity by on average 0.25 dex in halo mass, with more significant improvement around the most sensitive regions; (iii) the sensitivity to low-mass objects is largest for bright and complex lensed galaxies located inside the caustic curves and lensed into larger Einstein rings (i.e $r_{E}\geq1.0"$). We find that for the sensitive mock images considered in this work, the minimum mass that we can detect at the redshift of the lens lies between $1.5\times10^{8}$ and $3\times10^{9}M_{\odot}$. We derive analytic relations between $M_{\rm low}$, the SNR and resolution and discuss the impact of the lensing configuration and source structure. Our results start to fill the gap between approximate predictions and real data and demonstrate the challenging nature of calculating precise forecasts for gravitational imaging. In light of our findings, we discuss possible strategies for designing strong lensing surveys and the prospects for HST, Keck, ALMA, Euclid and other future observations.

We present the first data release of the MaNGA-ARO Survey of CO Targets (MASCOT), an ESO Public Spectroscopic Survey conducted at the Arizona Radio Observatory (ARO). We measure the CO(1-0) line emission in a sample of 187 nearby galaxies selected from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey that has obtained integral field unit (IFU) spectroscopy for a sample of ~ 10,000 galaxies at low redshift. The main goal of MASCOT is to probe the molecular gas content of star-forming galaxies with stellar masses > 10^9.5 M_solar and with associated MaNGA IFU observations and well-constrained quantities like stellar masses, star formation rates and metallicities. In this paper we present the first results of the MASCOT survey, providing integrated CO(1-0) measurements that cover several effective radii of the galaxy and present CO luminosities, CO kinematics, and estimated H2 gas masses. We observe that the decline of galaxy star formation rate with respect to the star formation main sequence (SFMS) increases with the decrease of molecular gas and with a reduced star formation efficiency, in agreement with results of other integrated studies. Relating the molecular gas mass fractions with the slope of the stellar age gradients inferred from the MaNGA observations, we find that galaxies with lower molecular gas mass fractions tend to show older stellar populations close to the galactic center, while the opposite is true for galaxies with higher molecular gas mass fractions, providing tentative evidence for inside-out quenching.

We investigate the presence of a central black hole (BH) in B023-G078, M31's most massive globular cluster. We present high-resolution, adaptive-optics assisted, integral-field spectroscopic kinematics from Gemini/NIFS that shows a strong rotation ($\sim$20 km/s) and a velocity dispersion rise towards the center (37 km/s). We combine the kinematic data with a mass model based on a two-component fit to $HST$ ACS/HRC data of the cluster to estimate the mass of a putative BH. Our dynamical modeling suggests a $>$3$\sigma$ detection of a BH component of 9.1$^{+2.6}_{-2.8}\times$10$^4$ M$_\odot$ (1$\sigma$ uncertainties). The inferred stellar mass of the cluster is 6.22$^{+0.03}_{-0.05}\times$10$^6$ M$_\odot$, consistent with previous estimates, thus the BH makes up 1.5% of its mass. We examine whether the observed kinematics are caused by a collection of stellar mass BHs by modeling an extended dark mass as a Plummer profile. The upper limit on the size scale of the extended mass is 0.56 pc (95% confidence), which does not rule out an extended mass. There is compelling evidence that B023-G078 is the tidally stripped nucleus of a galaxy with a stellar mass $>$10$^9$ M$_{\odot}$, including its high mass, two-component luminosity profile, color, metallicity gradient, and spread in metallicity. Given the emerging evidence that the central BH occupation fraction of $>$10$^9$ M$_{\odot}$ galaxies is high, the most plausible interpretation of the kinematic data is that B023-G078 hosts a central BH. This makes it the strongest BH detection in a lower mass ($<$10$^7$ M$_{\odot}$) stripped nucleus, and one of the few dynamically detected intermediate-mass BHs.

Galight is a Python-based open-source package that can be used to perform two-dimensional model fitting of optical and near-infrared images to characterize the light distribution of galaxies with components including a disk, bulge, bar, and quasar. The decomposition of stellar components has been demonstrated in published studies of inactive galaxies and quasar host galaxies observed by the Hubble Space Telescope and Subaru's Hyper Suprime-Cam. Galight utilizes the image modeling capabilities of lenstronomy while redesigning the user interface for the analysis of large samples of extragalactic sources. The package is user-friendly with some automatic features such as determining the cutout size of the modeling frame, searching for PSF-stars in field-of-view, estimating the noise map of the data, identifying all the objects to set the initial model, and associated parameters to fit them simultaneously. These features minimize the manpower and allow the automatic fitting tasks. The software is distributed under the MIT license. The source code, installation guidelines, and example notebooks code can be found at https://galight.readthedocs.io/en/latest/

The origin of hydrodynamical instability and turbulence in the Keplerian accretion disk is a long-standing puzzle. The flow therein is linearly stable. Here we explore the evolution of perturbation in this flow in the presence of an additional force. Such a force, which is expected to be stochastic in nature hence behaving as noise, could result from thermal fluctuations (however small be), grain-fluid interactions, feedback from outflows in astrophysical disks, etc. We essentially establish the evolution of nonlinear perturbation in the presence of Coriolis and external forces, which is the modified Landau equation. We obtain that even in the linear regime, under suitable forcing and Reynolds number, the otherwise least stable perturbation evolves to a very large saturated amplitude, leading to nonlinearity and plausible turbulence. Hence, forcing essentially leads a linear stable mode to unstable. We further show that nonlinear perturbation diverges at a shorter time-scale in the presence of force, leading to a fast transition to turbulence. Interestingly, the emergence of nonlinearity depends only on the force but not on the initial amplitude of perturbation, unlike the original Landau equation-based solution.

We present a novel method to infer the Dark Matter (DM) content and spatial distribution within galaxies, based on convolutional neural networks trained within state-of-the-art hydrodynamical simulations (Illustris TNG100). The framework we have developed is capable of inferring the DM mass distribution within galaxies of mass $~10^{11}-10^{13}M_{\odot}$ with very high performance from the gravitationally baryon dominated internal regions to the DM-rich, baryon-depleted outskirts of the galaxies. With respect to traditional methods, the one presented here also possesses the advantages of not relying on a pre-assigned shape for the DM distribution, to be applicable to galaxies not necessarily in isolation, and to perform very well even in the absence of spectroscopic observations

We present a study of optically-selected Type II AGN at 0.5 < z < 0.9 from the VIPERS and VVDS surveys, to investigate the connection between AGN activity and physical properties of their host galaxies. The host stellar mass is estimated through spectral energy distribution fitting with the CIGALE code, and star formation rates are derived from the [OII]$\lambda$3727 $\r{A}$ line luminosity. We find that 49% of the AGN host galaxies are on or above the main sequence (MS), 40% lie in the sub-MS locus, and 11% in the quiescent locus. Using the [OIII]$\lambda$5007 $\r{A}$ line luminosity as a proxy of the AGN power, we find that at fixed AGN power Type II AGN host galaxies show a bimodal behaviour: systems with host galaxy stellar mass <10$^{10}$ M$_{\odot}$, reside along the MS or in the starbursts locus (high-SF Type II AGN), while systems residing in massive host-galaxies (>10$^{10}$ M$_{\odot}$) show a lower level of star formation (low-SF Type II AGN). At all stellar masses, the offset from the MS is positively correlated with the AGN power. We interpret this correlation as evidence of co-evolution between the AGN and the host, possibly due to the availability of cold gas. In the most powerful AGN with host galaxies below the MS we find a hint, though weak, of asymmetry in the [OIII] line profile, likely due to outflowing gas, consistent with a scenario in which AGN feedback removes the available gas and halts the star formation in the most massive hosts.

We use Gaia EDR3 data to identify stars associated with six classical dwarf spheroidals (Draco, Ursa Minor, Sextans, Sculptor, Fornax, Carina) at their outermost radii, beyond their nominal King stellar limiting radius. For all of the dSphs examined, we find radial velocity matches with stars residing beyond the King limiting radius and with $> 50\%$ astrometric probability (four in Draco, two in Ursa Minor, eight in Sextans, two in Sculptor, twelve in Fornax, and five in Carina), indicating that these stars are associated with their respective dwarf spheroidals (dSphs) at high probability. We compare the positions of our candidate "extra-tidal" stars with the orbital tracks of the galaxies, and identify stars, both with and without radial velocity matches, that are consistent with lying along the orbital track of the satellites. Cross matching with publicly available catalogs of RR Lyrae, we find one RR Lyrae candidate with $> 50\%$ astrometric probability outside the limiting radius in each of Sculptor and Fornax, two such candidates in Draco, nine in Ursa Minor, seven in Sextans, and zero in Carina. Follow-up spectra on all of our candidates, including possible metallicity information, will help confirm association with their respective dSphs, and could represent evidence for extended stellar halos or tidal debris around these classical dSphs.

We present \textsc{Magrathea-Pathfinder}, a ray-tracing framework which accurately reconstructs the past light-cone of an observer in numerical simulations. Our code directly computes the 3D trajectory of light rays through the null geodesic equations, with the weak-field limit as its only approximation. Therefore, it does not rely on any other standard ray-tracing approximations such as plane-parallel, Born or multiple-lens. \textsc{Magrathea-Pathfinder} fully takes advantage of the small-scale clustering of matter by using adaptive integration steps and interpolation within an Adaptive-Mesh Refinement (AMR) structure to accurately account for the non-linear regime of structure formation. It uses MPI parallelization, C\texttt{++}11 \texttt{std::thread} multithreading, and is optimised for High-Performance Computing (HPC) as a post-processing tool for very large $N$-body simulations. In this paper, we describe how to produce realistic cosmological observables from numerical simulation using ray-tracing techniques, in particular the production of simulated catalogues and maps which accounts for all the effects at first order in metric perturbations (such as peculiar velocities, gravitational potential, Integrated Sachs-Wolfe, time delay, gravitational lensing, etc\ldots). We perform convergence tests of our gravitational lensing algorithms and conduct performance tests of the null geodesic integration procedures. \textsc{Magrathea-Pathfinder} provides sophisticated ray-tracing tools to make the link between real space ($N$-body simulations) and light-cone observables. This should be useful to refine existing cosmological probes and to build new ones beyond standard assumptions in order to prepare for next-generation large-scale structure surveys.

General Relativistic effects on the clustering of matter in the universe provide a sensitive probe of cosmology and gravity theories that can be tested with the upcoming generation of galaxy surveys. Here, we present a suite of large volume high-resolution N-body simulations specifically designed to generate light-cone data for the study of relativistic effects on lensing-matter observables. RayGalGroupSims (or in short RayGal) consists of two N-body simulations of $(2625\,h^{-1}\,{\rm Mpc})^3$ volume with $4096^3$ particles of a standard flat $\Lambda$CDM model and a non-standard $w$CDM phantom dark energy model. Light-cone data from the simulations have been generated using a parallel ray-tracing algorithm that has accurately solved billion geodesic equations. Catalogues and maps with relativistic weak-lensing which include post-Born effects, magnification bias (MB) and redshift space distortions (RSD) due to gravitational redshift, Doppler, transverse Doppler, Integrated Sachs-Wolfe/Rees-Sciama effects, are publicly released. Using this dataset, we are able to reproduce the linear and quasi-linear predictions from the Class relativistic code for the 10 (cross-)power spectra (3$\times$2 points) of the matter density fluctuation field and the gravitational convergence at $z=0.7$ and $z=1.8$. We find 1-30\% level contribution from both MB and RSD to the matter power spectrum, while the Fingers-of-God effect is visible at lower redshift in the non-linear regime. MB contributes at the $10-30\%$ level to the convergence power spectrum leading to a deviation between the shear power-spectrum and the convergence power-spectrum. MB also plays a significant role in the galaxy-galaxy lensing by decreasing the density-convergence spectra by $20\%$, while coupling non-trivial configurations (such as the one with the convergence at the same or even lower redshift than the density field).

We describe a pipeline to measure scintillation in fast radio bursts (FRBs) detected by CHIME/FRB in the 400-800 MHz band by analyzing the frequency structure of the FRB's spectrum. We use the pipeline to measure the characteristic frequency bandwidths of scintillation between $4-100$ kHz in 12 FRBs corresponding to timescales of $\sim$2-40 $\mu$s for 10 FRBs detected by CHIME/FRB. For the other two FRBs, we did not detect scintillation in the region our analysis is sensitive. We compared the measured scintillation timescales to the NE2001 predictions for the scintillation timescales from the Milky Way. We find a strong correlation to be an indication that in most instances, the observed scintillation of FRBs can be explained by the Milky Way.

Most of the currently known planets are small worlds with radii between that of the Earth and that of Neptune. The characterization of planets in this regime shows a large diversity in compositions and system architectures, with distributions hinting at a multitude of formation and evolution scenarios. Using photometry from the K2 satellite and radial velocities measured with the HARPS and CORALIE spectrographs, we searched for planets around the bright and slightly evolved Sun-like star HD 137496. We precisely estimated the stellar parameters, $M_*$ = 1.035 +/- 0.022 $M_\odot$, $R_*$ = 1.587 +/- 0.028 $R_\odot$, $T_\text{eff}$ = 5799 +/- 61 K, together with the chemical composition of the slightly evolved star. We detect two planets orbiting HD 137496. The inner planet, HD 137496 b, is a super-Mercury (an Earth-sized planet with the density of Mercury) with a mass of $M_b$ = 4.04 +/- 0.55 $M_\oplus$, a radius of $R_b = 1.31_{-0.05}^{+0.06} R_\oplus,$ and a density of $\rho_b = 10.49_{-1.82}^{+2.08}$ $\mathrm{g cm^{-3}}$. With an interior modeling analysis, we find that the planet is composed mainly of iron, with the core representing over 70% of the planet's mass ($M_{core}/M_{total} = 0.73^{+0.11}_{-0.12}$). The outer planet, HD 137496 c, is an eccentric ($e$ = 0.477 +/- 0.004), long period ($P$ = $479.9_{-1.1}^{+1.0}$ days) giant planet ($M_c\sin i_c$ = 7.66 +/- 0.11 $M_{Jup}$) for which we do not detect a transit. HD 137496 b is one of the few super-Mercuries detected to date. The accurate characterization reported here enhances its role as a key target to better understand the formation and evolution of planetary systems. The detection of an eccentric long period giant companion also reinforces the link between the presence of small transiting inner planets and long period gas giants.

Spheromak type flux ropes are increasingly used for modelling coronal mass ejections (CMEs). Many models aim in accurately reconstructing the magnetic field topology of CMEs, considering its importance in assessing their impact on modern technology and human activities in space and on ground. However, so far there is little discussion about how the details of the magnetic structure of a spheromak affect its evolution through the ambient field in the modelling domain, and what impact this has on the accuracy of magnetic field topology predictions. If the spheromak has its axis of symmetry (geometric axis) at an angle with respect to the direction of the ambient field, then the spheromak starts rotating so that its symmetry axis finally aligns with the ambient field. When using the spheromak in space weather forecasting models this tilting can happen already during insertion and significantly affects the results. In this paper we highlight this issue previously not examined in the field of space weather and we estimate the angle by which the spheromak rotates under different conditions. To do this we generated simple purely radial ambient magnetic field topologies (weak/strong positive/negative) and inserted spheromaks with varying initial speed and tilt, and magnetic helicity sign. We employ different physical and geometric criteria to locate the magnetic centre of mass and axis of symmetry of the spheromak. We confirm that spheromaks rotate in all investigated conditions and their direction and angle of rotation depend on the spheromak's initial properties and ambient magnetic field strength and orientation.

(Abridged) We present numerical simulations of the formation of the planetary companions to 47 UMa, rho CrB, and 51 Peg. They are assumed to have formed in situ according to the basic model that a core formed first by accretion of solid particles, then later it captured substantial amounts of gas from the protoplanetary disk. In most of the calculations we prescribe a constant accretion rate for the solid core. The evolution of the gaseous envelope assumes that: (1) it is in quasi-hydrostatic equilibrium, (2) the gas accretion rate is determined by the requirement that the outer radius of the planet is the place at which the thermal velocity of the gas allows it to reach the boundary of the planet's Hill sphere, (3) the gas accretion rate is limited, moreover, by the prescribed maximum rate at which the nebula can supply the gas, and (4) the growth of the planet stops once it obtains approximately the minimum mass determined from radial velocity measurements. Calculations are carried out through an initial phase during which solid accretion dominates, past the point of crossover when the masses of solid and gaseous material are equal, through the phase of rapid gas accretion, and into the final phase of contraction and cooling at constant mass. Alternative calculations are presented for the case of 47 UMa in which the solid accretion rate is calculated, not assumed, and the dissolution of planetesimals within the gaseous envelope is considered. In all cases there is a short phase of high luminosity (1e-3-1e-2 Lsun) associated with rapid gas accretion. The height and duration of this peak depend on uncertain model parameters. The conclusion is reached that in situ formation of all of these companions is possible under some conditions. However, it is more likely that orbital migration was an important component of the evolution, at least for the planets around rho CrB and 51 Peg.

Polluted white dwarfs that have accreted planetary material provide a unique opportunity to probe the geology of exoplanetary systems. However, the nature of the bodies which pollute white dwarfs is not well understood: are they small asteroids, minor planets, or even terrestrial planets? We present a novel method to infer pollutant masses from detections of Ni, Cr and Si. During core--mantle differentiation, these elements exhibit variable preference for metal and silicate at different pressures (i.e., object masses), affecting their abundances in the core and mantle. We model core--mantle differentiation self-consistently using data from metal--silicate partitioning experiments. We place statistical constraints on the differentiation pressures, and hence masses, of bodies which pollute white dwarfs by incorporating this calculation into a Bayesian framework. We show that Ni observations are best suited to constraining pressure when pollution is mantle-like, while Cr and Si are better for core-like pollution. We find 3 systems (WD0449-259, WD1350-162 and WD2105-820) whose abundances are best explained by the accretion of fragments of small parent bodies ($<0.2M_\oplus$). For 2 systems (GD61 and WD0446-255), the best model suggests the accretion of fragments of Earth-sized bodies, although the observed abundances remain consistent ($<3\sigma$) with the accretion of undifferentiated material. This suggests that polluted white dwarfs potentially accrete planetary bodies of a range of masses. However, our results are subject to inevitable degeneracies and limitations given current data. To constrain pressure more confidently, we require serendipitous observation of (nearly) pure core and/or mantle material.

Professional astronomy is historically not an environment of diverse identities. In recognizing that public outreach efforts affect career outcomes for young people, it is important to assess the demographics of those being reached and continually consider strategies for successfully engaging underrepresented groups. One such outreach event, the International Astronomical Youth Camp (IAYC), has a 50-year history and has reached ~1700 participants from around the world. We find that the IAYC is doing well in terms of gender (59% female, 4.7% non-binary at the most recent camp) and LGBT+ representation, whereas black and ethnic minorities are lacking. In this proceeding, we report the current landscape of demographics applying to and attending the IAYC; the efforts we are making to increase diversity amongst participants; the challenges we face; and our future plans to bridge these gaps, not only for the benefit of the camp but for society overall.

Based on recent Hubble Space Telescope (HST) observations in far-UV and groundbased observations in optical bands, Pavlov and colleagues have revealed a thermal component in the spectrum of the old pulsar B0950+08 (spin-down age 17.5 Myr) and estimated a neutron star (NS) surface temperature of $(1$--$3)\times 10^5$ K. Our new HST observations in the optical have allowed us to resolve the pulsar from a close-by galaxy and measure the optical fluxes more accurately. Using the newly measured fluxes and a new calibration of the HST's far-UV detector, we fit the optical-UV pulsar's spectrum with a model that consists of a nonthermal power-law ($f_\nu\propto \nu^\alpha$) and a thermal blackbody components. We obtained the spectral slope $\alpha=-0.3\pm 0.3$, considerably flatter than found from groundbased observations, and the best-fit temperature in the range of $(6$--$12)\times 10^4$ K (as seen by a distant observer), depending on interstellar extinction and NS radius. The temperature is lower than reported previously, but still much higher than predicted by NS passive cooling scenarios for such an old pulsar. This means that some heating mechanisms operate in NSs, e.g., caused by interaction of the faster rotating neutron superfluid with the slower rotating normal matter in the inner crust of the NS.

We report the results of a visual inspection of images of the Rapid ASKAP Continuum Survey (RACS) in search of extended radio galaxies (ERG) that reach or exceed linear sizes on the order of one Megaparsec. We searched a contiguous area of 1059deg$^2$ from RA$_{\rm J}$=20$^h$20$^m$ to 06$^h$20$^m$, and $-50^{\circ}<\rm{Dec}_J<-40^{\circ}$, which is covered by deep multi-band optical images of the Dark Energy Survey (DES), and in which previously only three ERGs larger than 1Mpc had been reported. For over 1800 radio galaxy candidates inspected, our search in optical and infrared images resulted in hosts for 1440 ERG, for which spectroscopic and photometric redshifts from various references were used to convert their largest angular size (LAS) to projected linear size (LLS). This resulted in 178 newly discovered giant radio sources (GRS) with LLS$>$1Mpc, of which 18 exceed 2Mpc and the largest one is 3.4Mpc. Their redshifts range from 0.02 to $\sim$2.0, but only 10 of the 178 new GRS have spectroscopic redshifts. For the 146 host galaxies the median $r$-band magnitude and redshift are 20.9 and 0.64, while for the 32 quasars or candidates these are 19.7 and 0.75. Merging the six most recent large compilations of GRS results in 458 GRS larger than 1Mpc, so we were able to increase this number by $\sim39\%$ to now 636.

Observations of molecular clouds, prestellar cores, and protoplanetary disks have established that the HNC/HCN ratio may be a potent diagnostic of molecular gas physical conditions. The processes that govern the relative abundances of these molecules nevertheless remain poorly understood. We seek to exploit the wide range of UV irradiation strengths within the 1 pc diameter Helix planetary nebula to explore the potential role of UV radiation in driving HNC/HCN. We performed IRAM 30 m and APEX 12 m radio line observations across six positions within the Helix Nebula, making use of radiative transfer and photodissociation modeling codes to interpret the results for line intensities and line ratios in terms of the molecular gas properties. We have obtained the first detections of the plasma-embedded Helix molecular knots (globules) in HCN, HNC, HCO+, and other trace molecules. Analysis of the HNC/HCN integrated line intensity ratio reveals an increase with radial distance from the Helix central star. In the context of molecular line ratios of other planetary nebulae from the literature, the HNC/HCN ratio appears to be anticorrelated with UV emission over four orders of magnitude in incident flux. Models of the photodissociation regions within the Helix using the RADEX and Meudon codes reveal strong constraints on column density of the molecular gas, as well as pressure and temperature. Analysis of the molecular ion HCO+ across the Helix indicates that X-ray irradiation is likely driving HCO+ production in the outer regions of planetary nebulae, where photodissociation is limited, yet cold gas and ionized molecules are abundant. Although the observational results clearly indicate that UV irradiation is important in determining the HNC/HCN ratio, our PDR modeling indicates that the UV flux gradient alone cannot reproduce the observed variation of HNC/HCN across the Helix.

Compact Symmetric Objects (CSOs), young jetted-AGN of overall projected size <1 kpc, are of great interest due to their youth and evolution. The classification was introduced to distinguish between ~95% of powerful compact extragalactic radio sources in flux density limited samples that are dominated by asymmetric emission due to relativistic beaming from jets aligned close to the line of sight, and ~5% of objects that are not. The original classification criteria were: (i) overall projected diameter smaller than ~1 kpc, (ii) identified center of activity, and (iii) symmetric jet structure about the center. There is confusion and erosion of the value of the CSO classification due to misclassifications. Many jets contain compact bright features outside core, resulting in a GPS total spectrum and a "compact double" appearance, and some objects with jet axes aligned close to the line of sight appear symmetric because the approaching jet is projected on both sides of the core. To eliminate the confusion, we propose adding (iv) slow radio variability and (v) low apparent velocity of bright features moving along the jets to the above CSO criteria. We are compiling a catalog of CSOs using these five criteria to eliminate the confusion of Doppler boosting.

J01020100-7122208 is a star whose origin and nature still challenges us. It was first believed to be a yellow super giant ejected from the Small Magellanic Cloud, but it was more recently claimed to be a red giant accelerated by the Milky Way's central black hole. In order to unveil its nature, we analysed photometric, astrometric and high resolution spectroscopic observations to estimate the orbit, age, and 16 elemental abundances. Our results show that this star has a retrograde and highly-eccentric orbit, $e=0.914_{-0.020}^{+0.016}$. Correspondingly, it likely crossed the Galactic disk at $550\;\mathrm{pc}$ from the Galactic centre. We obtained a spectroscopic mass and age of $1.09\pm0.10$ $M_\odot$ and $4.51\pm1.44$ Gyr respectively. Its chemical composition is similar to the abundance of other retrograde halo stars. We found that the star is enriched in europium, having [Eu/Fe] = 0.93 $\pm$ 0.24, and is more metal-poor than reported in the literature, with [Fe/H] = -1.30 $\pm$ 0.10. This information was used to conclude that J01020100-7122208 is likely not a star ejected from the central black of the Milky Way or from the Small Magellanic Cloud. Instead, we propose that it is simply a halo star which was likely accreted by the Milky Way in the distant past but its mass and age suggest it is probably an evolved blue straggler.

CHEOPS(CHaracterising ExOPlanet Satellite) is an ESA S-class mission that observes bright stars at high cadence from low-Earth orbit. The main aim of the mission is to characterize exoplanets that transit nearby stars using ultrahigh precision photometry. Here we report the analysis of transits observed by CHEOPS during its Early Science observing programme for four well-known exoplanets: GJ436b, HD106315b, HD97658b and GJ1132b. The analysis is done using pycheops, an open-source software package we have developed to easily and efficiently analyse CHEOPS light curve data using state-of-the-art techniques that are fully described herein. We show that the precision of the transit parameters measured using CHEOPS is comparable to that from larger space telescopes such as Spitzer Space Telescope and Kepler. We use the updated planet parameters from our analysis to derive new constraints on the internal structure of these four exoplanets.

We compare the evolution of Mercury's eccentricity in published Solar System integrations. These data sets are affected to different degrees by numerical chaos, because of how well they resolve Mercury's pericenter passage. We find statistically significant differences between two data sets, in the eccentricity distribution of Mercury over time. We also study pericenter resolution in a variety of symplectic maps used in the literature.

We explore the structure around shell-crossing time of cold dark matter protohaloes seeded by two or three crossed sine waves of various relative initial amplitudes, by comparing Lagrangian perturbation theory (LPT) up to 10th order to high-resolution cosmological simulations performed with the public Vlasov code ColDICE. Accurate analyses of the density, the velocity, and related quantities such as the vorticity are performed by exploiting the fact that ColDICE can follow locally the phase-space sheet at the quadratic level. To test LPT predictions beyond shell-crossing, we employ a ballistic approximation, which assumes that the velocity field is frozen just after shell-crossing. In the generic case, where the amplitudes of the sine waves are all different, high-order LPT predictions match very well the exact solution, even beyond collapse. As expected, convergence slows down when going from quasi-1D dynamics where one wave dominates over the two others, to the axial-symmetric configuration, where all the amplitudes of the waves are equal. It is also noticed that LPT convergence is slower when considering velocity related quantities. Additionally, the structure of the system at and beyond collapse given by LPT and the simulations agrees very well with singularity theory predictions, in particular with respect to the caustic and vorticity patterns that develop beyond collapse. Again, this does not apply to axial-symmetric configurations, that are still correct from the qualitative point of view, but where multiple foldings of the phase-space sheet produce very high density contrasts, hence a strong backreaction of the gravitational force.

We test the idea that bar pattern speeds decrease with time owing to angular momentum exchange with a dark matter halo. If this process actually occurs, then the locations of the corotation resonance and other resonances should generally increase with time. We therefore derive the angular velocity $\Omega$ and epicyclic frequency $\kappa$ as functions of galactocentric radius for 85 barred galaxies using photometric data. Mass maps are constructed by assuming a dynamical mass-to-light ratio and then solving the Poisson equation for the gravitatonal potential. The location of Lindblad resonances and the corotation resonance radius are then derived using the standard precession frequency curves in conjunction with bar pattern speeds recently estimated from the Tremaine-Weinberg method as applied to Integral Field Spectroscopy (IFS) data. Correlations between physical properties of bars and their host galaxies indicate that bar {\it length} and the corotation radius depend on the disk circular velocity while bar {\it strength} and pattern speed do not. As the bar pattern speed decreases, bar strength, length, and corotation radius incease, but when bars are subclassified into fast, medium, and slow domains, no significant change in bar length is found. Only a hint of an increase of bar strength from fast to slow bars is found. These results suggest that bar length in galaxies undergoes little evolution, being instead determined mainly by the size of their host galaxy.

A giant planet embedded in a protoplanetary disk opens a gap by tidal interaction, and properties of the gap strongly depend on the planetary mass and disk parameters. Many numerical simulations of this process have been conducted, but detailed simulations and analysis of gap formation by a super-Jupiter mass planet have not been thoroughly conducted. We performed two-dimensional numerical hydrodynamic simulations of the gap formation process by a super-Jupiter mass planet and examined the eccentricity of the gap. When the planet is massive, the radial motion of gas is excited, causing the eccentricity of the gap's outer edge to increase. Our simulations showed that the critical planetary mass for the eccentric gap was $\sim3~M_{\rm J}$ in a disk with $\alpha=4.0\times10^{-3}$ and $h/r=0.05$, a finding that was consistent with that reported in a previous work. The critical planetary mass for the eccentric gap depends on the viscosity and the disk scale height. We found that the critical mass could be described by considering a dimensionless parameter related to the gap depth. The onset of gap eccentricity enhanced the surface density inside the gap, shallowing the gap more than the empirical relation derived in previous studies for a planet heavier than the critical mass. Therefore, our results suggest that the mass accretion rate, which strongly depends on the gas surface density in the gap is also enhanced for super-Jupiter mass planets. These results may substantially impact the formation and evolution processes of super-Jupiter mass planets and population synthesis calculations.

The dynamical state of massive clumps is key to our understanding of the formation of massive stars. In this work, we study the kinematic properties of massive clumps using synthetic observations. We have previously compiled a very large catalog of synthetic dust-continuum compact sources from our 250 pc, SN-driven, star formation simulation. Here, we compute synthetic $\rm N_{2}H^{+}$ line profiles for a subsample of those sources and compare their properties with the observations and with those of the corresponding three-dimensional (3D) clumps in the simulation. We find that the velocity dispersion of the sources estimated from the $\rm N_{2}H^{+}$ line is a good estimate of that of the 3D clumps, although its correlation with the source size is weaker than the velocity-size correlation of the 3D clumps. The relation between the mass of the 3D clumps, $M_{\rm main}$, and that of the corresponding synthetic sources, $M_{\rm SED}$, has a large scatter and a slope of 0.5, $M_{\rm main} \propto M_{\rm SED}^{0.5}$, due to uncertainties arising from the observational band-merging procedure and from projection effects along the line of sight. As a result, the virial parameters of the 3D clumps are not correlated with the clump masses, even if a negative correlation is found for the compact sources, and the virial parameter of the most massive sources may significantly underestimate that of the associated clumps.

The Mars Reconnaissance Orbiter's (MRO's) Mars Color Imager (MARCI) has returned approximately daily, approximately global images of Mars since late 2006, in up to seven different colors, from ultraviolet through near-infrared. To-date, that is over 5300 Mars days of data, nearly eight full Mars Years, or more than 15 Earth years. The data are taken at up to nearly 500 meters per pixel, and the nearly circular orbit of MRO and its consistent early afternoon imaging provide an unprecedented baseline of data with which to study Mars' atmosphere and surface processes. Unfortunately, processing the MARCI data is difficult, fraught with exploding file sizes, issues that require workarounds in free software, and other problems that make this a severely under-utilized dataset. This paper discusses a workflow to process MARCI data to their fullest, including suggestions on how to work around issues unique to MARCI and its interaction with the current version of the free software ISIS (Integrated Software for Imagers and Spectrometers), and discussion of some trades that can be made to dramatically speed data processing are also described. Examples of processed MARCI images, mosaics, and color composites are shown, demonstrating the abilities of this workflow on global, regional, and local areas at the full, 96 pixels per degree scale afforded by MARCI.

Aims. We report the detection of long-period radial velocity (RV) variations in three giant stars, HD 19615, HD 150010, and HD 174205, using precise RV measurements. Methods. These detections are part of the Search for Exoplanets around Northern Circumpolar Stars (SENS) survey being conducted at the Bohyunsan Optical Astronomy Observatory (BOAO). The nature of the RV variations was investigated by analyzing the photometric and line shape variations. We found no variability with the RV period in these quantities and conclude that the RV variations are most likely caused by planetary companions. Results. Orbital solutions for the three stars yield orbital periods of 402 d, 562 d, and 582 d and minimum masses of 8.5 MJ , 2.4 MJ , and 4.2 MJ , respectively. These masses and periods are typical for planets around intermediate-mass stars, although some unclear interpretations and recent studies may being calling some planet convictions into question. Nevertheless, the SENS program is contributing to our knowledge of giant planets around intermediate-mass stars.

Plumes in a convective flow, whose flow structure is localised in space and time, are considered to be relevant to the turbulent transport in convection. The effective mass, momentum, and heat transports in the convective turbulence are investigated in the framework of time--space double averaging procedure, where a field quantity is decomposed into three parts: the spatiotemporal mean (spatial average of the time-averaged) field, the dispersion or coherent fluctuation (deviation from the spatiotemporal mean), and the random or incoherent fluctuation. With this double-averaging framework, turbulent correlations such as the Reynolds stress, turbulent mass flux, turbulent internal-energy flux, etc., in the mean-field equations are divided into the dispersion/coherent correlation part and the random/incoherent correlation part. The evolution equations of these two parts of the correlation show what are responsible for the conversion of the fluctuation energy between the coherent and incoherent components. By reckoning the plume as the coherent fluctuation, a transport model for the convective turbulence is constructed with the aid of the non-equilibrium effect along plume motions, and applied to a stellar convective flow. One of the prominent characteristics of a surface cooling-driven convection, the enhanced and localised turbulent mass flux below the surface layer, which cannot be reproduced at all by the usual eddy-diffusivity model with mixing length theory (MLT), is well reproduced by the present model with the non-equilibrium effect. Our results show that the incorporation of plume motion into turbulent transport model through the non-equilibrium effect is an important and very relevant extension of mean-field theory beyond the heuristic gradient transport model with MLT.

We report millimeter-VLBI results of low-luminosity active galactic nuclei (M 84 and M 87) up to 88 GHz with source-frequency phase-referencing observations. We detected the weak VLBI core and obtained the first image of M 84 at 88 GHz. The derived brightness temperature of M 84 core was about 7.2$\times$10$^9$ K, which could serve as a lower limit as the core down to 30 Schwarzschild radii was still un-resolved in our 88 GHz observations. We successfully determined the core-shifts of M 87 at 22-44 GHz and 44-88 GHz through source-frequency phase-referencing technique. The jet apex of M 87 could be deduced at about 46 $\mu$as upstream of the 43 GHz core from core-shift measurements. The estimated magnetic field strength of the 88 GHz core of M 87 is 4.8$\pm$2.4 G, which is at the same magnitude of 1-30 G near the event horizon probed by the Event Horizon Telescope.

We reexamine high-magnification microlensing events in the previous data collected by the KMTNet survey with the aim of finding planetary signals that were not noticed before. In this work, we report the planetary system KMT-2018-BLG-1988L that was found from this investigation. The planetary signal appears as a deviation with $\lesssim 0.2$~mag from a single-lens light curve and lasted for about 6 hours. The deviation exhibits a pattern of a dip surrounded by weak bumps on both sides of the dip. The analysis of the lensing light curve indicates that the signal is produced by a low mass-ratio ($q\sim 4\times 10^{-5}$) planetary companion located near the Einstein ring of the host star. The mass of the planet, $M_{\rm planet}=6.8^{+4.7}_{-3.5}~M_\oplus$ and $5.6^{+3.8}_{-2.8}~M_\oplus$ for the two possible solutions, estimated from the Bayesian analysis indicates that the planet is in the regime of a super-Earth. The host of the planet is a disk star with a mass of $M_{\rm host} = 0.47^{+0.33}_{-0.25}~M_\odot$ and a distance of $D_{\rm L}= 4.2^{+1.8}_{-.14}$~kpc. KMT-2018-BLG-1988Lb is the seventeenth microlensing planet with a mass below the upper limit of a super-Earth. The fact that 14 out of 17 microlensing planets with masses $\lesssim 10~M_\oplus$ were detected during the last 5 years since the full operation of the KMTNet survey indicates that the KMTNet database is an important reservoir of very low-mass planets.

When observed in quiet regions close to the solar limb, many strong resonance lines show conspicuous linear polarization signals, produced by scattering processes, with extended wing lobes. Recent studies indicate that, contrary to what was previously believed, the wing lobes are sensitive to the presence of relatively weak longitudinal magnetic fields through magneto-optical (MO) effects. We theoretically investigate the sensitivity of the scattering polarization wings of the Ca I 4227 {\AA} line to the MO effects, and we explore its diagnostic potential for inferring information on the longitudinal component of the photospheric magnetic field. We calculate the intensity and polarization profiles of the Ca I 4227 {\AA} line by numerically solving the problem of the generation and transfer of polarized radiation under non-local thermodynamic equilibrium conditions in one-dimensional semi-empirical models of the solar atmosphere, taking into account the joint action of the Hanle, Zeeman, and MO effects. We consider volume-filling magnetic fields as well as magnetic fields occupying a fraction of the resolution element. In contrast to the circular polarization signals produced by the Zeeman effect, we find that the linear polarization angle in the scattering polarization wings of Ca I 4227 presents a clear sensitivity, through MO effects, not only to the flux of the photospheric magnetic field, but also to the fraction of the resolution element that the magnetic field occupies. We identify the linear polarization angle in the wings of strong resonance lines as a valuable observable for diagnosing unresolved magnetic fields. Used in combination with observables that encode information on the magnetic flux and other properties of the observed atmospheric region, it can provide constraints on the filling factor of the magnetic field.

Solar flares are occasionally responsible for severe Space Weather events, which can affect space-borne and ground-based infrastructures, endangering anthropic technological activities and even human health and safety. Thus, an essential activity in the framework of Space Weather monitoring is devoted to the observation of the activity level of the Sun. In this context, the acquisition system of the Catania Solar Telescope has been recently upgraded in order to improve its contribution to the European Space Agency (ESA) - Space Weather Service Network through the ESA Portal, which represents the main asset for Space Weather in Europe. Here, we describe the hardware and software upgrades of the Catania Solar Telescope and the main data products provided by this facility, which include full-disc images of the photosphere and chromosphere, together with a detailed characterization of the sunspot groups. As a showcase of the observational capabilities of the revamped Catania Solar Telescope, we report the analysis of a B5.4 class flare occurred on 2020 December 7, simultaneously observed by the IRIS and SDO satellites.

This work is devoted to the study of the broadband 0.8-79 keV spectral and timing properties of the poorly studied X-ray pulsar XTE J1859+083 during its 2015 outburst based on the data from the NuSTAR and Swift observatories. We show that the source pulse profile has complex shape that depends on the energy band. Pulse fraction of XTE J1859+083 has constant value around 35% in the broad energy band, this behaviour is atypical for X-ray pulsars. At the same time its energy spectrum is typical of this class of objects and has a power-law shape with an exponential cutoff at high energies. No cyclotron absorption line was discovered in the source spectrum. On the basis of indirect method and the absence of a cyclotron line, an estimation was made for the magnetic field strength as less than $5\times10^{11}$ G or belonging to the interval from $5\times10^{12}$ to $2.0^{+0.9}_{-1.2}\times10^{13}$ G. Data from the NOT and SALT telescopes as well as optical and IR sky surveys allowed us also to study the nature of its optical companion. We have proposed and studied new possible candidates for the optical companion of XTE J1859+083 and the most likely candidate was identified. The results of the optical and IR photometry and spectroscopy of these possible companions showed that the system is a Be X-ray binary, showing Br$\gamma$, He I and strong H$\alpha$ spectral lines.

BL Lacs are one subclass of blazars with highly energetic $\gamma$-ray emission, which is strongly boosted by a relativistic beaming effect. The latest catalogue of the 10 years of {\it Fermi}/LAT data Abdollahi et al. (2020) and the $\gamma$-ray Doppler factors in Pei et al. (2020) provide us with a large BL Lac sample to study their jet emission morphologies and intrinsic properties. In this paper, we collected a sample of 294 {\it Fermi} BL Lacs and probed the correlations between the $\gamma$-ray emissions and luminosity distances. Our analyses give following conclusions: (1) the observed $\gamma$-ray emissions are really boosted by the $\gamma$-ray Doppler factor, and the intrinsic $\gamma$-ray emissions are closely correlated with luminosity distances. (2) the morphology of jet emissions for HBLs may be continuous, while that for IBLs may be the case of a moving sphere in the $\gamma$-ray bands.

Solar jets are impulsive, collimated plasma ejections that are triggered by magnetic reconnection. They are observed for many decades in various temperatures and wavelengths, therefore their kinematic characteristics, such as velocity and recurrence, have been extensively studied.Nevertheless, the high spatial resolution of the Interface Region Imaging Spectrograph (IRIS) launched in 2013 allowed us to make a step forward in the understanding of the relationship between surges and hot jets. In this paper we report on several results of recent studies of jets observed by IRIS. Cool and hot plasma have been detected with ejections of cool blobs having a speed reaching 300 km/s during the impulsive phase of jet formation and slow velocity surges surrounding hot jets after the reconnection phase. Plasma characteristics of solar jets, such as the emission measure, temperature, and density have been quantified. A multi-layer atmosphere at the reconnection site based on observed IRIS spectra has been proposed. IRIS evidenced bidirectional flows at reconnection sites, and tilt along the spectra which were interpreted as the signature of twist in jets. The search of possible sites for reconnection could be achieved by the analysis of magnetic topology. Combining Solar Dynamics Observatory/Helioseismic Magnetic Imager (SDO/HMI) vector magnetograms and IRIS observations, it was found that reconnection site could be located at null points in the corona as well as in bald patch regions low in the photosphere. In one case study a magnetic sketch could explain the initiation of a jet starting in a bald patch transformed to a current sheet in a dynamical way, and the transfer of twist from a flux rope to the jet during the magnetic reconnection process.

Supermassive black holes dominate the gravitational potential in galactic nuclei. In these dense environments, stars follow nearly Keplerian orbits and see their orbital planes relax through the potential fluctuations generated by the stellar cluster itself. For typical astrophysical galactic nuclei, the most likely outcome of this vector resonant relaxation (VRR) is that the orbital planes of the most massive stars spontaneously self-align within a narrow disc. We present a maximum entropy method to systematically determine this long-term distribution of orientations and use it for a wide range of stellar orbital parameters and initial conditions. The heaviest stellar objects are found to live within a thin equatorial disk. The thickness of this disk depends on the stars' initial mass function, and on the geometry of the initial cluster. This work highlights a possible (indirect) novel method to constrain the distribution of intermediate mass black holes in galactic nuclei.

In this work, we study the causal relations between a localised energy release and a remote prominence oscillation, where the prominence has a realistic thread-like structure. We used an open source magnetohydrodynamic (MHD) code known as MPI-AMRVAC to create a multithreaded prominence body. We introduced an additional energy source from which a shock wave originates, thereby inducing prominence oscillation. We studied two cases with different source amplitudes to analyze its effect on the oscillations. Our results show that the frequently used pendulum model does not suffice to fully estimate the period of the prominence oscillation, in addition to showing that the influence of the source and the thread-like prominence structure needs to be taken into account. Repeated reflections and transmissions of the initial shock wave occur at the specific locations of multiple high-temperature and high-density gradients in the domain. This includes the left and right transition region (TR) located at the footpoints of the magnetic arcade, as well as the various transition regions between the prominence and the corona (PCTR). This results in numerous interferences of compressional waves. They contribute to the restoring forces of the oscillation, causing the period to deviate from the expected pendulum model, in addition to leading to differences in attributed damping or even growth in amplitude between the various threads. Along with the global longitudinal motion that result from the shock impact, small-scale transverse oscillations are also evident. Multiple high-frequency oscillations represent the propagation of magnetoacoustic waves. The damping we see is linked to the conversion of energy and its exchange with the surrounding corona. Our simulations demonstrate the exchange of energy between different threads and their different modes of oscillation.

In this work, we collected a sample of BL Lacs, FR I and FR II(G) radio galaxies with available core and extended emissions from published works to discuss the unified schemes and estimate the Doppler factor for BL Lacs. Wilcoxon rank-sum test and Kolmogorov-Smirnov test both suggest that the probabilities for the distribution of the extended luminosity of BL Lacs and that of FR I and FR II(G) radio galaxies to be from the same parent distribution are $p_{\rm{WRS}}=0.779$ and $p_{\rm{K-S}}=0.326$, suggesting they are unified. Based on this unified schemes, we propose to estimate the Doppler factors for BL Lacs. Comparing the Doppler factor estimated by the fitting/regression method with those for the common sources in the literatures, we found a good linear correlation for common sources.

Prominence threads are very long and thin flux tubes which are partially filled with cold plasma. Observations have shown that transverse oscillations are frequent in these solar structures. The observations are usually interpreted as the fundamental kink mode, while the detection of the first harmonic remains elusive. Here, we aim to study how the density inhomogeneity in the longitudinal and radial directions modify the periods and damping times of kink oscillations, and how this effect would be reflected in observations. We solve the ideal magnetohydrodynamics equations through two different methods: a) performing 3D numerical simulations, and b) solving a 2D generalised eigenvalue problem. We study the dependence of the periods, damping times and amplitudes of transverse kink oscillations on the ratio between the densities at the centre and at the ends of the tube, and on the average density. We apply forward modelling on our 3D simulations to compute synthetic H$\alpha$ profiles. We confirm that the ratio of the period of the fundamental oscillation mode to the period of the first harmonic increases as the ratio of the central density to the footpoint density is increased or as the averaged density of the tube is decreased. We find that the damping times due to resonant absorption decrease as the central to footpoint density ratio increases. Contrary to the case of longitudinally homogeneous tubes, we find that the damping time to period ratio also increases as the density ratio is increased or the average density is reduced. We present snapshots and time-distance diagrams of the emission in the H$\alpha$ line. The results presented here have implications for the field of prominence seismology. While the H$\alpha$ emission can be used to detect the fundamental mode, the first harmonic is barely detectable in H$\alpha$. This may explain the lack of detections of the first harmonic.

Solar-type protostars harbor highly deuterated complex organics. While this degree of deuteration may provide important clues in studying the formation of these species, spectroscopic information on multiply deuterated isotopologs is often insufficient. In particular, searches for triply deuterated methanol, CD$_3$OH, are hampered by the lack of intensity information from a spectroscopic model. The aim of the present study is to develop such a model of CD$_3$OH in low-lying torsional states that is sufficiently accurate to facilitate further searches for CD$_3$OH in space. We performed a new measurement campaign for CD$_3$OH involving three spectroscopic laboratories that covers the 34 GHz-1.1 THz and the 20-900 cm$^{-1}$ ranges. The analysis was performed using the torsion-rotation Hamiltonian model based on the rho-axis method. We determined a model that describes the ground and first excited torsional states of CD$_3$OH, up to quantum numbers $J \leqslant 55$ and $K_a \leqslant 23$, and we derived a line list for radio-astronomical observations. This list is accurate up to at least 1.1 THz and should be sufficient for all types of radio-astronomical searches for this methanol isotopolog. It was used to search for CD$_3$OH in data from the Protostellar Interferometric Line Survey of IRAS 16293-2422 using ALMA. CD$_3$OH is securely detected in the data, with a large number of clearly separated and well-reproduced lines. We detected lines belonging to the ground and the first excited torsional states. The derived abundance of CD$_3$OH relative to non-deuterated isotopolog confirm the significant enhancement of this multiply deuterated variant. This finding is in line with other observations of multiply deuterated complex organic molecules and may serve as an important constraint on their formation models.

We present ALMA observations of a small but statistically complete sample of twelve 250 micron selected galaxies at $z=0.35$ designed to measure their dust submillimeter continuum emission as well as their CO(1-0) and atomic carbon [CI](3P1-3P0) spectral lines. This is the first sample of galaxies with global measures of all three $H_2$-mass tracers and which show star formation rates (4-26 Msun yr$^{-1}$) and infra-red luminosities ($1-6\times10^{11}$ Lsun) typical of star forming galaxies in their era. We find a surprising diversity of morphology and kinematic structure; one-third of the sample have evidence for interaction with nearby smaller galaxies, several sources have disjoint dust and gas morphology. Moreover two galaxies have very high $L_{CI}/L_{CO}$ ratios for their global molecular gas reservoirs; if confirmed, such extreme intensity ratios in a sample of dust selected, massive star forming galaxies presents a challenge to our understanding of ISM. Finally, we use the emission of the three molecular gas tracers, to determine the carbon abundance, $X_{ci}$, and CO-$\rm{H_2}$ conversion $\alpha_{co}$ in our sample, using a weak prior that the gas-to-dust ratio is similar to that of the Milky Way for these massive and metal rich galaxies. Using a likelihood method which simultaneously uses all three gas tracer measurements, we find mean values and errors on the mean of $\alpha_{co}=3.0\pm0.5\,\rm{Msun\,(K\,kms^{-1}\,pc^2)^{-1}}$ and $X_{ci}=1.6\pm0.1\times 10^{-5}$ (or $\alpha_{ci}=18.8\,K kms^{-1}\,pc^2 (Msun)^{-1}$) and $\delta_{GDR}=128\pm16$ (or $\alpha_{850}=5.9\times10^{12}\,\rm{W\,Hz^{-1}\, Msun^{-1}}$), where our starting assumption is that these metal rich galaxies have an average gas-to-dust ratio similar to that of the Milky Way centered on $\delta_{GDR}=135$.

Gas rich dusty circumstellar discs observed around young stellar objects are believed to be the birthplace of planets and planetary systems. Recent observations revealed that large-scale horseshoe-like brightness asymmetries are present in dozens of transitional protoplanetary discs. Theoretical studies suggest that these brightness asymmetries bf could be caused by large-scale anticyclonic vortices triggered by the Rossby Wave Instability (RWI), which can be excited at the edges of the accretionally inactive region, the dead zone edge. Since vortices may play a key role in planet formation, investigating the conditions of the onset of RWI and the long-term evolution of vortices is inevitable. The aim of our work was to explore the effect of disc geometry (the vertical thickness of the disc), viscosity, the width of the transition region at the dead zone edge, and the disc mass on the onset, lifetime, strength and evolution of vortices formed in the disc. We performed a parametric study assuming different properties for the disc and the viscosity transition by running 1980 2D hydrodynamic simulations in the locally isothermal assumption with disc self-gravity included. Our results revealed that long-lived, large-scale vortex formation favours a shallow surface density slope and low- or moderate disc masses with Toomre $Q \lesssim 1/h$, where $h$ is the geometric aspect ratio of the disc. In general, in low viscosity models, stronger vortices form. However, rapid vortex decay and re-formation is more widespread in these discs.

All 280 of the statistically-complete Palomar sample of nearby (<120 Mpc) galaxies dec > 20 degrees have been observed at 1.5 GHz as part of the LeMMINGs e-MERLIN legacy survey. Here, we present Chandra X-ray observations of the nuclei of 213 of these galaxies, including a statistically-complete sub-set of 113 galaxies in the declination range 40 degrees to 65 degrees. We observed galaxies of all optical spectral types, including 'active' galaxies (e.g., LINERs and Seyferts) and 'inactive' galaxies like HII galaxies and absorption line galaxies (ALG). The X-ray flux limit of our survey is 1.65$\times$10$^{-14}$~erg s$^{-1}$ cm$^{-2}$ (0.3$-$10 keV). We detect X-ray emission coincident within 2-arcsec of the nucleus in 150/213 galaxies, including 13/14 Seyferts, 68/77 LINERs, 13/22 ALGs and 56/100 HII galaxies, but cannot completely rule out contamination from non-AGN processes in sources with nuclear luminosities <10$^{39}$ erg s$^{-1}$. We construct an X-ray Luminosity function (XLF) and find that the local galaxy XLF, when including all AGN types, can be represented as a single power-law of slope $-0.54 \pm 0.06$. The Eddington ratio of the Seyferts is usually 2-4 decades higher than that of the LINERs, ALGs and HII galaxies, which are mostly detected with Eddington ratios <10$^{-3}$. Using [O III] line measurements and BH masses from the literature, we show that LINERs, HII galaxies and ALGs follow similar correlations to low luminosities, suggesting that some 'inactive' galaxies may harbour AGN.

The absorption of light by molecules in the atmosphere of Earth is a complication for ground-based observations of astrophysical objects. Comprehensive information on various molecular species is required to correct for this so called telluric absorption. We present a neural network autoencoder approach for extracting a telluric transmission spectrum from a large set of high-precision observed solar spectra from the HARPS-N radial velocity spectrograph. We accomplish this by reducing the data into a compressed representation, which allows us to unveil the underlying solar spectrum and simultaneously uncover the different modes of variation in the observed spectra relating to the absorption of $\mathrm{H_2O}$ and $\mathrm{O_2}$ in the atmosphere of Earth. We demonstrate how the extracted components can be used to remove $\mathrm{H_2O}$ and $\mathrm{O_2}$ tellurics in a validation observation with similar accuracy and at less computational expense than a synthetic approach with molecfit.

Plasma lensing is the refraction of low-frequency electromagnetic rays due to cold free electrons in the universe. For sources at a cosmological distance, there is observational evidence of elongated, complex plasma structures along the line of sight requiring a multi-lens-plane description. To investigate the limits of single-plane plasma lensing, we set up a double-plane lens with a projected Gaussian electron density in each lens plane. We compare double-plane scenarios with corresponding effective single-plane configurations. Our results show how double-plane lenses can be distinguished from single-plane lenses by observables, i.e. resolved multiple image positions, relative magnifications, time delays, and pulse shapes. For plasma lensing of fast radio bursts, the observed pulse shape may be dominated by the lensing effect, allowing us to neglect the intrinsic source pulse shape to distinguish different lensing configurations. The time-domain observables turn out to be the most salient features to tell multi- and single-plane lenses apart.

I give a summary of my research with Prof. Dr. Reinhard Schlickeiser. It mainly took place in the 1990s, when the Compton Gamma Ray Observatory was giving a startling new view of the gamma-ray sky. Our work focused on particle acceleration and radiation processes in the jets of active galactic nuclei (AGN). We pioneered the external Compton scattering model of blazars, where photons from outside the jet are intercepted and scattered to high energies by the radio-emitting electrons within the jet. Although we originally focused on external photons from the accretion disks of AGN, this process has been extended to include a range of external photon sources, and is now established as the primary source of gamma rays from flat spectrum radio quasars.

Thermal noise from the suspension fibres used in the mirror pendulums in current gravitational wave detectors is a critical noise source. Future detectors will require improved suspension performance with the specific ability to suspend much heavier masses to reduce radiation pressure noise, whilst retaining good thermal noise performance. In this letter, we propose and experimentally demonstrate a design for a large-scale fused silica suspension, demonstrating its suitability for holding an increased mass of 160 kg. We demonstrate the concepts for improving thermal noise via longer suspension fibres supporting a higher static stress. We present a full thermal noise analysis of our prototype, meeting requirements for conceptual 3rd generation detector designs such as the high frequency interferometer of the Einstein Telescope (ET-HF), and closely approaching that required for Cosmic Explorer (CE).

Long gamma-ray burst GRB 191016A was a bright and slow rising burst that was detected by the \textit{Swift} satellite and followed up by ground based Liverpool Telescope (LT). LT follow-up started $2411$-s after the \textit{Swift} Burst Alert Telescope (BAT) trigger using imager IO:O around the time of the late optical peak. From $3987-7687$-s, we used the LT polarimeter RINGO3 to make polarimetric and photometric observations of the GRB simultaneously in the $V,R$ and $I$ bands. The combined optical light curve shows an initial late peak followed by a decline until 6147-s, 6087-s, and 5247-s for $I,R$ and $V$ filters respectively followed by a flattening phase. There is evidence of polarization at all phases including polarization ($P = 14.6 \pm 7.2 \%$) which is coincident with the start of the flattening phase. The combination of the light curve morphology and polarization measurement favours an energy injection scenario where slower magnetised ejecta from the central engine catches up with the decelerating blast wave. We calculate the minimum energy injection to be $\Delta E / E>0.36$. At a later time combining the optical light curve from BOOTES (reported via GCN) and IO:O we see evidence of a jet break with jet opening angle 2 \degree.

As primary anchors of the distance scale, Cepheid stars play a crucial role in our understanding of the distance scale of the Universe because of their period-luminosity relation. Determining precise and consistent parameters (radius, temperature, color excess, and projection factor) of Cepheid pulsating stars is therefore very important. With the high-precision parallaxes delivered by the early third Gaia data release, we aim to derive various parameters of Cepheid stars in order to calibrate the period-luminosity and period-radius relations and to investigate the relation of period to p-factor. We applied an implementation of the parallax-of-pulsation method through the algorithm called Spectro-Photo-Interferometry of Pulsating Stars, which combines all types of available data for a variable star in a global modeling of its pulsation. We present the SPIPS modeling of a sample of 63 Galactic Cepheids. Adopting Gaia EDR3 parallaxes as an input associated with the best available dataset, we derive consistent values of parameters for these stars such as the radius, multiband apparent magnitudes, effective temperatures, color excesses, period changes, Fourier parameters, and the projection factor. We then derive new calibrations of the period-luminosity and period-radius relations. After investigating the dependences of the p-factor on the parameters of the stars, we find a high dispersion of its values and no evidence of its correlation with the period or with any other parameters. Statistically, the p-factor has an average value of p=1.26$\pm$0.07, but with an unsatisfactory agreement. In absence of any clear correlation between the p-factor and other quantities, the best agreement is obtained under the assumption that the p-factor can take any value in a band with a width of 0.15. This result highlights the need for a further examination of the physics behind the p-factor.

The detection of the radio signal from filaments in the cosmic web is crucial to distinguish possible magnetogenesis scenarios. We review the status of the different attempts to detect the cosmic web at radio wavelengths. This is put into the context of the advanced simulations of cosmic magnetism carried out in the last few years by our {\magcow} project. While first attempts of imaging the cosmic web with the MWA and LOFAR have been encouraging and could discard some magnetogenesis models, the complexity behind such observations makes a definitive answer still uncertain. A combination of total intensity and polarimetric data at low radio frequencies that the SKA and LOFAR2.0 will achieve is key to removing the existing uncertainties related to the contribution of many possible sources of signal along deep lines of sight. This will make it possible to isolate the contribution from filaments, and expose its deep physical connection with the origin of extragalactic magnetism.

We study the gas distribution and kinematics of the inner kpc of six moderately luminous (43.43$\leq$log$L_{\rm bol}\leq$44.83) nearby ($0.004 \leq z \leq 0.014$) Seyfert galaxies observed with the Near-infrared Integral Field Spectrograph (NIFS) in the J(1.25$\mu$m) and K(2.2$\mu$m) bands. We analyse the most intense emission lines detected on these spectral wavebands: [Fe II] 1.2570 $\mu$m and Pa$\beta$, which trace the ionised gas in the partially and fully ionised regions, and H$_2$2.1218$\mu$m, that traces the hot ($\sim 2000$ K) molecular gas. The dominant kinematic component is rotation in the disc of the galaxies, except for the ionised gas in NGC 5899 which shows only weak signatures of a disc component. We find ionised gas outflow in four galaxies, while signatures of H$_2$ outflows are seen in three galaxies. The ionised gas outflows display velocities of a few hundred km\,s$^{-1}$, and their mass-outflow rates are in the range $0.005 - 12.49$M$_{\odot}$yr$^{-1}$ . Their kinetic powers correspond to $0.005-0.7$ per cent of the AGN bolometric luminosities. Besides rotation and outflows signatures in some cases, the H 2 kinematics reveals also inflows in three galaxies. The inflow velocities are $50-80$km\,s$^{-1}$ and the mass inflow rates are in the range $1 - 9 \times 10^{-4}$M$_{\odot}$yr$^{-1}$ for hot molecular gas. These inflows might be only the hot skin of the total inflowing gas, which is expected to be dominated by colder gas. The mass inflow rates are lower than the current accretion rates to the AGN, and the ionised outflows are apparently disturbing the gas in the inner kpc.

We present a demonstration of the in-flight polarization angle calibration for the JAXA/ISAS second strategic large class mission, LiteBIRD, and estimate its impact on the measurement of the tensor-to-scalar ratio parameter, r, using simulated data. We generate a set of simulated sky maps with CMB and polarized foreground emission, and inject instrumental noise and polarization angle offsets to the 22 (partially overlapping) LiteBIRD frequency channels. Our in-flight angle calibration relies on nulling the EB cross correlation of the polarized signal in each channel. This calibration step has been carried out by two independent groups with a blind analysis, allowing an accuracy of the order of a few arc-minutes to be reached on the estimate of the angle offsets. Both the corrected and uncorrected multi-frequency maps are propagated through the foreground cleaning step, with the goal of computing clean CMB maps. We employ two component separation algorithms, the Bayesian-Separation of Components and Residuals Estimate Tool (B-SeCRET), and the Needlet Internal Linear Combination (NILC). We find that the recovered CMB maps obtained with algorithms that do not make any assumptions about the foreground properties, such as NILC, are only mildly affected by the angle miscalibration. However, polarization angle offsets strongly bias results obtained with the parametric fitting method. Once the miscalibration angles are corrected by EB nulling prior to the component separation, both component separation algorithms result in an unbiased estimation of the r parameter. While this work is motivated by the conceptual design study for LiteBIRD, its framework can be broadly applied to any CMB polarization experiment. In particular, the combination of simulation plus blind analysis provides a robust forecast by taking into account not only detector sensitivity but also systematic effects.

In dead zones of protoplanetary discs, it is assumed that micrometre-sized particles grow Brownian, sediment to the midplane and drift radially inward. When collisional compaction sets in, the growing aggregates collect slower and therefore dynamically smaller particles. This sedimentation and growth phase of highly porous ice and dust aggregates is simulated with laboratory experiments in which we obtained mm- to cm-sized ice aggregates with a porosity of 90\% as well as cm-sized dust agglomerates with a porosity of 85\%. We modelled the growth process during sedimentation in an analytical calculation to compute the agglomerate sizes when they reach the midplane of the protoplanetary disc. In the midplane, the dust particles form a thin dense layer and gain relative velocities by, e.g., the streaming instability or the onset of shear turbulence. To investigate also these collisions, we performed additional laboratory drop tower experiments with the high-porosity aggregates formed in the sedimentary-growth experiments and determined their mechanical parameters, including their sticking threshold velocity, which is important for their further collisional evolution on their way to form planetesimals. Finally, we developed a method to calculate the packing-density-dependent fundamental properties of our dust and ice agglomerates, the Young's modulus, the Poisson ratio, the shear viscosity and the bulk viscosity from compression measurements. With these parameters, it was possible to derive the coefficient of restitution which fits our measurements. In order to physically describe these outcomes, we applied a collision model. With this model, predictions about general dust-aggregate collisions are possible.

We study Luminous Blue Variable (LBV) candidate J004341.84+411112.0 in the Andromeda galaxy. We present optical spectra of the object obtained with the 6-m telescope of SAO RAS. The candidate shows typical LBV features in its spectra: broad and strong hydrogen lines and the He I lines with P Cigni profiles. Its remarkable spectral resemblance to the well known LBV P Cygni suggests a common nature of the objects and supports LBV classification of J004341.84+411112.0. We estimate the temperature, reddening, radius and luminosity of the star using its spectral energy distribution. Obtained bolometric luminosity of the candidate ($M_{\text{bol}}=-10.41\pm0.12$ mag) is quite similar to those of known LBV stars in the Andromeda galaxy. We analysed ten year light curve of the object in $R$ filter. The candidate demonstrates photometric variations of the order of 0.4 mag, with an overall brightness increasing trend $\Delta R > 0.1$ mag. Therewith, the corresponding colour variations of the object are fully consistent with LBV behavior when a star become cooler and brighter in the optical spectral range with a nearly constant bolometric luminosity. LBV-type variability of the object, similarity of its spectrum and estimated luminosity to those of known LBVs allows us to classify J004341.84+411112.0 as a LBV.

We study the variation of the gravitational Newton's constant on cosmological scales in scalar-tensor theories of gravity. We focus on the simplest models of scalar-tensor theories with a coupling to the Ricci scalar of the form $F(\sigma) = N_{pl}^2 + \xi\sigma^2$, such as extended Jordan-Brans-Dicke ($N_{pl}=0$), or a non-minimally coupled scalar field with $N_{pl}=M_{pl}$, which permits the gravitational constant to vary self-consistently in time and space. In addition, we allow the gravitational constant to differ from the Newton's constant $G$, i.e. $G_{\rm eff}(z=0) = G(1+\Delta)^2$. Combining the information from {\em Planck} 2018 CMB temperature, polarization and lensing, together with a compilation of BAO measurements from BOSS, we constrain the imbalance to $\Delta = -0.022 \pm 0.023$ (68% CL) and the coupling to $10^3\, \xi < 0.82$ (95% CL) for JBD and for a non-minimally coupled scalar field we constrain the imbalance to $\Delta > -0.018$ ($< 0.021$) and the coupling parameter to $\xi < 0.089$ ($\xi > - 0.041$) both at 95% CL. These constraints correspond to a variation of the gravitational constant now respect to the one in the radiation era to be smaller than 3% (95% CL) and to the ratio of the gravitational Newton's constant measured from cosmological scales and the one measured in a Cavendish-like experiment to be smaller than 4-15% (95% CL). With current data, we observe that the degeneracy between $\Delta$, the coupling $\xi$, and $H_0$ allows for a larger value of the Hubble constant increasing the agreement between the measurement of the Hubble constant by the SH0ES team and its value inferred by CMB data. Future data such as the combination of CMB anisotropies from LiteBIRD and CMB-S4, and large-scale structures galaxy clustering from DESI and galaxy shear from LSST will reduce the uncertainty to $\sigma(\Delta) = 0.004$.

Tidal disruption events (TDEs) are valuable probes of the demographics of supermassive black holes as well as the dynamics and population of stars in the centers of galaxies. In this Letter, we focus on studying how the TDE disk formation and circularization processes can impact the possibility of observing prompt TDE flares, based on recent theoretical developments. First, we investigate how the efficiency of disk formation is determined by the key parameters, namely, the black hole mass $M_{BH}$, the stellar mass $m_\star$, and the stellar orbital penetration parameter $\beta$. Then we apply the loss cone theory to calculate the differential TDE rate as a function of these three parameters. Combining these two results, we find that the rates of TDEs with prompt disk formation are significantly suppressed around lighter black holes, which provides a plausible explanation for why the observed TDE host black hole mass distribution peaks between $10^6$ and $10^7 M_\odot$. Therefore, the consideration of the disk formation efficiency is crucial for recovering the intrinsic black hole demographics from TDEs. Furthermore, we find that the efficiency of the disk formation process also impacts the distributions of both stellar orbital penetration parameter and stellar mass observed in TDEs.

Fluorescence emission lines are broadly applied in observation for diffuse medium in the universe. They are normally observed around strong pumping source, tracing the gas in circumstellar medium, reflection nebula, and H\,{\sc ii} regions, etc. They reside in UV/optical and infrared bands and hence could be directly observed with ground-base telescopes. In this letter, we demonstrate the polarization of fluorescence lines as a magnetic field tracer arising from ground state atomic alignment in diffuse medium, including our solar system, supernova remnants (SNRs), as well as quasi-stellar object (QSO) host galaxies. Two types of fluorescence emissions are considered: the primary fluorescence from the excited states; and the secondary fluorescence from the metastable state (forbidden lines). We find that the synergy of these lines could measure three-dimensional magnetic direction: the polarizations of the primary fluorescence lines could reveal the magnetic polar angle along the line-of-sight, whereas the polarization of forbidden lines traces the plane-of-sky magnetic direction. The expected degree of polarization is $P>10\%$. Polarizations of both types of fluorescence emissions have shown strong potential for observations, and are applicable to measure magnetic field within and beyond our galaxy.

Current theories of planetary evolution predict that infant giant planets have large radii and very low densities before they slowly contract to reach their final size after about several hundred million years. These theoretical expectations remain untested to date, despite the increasing number of exoplanetary discoveries, as the detection and characterisation of very young planets is extremely challenging due to the intense stellar activity of their host stars. However, the recent discoveries of young planetary transiting systems allow to place initial constraints on evolutionary models. With an estimated age of 20 million years, V1298\,Tau is one of the youngest solar-type stars known to host transiting planets: it harbours a multiple system composed of two Neptune-sized, one Saturn-sized, and one Jupiter-sized planets. Here we report the analysis of an intense radial velocity campaign, revealing the presence of two periodic signals compatible with the orbits of two of its planets. We find that planet b, with an orbital period of 24 days, has a mass of 0.64 Jupiter masses and a density similar to the giant planets of the Solar System and other known giant exoplanets with significantly older ages. Planet e, with an orbital period of 40 days, has a mass of 1.16 Jupiter masses and a density larger than most giant exoplanets. This is unexpected for planets at such a young age and suggests that some giant planets might evolve and contract faster than anticipated, thus challenging current models of planetary evolution.

Relativistic magnetohydrodynamic shocks are efficient particle accelerators, often invoked in the models of gamma-ray bursts (GRBs) and shock-powered fast radio bursts (FRBs). Most theoretical studies assume a perpendicular shock with an ordered magnetic field perpendicular to the shock normal. However, the degree of magnetization $\sigma$ and the magnetic field geometry in shock-powered GRB/FRB scenarios are still poorly constrained by observations. Analogous to the magnetization $\sigma$ associated with the total field strength, we define a tangential magnetization $\sigma_\perp$ associated with the tangential field component. We explore the jump conditions of magnetized relativistic shocks, either with an ordered field of arbitrary inclination angle or with a random field of arbitrary anisotropy. In either case, we find that the jump conditions of relativistic shocks are governed by the tangential magnetization $\sigma_\perp$ instead of the total magnetization $\sigma$, insensitive to the inclination angles or the anisotropy of the pre-shock magnetic field. The approximated analytical solution developed in this work could serve as a quick check for numerical simulations and apply to theoretical studies of GRBs/FRBs with a more general field geometry.

AMEGO-X, the All-sky Medium Energy Gamma-Ray Observatory eXplorer, is a proposed instrument designed to bridge the so-called "MeV gap" by surveying the sky with unprecedented sensitivity from ~100 keV to about one GeV. This energy band is of key importance for multi-messenger and multi-wavelength studies but it is nevertheless currently under-explored. AMEGO-X addresses this situation by proposing a design capable of detecting and imaging gamma rays via both Compton interactions and pair production processes. However, some of the objects that AMEGO-X will study, such as gamma-ray bursts progenitors and magnetars, extend to energies below ~100 keV where the dominant interaction becomes photoelectric absorption. These events deposit their energy in a single pixel of the detector. In this work we show how the ~3500 cm^2 effective area of the AMEGO-X tracker to events between ~25 keV to ~100 keV will be utilized to significantly improve its sensitivity and expand the energy range for transient phenomena. Although imaging is not possible for single-site events, we show how we will localize a transient source in the sky using their aggregate signal --within a 10deg radius for an event similar to GRB 170817A. This technique will triple the number of cosmological gamma-ray bursts seen by AMEGO-X, allow us to detect and resolve the pulsating tails of extragalactic magnetar giant flares, and increase the number of detected less-energetic magnetar bursts --some of which are associated with fast radio bursts. Overall, single-site events will increase the energy range of the sensitivity, expand the science program, and promptly alert the community of fainter transient events.

We present the Strong Lensing Science Collaboration's (SLSC) recommended observing targets for the science verification and science validation phases of commissioning. Our recommendations have been developed in collaboration with the Dark Energy Science Collaboration (DESC) Strong Lensing Topical Team. In summary, our key recommendations are as follows: (1) Prioritize fields that span the full range of declination observable from Cerro Pachon during the engineering focused Science Verification phase of commissioning, before concentrating on equatorial fields for the Science Validation surveys. (2) Observe quadruply lensed quasars as the ultimate test of the Active Optics system towards the end of the Science Verification phase of commissioning. These systems are the strongest tests known for delivered image quality. (3) Prioritize science validation survey fields (both single deep pointings and wide fields) that have been searched thoroughly by precursor surveys for strong lenses. (4) The optimal wide (~100 degree^2) science validation field would include the CFHT-LS W4 field, and overlap with the SDSS Stripe 82, DES-SN, KIDS and HSC-SSP fields. (5) The optimal single pointing science validation fields are the XMM-LSS and COSMOS Deep Drilling Fields, the equatorial Hubble Frontier Fields galaxy clusters, and strongly lensed quasars with measured time delays that are well-matched to commissioning timescales.

The reasons behind the Great Dimming and subsequent rising in the brightness of Betelgeuse between October 2019 and March 2020 still continue to baffle astronomers. It has been shown by George et. al. (2020) that critical slowing down preceded the dimming event. This suggested that the dimming was a result of the change in the nature of the nonlinear dynamics of the star. In this work we present additional evidence for dynamical changes in Betelgeuse prior to the Great Dimming event, using nonlinear time series analysis. We study the relations between the different bands in the photometry data collected from the Wing photometry (IR/near-IR) and Wasatonic observatory (V-band). We also analyse how the early warning signals studied previously changed during and after the Great Dimming.

In this work we investigate properties of intermediate-mass black holes (IMBHs) that escape from star clusters due to dynamical interactions. The studied models were simulated as part of the preliminary second survey carried out using the MOCCA code (MOCCA-SURVEY Database II), which is based on the Monte Carlo N-body method. We have found that IMBHs are more likely to be formed and ejected in models where both initial central density and central escape velocities have high values. Most of our studied objects escape in a binary with another black hole (BH) as their companion and they have masses between $100$ and $140\: M_{\odot}$. Escaping IMBHs tend to build-up mass most effectively through repeated mergers in a binary with BHs due to gravitational wave emission. Binaries play a key role in their ejection from the system as they allow these massive objects to gather energy needed for escape. The binaries in which IMBHs escape tend to have very high binding energy at the time of escape and the last interaction is strong but does not involve a massive intruder. These IMBHs gain energy needed to escape the cluster gradually in successive dynamical interactions. We present specific examples of the history of IMBH formation and escape from our star cluster models. We also discuss the observational implications of our findings and estimate that the merger rate of ejected IMBH binaries from star clusters is $\mathcal{R}=0.7\:\textrm{Gpc}^{-3}\:\textrm{yr}^{-1}$.

Mergers between galaxy clusters often drive shocks into the intra cluster medium (ICM), the effects of which are sometimes visible via temperature and density jumps in the X-ray, and via radio emission from relativistic particles energized by the shock's passage. Abell2108 was selected as a likely merger system through comparing the X-ray luminosity to the Planck Sunyaev-Zeldovich signal, where this cluster appeared highly X-ray underluminous. Follow up observations confirmed it to be a merging low mass cluster featuring two distinct subclusters, both with a highly disturbed X-ray morphology. Giant Metrewave Radio Telescope (GMRT) data covering 120-750MHz show an extended radio feature resembling a radio relic, near the location of a temperature discontinuity in the X-rays. We measure a Mach number from the X-ray temperature jump. Several characteristics of radio relics are found in Abell2108, making this cluster one of the few low mass mergers likely hosting a radio relic. The radio spectrum is exceptionally steep, and the radio power is very weak (P1.4GHz=1E22W/Hz). To account for the shock/relic offset, we propose a scenario in which the shock created the relic by re-accelerating a cloud of pre-existing relativistic electrons and then moved away, leaving behind a fading relic. The electron aging timescale derived from the high-frequency steepening in the relic spectrum is consistent with the shock travel time to the observed X-ray discontinuity. However, the lower flux in GMRT band 4 data causing the steepening could be due to instrumental limitations, and deeper radio data are needed to constrain the spectral slope of the relic at high frequencies.

The rotational evolution of a young stellar population can give informations about the rotation pattern of more evolved clusters. Combined with rotational period values of thousands of young stars and theoretical propositions about the redistribution and loss of stellar angular momentum, it allows us to trace the rotational history of stars according to their mass. We want to investigate how internal and environmental changes on single stars can change the rotational evolution of a young stellar population. We run Monte Carlo simulations of a young cluster composed by solar mass stars of 0.5, 0.8 and 1.0 M$_\odot$ from 1 to 550 Myr taking into account observational and theoretical parameters. In order to compare our results with the observations we run Kolmogorov-Smirnov tests. Our standard model is able to reproduce some clusters younger than h Per and marginally M37, which is 550 Myr old. Varying the disk fraction or the initial period distribution did not improve the results. However, when we run a model with a finer mass grid the Pleiades can be also reproduced. Changing the initial mass distribution to be similar to the empirical ONC mass function also gives good results. Modeling the evolution of a young synthetic cluster from pre-main sequence to early main sequence considering physical mechanisms of extraction and exchange of angular momentum can not be achieved successfully for all clusters for which we have enough rotational data. Clusters of about the same age present different rotational behaviors due perhaps to differences in their initial conditions.

We present a comprehensive multiwavelength investigation of a likely massive young cluster `IRAS 05100+3723' and its environment with the aim to understand its formation history and feedback effects. We find that IRAS 05100+3723 is a distant ($\sim$3.2 kpc), moderate mass ($\sim$500 \msun), young ($\sim$3 Myr) cluster with its most massive star being an O8.5V-type. From spectral modeling, we estimate the effective temperature and log $g$ of the star as $\sim$33,000 K and $\sim$3.8, respectively. Our radio continuum observations reveal that the star has ionized its environment forming an HII region of size $\sim$2.7 pc, temperature $\sim$5,700 K, and electron density $\sim$165 cm$^{-3}$. However, our large-scale dust maps reveal that it has heated the dust up to several parsecs ($\sim$10 pc) in the range 17$-$28 K and the morphology of warm dust emission resembles a bipolar HII region. From dust and $^{13}$CO gas analyses, we find evidences that the formation of the HII region has occurred at the very end of a long filamentary cloud around 3 Myr ago, likely due to edge collapse of the filament. We show that the HII region is currently compressing a clump of mass $\sim$2700 \msun at its western outskirts, at the junction of the HII region and filament. We observe several 70 $\mu$m point sources of intermediate-mass and class 0 nature within the clump. We attribute these sources as the second generation stars of the complex. We propose that the star formation in the clump is either induced or being facilitated by the compression of the expanding HII region onto the inflowing filamentary material.

We present new cooling models for carbon-oxygen white dwarfs with both H- and He-atmospheres, covering the whole relevant mass range, to extend our updated BaSTI (a Bag of Stellar Tracks and Isochrones) stellar evolution archive. They have been computed using core chemical stratifications obtained from new progenitor calculations, adopting a semiempirical initial-final mass relation. The physics inputs have been updated compared to our previous BaSTI calculations: ^{22}Ne diffusion in the core is now included, together with an updated CO phase diagram, and updated electron conduction opacities. We have calculated models with various different neon abundances in the core, suitable to study white dwarfs in populations with metallicities ranging from super-solar to metal poor, and have performed various tests/comparisons of the chemical stratification and cooling times of our models. Two complete sets of calculations are provided, for two different choices of the electron conduction opacities, to reflect the current uncertainty in the evaluation of the electron thermal conductivity in the transition regime between moderate and strong degeneracy, crucial for the H and He envelopes. We have also made a first, preliminary estimate of the effect -- that turns out to be generally small -- of Fe sedimentation on the cooling times of white dwarf models, following recent calculations of the phase diagrams of carbon-oxygen-iron mixtures. We make publicly available the evolutionary tracks from both sets of calculations, including cooling times and magnitudes in the Johnson-Cousins, Sloan, Pan-STARSS, Galex, Gaia-DR2, Gaia-eDR3, HST-ACS, HST-WFC3, and JWST photometric systems.

Perturbative analyses of planetary resonances commonly predict singularities and/or divergences of resonance widths at very low and very high eccentricities. We have recently re-examined the nature of these divergences using non-perturbative numerical analyses, making use of Poincar\'e sections but from a different perspective relative to previous implementations of this method. This perspective reveals fine structure of resonances which otherwise remains hidden in conventional approaches, including analytical, semi-analytical and numerical-averaging approaches based on the critical resonant angle. At low eccentricity, first order resonances do not have diverging widths but have two asymmetric branches leading away from the nominal resonance location. A sequence of structures called ``low-eccentricity resonant bridges" connecting neighboring resonances is revealed. At planet-grazing eccentricity, the true resonance width is non-divergent. At higher eccentricities, the new results reveal hitherto unknown resonant structures and show that these parameter regions have a loss of some -- though not necessarily entire -- resonance libration zones to chaos. The chaos at high eccentricities was previously attributed to the overlap of neighboring resonances. The new results reveal the additional role of bifurcations and co-existence of phase-shifted resonance zones at higher eccentricities. By employing a geometric point of view, we relate the high eccentricity phase space structures and their transitions to the shapes of resonant orbits in the rotating frame. We outline some directions for future research to advance understanding of the dynamics of mean motion resonances.

Recent observations of several peculiar over- and under-luminous type Ia supernovae infer indirect evidence for the violation of the Chandrasekhar mass-limit by suggesting the existence of super- and sub-Chandrasekhar limiting mass white dwarfs. In an attempt to explain these phenomena in the context of general relativistic extensions, we study these objects in Palatini $f(R)$ gravity. We obtain the super- and sub-Chandrasekhar limiting masses as well as the dynamical instability criteria for white dwarfs in the given gravitational theory. We further demonstrate that the conventional positivity condition $\partial{M}/\partial{\rho_\text{c}}>0$ with $M$ being the WD's mass with central density $\rho_\text{c}$, is also a valid criterion for stability in Palatini gravity.

We investigate the properties of a special class of singular solutions for a self-gravitating perfect fluid in general relativity: the singular isothermal sphere. For arbitrary constant equation-of-state parameter $w=p/\rho$, there exist static, spherically-symmetric solutions with density profile $\propto 1/r^2$, with the constant of proportionality fixed to be a special function of $w$. Like black holes, singular isothermal spheres possess a fixed mass-to-radius ratio independent of size, but no horizon cloaking the curvature singularity at $r=0$. For $w=1$, these solutions can be constructed from a homogeneous dilaton background, where the metric spontaneously breaks spatial homogeneity. We study the perturbative structure of these solutions, finding the radial modes and tidal Love numbers, and also find interesting properties in the geodesic structure of this geometry. Finally, connections are discussed between these geometries and dark matter profiles, the double copy, and holographic entropy, as well as how the swampland distance conjecture can obscure the naked singularity.

Neutrino telescopes provide strong sensitivity to sterile neutrino oscillations through matter-enhanced oscillation, occurring in the few TeV energy range for eV$^{2}$-scale neutrino mass-squared splittings. Prior searches have focused on $\nu_\mu$ disappearance, which has a particularly strong sensitivity to the mixing angle $\theta_{24}$ via $\nu_\mu\rightarrow\nu_s$ transitions. Nowadays, the $\nu_\mu\rightarrow\nu_e$ and $\nu_\mu\rightarrow\nu_\tau$ appearance channels have been considered less promising at neutrino telescopes, due in part to the much smaller target volume for cascades relative to tracks. This work explores the detectability of these signatures at neutrino telescopes given present constraints on sterile neutrino mixing, and as an example, forecasts the sensitivity of the IceCube Neutrino Observatory to the mixing angles $\theta_{14}$, $\theta_{24}$, and $\theta_{34}$ in the 3+1 sterile neutrino model using the cascade channel with ten years of data. We find that $\nu_\tau$ appearance signatures consistent with the existing IceCube $\nu_\mu$ disappearance best-fit point are discoverable for values of $\theta_{34}$ consistent with world constraints, and that the sterile neutrino parameters favored by the BEST and gallium anomalies are expected to be testable at the 95\% confidence level.

Time dependent magnetic fields can be sourced by spinning neutron stars, orbiting binaries and merging neutron stars. We consider electromagnetic radiation from axion condensates in the background of an alternating magnetic field. We find that a resonant peak in radiation can occur when the frequency of the alternating magnetic field is comparable with the axion mass scale. More interestingly, in situations where the frequency of the alternating magnetic field itself changes with time, as can be the case in binary mergers due to steady increase in orbital frequency, the resonant peak in radiation may occur for a range of axion mass scales scanned by the time-varying magnetic field frequency.

We explore the neutrinoless double beta ($0\nu \beta\beta$) decay induced by an ultralight dark matter field coupled to neutrinos. The effect on $0\nu\beta\beta$ decay is significant if the coupling violates the lepton number, for which the $\Delta L=2$ transition is directly driven by the dark matter field without further suppression of small neutrino masses. As the ultralight dark matter can be well described by a classical field, the effect features a periodic modulation pattern in decay events. However, we find that in the early Universe such coupling will be very likely to alter the standard cosmological results. In particular, the requirement of neutrino free-streaming before the matter-radiation equality severely constrains the parameter space, such that the future $0\nu \beta\beta$ decay experiments can hardly see any signal even with a meV sensitivity to the effective neutrino mass.

The symmetry breaking of grand unified gauge groups in the early Universe often leaves behind relic topological defects such as cosmic strings, domain walls, or monopoles. For some symmetry breaking chains, hybrid defects can form where cosmic strings attach to domain walls or monopoles attach to strings. In general, such hybrid defects are unstable, with one defect "eating" the other via the conversion of its rest mass into the other's kinetic energy and subsequently decaying via gravitational waves. In this work, we determine the gravitational wave spectrum from 1) the destruction of a cosmic string network by the nucleation of monopoles which cut up and "eat" the strings, 2) the collapse and decay of a monopole-string network by strings that "eat" the monopoles, 3) the destruction of a domain wall network by the nucleation of string-bounded holes on the wall that expand and "eat" the wall, and 4) the collapse and decay of a string-bounded wall network by walls that "eat" the strings. We call the gravitational wave signals produced from the "eating" of one topological defect by another gravitational wave gastronomy. We find that the four gravitational wave gastronomy signals considered yield unique spectra that can be used to narrow down the SO(10) symmetry breaking chain to the Standard Model and the scales of symmetry breaking associated with the consumed topological defects. Moreover, the systems we consider are unlikely to have a residual monopole or domain wall problem.

COSINE-100 is a direct detection dark matter experiment that aims to test DAMA/LIBRA's claim of dark matter discovery by searching for a dark matter-induced annual modulation signal with NaI(Tl) detectors. We present new constraints on the annual modulation signal from a dataset with a 2.82 yr livetime utilizing an active mass of 61.3 kg, for a total exposure of 173 kg$\cdot$yr. This new result features an improved event selection that allows for both lowering the energy threshold to 1 keV and a more precise time-dependent background model. In the 1-6 keV and 2-6 keV energy intervals, we observe best-fit values for the modulation amplitude of 0.0067$\pm$0.0042 and 0.0050$\pm$0.0047 counts/(day$\cdot$kg$\cdot$keV), respectively, with a phase fixed at 152.5 days.

The structured hadron-quark mixed phase, known as the pasta phase, is expected to appear in the core of massive neutron stars. Motivated by the recent advances in astrophysical observations, we explore the possibility of the appearance of quarks inside neutron stars and check its compatibility with current constraints. We investigate the properties of the hadron-quark pasta phases and their influences on the equation of state (EOS) for neutron stars. In this work, we extend the energy minimization (EM) method to describe the hadron-quark pasta phase, where the surface and Coulomb contributions are included in the minimization procedure. By allowing different electron densities in the hadronic and quark matter phases, the total electron chemical potential with the electric potential remains constant, and local ? equilibrium is achieved inside the Wigner-Seitz cell. The mixed phase described in the EM method shows the features lying between the Gibbs and Maxwell constructions, which is helpful for understanding the transition from the Gibbs construction (GC) to the Maxwell construction (MC) with increasing surface tension. We employ the relativistic mean-field model to describe the hadronic matter, while the quark matter is described by the MIT bag model with vector interactions. It is found that the vector interactions among quarks can significantly stiffen the EOS at high densities and help enhance the maximum mass of neutron stars. Other parameters like the bag constant can also affect the deconfinement phase transition in neutron stars. Our results show that hadron-quark pasta phases may appear in the core of massive neutron stars that can be compatible with current observational constraints.

We propose several extensions of the Starobinsky model of inflation, which obey all observational constraints on the inflationary parameters, by demanding that both the inflaton scalar potential in the Einstein frame and the $F(R)$ gravity function in the Jordan frame have the explicit dependence upon fields and parameters in terms of elementary functions. All our models are continuously connected to the original Starobinsky model via changing the parameters. We modify the Starobinsky $(R+R^2)$ model by adding an $R^3$-term, an $R^4$-term, and an $R^{3/2}$-term, respectively, and calculate the scalar potentials, the inflationary observables and the allowed limits on the deformation parameters in all these cases. We also deform the scalar potential of the Starobinsky model in the Einstein frame, in powers of $y=\exp\left(-\sqrt{\frac{2}{3}}\phi/M_{Pl}\right)$, where $\phi$ is the canonical inflaton (scalaron) field, calculate the corresponding $F(R)$ gravity functions in the two new cases, and find the restrictions on the deformation parameters in the lowest orders with respect to the variable $y$ that is physically small during slow-roll inflation.

Using the effective field theory framework for extended objects and the coset construction, we build the leading order effective action for the most general compact object allowed in an effective theory of gravity as general relativity. We derive all the covariant building blocks of the theory and construct the relevant terms in the action that are invariant under the symmetries of the problem. Our effective action contains a relativistic point mass term, and relativistic corrections due to spin, charge, size effects and dissipation. Each of the invariant operators that account for the internal structure of the compact object is accompanied by a coefficient that encapsulates its properties. We match the known coefficients of the effective action from the literature, and point out the unknown ones that are to be derived.

The thermal balance between the black hole and surrounding thermal radiation is an interesting topic. The primordial black hole is a type of black hole formed in the early universe, especially during the thermal radiation period of the early universe. Driven by gravity, thermal radiation will collapse to form a black hole. Because of the complexity of the actual situation, we use the simplest case to simulate: an isolated box full of thermal radiation, in what direction will it evolve? Whether to achieve a pure radiation thermal balance or collapse to form a black hole is limited by the difference in the initial conditions, and the outcome will also be different. If we consider the quantization of the black hole and the number of microscopic states of the radiation, is it consistent with the classic conclusion? These are the main points to be considered in this paper.

Second-order tensor modes induced by nonlinear gravity are a key component of the cosmological background of gravitational waves that have the potential to probe the primordial power spectrum at otherwise inaccessible scales. A widely-recognized issue is that these modes can contain gauge dependent degrees of freedom. This poses a problem since observable gravitational waves should not depend on any choice of gauge. The ambiguity arises from the non-covariant transverse-traceless condition attached to a specific gauge choice as well as the misidentification of gravitational waves after a gauge transformation. We show that the covariant transverse-traceless projection of the extrinsic curvature consistently extracts the kinetic energy of gravitational waves associated with a given hypersurface. This projection is gauge invariant up to second order in a perturbative expansion, providing that we identify the original hypersurface after the gauge transformation. We emphasize that the usual tensor modes in different gauges are physically different quantities and thus should not necessarily be directly compared.

The primordial Universe might be highly inhomogeneous. We perform the 3+1D Numerical Relativity simulation for the evolution of scalar field in an initial inhomogeneous expanding Universe, and investigate how it populates the landscape with both de Sitter (dS) and AdS vacua. The simulation results show that eventually either the field in different region separates into different vacua, so that the expanding dS or AdS bubbles (the bubble wall is expanding but the spacetime inside AdS bubbles is contracting) come into being with clear bounderies, or overall region is dS expanding with a few smaller AdS bubbles (which collapsed into black holes) or inhomogeneously collapsing.

Gravitational waves from the collision of binary neutron stars provide a unique opportunity to study the behaviour of supranuclear matter, the fundamental properties of gravity, and the cosmic history of our Universe. However, given the complexity of Einstein's Field Equations, theoretical models that enable source-property inference suffer from systematic uncertainties due to simplifying assumptions. We develop a hypermodel approach to compare and measure the uncertainty gravitational-wave approximants. Using state-of-the-art models, we apply this new technique to the binary neutron star observations GW170817 and GW190425 and the sub-threshold candidate GW200311_103121. Our analysis reveals subtle systematic differences between waveform models, and a frequency-dependence study suggests that this is due to the treatment of the tidal sector. This new technique provides a proving ground for model development, and a means to identify waveform-systematics in future observing runs where detector improvements will increase the number and clarity of binary neutron star collisions we observe.

We relate the known Oberth effect and the nonrelativistic analogue of the Penrose process. When a particle decays to two fragments, we derive the conditions on the angles under which debris can come out for such a process to occur. We also consider the decay and the Oberth effect in the relativistic case, when a particle moves in the background of the Schwarzschild black hole. This models the process when a rocket ejects fuel. Different scenarios are analyzed depending on what data are fixed. The efficiency of the process is found, in particular near the horizon and for a photon rocket (when the ejected particle is massless). We prove directly that the most efficient process occurs when fuel is ejected along the rocket trajectory. When this occurs on the horizon, the efficiency reaches 100 % for a photon rocket. We also consider briefly the scenario when a rocket hangs over a black hole due to continuous ejection of fuel. Then, the fuel mass decays exponentially with the proper time.

Recent analysis of the viability of solid state-based relic neutrino detectors has revealed the fundamental necessity for the use of heavy, $A>100$, $\beta$-decayers as neutrino targets. Of all heavy isotopes, $^{171}$Tm and $^{151}$Sm stand out for their sufficiently low decay energies, reasonable half-life times and stable daughter nuclei. However, the crucial bit of information, that is the soft neutrino capture cross-section is missing for both isotopes. The main reason for that is a particular type of $\beta$-decay, which precludes a simple link between the isotope's half-life time and the neutrino capture rate. Here we propose an experimental method to bypass this difficulty and obtain the capture cross-section of a soft neutrino by a given isotope from the isotope's $\beta$-spectrum.