Papers for Friday, Jul 16 2021

We present a new release of the RELTRANS model to fit the complex cross-spectrum of accreting black holes as a function of energy. The model accounts for continuum lags and reverberation lags self-consistently in order to consider the widest possible range of X-ray variability timescales. We introduce a more self-consistent treatment of the reverberation lags, accounting for how the time variations of the illuminating flux change the ionisation level of the accretion disc. This process varies the shape of the reflection spectrum in time causing an additional source of lags besides the light crossing delay. We also consider electron densities in the accretion disc up to $10^{20}$ cm$^{-3}$, which are found in most of the stellar mass black holes and in some AGN. These high densities increase the amplitude of the reverberation lags below $1$ keV since the reflection flux enhances in the same energy range. In addition, we investigate the properties of hard lags produced by variations in the power-law index of the continuum spectrum, which can be interpreted as due to roughly $3\%$ variability in the corona's optical depth and temperature. As a test case, we simultaneously fit the lag energy spectra in a wide range of Fourier frequency for the black hole candidate MAXI J1820+070 observed with NICER. The best fit shows how the reverberation lags contribute even at the longer timescales where the hard lags are important. This proves the importance of modelling these two lags together and self-consistently in order to constrain the parameters of the system.

Dark Matter (DM) halos with masses below $\sim10^{8}$ $M_{\odot}$, which would help to discriminate between DM models, may be detected through their gravitational effect on distant sources. The same applies to primordial black holes, considered as an alternative scenario to DM particle models. However, there is still no evidence for the existence of such objects. With the aim of finding compact objects in the mass range $\sim$ 10$^{6}$ -- 10$^{9}$$M_{\odot}$, we search for strong gravitational lenses on milli (mas)-arcseconds scales (< 150 mas). For our search, we used the Astrogeo VLBI FITS image database -- the largest publicly available database, containing multi-frequency VLBI data of 13828 individual sources. We used the citizen science approach to visually inspect all sources in all available frequencies in search for images with multiple compact components on mas-scales. At the final stage, sources were excluded based on the surface brightness preservation criterion. We obtained a sample of 40 sources that passed all steps and therefore are judged to be milli-arcsecond lens candidates. These sources are currently followed-up with on-going European VLBI Network (EVN) observations at 5 and 22 GHz. Based on spectral index measurements, we suggest that two of our candidates have a higher probability to be associated with gravitational lenses.

Direct collapse models for black hole (BH) formation predict massive ($\sim 10^5 M_{\odot}$) seeds, which hold great appeal as a means to rapidly grow the observed $\sim 10^9 M_{\odot}$ quasars by $z\gtrsim 7$; however, their formation requires fine-tuned conditions. In this work, we use cosmological zoom simulations to study systematically the impact of requiring: 1) low gas angular momentum, and 2) a minimum incident Lyman Werner (LW) flux radiation in order to form direct-collapse BH seeds. We start with a baseline model (introduced in Bhowmick et. al. 2021) that restricts black hole seed formation (with seed masses of $M_{\mathrm{seed}}=1.25\times10^{4},1\times10^{5} \& 8\times10^{5}M_{\odot}/h$) to occur only in haloes with a minimum total mass ($3000\times M_{\mathrm{seed}}$) and star forming, metal poor gas mass ($5\times M_{\mathrm{seed}}$). When seeding is further restricted to halos with low gas spins (i.e. smaller than the minimum value required for the gas disc to be gravitationally stable), the seeding frequency is suppressed by factors of $\sim6$ compared to the baseline model regardless of the mass threshold used. In contrast, imposing a minimum LW flux ($>10J_{21}$) disproportionately suppresses seed formation in $\lesssim10^9M_{\odot}/h$ halos, by factors of $\sim100$ compared to the baseline model. Very few BH merger events occur in the models with a LW flux criterion, and because early BH growth is dominated by mergers in our models, this results in only the most massive ($8\times10^{5}M_{\odot}/h$) seeds being able to grow to the supermassive regime ($\gtrsim10^6M_{\odot}/h$) by $z=7$. Our results therefore suggest that producing the bulk of the $z\gtrsim7$ BH population requires alternate seeding channels, early BH growth dominated by rapid or super-eddington accretion, massive seeding scenarios that do not depend on LW flux, or a combination of these possibilities.

Radial velocity (RV) surveys have discovered giant exoplanets on au-scale orbits with a broad distribution of eccentricities. Those with the most eccentric orbits are valuable laboratories for testing theories of high eccentricity migration. However, few such exoplanets transit their host stars thus removing the ability to apply constraints on formation from their bulk internal compositions. We report the discovery of Kepler-1704 b, a transiting 4.15 $M_{\rm J}$ giant planet on a 988.88 day orbit with the extreme eccentricity of $0.921^{+0.010}_{-0.015}$. Our decade-long RV baseline from the Keck I telescope allows us to measure the orbit and bulk heavy element composition of Kepler-1704 b and place limits on the existence of undiscovered companions. Kepler-1704 b is a failed hot Jupiter that was likely excited to high eccentricity by scattering events that possibly began during its gas accretion phase. Its final periastron distance was too large to allow for tidal circularization, so now it orbits it host from distances spanning 0.16 - 3.9 au. The maximum difference in planetary equilibrium temperature resulting from this elongated orbit is over 700 K. A simulation of the thermal phase curve of Kepler-1704 b during periastron passage demonstrates that it is a remarkable target for atmospheric characterization from the James Webb Space Telescope, which could potentially also measure the planet's rotational period as the hot spot from periastron rotates in and out of view. Continued characterization of the Kepler-1704 system promises to refine theories explaining the formation of hot Jupiters and cool giant planets like those in the solar system.

We investigate the extent to which the number of clusters of mass exceeding $10^{15}\,M_{\odot}\,h^{-1}$ within the local super-volume ($<135\mathrm{\,Mpc}h^{-1}$) is compatible with the standard $\Lambda$CDM cosmological model. Depending on the mass estimator used, we find that the observed number $N$ of such massive structures can vary between $0$ and $5$. Adopting $N=5$ yields $\Lambda$CDM likelihoods as low as $2.4\times 10^{-3}$ (with $\sigma_8=0.81$) or $3.8\times 10^{-5}$ (with $\sigma_8=0.74$). However, at the other extreme ($N=0$), the likelihood is of order unity. Thus, while potentially very powerful, this method is currently limited by systematic uncertainties in cluster mass estimates. This motivates efforts to reduce these systematics with additional observations and improved modelling.

We have used Galaxy Zoo DECaLS (GZD) to study strong and weak bars in disk galaxies. Out of the 314,000 galaxies in GZD, we created a volume-limited sample (0.01 < z < 0.05, Mr < -18.96) which contains 1,867 galaxies with reliable volunteer bar classifications in the ALFALFA footprint. In keeping with previous Galaxy Zoo surveys (such as GZ2), the morphological classifications from GZD agree well with previous morphological surveys. GZD considers galaxies to either have a strong bar (15.5%), a weak bar (28.1%) or no bar (56.4%), based on volunteer classifications on images obtained from the DECaLS survey. This places GZD in a unique position to assess differences between strong and weak bars. We find that the strong bar fraction is typically higher in quiescent galaxies than in star forming galaxies, while the weak bar fraction is similar. Moreover, we have found that strong bars facilitate the quenching process in star forming galaxies, finding higher fibre SFRs, lower gas masses and shorter depletion timescales in these galaxies compared to unbarred galaxies. However, we also found that any differences between strong and weak bars disappear when controlling for bar length. Based on this, we conclude that weak and strong bars are not fundamentally different phenomena. Instead, we propose that there is a continuum of bar types, which varies from 'weakest' to 'strongest'.

