Papers for Monday, May 10 2021

The statistical distributions of active galactic nuclei (AGN), i.e. accreting supermassive black holes (BHs), in mass, space and time, are controlled by a series of key properties, namely the BH-galaxy scaling relations, Eddington ratio distributions and fraction of active BHs (duty cycle). Shedding light on these properties yields strong constraints on the AGN triggering mechanisms whilst providing a clear baseline to create useful mock catalogues for the planning of large galaxy surveys. We here delineate a robust methodology to create mock AGN catalogs built on top of large N-body dark matter simulations via state-of-the-art semi-empirical models. We show that by using as independent tests the AGN clustering at fixed X-ray luminosity, galaxy stellar mass and BH mass, along with the fraction of AGN in groups and clusters, it is possible to significantly narrow down the choice in the relation between black hole mass and host galaxy stellar mass, the duty cycle, and the average Eddington ratio distribution, delivering well-suited constraints to guide cosmological models for the co-evolution of BHs and galaxies. Avoiding such a step-by-step methodology inevitably leads to strong degeneracies in the final mock catalogs, severely limiting their usefulness in understanding AGN evolution and in survey planning and testing.

We present accretion-disk structure measurements from UV-optical reverberation mapping observations of a sample of eight quasars at 0.24<z<0.85. Ultraviolet photometry comes from two cycles of Hubble Space Telescope monitoring, accompanied by multi-band optical monitoring by the Las Cumbres Observatory network and Liverpool Telescopes. The targets were selected from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project sample with reliable black-hole mass measurements from Hbeta reverberation mapping results. We measure significant lags between the UV and various optical griz bands using JAVELIN and CREAM methods. We use the significant lag results from both methods to fit the accretion-disk structure using a Markov chain Monte Carlo approach. We study the accretion disk as a function of disk normalization, temperature scaling, and efficiency. We find direct evidence for diffuse nebular emission from Balmer and FeII lines over discrete wavelength ranges. We also find that our best-fit disk color profile is broadly consistent with the Shakura \& Sunyaev disk model. We compare our UV-optical lags to the disk sizes inferred from optical-optical lags of the same quasars and find that our results are consistent with these quasars being drawn from a limited high-lag subset of the broader population. Our results are therefore broadly consistent with models that suggest longer disk lags in a subset of quasars, for example, due to a nonzero size of the ionizing corona and/or magnetic heating contributing to the disk response.

In this paper, we discuss the impact of the following laboratory experiments and astrophysical observation of neutron stars (NSs) on its equation of state (EoS): (a) The new measurement of neutron skin thickness of $\rm ^{208} \! Pb$, $R_{\rm skin}^{208} = 0.29 \pm 0.07$ fm by the PREX-II experiment. (b) The mass measurement of PSR J0740+6620 has been slightly revised down by including additional $\sim 1.5$ years of pulsar timing data. (c) A possible NICER observation giving the measurement of the radius of PSR J0740+6620 which probably has similar size as PSR J0030+0451. We combine these information using Bayesian statistics along with the previous LIGO/Virgo and NICER observations of NS using a hybrid nuclear+piecewise polytrope EoS parameterization. Our findings are as follows: (a). Adding PREX-II result yields the value of empirical parameter $L = 69^{+16}_{-16}$ MeV, $R_{\rm skin}^{208} = 0.20_{-0.04}^{+0.04}$ fm, and radius of a $1.4 M_{\odot}$ ($R_{1.4}) = 12.66_{- 0.47}^{+ 0.38}$ km at $1 \sigma$ confidence interval (CI). We find these inferred values are mostly dominated by the combined astrophysical observations as the measurement uncertainty in $R_{\rm skin}^{208}$ by PREX-II is much broader. However, a better measurement of $R_{\rm skin}^{208}$ could provide us a tighter constraint on $R_{1.4}$. (b) The revised mass measurement of PSR J0740+6620 has a very marginal effect on the NS EoS. (c) The possible NICER observation may help to estimate the $R_{1.4}$ within $\sim \pm 5\%$ accuracy at $90 \%$ CI which is pretty impressive.

Calculations of the evolution of cosmological perturbations generally involve solution of a large number of coupled differential equations to describe the evolution of the multipole moments of the distribution of photon intensities and polarization. However, this "Boltzmann hierarchy" communicates with the rest of the system of equations for the other perturbation variables only through the photon-intensity quadrupole moment. Here I develop an alternative formulation wherein this photon-intensity quadrupole is obtained via solution of two coupled integral equations -- one for the intensity quadrupole and another for the linear-polarization quadrupole -- rather than the full Boltzmann hierarchy. This alternative method of calculation provides some physical insight and a cross-check for the traditional approach. I describe a simple and efficient iterative numerical solution that converges fairly quickly. I surmise that this may allow current state-of-the-art cosmological-perturbation codes to be accelerated.

The non-resonant Kuiper belt objects (KBOs) between the 3:2 and 2:1 Neptunian mean motion resonances can be largely divided between a cold classical belt (CCB) and a hot classical belt (HCB). A notable difference between these two subpopulations is the prevalence of widely spaced, equal-mass binaries in the CCB and a much smaller but non-zero number in the HCB. The primary reason for this difference in binary rate remains unclear. Here using N-body simulations we examine whether close encounters with the giant planets during an early outer solar system instability may have disrupted primordial Kuiper Belt binaries that existed within the primordial Kuiper belt before they attained HCB orbits. We find that such encounters are very effective at disrupting binaries down to separations of ~1% of their Hill radius (as measured in the modern Kuiper belt), potentially explaining the paucity of widely spaced, equal mass binaries in the modern HCB. Moreover, we find that the widest binaries observed in the modern HCB are quite unlikely to survive planetary encounters, but these same planetary encounters can widen a small subset of tighter binaries to give rise to the small population of very wide binaries seen in today's HCB.

We assemble a sample of 17 low metallicity (7.45 < log(O/H)+12 < 8.12) galaxies with z < 0.1 found spectroscopically, without photometric pre-selection, in early data from the Hobby Eberly Telescope Dark Energy Experiment (HETDEX). Star forming galaxies that occupy the lowest mass and metallicity end of the mass-metallicity relation tend to be under sampled in continuum-based surveys as their spectra are typically dominated by emission from newly forming stars. We search for galaxies with high [OIII]$\lambda$5007 / [OII]$\lambda$3727, implying highly ionized nebular emission often indicative of low metallicity systems. With the Second Generation Low Resolution Spectrograph on the Hobby Eberly Telescope we acquired follow-up spectra, with higher resolution and broader wavelength coverage, of each low-metallicity candidate in order to confirm the redshift, measure the H$\alpha$ and [NII] line strengths and, in many cases, obtain deeper spectra of the blue lines. We find our galaxies are consistent with the mass-metallicity relation of typical low mass galaxies. However, galaxies in our sample tend to have similar specific star formation rates (sSFRs) as the incredibly rare "blueberry" galaxies found in (Yang et. al. 2017). We illustrate the power of spectroscopic surveys for finding low mass and metallicity galaxies and reveal that we find a sample of galaxies that are a hybrid between the properties of typical dwarf galaxies and the more extreme blueberry galaxies.

Temporal variability in flux and spectral shape is ubiquitous in the X-ray sky and carries crucial information about the nature and emission physics of the sources. The EPIC instrument on board the XMM-Newton observatory is the most powerful tool for studying variability even in faint sources. Each day, it collects a large amount of information about hundreds of new serendipitous sources, but the resulting huge (and growing) dataset is largely unexplored in the time domain. The project called Exploring the X-ray transient and variable sky (EXTraS) systematically extracted all temporal domain information in the XMM-Newton archive. This included a search and characterisation of variability, both periodic and aperiodic, in hundreds of thousands of sources spanning more than eight orders of magnitude in timescale and six orders of magnitude in flux, and a search for fast transients that were missed by standard image analysis. All results, products, and software tools have been released to the community in a public archive. A science gateway has also been implemented to allow users to run the EXTraS analysis remotely on recent XMM datasets. We give details on the new algorithms that were designed and implemented to perform all steps of EPIC data analysis, including data preparation, source and background modelling, generation of time series and power spectra, and search for and characterisation of different types of variabilities. We describe our results and products and give information about their basic statistical properties and advice on their usage. We also describe available online resources. The EXTraS database of results and its ancillary products is a rich resource for any kind of investigation in almost all fields of astrophysics. Algorithms and lessons learnt from our project are also a very useful reference for any current and future experiment in the time domain.

