Papers for Thursday, Apr 01 2021

Context. The birth environments of planetary systems are thought to influence planet formation and orbital evolution, through external photoevaporation and stellar flybys. Recent work has claimed observational support for this, in the form of a correlation between the properties of planetary systems and the local Galactic phase space density of the host star. In particular, Hot Jupiters are found overwhelmingly around stars in regions of high phase space density, which may reflect a formation environment with high stellar density. Aims. We instead investigate whether the high phase density may have a galactic kinematic origin: Hot Jupiter hosts may be biased towards being young and therefore kinematically cold, because tidal inspiral leads to the destruction of the planets on Gyr timescales, and the velocity dispersion of stars in the Galaxy increases on similar timescales. Methods. We use 6D positions and kinematics from Gaia for the Hot Jupiter hosts and their neighbours, and construct distributions of the phase space density. We investigate correlations between the stars' local phase space density and peculiar velocity. Results. We find a strong anticorrelation between the phase space density and the host star's peculiar velocity with respect to the Local Standard of Rest. Therefore, most stars in "high-density" regions are kinematically cold, which may be caused by the aforementioned bias towards detecting Hot Jupiters around young stars before the planets' tidal destruction. Conclusions. We do not find evidence in the data for Hot Jupiter hosts preferentially being in phase space overdensities compared to other stars, nor therefore for their originating in birth environments of high stellar density.

The complex structure of gas, metals, and dust in the interstellar and circumgalactic medium (ISM and CGM, respectively) in star-forming galaxies can be probed by Ly$\alpha$ emission and absorption, low-ionization interstellar (LIS) metal absorption, and dust reddening $E(B-V)$. We present a statistical analysis of the correlations among Ly$\alpha$ equivalent width (EW$_{Ly\alpha}$), LIS equivalent width (EW$_{LIS}$), and $E(B-V)$ in a sample of 159 star-forming galaxies at $z\sim2.3$. With measurements obtained from individual, deep rest-UV spectra and spectral-energy distribution (SED) modeling, we find that the tightest correlation exists between EW$_{LIS}$ and $E(B-V)$, although correlations among all three parameters are statistically significant. These results signal a direct connection between dust and metal-enriched HI gas, and that they are likely co-spatial. By comparing our results with the predictions of different ISM/CGM models, we favor a clumpy ISM/CGM model where dust resides in neutral gas clumps and Ly$\alpha$ photons escape through the low HI covering fraction/column density intra-clump medium. Finally, we investigate the factors that potentially contribute to the intrinsic scatter in the correlations studied in this work, including metallicity, outflow kinematics, Ly$\alpha$ production efficiency, and slit loss. Specifically, we find evidence that scatter in the relationship between EW$_{Ly\alpha}$ and $E(B-V)$ reflects the variation in metal-to-HI covering fraction ratio as a function of metallicity, and the effects of outflows on the porosity of the ISM/CGM. Future simulations incorporating star-formation feedback and the radiative transfer of Ly$\alpha$ photons will provide key constraints on the spatial distributions of neutral gas and dust in the ISM/CGM structure.

Planck data provide precise constraints on cosmological parameters when assuming the base $\Lambda$CDM model, including a 0.17\% measurement of the age of the Universe, $t_0=13.797 \pm 0.023\,{\rm Gyr}$. However, the persistence of the "Hubble tension" calls the base $\Lambda$CDM model's completeness into question and has spurred interest in models such as Early Dark Energy (EDE) that modify the assumed expansion history of the Universe. We investigate the effect of EDE on the redshift-time relation $z \leftrightarrow t$ and find that it differs from the base $\Lambda$CDM model by at least ${\approx} 4\%$ at all $t$ and $z$. As long as EDE remains observationally viable, any inferred $t \leftarrow z$ or $z \leftarrow t$ quoted to a higher level of precision do not reflect the current status of our understanding of cosmology. This uncertainty has important astrophysical implications: the reionization epoch - $10>z>6$ - corresponds to disjoint lookback time periods in the base $\Lambda$CDM and EDE models, and the EDE value of $t_0=13.25 \pm 0.17~{\rm Gyr}$ is in tension with published ages of some stars, star clusters, and ultra-faint dwarf galaxies. However, most published stellar ages do not include an uncertainty in accuracy (due to, e.g., uncertain distances and stellar microphysics) that is estimated to be $\sim7-10\%$, potentially reconciling stellar ages with $t_{0,\rm EDE}$. We discuss how the big data era for stars is providing extremely precise ages ($<1\%$) and how improved distances and treatment of stellar microphysics such as convection could result in ages accurate to $4-5\%$, comparable to the current accuracy of $t \leftrightarrow z$. Such precise and accurate stellar ages can provide detailed insight into the high-redshift Universe independent of a cosmological model.

It is well-established that the magnitude of the incident stellar flux is the single most important factor in determining the day-night temperature gradients and atmospheric chemistries of short-period gas giant planets. However it is likely that other factors, such as planet-to-planet variations in atmospheric metallicity, C/O ratio, and cloud properties, also contribute to the observed diversity of infrared spectra for this population of planets. In this study we present new 3.6 and 4.5 micron secondary eclipse measurements for five transiting gas giant planets: HAT-P-5b, HAT-P-38b, WASP-7b, WASP-72b, and WASP-127b. We detect eclipses in at least one bandpass for all five planets and confirm circular orbits for all planets except for WASP-7b, which shows evidence for a non-zero eccentricity. Building on the work of Garhart et al. (2020), we place these new planets into a broader context by comparing them with the sample of all planets with measured Spitzer secondary eclipses. We find that incident flux is the single most important factor for determining the atmospheric chemistry and circulation patterns of short-period gas giant planets. Although we might also expect surface gravity and host star metallicity to play a secondary role, we find no evidence for correlations with either of these two variables.

In order to study the circumgalactic medium (CGM) of galaxies we develop an automated pipeline to estimate the optical continuum of quasars and detect intervening metal absorption line systems with a matched kernel convolution technique and adaptive S/N criteria. We process $\sim$ one million quasars in the latest Data Release 16 (DR16) of the Sloan Digital Sky Survey (SDSS) and compile a large sample of $\sim$ 160,000 MgII absorbers, together with $\sim$ 70,000 FeII systems, in the redshift range $0.35<z_{abs}<2.3$. Combining these with the SDSS DR16 spectroscopy of $\sim1.1$ million luminous red galaxies (LRGs) and $\sim 200,000$ emission line galaxies (ELGs), we investigate the nature of cold gas absorption at $0.5<z<1$. These large samples allow us to characterize the scale dependence of MgII with greater accuracy than in previous work. We find that there is a strong enhancement of MgII absorption within $\sim 50$ kpc of ELGs, and the covering fraction within $0.5r_{\rm vir}$ of ELGs is 2-5 times higher than for LRGs. Beyond 50 kpc, there is a sharp decline in MgII for both kinds of galaxies, indicating a transition to the regime where the CGM is tightly linked with the dark matter halo. The MgII covering fraction correlates strongly with stellar mass for LRGs, but weakly for ELGs, where covering fractions increase with star formation rate. Our analysis implies that cool circumgalactic gas has a different physical origin for star forming versus quiescent galaxies.

Observational constraints on planet spin-axis has recently become possible, and revealed a system that favors a large spin-axis misalignment, a low stellar spin-orbit misalignment and a high eccentricity. To explain the origin of such systems, we propose a mechanism that could tilt the planet spin-axis during planet-planet scattering, which are natural outcomes of in-situ formation and disk migration. Specifically, we show that secular spin-orbit resonances could occur for a short time period during the scattering processes, and excite the misalignment of the planet spin-axis. This typically leads to planets with large spin-misalignment and a wide range of eccentricities and inclinations.