We have developed a new analytic method to compute the galaxy two-point correlation functions (TPCFs) efficiently and accurately. We have derived precise formulas of the normalized random-random pair counts $RR(r)$ as functions of the survey area dimensions. We have also suggested algorithms to compute the normalized data-random pair counts $DR(r)$ analytically for a given data set and survey area. This method is applicable to survey areas with rectangular, cuboidal, circular, or spherical shapes. With all edge corrections fully considered analytically, our method computes $RR(r)$ and $DR(r)$ with perfect accuracy and zero variance in $O(1)$ and $O(N_g)$ time, respectively. Together with the data-data pair counts $DD(r)$, they can be used to compute TPCFs with any estimator at a speed 4 to 5 orders of magnitude faster than the traditional brute-force Monte Carlo method, which has $O(N_r^2)$ scaling. This method is favored over the Monte Carlo method whenever applicable. It is also directly applicable to galaxy angular TPCFs when the survey area is close to Euclidean.

Aims: We present a statistical analysis of different characteristics of ringed spiral galaxies with the aim of assessing the effects of rings on disk galaxy properties. Methods: We built a catalog of ringed galaxies from the Sloan Digital Sky Survey Data Release 14 (SDSS-DR14). Via visual inspection of SDSS images, we classified the face-on spiral galaxies brighter than $g < 16.0$ mag into galaxies with: an inner ring, an outer ring, a nuclear ring, both an inner and an outer ring, and a pseudo-ring. In addition to rings, we recorded morphological types and the existence of bars, lenses, and galaxy pair companions with or without interaction. With the goal of providing an appropriate quantification of the influence of rings on galaxy properties, we also constructed a suitable control sample of non-ringed galaxies with similar redshift, magnitude, morphology, and local density environment distributions to those of ringed ones. Results: We found 1868 ringed galaxies, accounting for 22% of the full sample of spiral galaxies. In addition, within galaxies with ringed structures, 46% have an inner ring, 10% an outer ring, 20% both an inner and an outer ring, 6% a nuclear ring, and 18% a partial ring. Moreover, 64% of the ringed galaxies present bars. We also found that ringed galaxies have both a lower efficiency of star formation activity and older stellar populations (as derived with the $D_n(4000)$ spectral index) with respect to non-ringed disk objects from the control sample. Galaxies with ringed structures present an excess of high metallicity values compared to non-ringed ones, which show a $12 + \rm Log (\rm O /\rm H)$ distribution toward lower values. These findings seem to indicate that rings are peculiar structures that produce an accelerating galactic evolution, strongly altering the physical properties of their host galaxies.

We present a complete dynamical study of the intermediate polar and dwarf nova cataclysmic variable GK Per (Nova Persei 1901) based on a multi-site optical spectroscopy and $R$-band photometry campaign. The radial velocity curve of the evolved donor star has a semi-amplitude $K_2=126.4 \pm 0.9 \, \mathrm{km}\,\mathrm{s}^{-1}$ and an orbital period $P=1.996872 \pm 0.000009 \, \mathrm{d}$. We refine the projected rotational velocity of the donor star to $v_\mathrm{rot} \sin i = 52 \pm 2 \, \mathrm{km}\,\mathrm{s}^{-1}$ which, together with $K_2$, provides a donor star to white dwarf mass ratio $q=M_2/M_1=0.38 \pm 0.03$. We also determine the orbital inclination of the system by modelling the phase-folded ellipsoidal light curve and obtain $i=67^{\circ} \pm 5^{\circ}$. The resulting dynamical masses are $M_{1}=1.03^{+0.16}_{-0.11} \, \mathrm{M}_{\odot}$ and $M_2 = 0.39^{+0.07}_{-0.06} \, \mathrm{M}_{\odot}$ at $68$ per cent confidence level. The white dwarf dynamical mass is compared with estimates obtained by modelling the decline light curve of the $1901$ nova event and X-ray spectroscopy. The best matching mass estimates come from the nova light curve models and the X-ray data analysis by arXiv:1711.01727 that use the ratio between the Alfv\'en radius in quiescence and during dwarf nova outburst.

Using the multi-temperature observations from SDO/AIA on 30th December 2019, we provide a signature of prominence driven forced magnetic reconnection in the corona and associated plasma dynamics during 09:20 UT to 10:38 UT. A hot prominence segment erupts with a speed of 21 km/s and destabilises the entire prominence system. Thereafter, it rose upward in the north during 09:28 UT to 09:48 UT with a speed of 24 km/s. The eruptive prominence stretches overlying field lines upward with the speed of 27-28 km/s , which further undergo into the forced reconnection. The coronal plasma also flows in southward direction with the speed of 7 km/s, and both these inflows trigger the reconnection at 09:48 UT. Thereafter, the east and westward magnetic channels are developed and separated. The east-west reorganization of the magnetic fields starts creating bi-directional plasma outflows towards the limb with their respective speed of 28 km/s and 37 km/s. Their upper ends are diffused in the overlying corona, transporting another set of upflows with the speed of 22 km/s and 19 km/s. The multi-temperature plasma (Te=6.0-7.2) evolves and elongated upto a length of ~10^5 km on the reorganized fields. The hot plasma and remaining prominence threads move from reconnection region towards another segment of prominence in the eastward direction. The prominence-prominence/loop interaction and associated reconnection generate jet-like eruptions with the speed of 178-183 km/s. After the formation of jet, the overlying magnetic channel is disappeared in the corona.

Event Horizon Telescope (EHT) observations of the core of the galaxy M87 suggest an observational appearance dominated by a ring of approximately 40$\mu$as in diameter. The thickness of the ring is less certain: imaging efforts constrained it to be less than half the diameter (consistent with an imaging resolution of 20$\mu$as), while visibility-domain modeling suggested fractional widths of $10$-$20\%$. The fractional width is very interesting as it has the potential to discriminate between different astrophysical scenarios for the source; in fact, the $10$-$20\%$ range is so narrow as to be in tension with theoretical expectations. In the first of a series of papers on the width of the observed ring, we reproduce a subset of EHT visibility-domain modeling results and we explore whether alternative data analysis methods might favor thicker rings. We point out that the closure phase (and closure amplitude) likelihood function is not independent of residual station gain amplitudes, even at high signal-to-noise, and explore two approximations of practical interest: one standard in the field (and employed by the EHT collaboration), and a new one that we propose. Analyzing the public data, we find that the new likelihood approximation prefers somewhat thicker rings, more in line with theoretical expectations. Further analysis is needed, however, to determine which approximation is better for the EHT data.

Using the Keck Planet Imager and Characterizer (KPIC), we obtained high-resolution (R$\sim$35,000) $K$-band spectra of the four planets orbiting HR 8799. We clearly detected \water{} and CO in the atmospheres of HR 8799 c, d, and e, and tentatively detected a combination of CO and \water{} in b. These are the most challenging directly imaged exoplanets that have been observed at high spectral resolution to date when considering both their angular separations and flux ratios. We developed a forward modeling framework that allows us to jointly fit the spectra of the planets and the diffracted starlight simultaneously in a likelihood-based approach and obtained posterior probabilities on their effective temperatures, surface gravities, radial velocities, and spins. We measured $v\sin(i)$ values of $10.1^{+2.8}_{-2.7}$~km/s for HR 8799 d and $15.0^{+2.3}_{-2.6}$~km/s for HR 8799 e, and placed an upper limit of $< 14$~km/s of HR 8799 c. Under two different assumptions of their obliquities, we found tentative evidence that rotation velocity is anti-correlated with companion mass, which could indicate that magnetic braking with a circumplanetary disk at early times is less efficient at spinning down lower mass planets.

This list of contributions to the 37th International Cosmic Ray Conference in Berlin, Germany (12-23 July 2021) summarizes the latest results from the IceCube Neutrino Observatory. IceCube, completed 10 years ago at the geographic South Pole, comprises a surface detector designed to observe cosmic ray air showers, a cubic-kilometer array of optical sensors deployed deep in the ice sheet to observe TeV-PeV neutrinos, and a 15 Megaton deep-ice subdetector sensitive to >10 GeV neutrinos. Data from IceCube are used to investigate a broad set of key questions in physics and astrophysics, such as the origins of galactic and extragalactic cosmic rays, the fundamental properties of neutrinos, and searches for physics beyond the Standard Model. The papers in this index are grouped topically to highlight IceCube contributions related to neutrino and multi-messenger astrophysics, cosmic-ray physics, fundamental physics, education and public outreach, and research and development for next-generation neutrino observatories. Contributions related to IceCube-Gen2, the future extension of IceCube, are available in a separate collection.