The transfer function of the baryon power spectrum from redshift $z\approx 1100$ to today has recently been, for the first time, determined from data by Pardo and Spergel. We observe a remarkable coincidence between this function and the transport function of a cold ideal Fermi gas at different redshifts. Guided by this, we unveil an infinite set of critical temperatures of the relativistic ideal Fermi gas which depend on a very finely quantized long-distance cutoff. The sound horizon scale of Baryon Acoustic Oscillations (BAO) seems set such a cutoff, which dials a critical temperature that is subsequently reached during redshift. At the critical point the Fermi gas becomes scale invariant and may condense to subsequently undergo gravitational collapse, seeding small scale structure. We mention some profound implications including the apparent quantization of Fermi momentum conjugate to the cutoff and the corresponding "gapping" of temperature.

Reports of "cosmology in crisis" are in vogue, but as Mark Twain said, "the report of my death was an exaggeration". We explore what we might actually mean by the standard cosmological model, how tensions - or their apparent resolutions - might arise from too narrow a view, and why looking at the big picture is so essential. This is based on the seminar "All Cosmology, All the Time".

The exotic range of known planetary systems has provoked an equally exotic range of physical explanations for their diverse architectures. However, constraining formation processes requires mapping the observed exoplanet population to that which initially formed in the protoplanetary disc. Numerous results suggest that (internal or external) dynamical perturbation alters the architectures of some exoplanetary systems. Isolating planets that have evolved without any perturbation can help constrain formation processes. We consider the Kepler multiples, which have low mutual inclinations and are unlikely to have been dynamically perturbed. We apply a modelling approach similar to that of Mulders et al. (2018), additionally accounting for the two-dimensionality of the radius ($R =0.3-20\,R_\oplus$) and period ($P= 0.5-730$ days) distribution. We find that an upper limit in planet mass of the form $M_{\rm{lim}} \propto a^\beta \exp(-a_{\rm{in}}/a)$, for semi-major axis $a$ and a broad range of $a_{\rm{in}}$ and $\beta$, can reproduce a distribution of $P$, $R$ that is indistinguishable from the observed distribution by our comparison metric. The index is consistent with $\beta= 1.5$, expected if growth is limited by accretion within the Hill radius. This model is favoured over models assuming a separable PDF in $P$, $R$. The limit, extrapolated to longer periods, is coincident with the orbits of RV-discovered planets ($a>0.2$ au, $M>1\,M_{\rm{J}}$) around recently identified low density host stars, hinting at isolation mass limited growth. We discuss the necessary circumstances for a coincidental age-related bias as the origin of this result, concluding that such a bias is possible but unlikely. We conclude that, in light of the evidence that some planetary systems have been dynamically perturbed, simple models for planet growth during the formation stage are worth revisiting.

Feedback from an active galactic nucleus (AGN) is often implicated as a mechanism that leads to the quenching of galactic star formation. However, AGN-driven quenching is challenging to reconcile with observations that AGN hosts tend to harbour equal (or even excess) amounts of gas compared with inactive galaxies of similar stellar mass. In this paper, we investigate whether AGN feedback happens on sub-galactic (kpc) scales, an effect that might be difficult to detect with global gas measurements. Using kpc-scale measurements of molecular gas (Sigma_H2) and stellar mass (Sigma_*) surface densities taken from the EDGE-CALIFA survey, we show that the gas fractions of central AGN regions are typically a factor of ~2 lower than in star-forming regions. Based on four galaxies with the best spaxel statistics, the difference between AGN and star-forming gas fractions is seen even within a given galaxy, indicating that AGN feedback is able to deplete the molecular gas reservoir in the central few kpc.

We present high resolution millimeter continuum and CO line observations for the circumbinary disk around V892 Tau to constrain the stellar and disk properties. The total mass of the two near-equal-mass A stars is estimated to be $6.0\pm0.2\,M_{\odot}$ based on our models of the Keplerian-dominated gas disk rotation. The detection of strong ionized gas emission associated with the two stars at 8 mm, when combined with previous astrometric measurements in the near-infrared, provides an updated view of the binary orbit with $a=7.1\pm0.1$ au, $e=0.27\pm0.1$, and $P=7.7\pm0.2$ yr, which is about half of a previously reported orbital period. The binary orbital plane is proposed to be near coplanar to the circumbinary disk plane (with a mutual inclination of only $\Delta=8\pm4.2$ deg; another solution with $\Delta=113$ deg is less likely given the short re-alignment timescale). An asymmetric dust disk ring peaking at a radius of 0.''2 is detected at 1.3 mm and its fainter counterparts are also detected at the longer 8 and 9.8 mm. The CO gas disk, though dominated by Keplerian rotation, presents a mild inner and outer disk misalignment, such that the inner disk to the SW and outer disk to the NE appear brighter than their counterparts at the opposite disk sides. The radial extension of the disk, its asymmetric dust ring, and the presence of a disk warp could all be explained by the interaction between the eccentric binary and the circumbinary disk, which we assume were formed with non-zero mutual inclination. Some tentatively detected gas spirals in the outer disk are likely produced by interactions with the low mass tertiary component located 4'' to the northeast. Our analyses demonstrate the promising usage of V892 Tau as an excellent benchmark system to study the details of binary--disk interactions.

We study the properties of slow magneto-acoustic waves that are naturally excited due to turbulent convection and investigate their role in the energy balance of a plage region using three dimensional (3D) radiation-MHD simulations. We calculate the horizontally averaged (over the whole domain) frequency power spectra for both longitudinal and vertical (i.e. the component perpendicular to the surface) components of velocity. To compare our results with the observations we degrade the simulation data with Gaussian kernels having FWHM of 100 km and 200 km, and calculate horizontally averaged power spectra for the vertical component of velocity. The power spectra of the longitudinal component of velocity, averaged over field lines in the core of a kG magnetic flux concentration, reveal that the dominant period of oscillations shifts from around 6.5 minutes in the photosphere to around 4 minutes in the chromosphere. At the same time, the velocity power spectra, averaged horizontally over the whole domain, show that low frequency waves (approximately 6.5 minute period) may reach well into the chromosphere. Importantly, waves with frequencies above 5 mHz propagating along different field lines are found to be out of phase with each other even within a single magnetic concentration. The horizontally averaged power spectra of the vertical component of velocity at various effective resolutions show that the observed acoustic wave energy fluxes are underestimated, by a factor of three even if determined from observations carried out at a high spatial resolution of 200 km. Our results show that longitudinal waves carry (just) sufficient energy to heat the chromosphere in solar plage. We conjecture that current observations (with spatial resolution around 200 km) underestimate the energy flux by roughly a factor of three, or more if the observations have lower spatial resolution.