Imaging faint objects, such as exoplanets or disks, around nearby stars is extremely challenging because host star images are dominated by the telescope diffraction pattern. Using a coronagraph is an efficient solution for removing diffraction but requires an incoming wavefront with good quality to maximize starlight rejection. On the ground, the most advanced exoplanet imagers use extreme adaptive optics (ExAO) systems that are based on a deformable mirror (DM) with a large number of actuators to efficiently compensate for high-order aberrations and provide diffraction-limited images. While several exoplanet imagers with DMs using around 1500 actuators are now routinely operating on large telescopes to observe gas giant planets, future systems may require a tenfold increase in the number of degrees of freedom to look for rocky planets. In this paper, we explore wavefront correction with a secondary adaptive optics system that controls a very large number of degrees of freedom that are not corrected by the primary ExAO system. Using Marseille Imaging Testbed for High Contrast (MITHiC), we implement a second stage of adaptive optics with ZELDA, a Zernike wavefront sensor, and a spatial light modulator (SLM) to compensate for the phase aberrations of the bench downstream residual aberrations from adaptive optics. We demonstrate that their correction up to 137 cycles per pupil with nanometric accuracy is possible, provided there is a simple distortion calibration of the pupil and a moderate wavefront low-pass filtering. We also use ZELDA for a fast compensation of ExAO residuals, showing its promising implementation as a second-stage correction for the observation of rocky planets around nearby stars. Finally, we present images with a classical Lyot coronagraph on MITHiC and validate our ability to reach its theoretical performance with our calibration.

The process of momentum and energy transfer from a massive body moving through a background medium, known as dynamical friction (DF), is key to our understanding of many astrophysical systems. We present a series of high-resolution simulations of gaseous DF using Lagrangian hydrodynamics solvers, in the state-of-the-art multi-physics code, GIZMO. The numerical setup is chosen to allow direct comparison to analytic predictions for DF in the range of Mach 0.2<M<3. We investigate, in detail, the DF drag force, the radial structure of the wake, and their time evolution. The subsonic forces are shown to be well resolved, except at Mach numbers close to M=1. The supersonic cases, close to M=1, significantly under-shoot the predicted force. We find that for scenarios with 0.7<M<2, between 10%-50% of the expected DF force is missing. The origin of this deficit is mostly related to the wake structure close to the perturber, where the density profile of the Mach cone face shows significant smoothing, which does not improve with time. The spherical expanding perturbation of the medium is captured well by the hydro scheme, but it is the sharp density structure, at the transition from Mach cone to average density, that introduces the mismatch. However, we find a general improvement of the force deficit with time, though significant differencesremain, in agreement with other studies. This is due to (1) the structure of the far field wake being reproduced well, and (2) the fraction of total drag from the far field wake increasing with time. Dark matter sub-haloes, in typical cosmological simulations, occupy parameters similar to those tested here, suggesting that the DF which these sub-haloes experience is significantly underestimated, and hence their merger rate. Dynamical friction is a relevant benchmark and should be included as one of the standard hydro tests for astrophysical simulations.

Understanding the escape of ionizing (Lyman continuum) photons from galaxies is vital for determining how galaxies contributed to reionization in the early universe. While directly detecting Lyman continuum from high redshift galaxies is impossible due to the intergalactic medium, low redshift galaxies in principle offer this possibility, but requirie observations from space. The first local galaxy for which Lyman continuum escape was found is Haro11 , a luminous blue compact galaxy at z=0.02, where observations with the FUSE satellite revealed an escape fraction of 3.3 %. However the FUSE aperture covers the entire galaxy, and it is not clear from where the Lyman continuum is leaking out. Here we utilize HST/COS spectroscopy in the wavelength range 1100-1700 A of the three knots (A, B, and C) of Haro11 to study the presence of Ly-$\alpha$ emission and the properties of intervening gas. We find that all knots have bright Ly-$\alpha$ emission. UV absorption lines, originating in the neutral interstellar medium, as well as lines probing the ionized medium, are seen extending to blue shifted velocities of 500 km/s in all three knots, demonstrating the presence of an outflowing multiphase medium. We find that knots A and B have large covering fractions of neutral gas, making LyC escape along these sightlines improbable, while knot C has a much lower covering fraction ($\lesssim50$%). Knot C also has the the highest Ly-$\alpha$ escape fraction and we conclude that it is the most likely source of the escaping Lyman continuum detected in Haro11.

Variable stars in the compact elliptical galaxy M32 are identified, using three epochs of photometry from the Spitzer Space Telescope at 3.6 and 4.5 $\mu$m, separated by 32 to 381 days. We present a high-fidelity catalogue of sources detected in multiple epochs at both 3.6 and 4.5 $\mu$m, which we analysed for stellar variability using a joint probability error-weighted flux difference. Of these, 83 stars are identified as candidate large-amplitude, long-period variables, with 28 considered high-confidence variables. The majority of the variable stars are classified as asymptotic giant branch star candidates using colour-magnitude diagrams. We find no evidence supporting a younger, infrared-bright stellar population in our M32 field.

We have selected Abell 3266 to search for ram-pressure induced star formation as a global property of a merging cluster. Abell 3266 (z = 0.0594) is a high mass cluster that features a high velocity dispersion, an infalling subcluster near to the line of sight, and a strong shock front. These phenomena should all contribute to making Abell 3266 an optimum cluster to see the global effects of RPS induced star formation. Using archival X-ray observations and published optical data, we cross-correlate optical spectral properties ([OII, H$\beta$]), indicative of starburst and post starburst, respectively with ram-pressure, $\rho$v$^{2}$, calculated from the X-ray and optical data. We find that post-starburst galaxies, classified as E+A, occur at a higher frequency in this merging cluster than in the Coma cluster and at a comparable rate to intermediate redshift clusters. This is consistent with increased star formation due to the merger. However, both starburst and post-starburst galaxies are equally likely to be in a low or high ram pressure environment. From this result we infer that the duration of the starburst phase must be very brief so that: (1) at any time only a small fraction of the galaxies in a high ram pressure environment show this effect, and (2) most post-starburst galaxies are in an environment of low ram pressure due too their continued orbital motion in the cluster.