The electron density ($n_{e^{-}}$) plays an important role in setting the chemistry and physics of the interstellar medium. However, measurements of $n_{e^{-}}$ in neutral clouds have been directly obtained only toward a few lines of sight or they rely on indirect determinations. We use carbon radio recombination lines and the far-infrared lines of C$^{+}$ to directly measure $n_{e^{-}}$ and the gas temperature in the envelope of the integral shaped filament (ISF) in the Orion A molecular cloud. We observed the C$102\alpha$ and C$109\alpha$ carbon radio recombination lines (CRRLs) using the Effelsberg 100m telescope at ~2' resolution toward five positions in OMC-2 and OMC-3. Since the CRRLs have similar line properties, we averaged them to increase the signal-to-noise ratio of the spectra. We compared the intensities of the averaged CRRLs, and the 158 {\mu}m-[CII] and [$^{13}$CII] lines to the predictions of a homogeneous model for the C$^{+}$/C interface in the envelope of a molecular cloud and from this comparison we determined the electron density, temperature and C$^{+}$ column density of the gas. We detect the CRRLs toward four positions, where their velocity and widths (FWHM 2.3 km s$^{-1}$) confirms that they trace the envelope of the ISF. Toward two positions we detect the CRRLs, and the [CII] and [$^{13}$CII] lines with a signal-to-noise ratio >5, and we find $n_{e^{-}}=0.65\pm0.12$ cm$^{-3}$ and $0.95\pm0.02$ cm$^{-3}$, which corresponds to a gas density $n_{H}\approx5\times10^{3}$ cm$^{-3}$ and a thermal pressure of $p_{th}\approx4\times10^{5}$ K cm$^{-3}$. We also constrained the ionization fraction in the denser portions of the molecular cloud using the HCN(1-0) and C$_{2}$H(1-0) lines to $x(e^{-})<3\times10^{-6}$. The derived electron densities and ionization fraction imply that $x(e^{-})$ drops by a factor >100 between the C$^{+}$ layer and the regions probed by HCN(1-0).

IceCube-Gen2 is a planned extension of the IceCube Neutrino Observatory at the South Pole. The extension is optimized to search for sources of astrophysical neutrinos from TeV to EeV, and will improve the sensitivity of the observatory to neutrino point sources by a factor of five. The science case of IceCube-Gen2 is built on a successful decade of observations with IceCube. This index of contributions to the 37th International Cosmic Ray Conference in Berlin, Germany (12-23 July 2021) describes research and development efforts for IceCube-Gen2. Included are performance studies of next-generation optical sensors that will detect Cherenkov radiation from TeV-PeV cosmic rays and neutrinos; optimizations of the geometries of the surface and in-ice optical arrays; and estimates of the sensitivity of the proposed IceCube-Gen2 radio array to Askaryan emission from PeV-EeV neutrinos. Contributions related to the existing instrument, IceCube, are available in a separate collection.

Astrophysical plasma codes are built on atomic databases. In the current atomic databases, R-matrix electron-impact excitation data of O-like ions are limited. The accuracy of plasma diagnostics with O-like ions depends on the availability and accuracy of the atomic data. This is particularly relevant in the context of future observatories equipped with the next generation of high-resolution spectrometers. To obtain level-resolved effective collision strengths of O-like ions from \ion{Ne}{III} to \ion{Zn}{XXIII} (i.e., Ne$^{2+}$ to Zn$^{22+}$) over a wide range of temperatures. This includes transitions up to $nl=5d$ for each ion. We also aim to assess the accuracy of the new data, as well as their impact on solar atmosphere plasma diagnostics, compared to those available within the CHIANTI database. A large-scale R-matrix intermediate coupling frame transformation calculations were performed systematically for the O-like iso-electronic sequence. For each ion, 630 fine-structure levels were included in both the configuration interaction target and close-coupling collision expansions. The present results (energy levels, oscillator strengths, and effective collision strengths) of selected ions across the iso-electronic sequence are compared with those in archival databases and the literature. For selected ions across the iso-electronic sequence. We find general agreement with the few previous R-matrix calculations of collision strengths. We illustrate the improvements for a few solar plasma diagnostics over existing CHIANTI atomic models based on distorted wave data. The electron-impact excitation data are archived according to the Atomic Data and Analysis Structure (ADAS) data class adf04 and will be available in OPEN-ADAS.

As a result of various studies, it has been determined that several post-common envelope eclipsing binaries have variations in their orbital periods. These variations are thought to be caused by the existence of additional bodies in the system (hypothetical stars or planets) and/or other physical effects (such as angular momentum loss, magnetic activity) of the binary system. It is also known that the sdB+M eclipsing system NY Vir has shown such variations in the last decade, indicating additional objects and/or other physical effects. In this work, we present 51 new eclipse times for this system, which extend the time span of it is $O-C$ diagram by about three years, obtained between 2015 and 2021 using two different telescopes in Turkey. The data obtained in the last 3 years shows a new trend in the $O-C$ diagram differently from the predictions of the previous studies. Our model is consistent with the new $O-C$ diagram, which is statistically well fitted with the quadratic term and the additional two planets with masses of $M_3=2.74 \: M_\text{Jup}$ and $M_4=5.59 \: M_\text{Jup}$. However, the orbital period variation can also be related to magnetic activity. In order to better understand the mechanism causing the changes in the orbital period, new observation data is needed that will show at least one full cycle of the change in the O-C diagram.

Empirical trends in stellar X-ray and radio luminosities suggest that low mass ultracool dwarfs (UCDs) should not produce significant radio emission. Defying these expectations, strong non-thermal emission has been observed in a few UCDs in the 1-10 GHz range, with a variable component often attributed to global aurorae and a steady component attributed to other processes such as gyrosynchrotron emission. While both auroral and gyrosynchrotron emission peak near the critical frequency, only the latter radiation is expected to extend into millimeter wavelengths. We present ALMA 97.5 GHz and VLA 33 GHz observations of a small survey of 5 UCDs. LP 349-25, LSR J1835+3259, and NLTT 33370 were detected at 97.5 GHz, while LP 423-31 and LP 415-20 resulted in non-detections at 33 GHz. A significant flare was observed in NLTT 33370 that reached a peak flux of 4880 +/- 360 microJy, exceeding the quiescent flux by nearly an order of magnitude, and lasting 20 seconds. These ALMA observations show bright 97.5 GHz emission with spectral indices ranging from alpha = -0.76 to alpha = -0.29, suggestive of optically thin gyrosynchrotron emission. If such emission traces magnetic reconnection events, then this could have consequences for both UCD magnetic models and the atmospheric stability of planets in orbit around them. Overall, our results provide confirmation that gyrosynchrotron radiation in radio loud UCDs can remain detectable into the millimeter regime.

Using numerical simulations, we predict sky maps and light curves of gamma-ray emission from neutron stars in compact binaries, and in isolation. We briefly review some gamma-ray emission models, and reproduce sky maps from a standard isolated pulsar in the Separatrix Layer model. We then simulate isolated pulsars with several variations of a dipole magnetic field, including superpositions, and predict their gamma-ray emission. These simulations provide new heuristics on what can and cannot be inferred about the magnetic field configuration of pulsars from high-energy observations. We find that typical double-peak light curves can be produced by pulsars with significant multipole structure beyond a single dipole. We offer a simple approximation that is useful for rapid explorations of binary magnetic field structure. Finally, we predict the gamma-ray emission pattern from a compact black hole-neutron star binary moments before merger by applying the Separatrix Layer model to data simulated in full general relativity; we find that face-on observers receive little emission, equatorial observers see one broad peak, and more generic observers typically see two peaks.