Supernova remnants (SNRs) can be rich sources of information on the parent SN explosion. Thus investigating the transition from the phase of SN to that of SNR can be crucial to link these two phases of evolution. Here we aim to study the early development of SNR in more details, paying the major attention to the transition from the early-expansion stage to the Sedov stage and the role played by magnetic field in this transition. To this end, spherical magneto-hydrodynamic simulations of SNRs have been performed to study the evolution of magnetic field in young SNRs and explore a sequence of the SNR evolutionary stages in the pre-radiative epoch. Remnants of three supernova types are considered, namely, SNIa, SNIc and SNIIP, that covers a wide space of parameters relevant for SNRs. Changes in global characteristics and development of spatial distributions are analysed. It is shown that the radial component of magnetic field rapidly drops downstream of the forward shock. Therefore, the radially-aligned polarization patterns observed in few young SNRs cannot be reproduced in the one-dimensional MHD simulations. The period SNR takes for the transition from the earliest ejecta-driven phase to the Sedov phase is long enough, with its distinctive physical features, headed by the energy conversion from mostly kinetic one to a fixed ratio between the thermal and kinetic components. This transition worth to be distinguished as a phase in SNR evolutionary scheme. The updated sequence of stages in SNR evolution could be the free expansion (of gas) -- energy-conversion -- Sedov-Taylor -- post-adiabatic -- radiative.

Star-forming dwarf galaxies have properties similar to those expected in high-redshift galaxies. Hence, these local galaxies may provide insights into the evolution of the first galaxies, and the physical processes at work. We present a sample of eleven potential local analogs to high-$z$ (LAHz) galaxies. The sample consists of blue compact dwarf galaxies, selected to have spectral energy distributions that fit galaxies at $1.5<z<4$. We use SOFIA-HAWC+ observations combined with optical and near-infrared data to characterize the dust properties, star formation rate (SFR) and star formation histories (SFH) of the sample of LAHz. We employ Bayesian analysis to characterize the dust using two-component black-body models. Using the LIGHTNING package we fit the spectral energy distribution of the LAHz galaxies over the FUV-FIR wavelength range, and derive the SFH in five time-steps up to a look-back time of 13.3 Gyr. Of the eleven LAHz candidates, six galaxies have SFH consistent with no star formation activity at look-back times beyond 1 Gyr. The remaining galaxies show residual levels of star formation at ages $\gtrsim$1\,Gyr, making them less suitable as local analogs. The six young galaxies stand out in our sample by having the lowest gas-phase metallicities. They are characterized by warmer dust, having the highest specific SFR, and the highest gas mass fractions. The young age of these six galaxies suggests that merging is less important as a driver of the star formation activity. The six LAHz candidates are promising candidates for studies of the gas dynamics role in driving star formation.

We present broadband X-ray spectral-timing analysis of the new Galactic X-ray transient MAXI~J1348--630 using five simultaneous {\it AstroSat} and {\it NICER} observations. Spectral analysis using {\it AstroSat} data identify the source to be in the soft state for the first three observations and in a faint and bright hard state for the next two. Quasi-periodic oscillations at $\sim 0.9$ and $\sim 6.9$\,Hz, belonging to the type-C and type-A class are detected. In the soft state, the power density spectra are substantially lower (by a factor $> 5$) for the {\it NICER} (0.5--12 keV) band compared to the {\it AstroSat}/LAXPC (3--80 keV) one, confirming that the disk is significantly less variable than the Comptonization component. For the first time, energy-dependent fractional rms and time lag in the 0.5--80 keV energy band was measured at different Fourier frequencies, using the bright hard state observation. Hard time lag is detected for the bright hard state, while the faint one shows evidence for soft lag. A single-zone propagation model fits the LAXPC results in the energy band 3--80 keV with parameters similar to those obtained for Cygnus X--1 and MAXI J1820+070. Extending the model to lower energies, reveals qualitative similarities but having quantitative differences with the {\it NICER} results. These discrepancies could be because the {\it NICER} and {\it AstroSat} data are not strictly simultaneous and because the simple propagation model does not take into account disk emission. The results highlight the need for more joint coordinated observations of such systems by {\it NICER} and {\it AstroSat}.

We present results of a multi-line study of the filamentary infrared dark cloud G351.78-0.54 in the 1.3 and 0.8 mm wavelength bands. The lines of the three isotopologues of carbon monoxide CO, N$_2$H$^+$, CH$_3$CCH and HNCO were observed. The aim was to study the general structure of the filamentary cloud, its fragmentation and physical parameters with the emphasis on properties of dense clumps in this cloud. Several dense clumps are identified from the N$_2$H$^+$ (3-2) data, their masses and virial parameters are determined using the C$^{18}$O (2-1) line. Temperatures of some clumps are estimated from the CH$_3$CCH and HNCO data. Almost all clumps appear to be gravitationally unstable. The density estimates obtained from the C$^{18}$O (3-2)/(2-1) and N$_2$H$^+$ (3-2)/(1-0) intensity ratios are in the range $n \sim (0.3-3)\times 10^5$ cm$^{-2}$. The HNCO emission is detected exclusively toward the first clump which contains the luminous IR source IRAS 17233-3606, and indicates an even higher density. It is observed in the outflow, too. The velocity shift of the higher excitation HNCO lines may indicate a movement of the hot core relative the surrounding medium. In some clumps there is a velocity shift $\sim 1$ km s$^{-1}$ between N$_2$H$^+$ (3-2) and CO isotopologues. The large widths of the N$_2$H$^+$ (3-2) line in the clumps indicate an increase of the velocity dispersion in their dense interiors, which may be related to the star formation process. The N$_2$H$^+$ abundance drops toward the luminous IR source.

Based on the second Gaia data release (DR2) and APOGEE (DR16) spectroscopic surveys, wedefined two kinds of star sample: high-velocity thick disk (HVTD) with $v{\phi}>90km/s$ and metal-richstellar halo (MRSH) with $v{\phi}<90km/s$. Due to high resolution spectra data from APOGEE (DR16),we can analyze accurately the element abundance distribution of HVTD and MRSH. These elementsabundance constituted a multidimensional data space, and we introduced an algorithm method forprocessing multi-dimensional data to give the result of dimensionality reduction clustering. Accordingto chemical properties analysis, we derived that some HVTD stars could origin from the thin disk,and some MRSH stars from dwarf galaxies, but those stars which have similar chemical abundancecharacteristics in both sample may form in-situ.

We propose a new method to estimate ion escape from unmagnetized planets that combines observations and models. Assuming that upstream solar wind conditions are known, a computer model of the interaction between the solar wind and the planet is executed for different ionospheric ion production rates. This results in different amounts of mass loading of the solar wind. Then we obtain the ion escape rate from the model run that best fit observations of the bow shock location. As an example of the method we estimate the heavy ion escape from Mars on 2015-03-01 to be $2\cdot 10^{24}$ ions per second, using a hybrid plasma model and observations by MAVEN and Mars Express. This method enables studies of how escape depend on different parameters, and also escape rates during extreme solar wind conditions, applicable to studies of escape in the early solar system, and at exoplanets.

A significant fraction of barred spiral galaxies exhibits peanut/X-shaped structures in their central regions. Bars are known to rotate with a single pattern speed, and they eventually slow down over time due to the dynamical friction with the surrounding dark matter halo. However, the nature of the decay in pattern speed values and whether all peanut bars rotate with a single pattern speed remain to be investigated. Using N-body simulation of a collisionless stellar disc, we study the case of a long bar with a three-dimensional peanut structure prominent in both edge-on and face-on projections. We show that such a bar possesses three distinct peaks in the m=2 Fourier component. Using the Tremaine-Weinberg method, we measure the pattern speeds and demonstrate that the three regions associated with the three peaks rotate with different pattern speeds. The inner region, which is the core of the peanut, rotates slower than the outer regions. In addition, the pattern speed of the inner bar also decays faster than the outer bar with a decay timescale of 4.5 Gyr for the inner part and ~12.5 Gyr for the outer parts. This is manifested as a systematic offset in density and velocity dispersion maps between the inner and outer regions of the long peanut bar. We discuss the importance of our findings in the context of bar dynamics.