Planet-planet scattering best explains the eccentricity distribution of extrasolar giant planets. Past literature showed that the orbits of planets evolve due to planet-planet scattering. This work studies the spin evolution of planets in planet-planet scattering in 2-planet systems. Spin can evolve dramatically due to spin-orbit coupling made possible by the evolving spin and orbital precession during the planet-planet scattering phase. The main source of torque to planet spin is the stellar torque, and the total planet-plane torque contribution is negligible. As a consequence of the evolution of the spin, planets can end up with significant obliquity (the angle between a planet's own orbit normal and spin axis) like planets in our Solar System.

Context. Radiation-driven mass loss is key to our understanding of massive-star evolution. However, for low-luminosity O-type stars there are big discrepancies between theoretically predicted and empirically derived mass-loss rates (called the weak-wind problem). Aims. We compute radiation-line-driven wind models of a typical weak-wind star to determine its temperature structure and the corresponding impact on ultra-violet (UV) line formation. Methods. We carried out hydrodynamic simulations of the line-deshadowing instability (LDI) for a weak-wind star in the Galaxy. Subsequently, we used this LDI model as input in a short-characteristics radiative transfer code to compute synthetic UV line profiles. Results. We find that the line-driven weak wind is significantly shock heated to high temperatures and is unable to cool down effciently. This results in a complex temperature structure where more than half of the wind volume has temperatures significantly higher than the stellar effective temperature. Therefore, a substantial portion of the weak wind will be more ionised, resulting in a reduction of the UV line opacity and therefore in weaker line profiles for a given mass-loss rate. Quantifying this, we find that weak-wind mass-loss rates derived from unsaturated UV lines could be underestimated by a factor of between 10 and 100 if the high-temperature gas is not properly taken into account in the spectroscopic analysis. This offers a tentative basic explanation for the weak-wind problem: line-driven weak winds are not really weaker than theoretically expected, but rather a large portion of their wind volume is much hotter than the stellar effective temperature.

Kepler-30 is a unique target to study stellar activity and rotation in a young solar-like star accompanied by a compact planetary system. We use about 4 years of high-precision photometry collected by the Kepler mission to investigate the fluctuations caused by photospheric convection, stellar rotation, and starspot evolution as a function of the timescale. Our main goal is to apply methods for the analysis of timeseries to find the timescales of the phenomena that affect the light variations. We correlate those timescales with periodicities in the star as well as in the planetary system. We model the flux rotational modulation induced by active regions using spot modelling and apply the MFDMA in standard and multiscale versions for analysing the behaviour of variability and light fluctuations that can be associated with stellar convection and the evolution of magnetic fields on timescales ranging from less than 1 day up to about 35 days. The light fluctuations produced by stellar activity can be described by the multifractal Hurst index that provides a measure of their persistence. The spot modeling indicates a lower limit to the relative surface differential rotation of $\Delta \Omega/\Omega\sim 0.02\pm 0.01$ and suggests a short-term cyclic variation in the starspot area with a period of $\sim 34$ days, virtually close to the synodic period of 35.2 days of the planet Kepler-30b. By subtracting the two timeseries of the SAP and PDC Kepler pipelines, we reduce the rotational modulation and find a 23.1-day period close to the synodic period of Kepler-30c. This period also appears in the multifractal analysis as a crossover of the fluctuation functions associated with the characteristic evolutionary timescales of the active regions in Kepler-30 as confirmed by spot modelling. These procedures and methods may be greatly useful for analysing current TESS and future PLATO data.

In this paper, the formation of primordial black holes (PBHs) is reinvestigated using inflationary $\alpha$-attractors. Instead of using the conventional Press-Schechter theory to compute the abundance, the optimized peaks theory is used, which was developed in Ref. \cite{Yoo:2018kvb,Yoo:2020dkz}. This method takes into account how curvature perturbations play a r\^{o}le in modifying the mass of primordial black holes. Analyzing the model proposed in \cite{Mahbub:2019uhl} it is seen that the horizon mass of the collapsed Hubble patch is larger by $\mathcal{O}(10^2)$ compared to the usual computation. Moreover, PBHs can be formed from curvature power spectrum, $\mathcal{P}_{\zeta}(k)$, peaked at lower values using numerically favored threshold overdensities. As a result of the generally larger masses predicted, the peak of the power spectrum can be placed at $k>10^{15}\text{Mpc}^{-1}$ with which potential future constraints on the primordial power spectrum through gravitational waves (GWs) can be evaded.

On Titan, methane (CH4) and ethane (C2H6) are the dominant species found in the lakes and seas. In this study, we have combined laboratory work and modeling to refine the methane-ethane binary phase diagram at low temperatures and probe how the molecules interact at these conditions. We used visual inspection for the liquidus and Raman spectroscopy for the solidus. Through these methods we determined a eutectic point of 71.15$\pm$0.5 K at a composition of 0.644$\pm$0.018 methane - 0.356$\pm$0.018 ethane mole fraction from the liquidus data. Using the solidus data, we found a eutectic isotherm temperature of 72.2 K with a standard deviation of 0.4 K. In addition to mapping the binary system, we looked at the solid-solid transitions of pure ethane and found that, when cooling, the transition of solid I-III occurred at 89.45$\pm$0.2 K. The warming sequence showed transitions of solid III-II occurring at 89.85$\pm$0.2 K and solid II-I at 89.65$\pm$0.2 K. Ideal predictions were compared to molecular dynamics simulations to reveal that the methane-ethane system behaves almost ideally, and the largest deviations occur as the mixing ratio approaches the eutectic composition.

The third generation South Pole Telescope camera (SPT-3G) improves upon its predecessor (SPTpol) by an order of magnitude increase in detectors on the focal plane. The technology used to read out and control these detectors, digital frequency-domain multiplexing (DfMUX), is conceptually the same as used for SPTpol, but extended to accommodate more detectors. A nearly 5x expansion in the readout operating bandwidth has enabled the use of this large focal plane, and SPT-3G performance meets the forecasting targets relevant to its science objectives. However, the electrical dynamics of the higher-bandwidth readout differ from predictions based on models of the SPTpol system. To address this, we present an updated derivation for electrical crosstalk in higher-bandwidth DfMUX systems, and identify two previously uncharacterized contributions to readout noise. The updated crosstalk and noise models successfully describe the measured crosstalk and readout noise performance of SPT-3G, and suggest improvements to the readout system for future experiments using DfMUX, such as the LiteBIRD space telescope.

Most of our knowledge about the structure of the Milky Way has come from the study of variable stars. Among the variables, mimicking the periodic variation of pulsating stars, are the eclipsing binaries. These stars are important in astrophysics because they allow us to directly measure radii and masses of the components, as well as the distance to the system, thus being useful in studies of Galactic structure alongside pulsating RR Lyrae and Cepheids. Using the distinguishing features of their light curves, one can identify them using a semi-automated process. In this work, we present a strategy to search for eclipsing variables in the inner VVV bulge across an area of 13.4 sq. deg. within $1.68^{\rm o}<l<7.53^{\rm o}$ and $-3.73^{\rm o}<b<-1.44^{\rm o}$, corresponding to the VVV tiles b293 to b296 and b307 to b310. We accurately classify 212 previously unknown eclipsing binaries, including six very reddened sources. The preliminary analysis suggests these eclipsing binaries are located in the most obscured regions of the foreground disk and bulge of the Galaxy. This search is therefore complementary to other variable stars searches carried out at optical wavelengths.

Eccentricity has emerged as a potentially useful tool for helping to identify the origin of black hole mergers. However, owing to the large number of harmonics required to compute the amplitude of an eccentric signal, eccentric templates can be computationally very expensive, making statistical analyses to distinguish distributions from different formation channels very challenging. In this paper, we outline a method for estimating the signal-to-noise ratio for inspiraling binaries at lower frequencies such as those proposed for LISA and DECIGO. Our approximation can be useful more generally for any quasi-periodic sources. We argue that surprisingly, the signal-to-noise ratio evaluated at or near the peak frequency (of the power) is well approximated by using a constant noise curve, even if in reality the noise strain has power law dependence. We furthermore improve this initial estimate over our previous calculation to allow for frequency-dependence in the noise to expand the range of eccentricity and frequency over which our approximation applies. We show how to apply this method to get an answer accurate to within a factor of two over almost the entire projected observable frequency range. We emphasize this method is not a replacement for detailed signal processing. The utility lies chiefly in identifying theoretically useful discriminators among different populations and providing fairly accurate estimates for how well they should work. This approximation can furthermore be useful for narrowing down parameter ranges in a computationally economical way when events are observed. We furthermore show a distinctive way to identify events with extremely high eccentricity where the signal is enhanced relative to na\"ive expectations on the high frequency end.

In this paper, we present a model-independent approach to calibrate the largest quasar sample. Calibrating quasar samples is essentially constraining the parameters of the linear relation between the $\log$ of the ultraviolet (UV) and X-ray luminosities. This calibration allows quasars to be used as standardized candles. There is a strong correlation between the parameters characterizing the quasar luminosity relation and the cosmological distances inferred from using quasars as standardized candles. We break this degeneracy by using Gaussian process regression to model-independently reconstruct the expansion history of the Universe from the latest type Ia supernova observations. Using the calibrated quasar dataset, we further reconstruct the expansion history up to redshift of $z\sim 7.5$. Finally, we test the consistency between the calibrated quasar sample and the standard $\rm{\Lambda}CDM$ model based on the posterior probability distribution of the GP hyperparameters. Our results show that the quasar sample is in good agreement with the standard $\rm{\Lambda}CDM$ model in the redshift range of the supernova, despite of mildly significant deviations taking place at higher redshifts. Fitting the standard $\rm{\Lambda}CDM$ model to the calibrated quasar sample, we obtain a high value of the matter density parameter $\Omega_m = 0.382^{+0.045}_{-0.042}$, which is marginally consistent with the constraints from other cosmological observations.

So far quite a few ultraluminous X-ray (ULX) pulsars have been discovered. In this work, we construct a super-Eddington, magnetic accretion disk model to estimate the dipole magnetic field of eight ULX pulsars based on their observed spin-up variations and luminosities. We obtain two branches of dipole magnetic field solutions. They are distributed in the range of $B\sim (0.16-64.5)\times 10^{10}\,{\rm G}$ and $\sim (0.275-79.0)\times 10^{13}\,{\rm G}$ corresponding to the low- and high-$B$ solutions respectively. The low magnetic field solutions correspond to the state that the neutron stars are far away from the spin equilibrium, and the high magnetic field solutions are close to the spin equilibrium. The ultra-strong magnetic fields derived in Be-type ULX pulsars imply that the accretion mode in Be-type ULX pulsars could be more complicated than in the persistent ULX pulsars and may not be accounted for by the magnetized accretion disk model. We suggest that the transition between the accretor and the propeller regimes may be used to distinguish between the low- and high-$B$ magnetic field solutions in addition to the detection of the cyclotron resonance scattering features.

The recent report of an association of the gravitational-wave (GW) binary black hole (BBH) merger GW190521 with a flare in the Active Galactic Nuclei (AGN) J124942.3+344929 has generated tremendous excitement. However, GW190521 has one of the largest localization volumes amongst all of the GW events detected so far. The 90\% localization volume likely contains $7,400$ unobscured AGN brighter than $g \leq 20.5$ AB mag, and it results in a $\gtrsim 70\%$ probability of chance coincidence for an AGN flare consistent with the GW event. We present a Bayesian formalism to estimate the confidence of an AGN association by analyzing a population of BBH events with dedicated follow-up observations. Depending on the fraction of BBH arising from AGNs, counterpart searches of $\mathcal{O}(1)-\mathcal{O}(100)$ GW events are needed to establish a confident association, and more than an order of magnitude more for searches without followup (i.e, using only the locations of AGNs and GW events). Follow-up campaigns of the top $\sim 5\%$ (based on volume localization and binary mass) of BBH events with total rest frame mass $\ge 50~M_\odot$ are expected to establish a confident association during the next LIGO/Virgo/KAGRA observing run (O4), as long as the true value of the fraction of BBH giving rise to AGN flares is $>0.1$. Our formalism allows us to jointly infer cosmological parameters from a sample of BBH events that include chance coincidence flares. Until the confidence of AGN associations is established, the probability of chance coincidence must be taken into account to avoid biasing astrophysical and cosmological constraints.

Searching in the MaNGA IFU survey, I identify 9 galaxies that have strong Balmer absorption lines and weak nebular emission lines measured from the spectra integrated over the entire IFUs. The spectral features measured from the bulk of the stellar light make these galaxies local analogs of high-redshift spectroscopically-selected post-starburst galaxies, thus, proxies to understand the mechanisms producing post-starburst galaxies at high-redshifts. I present the distributions of absorption-line indices and emission-line strengths, as well as stellar kinematics of these local post-starburst galaxies. Almost all local post-starburst galaxies have central compact emission-line regions at the central $<1$ kpc, mostly powered by weak star-formation activities. The age-sensitive absorption line indices EW(H$\delta$) and Dn4000 indicate that the stellar populations at the outskirts are older. Toy stellar population synthesis models suggest that the entire galaxies are experiencing a rapid decline of star formation with residual star-formation activities at the centers. These features demand that most post-starburst galaxies are the aftermath of highly dissipative processes that drive gas into centers, invoke centrally-concentrated star formation, then quench the galaxies. Meanwhiles, when measurable, post-starburst galaxies have the directions of maximum stellar velocity gradients align with photometric major axes, which suggest against major mergers being the principal driving mechanism, while gas-rich minor mergers are plausible. While directly obtaining the same quality of spatially-resolved spectra of high-redshift post-starburst galaxies is very difficult, finding proper local counterparts provides an alternative to understand quenching processes in the distant Universe.

The Interface Region Imaging Spectrograph (IRIS) has been obtaining near- and far-ultraviolet images and spectra of the solar atmosphere since July 2013. The unique combination of near and far-ultraviolet spectra and images at subarcsecond resolution and high cadence allows the tracing of mass and energy through the critical interface between the solar surface and the corona or solar wind. IRIS has enabled research into the fundamental physical processes thought to play a role in the low solar atmosphere such as ion-neutral interactions, magnetic reconnection, the generation, propagation, and dissipation of various types of waves, the acceleration of non-thermal particles, and various small-scale instabilities. These new findings have helped provide novel insights into a wide range of phenomena including the discovery of non-thermal particles in coronal nanoflares, the formation and impact of spicules and other jets, resonant absorption and dissipation of Alfvenic waves, energy release associated with braiding of magnetic field lines, the thermal instability in the chromosphere-corona mass and energy cycle, the contribution of waves, turbulence, and non-thermal particles in the energy deposition during flares and smaller-scale events such as UV bursts, and the role of flux ropes and other mechanisms in triggering CMEs. IRIS observations have also advanced studies of the connections between solar and stellar physics. Advances in numerical modeling, inversion codes, and machine learning techniques have played a key role in driving these new insights. With the advent of exciting new instrumentation both on the ground (e.g., DKIST, ALMA) and space-based (e.g., Parker Solar Probe, Solar Orbiter), we aim to review new insights based on IRIS observations or related modeling, and highlight some of the outstanding challenges that have been brought to the fore.

The Martian Moons Exploration (MMX) spacecraft is a JAXA mission to Mars and its moons Phobos and Deimos. MMX will carry the Circum-Martian Dust Monitor (CMDM) which is a newly developed light-weight ($\mathrm{650\,g}$) large area ($\mathrm{1\,m^2}$) dust impact detector. Cometary meteoroid streams (also referred to as trails) exist along the orbits of comets, forming fine structures of the interplanetary dust cloud. The streams consist predominantly of the largest cometary particles (with sizes of approximately $\mathrm{100\,\mu m}$ to 1~cm) which are ejected at low speeds and remain very close to the comet orbit for several revolutions around the Sun. The Interplanetary Meteoroid Environment for eXploration (IMEX) dust streams in space model is a new and recently published universal model for cometary meteoroid streams in the inner Solar System. We use IMEX to study the detection conditions of cometary dust stream particles with CMDM during the MMX mission in the time period 2024 to 2028. The model predicts traverses of 12 cometary meteoroid streams with fluxes of $\mathrm{100\,\mu m}$ and bigger particles of at least $\mathrm{10^{-3}\,m^{-2}\,day^{-1}}$ during a total time period of approximately 90~days. The highest flux of $\mathrm{0.15\,m^{-2}\,day^{-1}}$ is predicted for comet 114P/Wiseman-Skiff in October 2026. With its large detection area and high sensitivity CMDM will be able to detect cometary meteoroid streams en route to Phobos. Our simulation results for the Mars orbital phase of MMX also predict the occurrence of meteor showers in the Martian atmosphere which may be observable from the Martian surface with cameras on board landers or rovers. Finally, the IMEX model can be used to study the impact hazards imposed by meteoroid impacts on to large-area spacecraft structures that will be particularly necessary for crewed deep space missions.

Measurement of magnetic field strengths in a molecular cloud is essential for determining the criticality of magnetic support against gravitational collapse. In this paper, as part of the JCMT BISTRO survey, we suggest a new application of the Davis-Chandrasekhar-Fermi (DCF) method to estimate the distribution of magnetic field strengths in the OMC-1 region. We use observations of dust polarization emission at 450 $\mu$m and 850 $\mu$m, and C$^{18}$O (3-2) spectral line data obtained with the JCMT. We estimate the volume density, the velocity dispersion and the polarization angle dispersion in a box, 40$''$ $\times$ 40$''$ (5$\times$5 pixels), which moves over the OMC-1 region. By substituting three quantities in each box to the DCF method, we get magnetic field strengths over the OMC-1 region. We note that there are very large uncertainties in inferred field strengths, as discussed in detail in this paper. The field strengths vary from 0.8 to 26.4 mG and their mean value is about 6 mG. Additionally, we obtain maps of the mass-to-flux ratio in units of a critical value and the Alfv$\acute{e}$n mach number. The central parts of the BN-KL and South (S) clumps in the OMC-1 region are magnetically supercritical, so the magnetic field cannot support the clumps against gravitational collapse. However, the outer parts of the region are magnetically subcritical. The mean Alfv$\acute{e}$n mach number is about 0.4 over the region, which implies that the magnetic pressure exceeds the turbulent pressure in the OMC 1 region.

Cataclysmic variable SDSS J143317.78+101123.3 with P=0.054241 d (78.1-min) was suspected to be a possible dwarf nova of WZ Sge type with the sub-stellar donor, but without detected outbursts so far. Checking the newly available data from ATLAS survey has revealed the outburst by at least 6 magnitudes in September 2020, thus confirming the dwarf nova nature of this object with the brown dwarf secondary. Other projects and individual observers have stopped their monitoring of this target several days before the outburst. This finding strengthens the value of observing the twilight zone by the professional surveys and amateurs.

Context. Recent observations have shown that magnetic flux cancellation occurs at the photosphere more frequently than previously thought. Aims. In order to understand the energy release by reconnection driven by flux cancellation, we previously studied a simple model of two cancelling polarities of equal flux. Here, we further develop our analysis to achieve a more general setup where the two cancelling polarities have unequal magnetic fluxes and where many new features are revealed. Methods. We carried out an analytical study of the cancellation of two magnetic fragments of unequal and opposite flux that approach one another and are located in an overlying horizontal magnetic field. Results. The energy release as microflares and nanoflares occurs in two main phases. During phase 1a, a separator is formed and reconnection is driven at it as it rises to a maximum height and then moves back down to the photosphere, heating the plasma and accelerating plasma jets in the process. During phase 1b, once the separator moves back to the photosphere, it bifurcates into two null points. Reconnection is no longer driven at the separator and an isolated magnetic domain connecting the two polarities is formed. During phase 2, the polarities cancel out at the photosphere as magnetic flux submerges below the photosphere and as reconnection occurs at and above the photosphere and plasma jets and a mini-filament eruption can be produced.

The recent launch of Solar Orbiter and BepiColombo opened a brief window in which these two spacecraft were positioned in a constellation that allows for the detailed sampling of any Earth-directed CMEs. Fortunately, two such events occurred with in situ detections of an ICME by Solar Orbiter on the 19th of April and the 28th of May 2020. These two events were subsequently also observed in situ by BepiColombo and Wind around a day later. We attempt to reconstruct the observed in situ magnetic field measurements for all three spacecraft simultaneously using an empirical magnetic flux rope model. This allows us to test the validity of our flux rope model on a larger and more global scale and allows for cross-validation of the analysis with different spacecraft combinations. Finally, we can also compare the results from the in situ modeling to remote observations obtained from the STEREO-A heliospheric imagers. We make use of the 3D coronal rope ejection model in order to simulate the ICME evolution. We adapt a previously developed ABC-SMC fitting algorithm for the application to multi point scenarios. We show that we are able to generally reconstruct the flux ropes signatures at three different spacecraft positions simultaneously using our model in combination with the flux rope fitting algorithm. For the well-behaved 19th of April ICME our approach works very well. The 28th of May ICME, on the other hand, shows the limitations of our approach. Unfortunately, the usage of multi point observations for these events does not appear to solve inherent issues, such as the estimation of the magnetic field twist or flux rope aspect-ratios due to the specific constellation of the spacecraft positions. As our general approach can be used for any fast forward simulation based model we give a blueprint for future studies using more advanced ICME models.

In the direct detection of the galactic dark matter, experiments using cryogenic solid-state detectors or noble liquids play for years a very relevant role, with increasing target mass and more and more complex detection systems. But smaller projects, based on very sensitive, advanced detectors following new technologies, could help in the exploration of the different proposed dark matter scenarios too. There are experiments focused on the observation of distinctive signatures of dark matter, like an annual modulation of the interaction rates or the directionality of the signal; other ones are intended to specifically investigate low mass dark matter candidates or particular interactions. For this kind of dark matter experiments at small scale, the physics case will be discussed and selected projects will be described, summarizing the basics of their detection methods and presenting their present status, recent results and prospects.

In Budai et al. (2020) we argued that angular non-stationarities of gamma-ray burst (GRB) jets can result in a statistical connection between the angle values deduced from jet break times and the variabilities of prompt light curves. The connection should be an anti-correlation if luminosity densities of jets follow a power-law or a uniform profile, and a correlation if they have a Gaussian profile. In this follow-up paper, we search for the connection by measuring Spearman's rank correlation coefficient in a sample of 19 long GRBs observed by the Swift satellite. Using 16 of the GRBs with well-defined angle measurements, we find $\rho = -0.20_{-0.09}^{+0.09}$ and $p = 0.46_{-0.19}^{+0.23}$. Adding three more GRBs to the sample, each with a pair of equally possible angle values, can strengthen the anti-correlation to $\rho=-0.31_{-0.08}^{+0.07}$ and $p=0.19_{-0.09}^{+0.14}$. We show that these results are incompatible with non-stationary jets having Gaussian profiles, and that $\gtrsim\!100$ GRBs with observed afterglows would be needed to confirm the potential existence of the angle-variability anti-correlation with $3\sigma$ significance. If the connection is real, GRB jet angles would be constrainable from prompt gamma light curves, without the need of afterglow observations.

Cosmological inference from cluster number counts is systematically limited by the accuracy of the mass calibration, i.e. the empirical determination of the mapping between cluster selection observables and halo mass. In this work we demonstrate a method to quantitatively determine the bias and uncertainties in weak-lensing mass calibration. To this end, we extract a library of projected matter density profiles from hydrodynamical simulations. Accounting for shear bias and noise, photometric redshift uncertainties, mis-centering, cluster member contamination, cluster morphological diversity, and line-of-sight projections, we produce a library of shear profiles. Fitting a one-parameter model to these profiles, we extract the so-called \emph{weak lensing mass} $M_\text{WL}$. Relating the weak-lensing mass to the halo mass from gravity-only simulations with the same initial conditions as the hydrodynamical simulations allows us to estimate the impact of hydrodynamical effects on cluster number counts experiments. Creating new shear libraries for $\sim$1000 different realizations of the systematics, provides a distribution of the parameters of the weak-lensing to halo mass relation, reflecting their systematic uncertainty. This result can be used as a prior for cosmological inference. We also discuss the impact of the inner fitting radius on the accuracy, and determine the outer fitting radius necessary to exclude the signal from neighboring structures. Our method is currently being applied to different Stage~III lensing surveys, and can easily be extended to Stage~IV lensing surveys.

Misalignments within protoplanetary discs are now commonly observed, and features such as shadows in scattered light images indicate departure from a co-planar geometry. VLT/SPHERE observations of the disc around HD 143006 show a large-scale asymmetry, and two narrow dark lanes which are indicative of shadowing. ALMA observations also reveal the presence of rings and gaps in the disc, along with a bright arc at large radii. We present new hydrodynamic simulations of HD 143006, and show that a configuration with both a strongly inclined binary and an outer planetary companion is the most plausible to explain the observed morphological features. We compute synthetic observations from our simulations, and successfully reproduce both the narrow shadows and the brightness asymmetry seen in IR scattered light. Additionally, we reproduce the large dust observed in the mm continuum, due to a 10 Jupiter mass planet detected in the CO kinematics. Our simulations also show the formation of a circumplanetary disc, which is misaligned with respect to the outer disc. The narrow shadows cast by the inner disc and the planet-induced "kink" in the disc kinematics are both expected to move on a time-scale of $\sim$ 5-10 years, presenting a potentially observable test of our model. If confirmed, HD 143006 would be the first known example of a circumbinary planet on a strongly misaligned orbit.

Collision and disruption processes of proto-planetary bodies in the early solar system are key to understanding the genesis of diverse types of main-belt asteroids. Mesosiderites are stony-iron meteorites that formed by mixing of howardite-eucrite-diogenite-like crust and molten core materials and provide unique insights into the catastrophic break-up of differentiated asteroids. However, the enigmatic formation process and the poorly constrained timing of metal-silicate mixing complicate the assignment to potential parent bodies. Here we report high-precision uranium-lead dating of mesosiderite zircons by isotope dilution thermal ionization mass spectrometry, revealing initial crust formation 4,558.5 +/- 2.1 million years ago and metal-silicate mixing at 4,525.39 +/- 0.85 million years. The two distinct ages coincide with the timing of crust formation and a large-scale reheating event on the eucrite parent body, likely the asteroid Vesta. This chronological coincidence corroborates that Vesta is the parent body of mesosiderite silicates. Mesosiderite formation on Vesta can be explained by a hit-and-run collision 4,525.4 million years ago that caused the thick crust observed by NASA's Dawn mission and explains the missing olivine in mesosiderites, howardite-eucrite-diogenite meteorites, and vestoids.

We study the decay of cosmic string loops in the Abelian-Higgs model. We confirm earlier results that loops formed by intersections of infinite strings formed from random-field initial conditions disappear quickly, with lifetime proportional to their initial rest-frame length $\ell_\text{init}$. We study a population with $\ell_\text{init}$ up to $6000$ inverse mass units, and measure the proportionality constant to be $0.14\pm0.04$, independently of the initial lengths. We propose a new method to construct oscillating non-self intersecting loops from initially stationary strings, and show that by contrast these loops have lifetimes scaling approximately as $\ell_\text{init}^2$, in line with previous works on artificially created string configurations. We show that the oscillating strings have mean-square velocity $\bar{v}^2 \simeq 0.500 \pm 0.004$, consistent with the Nambu-Goto value of $1/2$, while the network loops have $\bar{v}^2 \simeq 0.40 \pm 0.04$. We argue that whatever the mechanism behind the network loop decay is, it is non-linear, can only be suppressed by careful tuning of initial conditions, and is much stronger than gravitational radiation. An implication is that one cannot use the Nambu-Goto model to derive robust constraints on the tension of field theory strings. We advocate parametrising the uncertainty as the fraction $f_\text{NG}$ of Nambu-Goto-like loops surviving to radiate gravitationally. None of the 31 large network loops created survived longer than 0.25 of their initial length, so one can estimate that $f_\text{NG}<0.1$ at $95$% confidence level. If the recently reported NANOgrav signal is due to cosmic strings, $f_\text{NG}$ must be greater than $10^{-3}$ in order not to violate bounds from the Cosmic Microwave Background.

New radio telescope arrays offer unique opportunities for large-scale commensal SETI surveys. Ethernet-based architectures are allowing multiple users to access telescope data simultaneously by means of multicast Ethernet subscriptions. Breakthrough Listen will take advantage of this by conducting a commensal SETI survey on the MeerKAT radio telescope in South Africa. By subscribing to raw voltage data streams, Breakthrough Listen will be able to beamform commensally anywhere within the field of view during primary science observations. The survey will be conducted with unprecedented speed by forming and processing 64 coherent beams simultaneously, allowing the observation of several million objects within a few years. Both coherent and incoherent observing modes are planned. We present the list of desired sources for observation and explain how these sources were selected from the Gaia DR2 catalog. Given observations planned by MeerKAT's primary telescope users, we discuss their effects on the commensal survey and propose a commensal observing strategy in response. Finally, we outline our proposed approach towards observing one million nearby stars and analyse expected observing progress in the coming years.

We present here a self-consistent cosmological zoom-in simulation of a triple supermassive black hole (SMBH) system forming in a complex multiple galaxy merger. The simulation is run with an updated version of our code KETJU, which is able to follow the motion of SMBHs down to separations of tens of Schwarzschild radii while simultaneously modeling the large-scale astrophysical processes in the surrounding galaxies, such as gas cooling, star formation, and stellar and AGN feedback. Our simulation produces initially a SMBH binary system for which the hardening process is interrupted by the late arrival of a third SMBH. The KETJU code is able to accurately model the complex behavior occurring in such a triple SMBH system, including the ejection of one SMBH to a kiloparsec-scale orbit in the galaxy due to strong three-body interactions as well as Lidov-Kozai oscillations suppressed by relativistic precession when the SMBHs are in a hierarchical configuration. One pair of SMBHs merges $\sim 3\,\mathrm{Gyr}$ after the initial galaxy merger, while the remaining binary is at a parsec-scale separation when the simulation ends at redshift $z=0$. We also show that KETJU can capture the effects of the SMBH binaries and triplets on the surrounding stellar population, which can affect the binary merger timescales as the stellar density in the system evolves. Our results demonstrate the importance of dynamically resolving the complex behavior of multiple SMBHs in galactic mergers, as such systems cannot be readily modeled using simplified semi-analytic models.

The habitable zone is the region around a star where standing bodies of liquid water can be stable on a planetary surface. Its width is often assumed to be dictated by the efficiency of the carbonate-silicate cycle, which has maintained habitable surface conditions on our planet for billions of years. This cycle may be inhibited by surface condensation of significant amounts of $CO_{\mathrm{2}}$ ice, which is likely to occur on distant planets containing high enough levels of atmospheric $CO_{\mathrm{2}}$. Such a process could permanently trap $CO_{\mathrm{2}}$ ice within the planet, threatening its long-term habitability. Recent work has modeled this scenario for initially cold and icy planetary bodies orbiting the Sun. Here, we use an advanced energy balance model to consider both initially warm and cold rapidly-rotating planets orbiting F - K stars. We show that the range of orbital distances where significant surface $CO_{\mathrm{2}}$ ice condensation occurs is significantly reduced for warm start planets. Star type does not affect this conclusion, although surface $CO_{\mathrm{2}}$ ice condenses over a larger fraction of the habitable zone around hotter stars. The warm start simulations are thus consistent with 1-D model predictions, suggesting that the classical habitable zone limits in those earlier models are still valid. We also find that the cold start simulations exhibit trends that are consistent with those of previous work for the Sun although we now extend the analysis to other star types.

There have recently been detections of radio emission from low-mass stars, some of which are indicative of star-planet interactions. Motivated by these exciting new results, in this paper we present Alfv\'en wave-driven stellar wind models of the two active planet-hosting M dwarfs Prox Cen and AU Mic. Our models incorporate large-scale photospheric magnetic field maps reconstructed using the Zeeman-Doppler Imaging method. We obtain a mass-loss rate of $0.25~\dot{M}_{\odot}$ for the wind of Prox Cen. For the young dwarf AU Mic, we explore two cases: a low and high mass-loss rate. Depending on the properties of the Alfv\'en waves which heat the corona in our wind models, we obtain mass-loss rates of $27$ and $590~\dot{M}_{\odot}$ for AU Mic. We use our stellar wind models to assess the generation of electron cyclotron maser instability emission in both systems, through a mechanism analogous to the sub-Alfv\'enic Jupiter-Io interaction. For Prox Cen we do not find any feasible scenario where the planet can induce radio emission in the star's corona, as the planet orbits too far from the star in the super-Alfv\'enic regime. However, in the case that AU Mic has a stellar wind mass-loss rate of $27~\dot{M}_{\odot}$, we find that both planets b and c in the system can induce radio emission from $\sim10$ MHz to 3 GHz in the corona of the host star for the majority of their orbits, with peak flux densities of $\sim10$ mJy. Detection of such radio emission would allow us to place an upper limit on the mass-loss rate of the star.

We demonstrate that certain astrophysical distributions can be modelled with the truncated Weibull distribution, which can lead to some insights: in particular, we report the average value, the $r$th moment, the variance, the median, the mode, the generation of random numbers, and the evaluation of the two parameters with maximum likelihood estimators. The first application of the Weibull distribution is to the initial mass function for stars. The magnitude version of the Weibull distribution is applied to the luminosity function for the Sloan Digital Sky Survey (SDSS) galaxies and to the photometric maximum of the 2MASS Redshift Survey (2MRS) galaxies. The truncated Weibull luminosity function allows us to model the average value of the absolute magnitude as a function of the redshift for the 2MRS galaxies.

The current event rate estimates of long gamma-ray bursts based on distinct methods or samples especially at lower redshift are largely debated, which motivates us to re-study the dependence of luminosity function and event rates for different burst samples on the criteria of sample selection and threshold effect in this letter. To ensure the sample completeness as possible, we have chosen two samples including 88 and 118 long bright bursts with known redshift and peak flux over 2.6 ph cm$^{-2}$ s$^{-1}$. It is found that the evolution of luminosity with redshift can be expressed by $L\propto(1+z)^k$ with a diverse $k$ relied more on the sample selection. Interestingly, the cumulative distributions of either non-evolving luminosities or redshifts are found to be also determined by the sample selection rather the instrumental sensitivity. Nevertheless, the non-evolving luminosities of our samples are similarly distributed with a comparable break luminosity of $L_0\sim10^{51}$ erg s$^{-1}$. Importantly, we verify with a K-S test that three cases of event rates for the two burst samples evolve with redshift similarly except a small discrepancy due to sampling differences at low-redshift of $z<1$, in which all event rates show an excess of gaussian profile instead of monotonous decline. Most importantly, it is found that the low-redshift burst event rates violate the star formation rates, while both of them are good in agreement with each other in the higher-redshift regions as many authors discovered previously. Consequently, we predict that two types of long gamma-ray bursts should be expected on the basis of whether they match the star formation or not.

Galaxy lenses are frequently modeled as an elliptical mass distribution with external shear and isothermal spheres to account for secondary and line-of-sight galaxies. There is statistical evidence that some fraction of observed quads are inconsistent with these assumptions, and require a dipole-like contribution to the mass with respect to the light. Simplifying assumptions about the shape of mass distributions can lead to the incorrect recovery of parameters such as $H_0$. We create several tests of synthetic quad populations with different deviations from an elliptical shape, then fit them with an ellipse+shear model, and measure the recovered values of $H_0$. Kinematic constraints are not included. We perform two types of fittings -- one with a single point source and one with an array of sources emulating an extended source. We carry out two model-free comparisons between our mock quads and the observed population. One result of these comparisons is a statistical inconsistency not yet mentioned in the literature: the image distance ratios with respect to the lens center of observed quads appear to span a much wider range than those of synthetic or simulated quads. Bearing this discrepancy in mind, our mock populations can result in biases on $H_0$ $\sim10\%$.

The merger rate of primordial black holes depends on their initial clustering. In the absence of primordial non-Gaussianity correlating short and large-scales, primordial black holes are distributed \`a la Poisson at the time of their formation. However, primordial non-Gaussianity of the local-type may correlate primordial black holes on large-scales. We show that future experiments looking for CMB $\mu$-distortion will test the hypothesis of initial primordial black hole clustering induced by local non-Gaussianity.

The very first detection of gravitational waves from a neutron star binary merger, GW170817, exceeded all expectations. The event was relatively nearby, which may point to a relatively high merger rate. It was possible to extract finite-size effects from the gravitational-wave signal, which constrains the nuclear equation of state. Also, an electromagnetic counterpart was detected at many wavebands from radio to gamma rays marking the begin of a new multi-messenger era involving gravitational waves. We describe how multi-messenger observations of GW170817 are employed to constrain the nuclear equation of state. Combining the information from the optical emission and the mass measurement through gravitational waves leads to a lower limit on neutron star radii. According to this conservative analysis, which employs a minimum set of assumptions, the radii of neutron stars with typical masses should be larger than about 10.7~km. This implies a lower limit on the tidal deformability of about 210, while much stronger lower bounds are not supported by the data of GW170817. The multi-messenger interpretation of GW170817 rules out very soft nuclear matter and complements the upper bounds on NS radii which are derived from the measurement of finite-size effects during the pre-merger phase. We highlight the future potential of multi-messenger observations and of GW measurements of the postmerger phase for constraining the nuclear equation of state. Finally, we propose an observing strategy to maximize the scientific yield of future multi-messenger observations.

Classical T Tauri stars are young low-mass systems still accreting material from their disks. These systems are dynamic on timescales of hours to years. The observed variability can help us infer the physical processes that occur in the circumstellar environment. We aim at understanding the dynamics of the magnetic interaction between the star and the inner accretion disk in young stellar objects. We present the case of the young stellar system V2129 Oph, which is a well-known T Tauri star. We performed a time series analysis of this star using high-resolution spectroscopic data at optical and infrared wavelengths from CFHT/ESPaDOnS, ESO/HARPS and CFHT/SPIRou. The new data sets allowed us to characterize the accretion-ejection structure in this system and to investigate its evolution over a timescale of a decade via comparisons to previous observational data. We measure radial velocity variations and recover a stellar rotation period of 6.53d. However, we do not recover the stellar rotation period in the variability of various circumstellar lines, such as H$\alpha$ and H$\beta$ in the optical or HeI 1083nm and Pa$\beta$ in the infrared. Instead, we show that the optical and infrared line profile variations are consistent with a magnetospheric accretion scenario that shows variability with a period of about 6.0d, shorter than the stellar rotation period. Additionally, we find a period of 8.5d in H$\alpha$ and H$\beta$ lines, probably due to a structure located beyond the corotation radius, at a distance of 0.09au. We investigate whether this could be accounted for by a wind component, twisted or multiple accretion funnel flows, or an external disturbance in the inner disk. We conclude that the dynamics of the accretion-ejection process can vary significantly on a timescale of just a few years, presumably reflecting the evolving magnetic field topology at the stellar surface.

Continuum gamma-ray emission produced by interactions of cosmic rays with interstellar matter and radiation fields is a probe of non-thermal particle populations in galaxies. After decades of continuous improvements in experimental techniques and an ever-increasing sky and energy coverage, gamma-ray observations reveal in unprecedented detail the properties of galactic cosmic rays. A variety of scales and environments are now accessible to us, from the local interstellar medium near the Sun and the vicinity of cosmic-ray accelerators, out to the Milky Way at large and beyond, with a growing number of gamma-ray emitting star-forming galaxies. Gamma-ray observations have been pushing forward our understanding of the life cycle of cosmic rays in galaxies and, combined with advances in related domains, they have been challenging standard assumptions in the field and have spurred new developments in modelling approaches and data analysis methods. We provide a review of the status of the subject and discuss perspectives on future progress.

Periodic quasars have been suggested to host supermassive binary black holes (BBHs) in their centers, and their optical/UV periodicities are interpreted as caused by either the Doppler-boosting (DB) effect of continuum emission from the disk around the secondary black hole (BH) or intrinsic accretion rate variation. However, no other definitive evidence has been found to confirm such a BBH interpretation(s). In this paper, we investigate the responses of broad emission lines (BELs) to the continuum variations for these quasars under two BBH scenarios, and check whether they can be distinguished from each other and from that of a single BH system. We assume a simple circumbinary broad-line region (BLR) model, compatible with BLR size estimates, with a standard $\Gamma$ distribution of BLR clouds. We find that BELs may change significantly and periodically under the BBH scenarios due to (1) the position variation of the secondary BH and (2) the DB effect, if significant, and/or intrinsic variation, which is significantly different from the case of a single BH system. For the two BBH scenarios, the responses of BELs to (apparent) continuum variations, caused by the DB effect or intrinsic rate variation, are also significantly different from each other, mainly because the DB effect has a preferred direction along the direction of motion of the secondary BH, while that due to intrinsic variation does not. Such differences in the responses of BELs from different scenarios may offer a robust way to distinguish different interpretations of periodic quasars and to identify BBHs, if any, in these systems.

We present ALMA observations of 101 protoplanetary disks within the star-forming region Lynds 1641 in the Orion Molecular Cloud A. Our observations include 1.33 mm continuum emission and spectral windows covering the J=2-1 transition of $^{12}$CO, $^{13}$CO, and C$^{18}$O. We detect 89 protoplanetary disks in the dust continuum at the 4$\sigma$ level ($\sim$88% detection rate) and 31 in $^{12}$CO, 13 in $^{13}$CO, and 4 in C$^{18}$O. Our sample contains 23 transitional disks, 20 of which are detected in the continuum. We target infrared-bright Class II objects, which biases our sample towards massive disks. We determine dust masses or upper limits for all sources in our sample and compare our sample to protostars in this region. We find a decrease in dust mass with evolutionary state. We also compare this sample to other regions surveyed in the (sub-)millimeter and find that Lynds 1641 has a relatively massive dust disk population compared to regions of similar and older ages, with a median dust mass of 11.1$^{+32.9}_{-4.6}$ $M_\oplus$ and 27% with dust masses equal to or greater than the minimum solar nebula dust mass value of $\sim$30 $M_\oplus$. We analyze the disk mass-accretion rate relationship in this sample and find that the viscous disk lifetimes are similar to the age of the region, however with a large spread. One object, [MGM2012] 512, shows large-scale ($>$5000 AU) structure in both the dust continuum and the three gas lines. We discuss potential origins for this emission, including an accretion streamer with large dust grains.

The present paper gives new elements for the light variations of the EA variable star AY Peg, on the basis of new times of minimum performed visually and with ccd by members of GEOS between 1985 and 2018, and the ASAS-SN set of data available. On one hand, we can establish a new ephemeris with a possible quadratic term, and on the other hand, the amplitude of the primary minimum appears much deeper than the one given in GCVS. AY Peg varies between 13.1 and 15.6 magnitude at its primary eclipse.

Radio and $\gamma$-ray loud narrow-line Seyfert 1 galaxies (NLS1s) are unique objects to study the formation and evolution of relativistic jets, as they are believed to have high accretion rates and powered by low mass black holes contrary to that known for blazars. However, only about a dozen $\gamma$-ray detected NLS1s ($\gamma$-NLS1s) are known to date and all of them are at $z\le1$. Here, we report the identification of a new $\gamma$-ray emitting NLS1 TXS 1206+549 at $z=1.344$. A near-infrared spectrum taken with the Subaru telescope showed H$\beta$ emission line with FWHM of $1194\pm77$ km s$^{-1}$ and weak [O III] emission line but no optical Fe II lines, due to the limited wavelength coverage and poor signal-to-noise ratio. However, UV Fe II lines are present in the SDSS optical spectrum. The source is very radio-loud, unresolved, and has a flat radio spectrum. The broadband SED of the source has the typical two hump structure shown by blazars and other $\gamma$-NLS1s. The source exhibits strong variability at all wavelengths such as the optical, infrared, and $\gamma$-ray bands. All these observed characteristics show that TXS 1206+549 is the most distant $\gamma$-NLS1 known to date.

We determine the starspot detection rate in exoplanetary transit light curves for M and K dwarf stars observed by the Transiting Exoplanet Survey Satellite (TESS) using various starspot filling factors and starspot distributions. We used $3.6\times10^9$ simulations of planetary transits around spotted stars using the transit-starspot model \prism. The simulations cover a range of starspot filling factors using one of three distributions: uniform, polar-biased, and mid-latitude. After construction of the stellar disc and starspots, we checked the transit cord for starspots and examined the change in flux of each starspot to determine whether or not a starspot anomaly would be detected. The results were then compared to predicted planetary detections for TESS. The results show that for the case of a uniform starspot distribution, $64\pm9$ M dwarf and $23\pm4$ K dwarf transit light curves observed by TESS will contain a starspot anomaly. This reduces to $37\pm6$ M dwarf and $12\pm2$ K dwarf light curves for a polar-biased distribution and $47\pm7$ M dwarf and $21\pm4$ K dwarf light curves for a mid-latitude distribution. Currently there are only 17 M dwarf and 10 K dwarf confirmed planetary systems from TESS, none of which are confirmed as showing starspot anomalies. All three starspot distributions can explain the current trend. However, with such a small sample, a firm conclusion can not be made at present. In the coming years when more TESS M and K dwarf exoplanetary systems have been detected and characterised, it will be possible to determine the dominant starspot distribution.

In linear stability analysis of field quantities described by partial differential equations, the well-established classical theory is all but impossible to apply to concrete problems in its entirety even for uniform backgrounds when the spatial dimension is more than 1. In this study, using the Lefschetz thimble method, we develop a new formalism to give an explicit expression to the asymptotic behavior of linear perturbations. It is not only more mathematically rigorous than the previous theory but also useful practically in its applications to realistic problems, and, as such, has an impact on broad subjects in physics.

We explore color-kinematic duality for tree-level AdS/CFT correlators in momentum space. We start by studying the bi-adjoint scalar in AdS at tree-level as an illustrative example. We follow this by investigating two forms of color-kinematic duality in Yang-Mills theory, the first for the integrated correlator in AdS$_4$ and the second for the integrand in general AdS$_{d+1}$. For the integrated correlator, we find color-kinematics does not yield additional relations among $n$-point, color-ordered correlators. To study color-kinematics for the AdS$_{d+1}$ Yang-Mills integrand, we use a spectral representation of the bulk-to-bulk propagator so that AdS diagrams are similar in structure to their flat space counterparts. Finally, we study color KLT relations for the integrated correlator and double-copy relations for the AdS integrand. We find that double-copy in AdS naturally relates the bi-adjoint theory in AdS$_{d+3}$ to Yang-Mills in AdS$_{d+1}$. We also find a double-copy relation at three-points between Yang-Mills in AdS$_{d+1}$ and gravity in AdS$_{d-1}$ and comment on the higher-point generalization. By analytic continuation, these results on AdS/CFT correlators can be translated into statements about the wave function of the universe in de Sitter.

We explore the prospects for direct detection of dark energy by current and upcoming terrestrial dark matter direct detection experiments. If dark energy is driven by a new light degree of freedom coupled to matter and photons then dark energy quanta are predicted to be produced in the Sun. These quanta free-stream towards Earth where they can interact with Standard Model particles in the detection chambers of direct detection experiments, presenting the possibility that these experiments could be used to test dark energy. Screening mechanisms, which suppress fifth forces associated with new light particles, and are a necessary feature of many dark energy models, prevent production processes from occurring in the core of the Sun, and similarly, in the cores of red giant, horizontal branch, and white dwarf stars. Instead, the coupling of dark energy to photons leads to production in the strong magnetic field of the solar tachocline via a mechanism analogous to the Primakoff process. This then allows for detectable signals on Earth while evading the strong constraints that would typically result from stellar probes of new light particles. As an example, we examine whether the electron recoil excess recently reported by the XENON1T collaboration can be explained by chameleon-screened dark energy, and find that such a model is preferred over the background-only hypothesis at the $2.0\sigma$ level, in a large range of parameter space not excluded by stellar (or other) probes. This raises the tantalizing possibility that XENON1T may have achieved the first direct detection of dark energy. Finally, we study the prospects for confirming this scenario using planned future detectors such as XENONnT, PandaX-4T, and LUX-ZEPLIN.

Large-scale bulk flows are commonplace in the universe. This means that observers living in typical galaxies, like our Milky Way, do not follow the mean universal expansion, but have peculiar velocities relative to it. Using relativistic linear cosmological perturbation theory, we show that bulk peculiar motions introduce a characteristic length scale. The latter is analogous to the familiar Jeans length and marks the threshold below which the linear kinematics, as seen by observers living inside the bulk flow, are dominated by relative-motion effects. On these scales cosmological measurements can vary considerably between the rest-frame of the bulk flow and that of the smooth Hubble expansion, due to relative-motion effects alone. In this work, we look into and compare the measurements of the deceleration parameter. We find that the associated "peculiar Jeans length" has typical size of the order of few hundred Mpc. On smaller scales, the deceleration parameter measured by the bulk-flow observers can be considerably larger (or lower) than its Hubble-frame counterpart. This, in turn, depends on whether the peculiar motion is slightly expanding (or contracting) relative to the background universal expansion. Then, assuming that expanding and contracting bulk flows are randomly distributed, nearly half of the observers in the universe may be misled to think that their cosmos is over-decelerated. The rest of them, on the other hand, may come to believe that the universe is under-decelerated, or even accelerated in some cases.

Quantum mechanics allows for states in macroscopic superpositions, but they ordinarily undergo rapid decoherence due to interactions with their environment. A system that only interacts gravitationally, such as an arrangement of dark matter (DM), may exhibit slow decoherence. In this letter, we compute the decoherence rate of a quantum object in general relativity. For axion DM in a superposition of the field's phase, we find that DM in the Milky Way is robust against decoherence, while a spatial superposition is not. This novel phase behavior may be relevant to direct detection experiments.

Besides the relation between the wave vector $\bm k$ and the complex frequency $\omega$, wave polarization is useful for characterizing the properties of a plasma wave. The polarization of the electromagnetic fields, $\delta \bm E$ and $\delta \bm B$, have been widely used in plasma physics research. Here, we derive equations for the density and velocity perturbations, $\delta n_s$ and $\delta{\bm v}_s$, respectively, of each species in the electromagnetic kinetic plasma dispersion relation by using their relation to the species current density perturbation $\delta {\bm J}_s$. Then we compare results with those of another commonly used plasma dispersion code (WHAMP) and with those of a multi-fluid plasma dispersion relation. We also summarize a number of useful polarization quantities, such as magnetic ellipticity, orientation of the major axis of the magnetic ellipse, various ratios of field energies and kinetic energies, species compressibility, parallel phase ratio, Alfv\'en-ratio, etc., which are useful for plasma physics research, especially for space plasma studies. This work represents an extension of the BO electromagnetic dispersion code [H.S. Xie, Comput. Phys. Comm. 244 (2019) 343-371] to enhance its calculation of polarization and to include the capability of solving the electromagnetic magnetized multi-fluid plasma dispersion relation.

Since the initial detection of Gravitational Waves in 2015, 50 candidate events have been reported by the LIGO-Virgo-KAGRA collaboration. As the current generation of detectors move towards their design sensitivity the rate of these detections will increase. The next generation of detectors are likely to have high enough sensitivities that multiple merging binaries will be visible at the same time. In this paper we show that this is likely to happen before the end of the decade, with the move to the LIGO-Voyager detector. We investigate the situation of overlapping Binary-Black-Hole mergers in these detectors. We find that current parameter estimation techniques are capable of distinguishing the louder of two merging BBH events, without significant bias, when their merger times are not less than $\sim0.1$ seconds apart and when the ratio of the signal-to-noise ratios of the systems is uneven. This region of overlapping parameter space is dependent upon the sky locations of the signals and the relation of those locations to the light travel time between detectors. We also find that, if two signals are highly overlapping, then the recovered set of parameters often show strong evidence of precession. Finally we show that bias can occur even when the signal causing the bias is below the detection threshold.

Partially recorded data are frequently encountered in many applications. In practice, such datasets are usually clustered by removing incomplete cases or features with missing values, or by imputing missing values, followed by application of a clustering algorithm to the resulting altered data set. Here, we develop clustering methodology through a model-based approach using the marginal density for the observed values, using a finite mixture model of multivariate $t$ distributions. We compare our algorithm to the corresponding full expectation-maximization (EM) approach that considers the missing values in the incomplete data set and makes a missing at random (MAR) assumption, as well as case deletion and imputation. Since only the observed values are utilized, our approach is computationally more efficient than imputation or full EM. Simulation studies demonstrate that our approach has favorable recovery of the true cluster partition compared to case deletion and imputation under various missingness mechanisms, and is more robust to extreme MAR violations than the full EM approach since it does not use the observed values to inform those that are missing. Our methodology is demonstrated on a problem of clustering gamma-ray bursts and is implemented in the https://github.com/emilygoren/MixtClust R package.