Super-puffs -- low-mass exoplanets with extremely low bulk density -- are attractive targets for exploring their atmospheres and formation processes. Recent studies suggested that the large radii of super-puffs may be caused by atmospheric dust entrained in the escaping atmospheres. In this study, we investigate how the dust grows in escaping atmospheres and influence the transit radii using a microphysical model of grain growth. Collision growth is efficient in many cases, leading to hinder the upward transport of dust via enhanced gravitational settling. We find that dust abundance in the outflow hardly exceeds the Mach number at the dust production region. Thus, dust formed at upper atmospheres, say $P\lesssim{10}^{-5}$ bar, are needed to launch a dusty outflow with high dust abundance. With sufficiently high dust production altitudes and rates, the dusty outflow can enhance the observable radius by a factor of $\sim$2 or even more. We suggest that photochemical haze is a promising candidate of high-altitude dust that can be entrained in the outflow. We also compute the synthetic transmission spectra of super-puff atmospheres and demonstrate that the dusty outflow produces a broad spectral slope and obscures molecular features, in agreement with recently reported featureless spectra. Lastly, using an interior structure model, we suggest that the atmospheric dust could drastically enhance the observable radius only for planets in a narrow mass range of $\sim2$--$5M_{\rm \oplus}$, in which the boil-off tends to cause total atmospheric loss. This may explain why super-puffs are uncommon despite the suggested universality of photochemical hazes.

Using a filter in the GROWTH Marshal based on color and the amplitude and the timescale of variability, we have identified 372 objects as known or candidate cataclysmic variables (CVs) during the second year of operation of the Zwicky Transient Facility (ZTF). From the available difference imaging data, we found that 93 are previously confirmed CVs, and 279 are strong candidates. Spectra of four of the candidates confirm them as CVs by the presence of Balmer emission lines, while one of the four has prominent HeII lines indicative of containing a magnetic white dwarf. Gaia EDR3 parallaxes are available for 154 of these systems, resulting in distances from 108-2096 pc and absolute magnitudes in the range of 7.5-15.0, with the largest number of candidates between 10.5-12.5. The total numbers are 21% higher than from the previous year of the survey with a greater number of distances available but a smaller percentage of systems close to the Galactic plane. Comparison of these findings with a machine learning method of searching all the light curves reveals large differences in each dataset related to the parameters involved in the search process.

We found that the 1899 nova V606 Aql currently shows dwarf nova outbursts with a typical cycle length of 270 d and amplitudes of ~1.5 mag using Public Data Release of Zwicky Transient Facility observations. The low mass-transfer rate in quiescence has been suggested to explain the large eruption amplitude (Tappert et al., 2016), and the present detection of dwarf nova outbursts supports this interpretation. The transition to the dwarf nova state more than 100 yr after the nova eruption gives credence to the hibernation scenario. The absolute magnitude estimated from dwarf nova outbursts suggests that V606 Aql should have been a fast nova and the presence of high excitation lines in quiescence would be explained by the presence of a massive white dwarf.

The recent increase in well-localised fast radio bursts (FRBs) has facilitated in-depth studies of global FRB host properties, the source circumburst medium, and the potential impacts of these environments on the burst properties. The Australian Square Kilometre Array Pathfinder (ASKAP) has localised 11 FRBs with sub-arcsecond to arcsecond precision, leading to sub-galaxy localisation regions in some cases and those covering much of the host galaxy in others. The method used to astrometrically register the FRB image frame for ASKAP, in order to align it with images taken at other wavelengths, is currently limited by the brightness of continuum sources detected in the short-duration ('snapshot') voltage data captured by the Commensal Real-Time ASKAP Fast Transients (CRAFT) software correlator, which are used to correct for any frame offsets due to imperfect calibration solutions and estimate the accuracy of any required correction. In this paper, we use dedicated observations of bright, compact radio sources in ASKAP's low- and mid-frequency bands to investigate the typical astrometric accuracy of the positions obtained using this so-called 'snapshot' technique. Having captured these data with both the CRAFT software and ASKAP hardware correlators, we also compare the offset distributions obtained from both data products to estimate a typical offset between the image frames resulting from the differing processing paths, laying the groundwork for future use of the longer-duration, higher signal-to-noise ratio data recorded by the hardware correlator. We find typical offsets between the two frames of $\sim 0.6$ and $\sim 0.3$ arcsec in the low- and mid-band data, respectively, for both RA and Dec. We also find reasonable agreement between our offset distributions and those of the published FRBs. <Abridged>

Photochemical hazes are important opacity sources in temperate exoplanet atmospheres, hindering current observations from characterizing exoplanet atmospheric compositions. The haziness of an atmosphere is determined by the balance between haze production and removal. However, the material-dependent removal physics of the haze particles is currently unknown under exoplanetary conditions. Here we provide experimentally-measured surface energies for a grid of temperate exoplanet hazes to characterize haze removal in exoplanetary atmospheres. We found large variations of surface energies for hazes produced under different energy sources, atmospheric compositions, and temperatures. The surface energies of the hazes were found to be the lowest around 400 K for the cold plasma samples, leading to the lowest removal rates. We show a suggestive correlation between haze surface energy and atmospheric haziness with planetary equilibrium temperature. We hypothesize that habitable zone exoplanets could be less hazy, as they would possess high-surface-energy hazes which can be removed efficiently.

We studied Zwicky Transient Facility (ZTF) light curves of 34 dwarf nova candidates discovered by All-Sky Automated Survey for Supernovae (ASAS-SN) between 2020 May 12 and September 9 and found 6 AM CVn-type candidates. All objects showed short outbursts (post-superoutburst rebrightenings) on the fading tail. Two objects (ASASSN-20eq, ASASSN-20la) showed double superoutbursts. Three objects (ASASSN-20jt, ASASSN-20ke, and ASASSN-20lr) showed short superoutbursts (5-6 d). These features in the light curve can be used in discriminating AM CVn-type candidates from hydrogen-rich systems. In contrast to hydrogen-rich systems, some object did not show red color excess during the rebrightening or fading tail phase. We interpret that this is due to the higher ionization temperature in helium disks. Two objects had long (likely) supercycles: ASASSN-20gx (8.5 yr) and ASASSN-20lr (7 yr). We provide a scheme for identifying AM CVn-type candidates based on the light curve characteristics.

Ultramassive white dwarfs are extreme endpoints of stellar evolution. Recent findings, such as a missing multi-Gyr cooling delay for a number of ultramassive white dwarfs and a white dwarf with a quasi-Chandrasekhar mass, motivate a better understanding of their evolution. A key process still subject to important uncertainties is the crystallization of their dense cores, which are generally assumed to be constituted of $^{16}$O, $^{20}$Ne, and a mixture of several trace elements (most notably $^{23}$Na and $^{24}$Mg). In this work, we use our recently developed Clapeyron integration technique to compute accurate phase diagrams of three-component mixtures relevant to the modeling of O/Ne ultramassive white dwarfs. We show that, unlike the phase separation of $^{22}$Ne impurities in C/O cores, the phase separation of $^{23}$Na impurities in O/Ne white dwarfs cannot lead to the enrichment of their cores in $^{23}$Na via a distillation process. This severely limits the prospect of transporting large quantities of $^{23}$Na toward the center of the star, as needed in the white dwarf core collapse mechanism recently proposed by Caiazzo et al. We also show that despite representing $\approx 10\%$ of the ionic mixture, $^{23}$Na and $^{24}$Mg impurities only have a negligible impact on the O/Ne phase diagram, and the two-component O/Ne phase diagram can be safely used in white dwarf evolution codes. We provide analytic fits to our high-accuracy O/Ne phase diagram for implementation in white dwarf models.