The final aim of this paper is to expand the sparse group of X-ray binaries with gamma-ray counterparts as laboratories to study high-energy processes under physical conditions that periodically repeat. A follow-up of a candidate system has been carried out. We applied both photometric and spectroscopic techniques in the optical domain together with a period analysis using the phase dispersion minimization and CLEAN methods. A tentative period search was also conducted in the gamma-ray domain. Our main result is having established the binary nature of the optical star and X-ray emitter HD 3191 towards the Fermi gamma-ray source 4FGL J0035.8+6131, the last one proposed to be associated with a blazar of an unknown type. An orbital period close to 16 days is reported for HD 3191 together with a likely rotation, or pulsation, period of about 0.6 d. Although no convincing evidence for the orbital cycle has been found in the Fermi light curve up to now, the confirmed presence of a high-mass X-ray binary towards 4FGL J0035.8+6131 now strengthens the need for caution about its true nature.

Using recent observational data, we construct a set of multi-component equilibrium models of the disk of a Milky Way-like galaxy. The disk dynamics are studied using collisionless-gaseous numerical simulations, based on the joined integration of the equations of motion for the collision-less particles using direct integration of gravitational interaction and the gaseous SPH-particles. We find that after approximately one Gyr, a prominent central bar is formed having a semi-axis length of about three kpc, together with a multi-armed spiral pattern represented by a superposition of $m=$ 2-, 3-, and 4-armed spirals. The spiral structure and the bar exist for at least 3 Gyr in our simulations. The existence of the Milky Way bar imposes limitations on the density distributions in the subsystems of the Milky Way galaxy. We find that a bar does not form if the radial scale length of the density distribution in the disk exceeds 2.6 kpc. As expected, the bar formation is also suppressed by a compact massive stellar bulge. We also demonstrate that the maximum value in the rotation curve of the disk of the Milky Way galaxy, as found in its central regions, is explained by non-circular motion due to the presence of a bar and its orientation relative to an observer.

Coronal loops are building blocks of solar active regions. However, their formation mechanism is still not well understood. Here we present direct observational evidence for the formation of coronal loops through magnetic reconnection as new magnetic fluxes emerge into the solar atmosphere. Extreme-ultraviolet observations of the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) clearly show the newly formed loops following magnetic reconnection within a plasma sheet. Formation of the loops is also seen in the h{\alpha} line-core images taken by the New Vacuum Solar Telescope. Observations from the Helioseismic and Magnetic Imager onboard SDO show that a positive-polarity flux concentration moves towards a negative-polarity one with a speed of ~0.4 km/s, before the formation of coronal loops. During the loop formation process, we found signatures of flux cancellation and subsequent enhancement of the transverse field between the two polarities. The three-dimensional magnetic field structure reconstructed through a magnetohydrostatic model shows field lines consistent with the loops in AIA images. Numerous bright blobs with an average width of 1.37 Mm appear intermittently in the plasma sheet and move upward with a projected velocity of ~114 km/s. The temperature, emission measure and density of these blobs are about 3 MK, 2.0x10^(28) cm^(-5) and 1.2x10^(10) cm^(-3), respectively. A power spectral analysis of these blobs indicates that the observed reconnection is likely not dominated by a turbulent process. We have also identified flows with a velocity of 20 to 50 km/s towards the footpoints of the newly formed coronal loops.

The Central Molecular Zone (CMZ), a $\sim$200 pc sized region around the Galactic Centre, is peculiar in that it shows a star formation rate (SFR) that is suppressed with respect to the available dense gas. To study the SFR in the CMZ, young stellar objects (YSOs) can be investigated. Here we present radio observations of 334 2.2 $\mu$m infrared sources that have been identified as YSO candidates. Our goal is to investigate the presence of centimetre wavelength radio continuum counterparts to this sample of YSO candidates which we use to constrain the current SFR in the CMZ. As part of the GLOSTAR survey, D-configuration VLA data was obtained for the Galactic Centre, covering -2$^{\circ}<l<$2$^{\circ}$ and -1$^{\circ}<b<$1$^{\circ}$, with a frequency coverage of 4-8 GHz. We matched YSOs with radio continuum sources based on selection criteria and classified these radio sources as potential HII regions and determined their physical properties. Of the 334 YSO candidates, we found 35 with radio continuum counterparts. We find that 94 YSOs are associated with dense dust condensations identified in the 870 $\mu$m ATLASGAL survey, of which 14 have a GLOSTAR counterpart. Of the 35 YSOs with radio counterparts, 11 are confirmed as HII regions, based on their spectral indices and the literature. We estimated their Lyman continuum photon flux in order to estimate the mass of the ionising star. Combining these with known sources, the present-day SFR in the CMZ is calculated to be $\sim$0.068 M$_{\odot}$ yr$^{-1}$, which is $\sim$6.8$\%$ of the Galactic SFR. Candidate YSOs that lack radio counterparts may not have yet evolved to the stage of exhibiting an HII region or, conversely, are older and have dispersed their natal clouds. Since many lack dust emission, the latter is more likely. Our SFR estimate in the CMZ is in agreement with previous estimates in the literature.

We use a combination of Hipparcos space mission data with the USNO dedicated ground-based astrometric program URAT-Bright designed to complement and verify Gaia results for the brightest stars in the south to estimate the small perturbations of observed proper motions caused by exoplanets. One of the 1423 bright stars in the program, $\delta$ Pav, stands out with a small proper motion difference between our long-term estimate and Gaia EDR3 value, which corresponds to a projected velocity of $(-17,+13)$ m s$^{-1}$. This difference is significant at a 0.994 confidence in the RA component, owing to the proximity of the star and the impressive precision of proper motions. The effect is confirmed by a comparison of long-term EDR3-Hipparcos and short-term Gaia EDR3 proper motions at a smaller velocity, but with formally absolute confidence. We surmise that the close Solar analog $\delta$ Pav harbors a long-period exoplanet similar to Jupiter.

Future experiments based on the observation of Earth's atmosphere from sub-orbital and orbital altitudes plan to include optical Cherenkov cameras to observe extensive air showers produced by high-energy cosmic radiation via its interaction with both the Earth and its atmosphere. As discussed elsewhere, particularly relevant is the case of upward-moving showers initiated by astrophysical neutrinos skimming and interacting in the Earth. The Cherenkov cameras, by looking above Earth's limb, can also detect cosmic rays with energies starting from less than a PeV up to the highest energies (tens of EeV). Using a customized computation scheme to determine the expected optical Cherenkov signal from these high-energy cosmic rays, we estimate the sensitivity and event rate for balloon-borne and satellite-based instruments, focusing our analysis on the Extreme Universe Space Observatory aboard a Super Pressure Balloon 2 (EUSO-SPB2) and the Probe of Extreme Multi-Messenger Astrophysics (POEMMA) experiments. We find the expected event rates to be larger than hundreds of events per hour of experimental live time, enabling a promising overall test of the Cherenkov detection technique from sub-orbital and orbital altitudes as well as a guaranteed signal that can be used for understanding the response of the instrument.

We present $^{56}$Ni mass estimates for 110 normal Type II supernovae (SNe II), computed here from their luminosity in the radioactive tail. This sample consists of SNe from the literature, with at least three photometric measurements in a single optical band within 95-320 d since explosion. To convert apparent magnitudes to bolometric ones, we compute bolometric corrections (BCs) using 15 SNe in our sample having optical and near-IR photometry, along with three sets of SN II atmosphere models to account for the unobserved flux. We find that the $I$- and $i$-band are best suited to estimate luminosities through the BC technique. The $^{56}$Ni mass distribution of our SN sample has a minimum and maximum of 0.005 and 0.177 M$_{\odot}$, respectively, and a selection-bias-corrected average of $0.037\pm0.005$ M$_{\odot}$. Using the latter value together with iron isotope ratios of two sets of core-collapse (CC) nucleosynthesis models, we calculate a mean iron yield of $0.040\pm0.005$ M$_{\odot}$ for normal SNe II. Combining this result with recent mean $^{56}$Ni mass measurements for other CC SN subtypes, we estimate a mean iron yield $<$0.068 M$_{\odot}$ for CC SNe, where the contribution of normal SNe II is $>$36 per cent. We also find that the empirical relation between $^{56}$Ni mass and steepness parameter ($S$) is poorly suited to measure the $^{56}$Ni mass of normal SNe II. Instead, we present a correlation between $^{56}$Ni mass, $S$, and absolute magnitude at 50 d since explosion. The latter allows to measure $^{56}$Ni masses of normal SNe II with a precision around 30 per cent.

From an on-going survey of the Galactic bulge, we have discovered a number of compact, steep spectrum radio sources. In this present study we have carried out more detailed observations for two of these sources, located 43 arcmin and 12.7 deg from the Galactic Center. Both sources have a very steep spectrum (alpha ~ -3) and are compact, with upper limits on the angular size of 1-2 arcsec. Their flux densities appear to be relatively steady on timescales of years, months, and hours, with no indications of rapid variability or transient behavior. We detect significant circularly polarized emission from both sources, but only weak or upper limits on linear polarization. Neither source has a counterpart at other wavelengths and deep, high-frequency searches fail to find pulsations. We compare their source properties with other known compact, non-thermal source populations in the bulge (e.g. X-ray binaries, magnetars, the Burper, cataclysmic variables). Our existing data support the hypothesis that they are scatter broadened millisecond or recycled pulsars, either at the bulge or along the line of sight. We also consider the possibility that they may be a new population of Galactic radio sources which share similar properties as pulsars but lack pulsations; a hypothesis that can be tested by future large-scale synoptic surveys.

Data assimilation is an increasingly popular technique in Mars atmospheric science, but its effect on the mean states of the underlying atmosphere models has not been thoroughly examined. The robustness of results to the choice of model and assimilation algorithm also warrants further study. We investigate these issues using two Mars general circulation models (MGCMs), with particular emphasis on zonal wind and temperature fields. When temperature retrievals from the Mars Global Surveyor Thermal Emission Spectrometer (TES) are assimilated into the U.K.-Laboratoire de M\'et\'eorologie Dynamique (UK-LMD) MGCM to create the Mars Analysis Correction Data Assimilation (MACDA) reanalysis, low-level zonal jets in the winter northern hemisphere shift equatorward and weaken relative to a free-running control simulation from the same MGCM. The Ensemble Mars Atmosphere Reanalysis System (EMARS) reanalysis, which is also based on TES temperature retrievals, also shows jet weakening (but less if any shifting) relative to a control simulation performed with the underlying Geophysical Fluid Dynamics Laboratory (GFDL) MGCM. Examining higher levels of the atmosphere, monthly mean three-dimensional temperature and zonal wind fields are in generally better agreement between the two reanalyses than between the two control simulations. In conjunction with information about the MGCMs' physical parametrizations, intercomparisons between the various reanalyses and control simulations suggest that overall the EMARS control run is plausibly less biased (relative to the true state of the Martian atmosphere) than the MACDA control run. Implications for future observational studies are discussed.

The proposition that life can spread from one planetary system to another (interstellar panspermia) has a long history, but this hypothesis is difficult to test through observations. We develop a mathematical model that takes parameters such as the microbial survival lifetime, the stellar velocity dispersion, and the dispersion of ejecta into account in order to assess the prospects for detecting interstellar panspermia. We show that the correlations between pairs of life-bearing planetary systems (embodied in the pair-distribution function from statistics) may serve as an effective diagnostic of interstellar panspermia, provided that the velocity dispersion of ejecta is greater than the stellar dispersion. We provide heuristic estimates of the model parameters for various astrophysical environments, and conclude that open clusters and globular clusters appear to represent the best targets for assessing the viability of interstellar panspermia.