Giant radio galaxies provide important clues into the life cycles and triggering mechanisms of radio jets. With large-scale jets spanning 1.8 Mpc, ESO 422-G028 ($z = 0.038$) is a giant radio galaxy that also exhibits signs of restarted jet activity in the form of pc-scale jets. We present a study of the spatially-resolved stellar and gas properties of ESO 422-G028 using optical integral field spectroscopy from the WiFeS spectrograph. In addition to the majority $\sim 13\,\rm Gyr$ old stellar population, ESO 422-G028 exhibits a much younger ($\lesssim 10\,\rm Myr$ old) component with an estimated mass of $ 10^{7.6}\,\rm M_\odot$ which is predominantly located in the North-West region of the galaxy. Unusually, the ionised gas kinematics reveal two distinct disks traced by narrow ($\sigma_{\rm H\alpha} < 100 \,\rm km\,s^{-1}$) and broad ($\sigma_{\rm H\alpha} > 150 \,\rm km\,s^{-1}$) H$\alpha$ emission respectively. Both ionised gas disks are misaligned with the axis of stellar rotation, suggesting an external origin. This is consistent with the prominent interstellar Na D absorption, which traces a $1 - 3 \,\rm M_\odot \, yr^{-1}$ inflow of neutral gas from the North. We posit that an inflow of gas - either from an accretion event or a gas-rich merger - has triggered both the starburst and the restarted jet activity, and that ESO 422-G028 is potentially on the brink of an epoch of powerful AGN activity.

The Hubble tension might be resolved by injecting a new energy component, called Early Dark Energy (EDE), prior to recombination. An Anti-de Sitter (AdS) phase around recombination can make the injected energy decay faster, which thus allows a higher EDE fraction (so larger $H_0$) while prevents degrading the CMB fit. In this work, we test the AdS-EDE model with CMB and Large-Scale Structure (LSS) data. Our CMB dataset consists of low-$\ell$ part of Planck TT spectrum and SPTpol polarization and lensing measurements, since this dataset predicts the CMB lensing effect consistent with $\Lambda$CDM expectation. Combining it with BAO and Pantheon data, we find the bestfit values $H_0=71.92$ km/s/Mpc and $H_0=73.29$ km/s/Mpc without and with the SH0ES prior, respectively. Including cosmic shear and galaxy clusters data, we have $H_0=71.87$ km/s/Mpc and $S_8=0.785$, i.e. only $1.3\sigma$ discrepancy with direct $S_8$ measurement.

We present ${\tt nimbus}$ : a hierarchical Bayesian framework to infer the intrinsic luminosity parameters of kilonovae (KNe) associated with gravitational-wave (GW) events, based purely on non-detections. This framework makes use of GW 3-D distance information and electromagnetic upper limits from a given survey for multiple events, and self-consistently accounts for finite sky-coverage and probability of astrophysical origin. The framework is agnostic to the brightness evolution assumed and can account for multiple electromagnetic passbands simultaneously. Our analyses highlight the importance of accounting for model selection effects, especially in the context of non-detections. We show our methodology using a simple, two-parameter linear brightness model, taking the follow-up of GW190425 with the Zwicky Transient Facility (ZTF) as a single-event test case for two different prior choices of model parameters -- (i) uniform/uninformative priors and (ii) astrophysical priors based on surrogate models of Monte Carlo radiative transfer simulations of KNe. We present results under the assumption that the KN is within the searched region to demonstrate functionality and the importance of prior choice. Our results show consistency with ${\tt simsurvey}$ -- an astronomical survey simulation tool used previously in the literature to constrain the population of KNe. While our results based on uniform priors strongly constrain the parameter space, those based on astrophysical priors are largely uninformative, highlighting the need for deeper constraints. Future studies with multiple events having electromagnetic follow-up from multiple surveys should make it possible to constrain the KN population further.

Galaxy morphologies, kinematics, and stellar populations are thought to link to each other. However, both simulations and observations have pointed out mismatches therein. In this work, we study the nature and origin of the present-day quenched, bulge-dominated, but dynamically cold galaxies within a stellar mass range of $10.3<\log\,M_{\ast}/\mathrm{M_{\odot}}<11.2$ in the IllustrisTNG-100 Simulation, as a companion paper of Lu et al.(2021), which aimed at the star-forming but dynamically hot disc galaxies within a lower stellar mass range of $9.7<\log\,M_{\ast}/\mathrm{M_{\odot}}<10.3$. We compare cold quenched population with a population of normal star-forming dynamically cold disc galaxies and a population of normal quenched dynamically hot elliptical galaxies within the same mass range. The populations of the present-day quenched and bulge-dominated galaxies (both being dynamically cold and hot) used to have significantly higher star-formation rates and thinner morphologies at redshifts of z~2. They have experienced more frequent larger mass-ratio mergers below z~0.7 in comparison to their star-forming disc counterparts, which is responsible for the formation of their bulge-dominated morphologies. The dynamically cold populations (both being star-forming and quenched) have experienced more frequent prograde and tangential mergers especially below z~1, in contrast to the dynamically hot ellipticals, which have had more retrograde and radial mergers. Such different merging histories can well explain the differences on the cold and hot dynamical status among these galaxies. We point out that the real-world counterparts of these dynamically cold and hot bulge-dominated quenched populations are the fast- and slow-rotating early-type galaxies, respectively, as seen in observations and hence reveal the different evolution paths of these two distinct populations of early-type galaxies.

Smoothed particle hydrodynamics (SPH) is positioned as having ideal conservation properties. When properly implemented, conservation of total mass, energy, and both linear and angular momentum is guaranteed exactly, up to machine precision. This is particularly important for some applications in computational astrophysics, such as binary dynamics, mergers, and accretion of compact objects (neutron stars, black holes, and white dwarfs). However, in astrophysical applications that require the inclusion of gravity, calculating pairwise particle interactions becomes prohibitively expensive. In the Fast Multipole Method (FMM), they are, therefore, replaced with symmetric interactions between distant clusters of particles (contained in the tree nodes) Although such an algorithm is linear momentum-conserving, it introduces spurious torques that violate conservation of angular momentum. We present a modification of FMM that is free of spurious torques and conserves angular momentum explicitly. The new method has practically no computational overhead compared to the standard FMM.

A new mission about twenty years after Gaia with similar astrometric performance would be important for all branches of astronomy. The two missions together would, e.g., give much more accurate motions of the common objects due to the large epoch difference. By adding a Near-InfraRed (NIR) capability to the new mission we will be able to peer into the obscured regions of the Galaxy and measure up to 10 or 12 billion new objects and reveal many new sciences in the process. ESA has now ranked the development of this mission so high that a launch about 2045 is quite probable. A brief history of the project is included.

In this paper, we carry out multiwavelength observations of three recurring jets on 2014 November 7. The jets originated from the same region at the edge of AR 12205 and propagated along the same coronal loop. The eruptions were generated by magnetic reconnection, which is evidenced by continuous magnetic cancellation at the jet base. The projected initial velocity of the jet2 is 402 km s. The accelerations in the ascending and descending phases of jet2 are not consistent, the former is considerably larger than the value of solar gravitational acceleration at the solar surface, while the latter is lower than solar gravitational acceleration. There are two possible candidates of extra forces acting on jet2 during its propagation. One is the downward gas pressure from jet1 when it falls back and meets with jet2. The other is the viscous drag from the surrounding plasma during the fast propagation of jet2. As a contrast, the accelerations of jet3 in the rising and falling phases are constant, implying that the propagation of jet3 is not significantly influenced byextra forces.

Binary star systems represent a significant proportion of the Galactic stellar population, with X-ray binaries being an important subset of these for high energy astrophysics. Although hundreds of X-ray binaries are detected in the Milky Way and beyond, only 12 of these systems are listed in the 4FGL-DR2, the latest Fermi-LAT point source catalogue. With such a small number detectable by Fermi-LAT, much is still unknown about the mechanisms by which these systems emit $\gamma$-rays. We present the method and current status of our large-scale survey of the X-ray binary population using over 12 years of Fermi-LAT data, and current catalogues and background models.

Helioseismology is the study of the solar interior using observations of oscillations at the surface. It suffers from systematic errors, such as a center-to-limb error in travel-time measurements. Understanding these errors requires a good understanding of the nontrivial relationship between wave displacement and helioseismic observables. The wave displacement causes perturbations in the atmospheric thermodynamical quantities which perturb the opacity, the optical depth, the source function, and the local ray geometry, thus affecting the emergent intensity. We aim to establish the most complete relationship up to now between the displacement and the intensity perturbation by solving the radiative transfer problem in the atmosphere. We derive an expression for the intensity perturbation caused by acoustic oscillations at any point on the solar disk by applying the first-order perturbation theory. As input, we consider adiabatic modes of oscillation of different degrees. The background and the perturbed intensities are computed considering the main sources of opacity in the continuum. We find that, for all modes, the perturbations to the thermodynamical quantities are not sufficient to model the intensity. In addition, the geometrical effects due to the displacement must be taken into account as they lead to a difference in amplitude and a phase shift between the temperature at the surface and intensity perturbations. The closer to the limb, the larger the differences. This work presents improvements for the computation of the intensity perturbations, in particular for high-degree modes, and explains differences in intensity computations in earlier works. The phase shifts and amplitude differences between the temperature and intensity perturbations increase towards the limb. This should help to interpret some of the systematic center-to-limb effects observed in local helioseismology.

We present a new and very fast method for producing microlensing magnification maps at high optical depths. It is based on the combination of two approaches: (a) the two-dimensional Poisson solver for a deflection potential and (b) inverse polygon mapping. With our method we extremely reduce the computing time for the generation of magnification patterns and avoid the use of highly demanding computer resources. For example, the generation of a magnification map of size 2000 x 2000 pixels, covering a region of 20 Einstein radii, takes a few seconds on a state-of-the-art laptop. The method presented here will facilitate the massive production of magnification maps for extragalactic microlensing studies within the forthcoming surveys without the need for large computer clusters. The modest demand of computer power and a fast execution time allow the code developed here to be placed on a standard server and thus provide the public online access through a web-based interface.

Very precise observational data are needed for studying the stellar cluster parameters (distance, reddening, age, metallicity) and cluster internal kinematics. In turn, these give us an insight into the properties of our Galaxy, for example, by giving us the ability to trace Galactic spiral structure, star formation rates and metallicity gradients. We investigated the available Gaia DR2 catalogue of 1229 open clusters and studied cluster distances, sizes and membership distributions in the 3D space. An appropriate analysis of the parallaxto-distance transformation problem is presented in the context of getting distances toward open clusters and estimating their sizes. Based on our investigation of the Gaia DR2 data we argue that, within 2 kpc, the inverse-parallax method gives comparable results (distances and sizes) as the Bayesian approach based on the exponentially decreasing volume density prior. Both of these methods show very similar dependence of the line-of-sight elongation of clusters (needle-like shapes resulting from the parallax uncertainties) on the distance. We also looked at a measure of elongations of the studied clusters and find the maximum distance of 665 pc at which a spherical fit still contains about half of the stellar population of a cluster. It follows from these results that the 3D structure of an open cluster cannot be properly studied beyond about 500 pc when using any of mentioned standard transformations of parallaxes to distances.

Recent advances in techniques of critical close reading of historical texts can now be applied to records of pre-telescopic celestial observations - allowing significant progress for analyzing and solving orbits of past comets: we exemplify our method by solving the orbit of the comet in AD 760 only with historical observations and then identify it with 1P/Halley. A detailed eyewitness record with drawing of a comet in AD 760 in the Syriac Chronicle of Zuqnin (finished AD 775/6) was not yet included in the study of its orbit - the Chinese reports alone do not yield a sufficient number of dated positions. We analyze the Syriac and Chinese sources with critical methods for quantitative astronomical usage, we also consider a few further records from the Mediterranean and West Asian area. With our conservatively derived dated positions we can determine the best fitting Keplerian orbital solution by least squares fitting yielding the orbital elements; the parameter ranges for non-periodic solutions and highly eccentric periodic solutions are consistent with each other. The allowed parameter ranges for perihelion distance and inclination are sufficiently small to identify the comet with 1P/Halley. Although 1P/Halley is the only comet, where the telescopic orbit is credibly linked to pre-telescopic returns, e.g. to AD 760, our identification confirms claims from extrapolating telescopic observations backward in time - here independently based on historical data. In particular, we obtained a precise perihelion time (760 May 19.1 \pm 1.7). The inferior conjunction between comet and Sun as on the previously published orbit (760 May 31.9) is shifted by about one day compared to our new orbit (June 1.8), only the new one is consistent with the last observation (June 1.0) before conjunction as reported in the Chronicle of Zuqnin. We also study the comet's brightness evolution ...

Studies of planet migration derived from disc planet interactions began before the discovery of exoplanets. The potential importance of migration for determining orbital architectures being realised, the field received greater attention soon after the initial discoveries of exoplanets. Early studies based on very simple disc models indicated very fast migration times for low mass planets that raised questions about its relevance. However, more recent studies, made possible with improving resources, that considered improved physics and disc models revealed processes that could halt or reverse this migration. That in turn led to a focus on special regions in the disc where migration could be halted. In this way the migration of low mass planets could be reconciled with formation theories. In the case of giant planets which have a nonlinear interaction with the disc, the migration should be slower and coupled to the evolution of the disc. The latter needs to be considered more fully to make future progress in all cases. Here we are primarily concerned with processes where migration is connected with the presence of the protopolanetary disk. Migration may also be induced by disc-free gravitational interactions amongst planets or with binary companions. This is only briefly discussed here.

As a guide for astronomers new to the field of technosignature search (i.e. SETI), I present an overview of some of its observational and theoretical approaches. I review some of the various observational search strategies for SETI, focusing not on the variety of technosignatures that have been proposed or which are most likely to be found, but on the underlying philosophies that motivate searches for them. I cover passive versus active searches, ambiguous versus dispositive kinds of technosignatures, commensal or archival searches versus dedicated ones, communicative signals versus "artifacts", "active" versus derelict technologies, searches for beacons versus eavesdropping, and model-based versus anomaly-based searches. I also attempt to roughly map the landscape of technosignatures by kind and the scale over which they appear. I also discuss the importance of setting upper limits in SETI, and offer a heuristic for how to do so in a generic SETI search. I mention and attempt to dispel several misconceptions about the field. I conclude with some personal observations and recommendations for how to practice SETI, including how to choose good theory projects, how to work with experts and skeptics to improve one's search, and how to plan for success.

A cosmic string wake produces a distinct non-Gaussian signal in 21-cm intensity maps at redshifts above that of reionization. While the string signal is (locally) larger in amplitude than the signal of the Gaussian fluctuations of the $\Lambda$CDM model, they are overwhelmed (even locally in position space) by astrophysical and instrumental foregrounds. Here, we study to what extent the signal can be extracted from noisy interferometric data. The narrowness of the string-induced feature in redshift direction allows for a subtraction of astrophysical and instrumental foregrounds. Based on the specific geometry of the string signal we identify a particular three-point statistic which is promising in order to extract the signal, and we find that, having in mind a telescope of specifications similar to that of the MWA instrument, the string signal can be successfully extracted for a value of the string tension of $G\mu = 3 \times 10^{-7}$. Prospects for further improvements of the analysis are discussed.

Local variations in the intergalactic medium (IGM) neutral hydrogen fraction will affect the Ly-$\alpha$ absorption signature of protoclusters identified in tomographic surveys. Using the IllustrisTNG simulations, we investigate how the AGN proximity effect and hot, collisionally ionised gas arising from gravitational infall and black hole feedback changes the Ly-$\alpha$ absorption associated with $M_{z=0}\simeq10^{14}\,M_\odot$ protoclusters at $z\simeq2.4$. We find that protocluster galaxy overdensities exhibit a weak anti-correlation with Ly-$\alpha$ transmission in IGM transmission maps, but local HI ionisation enhancements due to hot $T>10^{6}\rm\,K$ gas or nearby AGN can disrupt this relationship within individual protoclusters. On average, however, we find that strong reductions in the IGM neutral fraction are limited to within $\lesssim 5h^{-1}\,\textrm{cMpc}$ of the dark matter haloes. Local ionisation enhancements will therefore have a minimal impact on the completeness of protocluster identification in tomographic surveys if smoothing Ly-$\alpha$ transmission maps over scales of $\sim4 h^{-1}\,\textrm{cMpc}$, as is typically done in observations. However, if calibrating the relationship between the matter density and Ly-$\alpha$ transmission in tomographic maps using simple analytical models for the Ly-$\alpha$ forest opacity, the presence of hot gas around haloes can still result in systematically lower estimates of $M_{z=0}$ for the most massive protoclusters.