Galaxy mergers occur frequently in the early universe and bring multiple supermassive black holes (SMBHs) into the nucleus, where they may eventually coalesce. Identifying post-merger-scale (i.e., <~a few kpc) dual SMBHs is a critical pathway to understanding their dynamical evolution and successive mergers. While serendipitously discovering kpc-scale dual SMBHs at z<1 is possible, such systems are elusive at z>2, but critical to constraining the progenitors of SMBH mergers. The redshift z~2 also marks the epoch of peak activity of luminous quasars, hence probing this spatial regime at high redshift is of particular significance in understanding the evolution of quasars. However, given stringent resolution requirements, there is currently no confirmed <10 kpc physical SMBH pair at z>2. Here we report two sub-arcsec double quasars at z>2 discovered from a targeted search with a novel astrometric technique, demonstrating a high success rate (~50%) in this systematic approach. These high-redshift double quasars could be the long-sought kpc-scale dual SMBHs, or sub-arcsec gravitationally-lensed quasar images. One of these double quasars (at z=2.95) was spatially resolved with optical spectroscopy, and slightly favors the scenario of a physical quasar pair with a projected separation of 3.5 kpc (0.46"). Follow-up observations of double quasars discovered by this targeted approach will be able to provide the first observational constraints on kpc-scale dual SMBHs at z>2.

The robustness of multifield inflation to the physics of reheating is investigated. In order to carry out this study, reheating is described in detail by means of a formalism which tracks the evolution of scalar fields and perfect fluids in interaction (the inflatons and their decay products). This framework is then used to establish the general equations of motion of the background and perturbative quantities controlling the evolution of the system during reheating. Next, these equations are solved exactly by means of a new numerical code and new analytical techniques, allowing us to interpret and approximate these solutions, are developed. As an illustration of a physical prediction that could be affected by the micro-physics of reheating, the amplitude of non-adiabatic perturbations in double inflation is considered. It is found that ignoring the fine-structure of reheating, as usually done in the standard approach, can lead to differences as big as $\sim 50\%$, while our semi-analytic estimates can reduce this error to $\sim 10\%$. We conclude that, in multifield inflation, tracking the perturbations through the details of the reheating process is important and, to achieve good precision, requires the use of numerical calculations.

We discuss how to efficiently and reliably estimate the level of agreement and disagreement on parameter determinations from different experiments, fully taking into account non-Gaussianities in the parameter posteriors. We develop two families of scalable algorithms that allow us to perform this type of calculations in increasing number of dimensions and for different levels of tensions. One family of algorithms rely on kernel density estimates of posterior distributions while the other relies on machine learning modeling of the posterior distribution with normalizing flows. We showcase their effectiveness and accuracy with a set of benchmark examples and find both methods agree with each other and the true tension within $0.5\sigma$ or better. This allows us to study the level of internal agreement between different measurements of the clustering of cosmological structures from the Dark Energy Survey and their agreement with measurements of the Cosmic Microwave Background from the Planck satellite.

Tomographic three-dimensional 21 cm images from the epoch of reionization contain a wealth of information about the reionization of the intergalactic medium by astrophysical sources. Conventional power spectrum analysis cannot exploit the full information in the 21 cm data because the 21 cm signal is highly non-Gaussian due to reionization patchiness. We perform a Bayesian inference of the reionization parameters where the likelihood is implicitly defined through forward simulations using density estimation likelihood-free inference (DELFI). We adopt a trained 3D Convolutional Neural Network (CNN) to compress the 3D image data into informative summaries (DELFI-3D CNN). We show that this method recovers accurate posterior distributions for the reionization parameters. Our approach outperforms earlier analysis based on two-dimensional 21 cm images. In contrast, an MCMC analysis of the 3D lightcone-based 21 cm power spectrum alone and using a standard explicit likelihood approximation results in inaccurate credible parameter regions both in terms of the location and shape of the contours. Our proof-of-concept study implies that the DELFI-3D CNN can effectively exploit more information in the 3D 21 cm images than a 2D CNN or power spectrum analysis. This technique can be readily extended to include realistic effects and is therefore a promising approach for the scientific interpretation of future 21 cm observation data.

The chemical pathways linking the small organic molecules commonly observed in molecular clouds to the large, complex, polycyclic species long-suspected to be carriers of the ubiquitous unidentified infrared emission bands remain unclear. To investigate whether the formation of mono- and poly-cyclic molecules observed in cold cores could form via the bottom-up reaction of ubiquitous carbon-chain species with, e.g. atomic hydrogen, a search is made for possible intermediates in data taken as part of the GOTHAM (GBT Observations of TMC-1 Hunting for Aromatic Molecules) project. Markov-Chain Monte Carlo (MCMC) Source Models were run to obtain column densities and excitation temperatures. Astrochemical models were run to examine possible formation routes, including a novel grain-surface pathway involving the hydrogenation of C$_6$N and HC$_6$N, as well as purely gas-phase reactions between C$_3$N and both propyne (CH$_3$CCH) and allene (CH$_2$CCH$_2$), as well as via the reaction CN + H$_2$CCCHCCH. We report the first detection of cyanoacetyleneallene (H$_2$CCCHC$_3$N) in space toward the TMC-1 cold cloud using the Robert C. Byrd 100 m Green Bank Telescope (GBT). Cyanoacetyleneallene may represent an intermediate between less-saturated carbon-chains, such as the cyanopolyynes, that are characteristic of cold cores and the more recently-discovered cyclic species like cyanocyclopentadiene. Results from our models show that the gas-phase allene-based formation route in particular produces abundances of H$_2$CCCHC$_3$N that match the column density of $2\times10^{11}$ cm$^{-2}$ obtained from the MCMC Source Model, and that the grain-surface route yields large abundances on ices that could potentially be important as precursors for cyclic molecules.

Radiative feedback can influence subsequent star formation. We quantify the heating from OB stars in the local star-forming regions in the JCMT Gould Belt survey. Dust temperatures are calculated from 450/850 micron flux ratios from SCUBA-2 observations at the JCMT assuming a fixed dust opacity spectral index $\beta=1.8$. Mean dust temperatures are calculated for each submillimetre clump along with projected distances from the main OB star in the region. Temperature vs. distance is fit with a simple model of dust heating by the OB star radiation plus the interstellar radiation field and dust cooling through optically thin radiation. Classifying the heating sources by spectral type, O-type stars produce the greatest clump average temperature rises and largest heating extent, with temperatures over 40 K and significant heating out to at least 2.4 pc. Early-type B stars (B4 and above) produce temperatures of over 20 K and significant heating over 0.4 pc. Late-type B stars show a marginal heating effect within 0.2 pc. For a given projected distance, there is a significant scatter in clump temperatures that is due to local heating by other luminous stars in the region, projection effects, or shadowing effects. Even in these local, `low-mass' star-forming regions, radiative feedback is having an effect on parsec scales, with 24% of the clumps heated to at least 3 K above the 15 K base temperature expected from heating by only the interstellar radiation field, and a mean dust temperature for heated clumps of 24 K.

We use the Fundamental Plane of Elliptical Galaxies to constrain the so-called Hybrid Gravity, a modified theory of gravity where General Relativity is improved by further degrees of freedom of metric-affine Palatini formalism of $f(\cal R)$ gravity. Because the Fundamental Plane is connected to the global properties of elliptical galaxies, it is possible to obtain observational constraints on the parameters of Hybrid Gravity in the weak field limit. We analyze also the velocity distribution of elliptical galaxies comparing our theoretical results in the case of Hybrid Gravity with astronomical data for elliptical galaxies. In this way, we are able to constrain the Hybrid Gravity parameters $m_\phi$ and $\phi_0$. We show that the Fundamental Plane, i.e. $v_c/\sigma$ relations, can be used as a standard tool to probe different theories of gravity in the weak field limit. We conclude that Hybrid Gravity is able to explain elliptical galaxies with different stellar kinematics without the dark matter hypothesis.

We provide a new analysis technique to measure the effect of the isotropic polarization rotation, induced by e.g. the isotropic cosmic birefringence from axion-like particles and a miscalibration of CMB polarization angle, via mode coupling in the cosmic microwave background (CMB). Several secondary effects such as gravitational lensing and CMB optical-depth anisotropies lead to mode coupling in observed CMB anisotropies, i.e., non-zero off-diagonal elements in the observed CMB covariance. To derive the mode coupling, however, we usually assume no parity violation in the observed CMB anisotropies. We first derive a new contribution to the CMB mode coupling arising from parity violation in observed CMB. Since the isotropic polarization rotation leads to parity violation in the observed CMB anisotropies, we then discuss the use of the new mode coupling for constraining the isotropic polarization angle. We find that constraints on the isotropic polarization angle by measuring the new mode-coupling contribution are comparable to that using the $EB$ cross-power spectrum in future high-sensitivity polarization experiments such as CMB-S4 and LiteBIRD. Thus, this technique can be used to cross-check results obtained by the use of the $EB$ cross-power spectrum.

Wu & Peek (2020) predict SDSS-quality spectra based on Pan-STARRS broad-band \textit{grizy} images using machine learning (ML). In this letter, we test their prediction for a unique object, UGC 2885 ("Rubin's galaxy"), the largest and most massive, isolated disk galaxy in the local Universe ($D<100$ Mpc). After obtaining the ML predicted spectrum, we compare it to all existing spectroscopic information that is comparable to an SDSS spectrum of the central region: two archival spectra, one extracted from the VIRUS-P observations of this galaxy, and a new, targeted MMT/Binospec observation. Agreement is qualitatively good, though the ML prediction prefers line ratios slightly more towards those of an active galactic nucleus (AGN), compared to archival and VIRUS-P observed values. The MMT/Binospec nuclear spectrum unequivocally shows strong emission lines except H$\beta$, the ratios of which are consistent with AGN activity. The ML approach to galaxy spectra may be a viable way to identify AGN supplementing NIR colors. How such a massive disk galaxy ($M^* = 10^{11}$ M$_\odot$), which uncharacteristically shows no sign of interaction or mergers, manages to fuel its central AGN remains to be investigated.

In this work, we study the effects of temperature on magnetic white dwarfs. We model their interior as a nuclei lattice surrounded by a relativistic free Fermi gas of electrons, accounting for effects from temperature, Landau levels and anomalous magnetic moment. We find that, at low densities (corresponding to the outer regions of star), both temperature and magnetic field effects play an important role in the calculation of microscopic thermodynamical quantities. To study macroscopic stellar structures within a general-relativistic approach, we solve numerically the coupled Einstein's-Maxwell's equations for fixed entropy per particle configurations and discuss how temperature affects stellar magnetic field profiles, masses and radii.

In this study we investigate possible applications of observed S2 orbit around Galactic Center for constraining the Yukawa gravity at scales in the range between several tens and several thousands astronomical units (AU) to obtain graviton mass constraints. In our model we suppose that bulk distribution of matter (includes stellar cluster, interstellar gas distribution and dark matter) exists near Supermassive Black Hole (SMBH) in our Galactic Center. We obtain the values of orbital precession angle for different values of mass density of matter and we require that the value of orbital precession is the same like in General Relativity (GR). From that request we determine gravity parameter $\lambda$ and the upper value for graviton mass. We found that in the cases where the density of extended mass is higher, the maximum allowed value for parameter $\lambda$ is smaller and the upper limit for graviton mass is higher. It is due to the fact that the extended mass causes the retrograde orbital precession. We believe that this study is a very efficient tool to evaluate a gravitational potential at the Galactic Center, parameter $\lambda$ of the Yukawa gravity model, and to constrain the graviton mass.

A number of recent, low-redshift, lensing measurements hint at a universe in which the amplitude of lensing is lower than that predicted from the $\Lambda$CDM model fit to the data of the Planck CMB mission. Here we use the auto- and cross-correlation signal of unWISE galaxies and Planck CMB lensing maps to infer cosmological parameters at low redshift. In particular, we consider three unWISE samples (denoted as "blue," "green" and "red") at median redshifts z~0.6, 1.1, and 1.5, which fully cover the Dark Energy dominated era. Our cross-correlation measurements, with combined significance S/N ~ 80, are used to infer the amplitude of low-redshift fluctuations, $\sigma_8$, the fraction of matter in the Universe, $\Omega_m$,and the combination $S_8 = \sigma_8 (\Omega_m/0.3)^{0.5}$ to which these low-redshift lensing measurements are most sensitive. The combination of blue and green samples gives a value $S_8 =0.776\pm0.017$, that is fully consistent with other low-redshift lensing measurements and in $2.6\sigma$ tension with the CMB predictions from Planck. This is noteworthy, because CMB lensing probes the same physics as previous galaxy lensing measurements, but with very different systematics, thus providing an excellent complement to previous measurements.

The Event Horizon Telescope recently captured images of the supermassive black hole in the center of the M87 galaxy, which show a ring-like emission structure with the South side only slightly brighter than the North side. This relatively weak asymmetry in the brightness profile along the ring has been interpreted as a consequence of the low inclination of the observer (around 17 deg for M87), which suppresses the Doppler beaming and boosting effects that might otherwise be expected due to the nearly relativistic velocities of the orbiting plasma. In this work, we use a large suite of general relativistic magnetohydrodynamic simulations to reassess the validity of this argument. By constructing explicit counter examples, we show that low-inclination is a sufficient but not necessary condition for images to have low brightness asymmetry. Accretion flow models with high accumulated magnetic flux close to the black hole horizon (the so-called magnetically arrested disks) and low black-hole spins have angular velocities that are substantially smaller than the orbital velocities of test particles at the same location. As a result, such models can produce images with low brightness asymmetry even when viewed edge on.

Recent theoretical studies have suggested that a magnetic field may play a crucial role in the first star formation in the universe. However, the influence of the magnetic field on the first star formation has yet to be understood well. In this study, we perform three-dimensional magnetohydrodynamic simulations taking into account all the relevant cooling processes and non-equilibrium chemical reactions up to the protostar density, in order to study the collapse of magnetized primordial gas cores with self-consistent thermal evolution. Our results show that the thermal evolution of the central core is hardly affected by a magnetic field, because magnetic forces do not prevent the contraction along the fields lines. We also find that the magnetic braking extracts the angular momentum from the core and suppresses fragmentation depending on the initial strength of the magnetic field. The angular momentum transport by the magnetic outflows is less effective than that by the magnetic braking because the outflows are launched only in a late phase of the collapse. Our results indicate that the magnetic effects become important for the field strength $B> 10^{-8}(n_{\rm H}/1\ \rm cm^{-3})^{2/3}\ \rm G$, where $n_{\rm H}$ is the number density, during the collapse phase. Finally, we compare our results with simulations using a barotropic approximation and confirm that this approximation is reasonable at least for the collapse phase. Nevertheless, self-consistent treatment of the thermal and chemical processes is essential for extending simulations to the accretion phase, in which radiative feedback by protostars plays a crucial role.

Many non-linear scalar field theories possess a screening mechanism that can suppress any associated fifth force in dense environments. As a result, these theories can evade local experimental tests of new forces. Chameleon-like screening, which occurs because of non-linearities in the scalar potential or the coupling to matter, is well understood around extended objects. However, many experimental tests of these theories involve objects with spatial extent much smaller than the scalar field's Compton wavelength, and which could therefore be considered point-like. In this work, we determine how the fifth forces are screened in the limit that the source objects become extremely compact.

The LIGO and Virgo observatories have reported 39 new gravitational-wave detections during the first part of the third observation run, bringing the total to 50. Most of these new detections are consistent with binary black-hole coalescences, making them suitable targets to search for gravitational-wave memory, a non-linear effect of general relativity. We extend a method developed in previous publications to analyse these events to determine a Bayes factor comparing the memory hypothesis to the no-memory hypothesis. Specifically, we calculate Bayes factors using two waveform models with higher-order modes that allow us to analyse events with extreme mass ratios and precessing spins, both of which have not been possible before. Depending on the waveform model we find a combined $\ln \mathrm{BF}_{\mathrm{mem}} = 0.024$ or $\ln \mathrm{BF}_{\mathrm{mem}} = 0.049$ in favour of memory. This result is consistent with recent predictions that indicate $\mathcal{O}(2000)$ binary black-hole detections will be required to confidently establish the presence or absence of memory.

A single space-based gravitational wave detector will push the boundaries of astronomy and fundamental physics. Having a network of two or more detectors would significantly improve source localization. Here we consider how dual networks of space-based detectors would improve parameter estimation of massive black hole binaries. We consider two scenarios: a network comprised of the Laser Interferometer Space Antenna (LISA) and an additional LISA-like heliocentric detector (e.g. Taiji); and a network comprised of LISA with an an additional geocentric detector (e.g. TianQin). We use Markov chain Monte Carlo techniques and Fisher matrix estimates to explore the impact of a two detector network on sky localization and distance determination. The impact on other source parameters is also studied. With the addition of a Taiji or TianQin, we find orders of magnitude improvements in sky localization for the more massive MBHBs, while also seeing improvements for lower mass systems, and for other source parameters.

The task of morphological classification is complex for simple parameterization, but important for research in the galaxy evolution field. Future galaxy surveys (e.g. EUCLID) will collect data about more than a $10^9$ galaxies. To obtain morphological information one needs to involve people to mark up galaxy images, which requires either a considerable amount of money or a huge number of volunteers. We propose an effective semi-supervised approach for galaxy morphology classification task, based on active learning of adversarial autoencoder (AAE) model. For a binary classification problem (top level question of Galaxy Zoo 2 decision tree) we achieved accuracy 93.1% on the test part with only 0.86 millions markup actions, this model can easily scale up on any number of images. Our best model with additional markup achieves accuracy of 95.5%. To the best of our knowledge it is a first time AAE semi-supervised learning model used in astronomy.

Amongst all the renewable energy resources (RES), solar is the most popular form of energy source and is of particular interest for its widely integration into the power grid. However, due to the intermittent nature of solar source, it is of the greatest significance to forecast solar irradiance to ensure uninterrupted and reliable power supply to serve the energy demand. There are several approaches to perform solar irradiance forecasting, for instance satellite-based methods, sky image-based methods, machine learning-based methods, and numerical weather prediction-based methods. In this paper, we present a review on short-term intra-hour solar prediction techniques known as nowcasting methods using sky images. Along with this, we also report and discuss which sky image features are significant for the nowcasting methods.

Accurate extractions of gravitational wave signal waveforms are essential to validate a detection and to probe the astrophysics behind the sources producing the gravitational waves. This however, could be difficult in realistic scenarios where the signals detected by the gravitational wave detectors could be contanimnated with non-stationary and non-Gaussian noise. In this paper, we demonstrate for the first time that a deep learning architecture, consisting of Convolutional Neural Network and bi-directional Long Short-Term Memory components can be used to extract all ten detected binary black hole gravitational wave waveforms from the detector data of LIGO-Virgo's first and second science runs with a high accuracy of 0.97 overlap compared to published waveforms.

The production of heavy-mass elements due to the rapid neutron-capture mechanism (r-process) is associated with astrophysical scenarios, such as supernovae and neutron-star mergers. In the r-process the capture of neutrons is followed by $\beta$-decays until nuclear stability is reached. A key element in the chain of nuclear weak-decays leading to the production of isotopes may be the change of the parameters controlling the neutrino sector, due to the mixing of active and sterile species. In this work we have addressed this question and calculated $\beta$-decay rates for the nuclei involved in the r-process chains as a function of the neutrino mixing parameters. These rates were then used in the calculation of the abundance of the heavy elements produced in core-collapse supernova and in neutron-star mergers, starting from different initial mass-fraction distributions. The analysis shows that the core-collapse supernova environment contributes with approximately $30\%$ of the total heavy nuclei abundance while the neutron-star merger contributes with about $70\%$ of it. Using available experimental data we have performed a statistical analysis to set limits on the active-sterile neutrino mixing angle and found a best-fit value $\sin^2 2\theta_{14}=0.22$, a value comparable with those found in other studies reported in the literature.

The dynamics of accretion disk is considered taking into account the Lense-Thirring precession in the presence of cosmological constant $\Lambda$ of Schwarzschild-de Sitter. The nodal and apsidal frequencies are obtained, and the role of $\Lambda$ is revealed in their properties, including the consequences for the Bardeen-Petterson effect.

Observation of high energy cosmic neutrinos by ICECUBE has ushered in a new era in exploring both cosmos and new physics beyond the Standard Model (SM). In the standard picture, although mostly $\nu_\mu$ and $\nu_e$ are produced in the source, oscillation will produce $\nu_\tau$ {\it en route}. Certain beyond SM scenarios, like interaction with ultralight DM can alter this picture. Thus, the flavor composition of the cosmic neutrino flux can open up the possibility of exploring certain beyond the SM scenarios that are inaccessible otherwise. We show that the $\tau$ flavor holds a special place among the neutrino flavors in elucidating new physics. Interpreting the two anomalous events observed by ANITA as $\nu_\tau$ events makes the tau flavor even more intriguing. We study how the detection of the two tau events by ICECUBE constrains the interaction of the neutrinos with ultralight dark matter and discuss the implications of this interaction for even higher energy cosmic neutrinos detectable by future radio telescopes such as ARA, ARIANNA and GRAND. We also revisit the $3+1$ neutrino scheme as a solution to the two anomalous ANITA events and clarify a misconception that exists in the literature about the evolution of high energy neutrinos in matter within the $3+1$ scheme with a possibility of scattering off nuclei. We show that the existing bounds on the flux of $\nu_\tau$ with energy of EeV rules out this solution for the ANITA events. We show that the $3+1$ solution can be saved from both this bound and from the bound on the extra relativistic degrees of freedom in the early universe by turning on the interaction of neutrinos with ultralight dark matter.

Baryonic matter close to the saturation density is very likely to present complex inhomogeneous structures collectively known under the name of pasta phase. At finite temperature, the different geometric structures are expected to coexist, with potential consequences on the neutron star crust conductivity and neutrino transport in supernova matter. In the framework of a statistical multi-component approach, we calculate the composition of matter in the pasta phase considering density, proton fraction, and geometry fluctuations. Using a realistic energy functional from relativistic mean field theory and a temperature and isospin dependent surface tension fitted from Thomas-Fermi calculations, we show that different geometries can coexist in a large fraction of the pasta phase, down to temperatures of the order of the crystallization temperature of the neutron star crust. Quantitative estimates of the charge fluctuations are given.

This is the translation from Latin of E547 'Determinatio facilis orbitae cometae, cuius transitum per eclipticam bis observare licuit', in which Euler addresses the determination of a comet's parabolic orbit, with the Sun at the focus, from two astronomical observations from the earth, when the comet crosses the ecliptic at the ascending and descending nodes. The key point of the calculation is the solution of a fourth degree polynomial, from which the determination of the orbital parameters are determined from one of its roots.

This is a translation from Latin of E840 'De motu cometarum in orbitis parabolicis, solem in foco habentibus', in which Euler addresses six problems related to comets in heliocentric parabolic orbits. Problem 1: Find the true anomaly of a heliocentric comet from the latus rectum of the orbit and the medium Earth to Sun distance. Problem 2: Find the orbit of a heliocentric comet from three given positions. Problem 3: Knowing the orbit of a comet, and the instant in time in which it dwells in the perihelion, define its longitude and latitude at any time. Problem 4: From two locations of a heliocentric comet, find the inclination of the comet's orbit in relation to the ecliptic, and the positions of the nodes. Problem 5: From the time before or after the comet had reached the perihelion, and from the comet's distance to the perihelion as seen from the Sun, find the same distance in another time before or after it had appeared in the perihelion. Problem 6: Find the orbit of a comet from three given heliocentric longitudes and latitudes. From these problems, several corollaries and scholia are derived.

With approximately $50$ binary black hole events detected by LIGO/Virgo to date and many more expected in the next few years, gravitational-wave astronomy is shifting from individual-event analyses to population studies aimed at understanding the formation scenarios of these sources. There is strong evidence that the black hole mergers detected so far belong to multiple formation channels. We perform a hierarchical Bayesian analysis on the GWTC-2 catalog using a combination of ab-initio astrophysical formation models (including common envelope, globular clusters, and nuclear star clusters) as well as a realistic population of primordial black holes formed in the early universe. The evidence for a primordial population is decisively favored compared to the null hypothesis and the inferred fraction of primordial black holes in the current data is estimated at $0.27^{+0.28}_{-0.24}$ ($90\%$ credible interval), a figure which is robust against different assumptions on the astrophysical populations. The primordial formation channel can explain events in the upper mass gap such as GW190521, which are in tension with astrophysical formation scenarios. Our results suggest the tantalizing possibility that LIGO/Virgo may have already detected black holes formed after inflation. This conclusion can ultimately be confirmed in the era of third-generation interferometers.

We propose a pseudo-Goldstone boson dark matter (pGDM) particle in $SO(10)$ grand unified theory (GUT). Due to its Goldstone nature, this pGDM evades the direct DM detection experiments which, otherwise, severely constrain the parameter space of DM models. In $SO(10)$, the pGDM is embedded as a linear combination of the Standard Model (SM) singlet scalars in ${\bf 16_H}$ and ${\bf 126_H}$ representations. We consider two scenarios for the intermediate route of $SO(10)$ symmetry breaking (SB) to the SM: $SU(5) \times U(1)_X$ and Pati-Salam the $SU(4)_c \times SU(2)_L \times SU(2)_R$ (4-2-2) gauge groups. The vacuum expectation value of ${\bf 126_H}$, which triggers the breaking of $U(1)_X$ and 4-2-2 symmetry in the two scenarios, respectively, determines the pGDM lifetime whose astrophysical lower bound provides one of the most stringent constraints. For the 4-2-2 route to $SO(10)$, the successful SM gauge coupling unification requires the 4-2-2 breaking scale to be ${\cal O} (10^{11})$ GeV, and most of the parameter space is excluded. For the $SU(5) \times U(1)_X$ route, on the other hand, the $U(1)_X$ breaking scale can be significantly higher, and a wide range of the parameter space is allowed. Furthermore, the proton lifetime in the $SU(5)$ case is predicted to be $4.53 \times 10^{34}$ years, which lies well within the sensitivity reach of the Hyper-Kamiokande experiment. We also examine the constraints on the model parameter space from the Large Hadron Collider and the indirect DM search by Fermi-LAT and MAGIC experiments.

We argue that, for generic string compactifications, dark energy is likely to signal the beginning of the end of our universe, perhaps even through decompactification, with possible implications for the cosmological coincidence problem. Thanks to the scarcity (absence?) of stable de Sitter vacua, dark energy in string theory is expected to take the form of a quintessence field in slow roll. As it rolls, a tower of heavy states will generically descend, triggering an apocalyptic phase transition in the low energy cosmological dynamics after at most a few hundred Hubble times. As a result, dark energy domination cannot continue indefinitely and there is at least a percentage chance that we find ourselves in the first Hubble epoch. We use a toy model of quintessence coupled to a tower of heavy states to explicitly demonstrate the breakdown in the cosmological dynamics as the tower becomes light. This occurs through a large number of corresponding particles being produced after a certain time, overwhelming quintessence. We also discuss some implications for early universe inflation.