In light of our previous work \cite{Liu:2019xhn}, we investigate the possibility of formation for primordial black-hole during preheating period, in which we have implemented the instability of the Mathieu equation. For generating sufficient enough enhanced power spectrum, we choose some proper parameters belonging to the narrow resonance. To characterize the full power spectrum, the enhanced part of the power spectrum is depicted by the $\delta$ function at some specific scales, which is highly relevant with the mass of inflaton due to the explicit coupling between the curvaton and inflaton. After the inflationary period (including the preheating period), there is only one condition satisfying with the COBE normalization upper limit. Thanks to the huge choices for this mass parameter, we can simulate the value of abundance of primordial black holes nearly covering all of the mass ranges, in which we have given three special cases. One case could account for the dark matter in some sense since the abundance of a primordial black hole is about $75\%$. At late times, the relic of exponential potential could be approximated to a constant of the order of cosmological constant dubbed as a role of dark energy. Thus, our model could unify dark energy and dark matter from the perspective of phenomenology. Finally, it sheds new light for exploring Higgs physics.

Measurements of the gravitational-wave signals from neutron star mergers allow scientists to learn about the interior of neutron stars and the properties of dense nuclear matter. The study of neutron star mergers is usually performed with computational fluid dynamics codes, mostly in Eulerian but also in Lagrangian formulation such as smoothed particle hydrodynamics (SPH). Codes include our best knowledge of nuclear matter in the form of an equation of state as well as effects of general relativity (GR). However, one important aspect of neutron stars is usually ignored: the solid nature of their crust. The solid matter in the crust is the strongest material known in nature which could lead to a multitude of possible observational effects that have not been studied in dynamical simulations yet. The crust could change the way a neutron star deforms during a merger, leaving an imprint in the gravitational wave signal. It could even shatter during the inspiral, producing a potentially observable electromagnetic signal. Here, we present a first study of the dynamical behavior of neutron stars with a solid crust and fixed GR background with FleCSPH. FleCSPH is a general-purpose SPH code, developed at Los Alamos National Laboratory. It features an efficient algorithm for gravitational interactions via the Fast Multipole Method, which, together with the implemented nuclear equation of state, makes it appropriate for astrophysical applications. The solid material dynamics is described via the elastic-perfectly plastic model with maximum-strain breaking. Despite its simplicity, the model reproduces the stress-strain curve of crustal material as extracted from microphysical simulations very well. We present first tests of our implementation via simulations of neutron star oscillations and give an outlook on our study of the dynamical behavior of the solid crust in neutron star merger events.

It was suggested that the $\gamma$-ray halo around Geminga might not be interpreted by slow-diffusion. If the ballistic regime of electron/positron propagation is considered, the Geminga halo may be explained even with a large diffusion coefficient. In this work, we examine this effect by taking the generalized J\"uttner propagator as the approximate relativistic Green's function for diffusion and find that the morphology of the Geminga halo can be marginally fitted in the fast-diffusion scenario. However, the recently discovered $\gamma$-ray halo around PSR J0622$+$3749 at LHAASO cannot be explained by the same effect and slow diffusion is the only solution. Furthermore, both the two pulsar halos require a conversion efficiency from the pulsar spin-down energy to the high energy electrons/positrons much larger than 100\%, if they are interpreted by this ballistic transport effect. Therefore, we conclude that slow diffusion is necessary to account for the $\gamma$-ray halos around pulsars.

Using photometric ULTRACAM observations of three new short period cataclysmic variables, we model the primary eclipse lightcurves to extract the orbital separation, masses, and radii of their component stars. We find donor masses of 0.060 +/- 0.008 solar masses, 0.042 +/- 0.001 solar masses, and 0.042 +/- 0.004 solar masses, two being very low-mass sub-stellar donors, and one within 2 sigma of the hydrogen burning limit. All three of the new systems lie close to the modified, "optimal" model evolutionary sequence of Knigge et al. (2011). We briefly re-evaluate the long-standing discrepancy between observed donor mass and radius data, and theoretical CV evolutionary tracks. By looking at the difference in the observed period at each mass and the period predicted by the Knigge et al. (2011) evolutionary sequence, we qualitatively examine the form of excess angular momentum loss that is missing from the models below the period gap. We show indications that the excess angular momentum loss missing from CV models grows in importance relative to gravitational losses as the period decreases. Detailed CV evolutionary models are necessary to draw more quantitative conclusions in the future.

We present a comparison of simulation-based inference to full, field-based analytical inference in cosmological data analysis. To do so, we explore parameter inference for two cases where the information content is calculable analytically: Gaussian random fields whose covariance depends on parameters through the power spectrum; and correlated lognormal fields with cosmological power spectra. We compare two inference techniques: i) explicit field-level inference using the known likelihood and ii) implicit likelihood inference with maximally informative summary statistics compressed via Information Maximising Neural Networks (IMNNs). We find that a) summaries obtained from convolutional neural network compression do not lose information and therefore saturate the known field information content, both for the Gaussian covariance and the lognormal cases, b) simulation-based inference using these maximally informative nonlinear summaries recovers nearly losslessly the exact posteriors of field-level inference, bypassing the need to evaluate expensive likelihoods or invert covariance matrices, and c) even for this simple example, implicit, simulation-based likelihood incurs a much smaller computational cost than inference with an explicit likelihood. This work uses a new IMNNs implementation in JAX that can take advantage of fully-differentiable simulation and inference pipeline. We also demonstrate that a single retraining of the IMNN summaries effectively achieves the theoretically maximal information, enhancing the robustness to the choice of fiducial model where the IMNN is trained.

We present a new generation of substellar atmosphere and evolution models, appropriate for application to studies of L, T, and Y-type brown dwarfs and self-luminous extrasolar planets. The atmosphere models describe the expected temperature-pressure profiles and emergent spectra of atmospheres in radiative-convective equilibrium with effective temperatures and gravities within the ranges $200\le T_{\rm eff}\le2400\,\rm K$ and $2.5\le \log g \le 5.5$. These ranges encompass masses from about 0.5 to 85 Jupiter masses for a set of metallicities ($[{\rm M/H}] = -0.5$ to $+0.5$), C/O ratios (from 0.5 to 1.5 times that of solar), and ages. The evolution tables describe the cooling of these substellar objects through time. These models expand the diversity of model atmospheres currently available, notably to cooler effective temperatures and greater ranges in C/O. Notable improvements from past such models include updated opacities and atmospheric chemistry. Here we describe our modeling approach and present our initial tranche of models for cloudless, chemical equilibrium atmospheres. We compare the modeled spectra, photometry, and evolution to various datasets.

Surprisingly strong CO emission has been observed from more than a dozen debris disks around nearby main-sequence stars. The origin of this CO is unclear, in particular whether it is left over from the protoplanetary disk phase or is second-generation material released from collisions between icy bodies like debris dust. The primary unexplored avenue for distinguishing the origin of the material is understanding its molecular composition. Here we present a deep search for five molecules (CN, HCN, HCO+, SiO, and CH3OH) in the debris disk around 49 Ceti. We take advantage of the high sensitivity of the Atacama Large Millimeter/submillimeter Array (ALMA) at Band 7 to integrate for 3.2 hours at modest spatial (1") and spectral (0.8 km/s) resolution. Our search yields stringent upper limits on the flux of all surveyed molecular lines, which imply abundances relative to CO that are orders of magnitude lower than those observed in protoplanetary disks and Solar System comets, and also those predicted in outgassing models of second-generation material. However, if CI shielding is responsible for extending the lifetime of any CO produced in second-generation collisions, as proposed by Kral et al. (2018), then the line ratios do not reflect true ice phase chemical abundances, but rather imply that CO is shielded by its own photodissociation product, CI, but other molecules are rapidly photodissociated by the stellar and interstellar radiation field.

Asymptotic Giant Branch (AGB) stars are major contributors of cosmic dust to the universe. Typically, dust around AGB stars is investigated via radiative transfer (RT) modeling, or via simple deconstruction of observed spectra. However, methodologies applied vary. Using archival spectroscopic, photometric, and temporal data for the archetypal dusty star, Mira, we identify its circumstellar silicate dust grains. This is achieved by matching the positions and widths of observed spectral features with laboratory data. To do this comparison properly, it is necessary to account for the continuum emission. Here we investigate various ways in which a continuum is eliminated from observational spectra and how it affects the interpretation of spectral features. We find that while the precise continuum shapes and temperatures do not have a critical impact on the positions and shapes of dust spectral features, it is important to eliminate continua in a specific way. It is important to understand what contributes to the spectrum in order to remove the continuum in a way that allows comparison with laboratory spectra of candidate dust species. Our methodologies are applicable to optically thin systems, like that of Mira. Higher optical depths will require RT modeling, which cannot include many different potential astrominerals because there is a lack of complex refractive indices. Finally, we found that the classic silicate feature exhibited by Mira is not consistent with a real amorphous silicate alone but may be best explained with a small alumina contribution to match the observed FWHM of the 10 micron feature.

We present diameters and albedos computed for the near-Earth and Main Belt asteroids observed by the Near-Earth Object Wide-field Infrared Survey Explorer (NEOWISE) spacecraft during the sixth and seventh years of its Reactivation mission. These diameters and albedos are calculated from fitting thermal models to NEOWISE observations of $199$ NEOs and $5851$ MBAs detected during the sixth year of the survey, and $175$ NEOs and $5861$ MBAs from the seventh year. Comparisons of the near-Earth object diameters derived from Reactivation data with those derived from the WISE cryogenic mission data show a $\sim30\%$ relative uncertainty. This larger uncertainty compared to data from the cryogenic mission is due to the need to assume a beaming parameter for the fits to the shorter wavelength data that the Reactivation mission is limited to. We also present an analysis of the orbital parameters of the Main Belt asteroids that have been discovered by NEOWISE during Reactivation, finding that these objects tend to be on orbits that result in their perihelia being far from the ecliptic, and thus missed by other surveys. To date, the NEOWISE Reactivation survey has provided thermal fits of $1415$ unique NEOs. Including the mission phases before spacecraft hibernation increases the count of unique NEOs characterized to $1845$ from WISE's launch to the present.

Although the paradigm of inflation has been extensively studied to demonstrate how macroscopic inhomogeneities in our universe originate from quantum fluctuations, most of the established literature ignores the crucial role that entanglement between the modes of the fluctuating field plays in its observable predictions. In this paper, we import techniques from quantum information theory to reveal hitherto undiscovered predictions for inflation which, in turn, signals how quantum entanglement across cosmological scales can affect large scale structure. Our key insight is that observable long-wavelength modes must be part of an open quantum system, so that the quantum fluctuations can decohere in the presence of an environment of short-wavelength modes. By assuming the simplest model of single-field inflation, and considering the leading order interaction term from the gravitational action, we derive a universal lower bound on the observable effect of such inescapable entanglement.

We identify potential sources of decoherence for $U(1)$ gauge bosons from a cosmological standpoint. Besides interactions with different species in the cosmological medium, we also consider effects due to the expansion of the Universe, which can produce particles (especially scalars) than can potentially interact with the photon in a quantum state. We look in particular at the case of axion-like particles and their predicted decay channels in our analysis. These interactions are shown to have a negligible effect as far as decoherence goes. Interaction rates with CMB radiation or through Thomson scattering are small, so that the interstellar medium remains the biggest decoherence factor. Thus, quantum teleportation experiments with photon energies in the range $1$-$10$ keV should be feasible at cosmological distances up to the galaxy formation epoch or beyond ($z \sim 100$).

We derive gravitational waves in a theory with non-local curvature corrections to the Hilbert-Einstein Lagrangian. In addition to the standard two massless tensor modes, with plus and cross polarizations, helicity 2 and angular frequency $\omega_{1}$, we obtain a further scalar massive mode with helicity 0 and angular frequency $\omega_{2}$, whose polarization is transverse. It is a breathing mode, which, at the lowest order of an effective parameter $\gamma$, presents a speed difference between nearly null and null plane waves. Finally, the quasi-Lorentz $E(2)$-invariant class for the non-local gravity is type $N_{3}$, according to the Petrov classification. This means that the presence (or absence) of gravitational wave modes is observer-independent.

The vacuum transition probabilities for anisotropic universes in the presence of a scalar field potential in the WKB approximation are studied. We follow the work by Cespedes et al [arXiv:2011.13936 [hep-th]], which discuss these transitions in the isotropic context using the Wheeler-DeWitt equation, the Lorentzian Hamiltonian approach and the thin wall limit. First, we propose a general procedure to adapt their formalism to compute the decay rates for any superspace model. Then we apply it to compute the transition probabilities of an FLRW metric with both positive and zero curvature, reproducing in this way one of the results obtained at Cespedes et al. We then proceed to apply the formalism to three anisotropic metrics, namely, Kantowski-Sachs, Bianchi III and biaxial Bianchi IX to compute the rate decays for these three cases. In the process we find that this method involves some conditions which relates the effective number of independent degrees of freedom resulting on all probabilities being described with only two independent variables. For the Bianchi III metric, we find that a general effect of anisotropy is to decrease the transition probability as the degree of anisotropy is increased, having as the isotropic limit the flat FLRW result.

We study the dynamics of a homogeneous, isotropic, and positively curved universe in the presence of a SU(2) gauge field or a triplet of mutually orthogonal vector fields. In the SU(2) case we use the previously known ansatz for the gauge-field configuration, but the case without non-abelian symmetries is more nontrivial and we develop a new ansatz. We in particular consider axion-SU(2) inflation and inflation with vector fields having U(1)$\times$U(1)$\times$U(1) symmetry, and analyze their dynamics in detail numerically. Novel effects of the spatial curvature come into play through vector fields, which causes unconventional pre-inflationary dynamics. It is found that the closed universe with vector fields is slightly more stable against collapse than that filled solely with an inflaton field.

We show how a relativistic Langevin equation can be derived from a Lorentz-covariant version of the Caldeira-Leggett particle-bath Lagrangian. In one of its limits, we identify the obtained equation with the Langevin equation used in contemporary extensions of statistical mechanics to the near-light-speed motion of a Brownian particle in non-relativistic dissipative fluids. The proposed framework provides a more rigorous and first-principles form of the Langevin equation often quoted or postulated as ansatz in previous works. We then refine the aforementioned results by considering more terms in the particle-bath coupling, which improves the precision of the approximation for fully relativistic settings where not only the tagged particle but also the thermal bath motion is relativistic. We discuss the implications of the apparent breaking of space-time translation and parity invariance, showing that these effects are not necessarily in contradiction with the assumptions of statistical mechanics. The intrinsically non-Markovian character of the fully relativistic generalized Langevin equation derived here, and of the associated fluctuation-dissipation theorem, is also discussed.

Future gravitational wave detectors have been projected to be able to probe the nature of compact objects in great detail. In this work, we study the potential observability of the Planck length physics with the tidal deformability of the compact objects in an inspiraling binary. We find that despite the error in the Planck scale distance resolution being exponentially sensitive to errors in the Love number, it is possible to probe them with extreme mass ratio inspirals. We also consider the consequences of the Love number not being sensitive to the Planck scale logarithmically. We discuss how the quantum effects can affect the gravitational wave observables in that scenario.

Axion dark matter can resonantly convert to photons in the magnetosphere of neutron stars, possibly giving rise to radio signals observable on Earth. This method for the indirect detection of axion dark matter has recently received significant attention in the literature. The calculation of the radio signal is complicated by a number of effects; most importantly, the gravitational infall of the axions onto the neutron star accelerates them to semi-relativistic speed, and the neutron star magnetosphere is highly anisotropic. Both of these factors complicate the calculation of the conversion of axions to photons. In this work, we present the first fully three-dimensional calculation of the axion-photon conversion in highly magnetised anisotropic media. Depending on the axion trajectory, this calculation leads to orders-of-magnitude differences in the conversion compared to the simplified one-dimensional calculation used so far in the literature, altering the directionality of the produced photons. Our results will have important implications for the radio signal one would observe in a telescope.