Papers for Tuesday, Feb 02 2021

We show that the perturbation at the peak of the light curve of microlensing event KMT-2019-BLG-0371 is explained by a model with a mass ratio between the host star and planet of $q \sim 0.08$. Due to the short event duration ($t_{\rm E} \sim 6.5\ $ days), the secondary object in this system could potentially be a massive giant planet. A Bayesian analysis shows that the system most likely consists of a host star with a mass $M_{\rm h} = 0.09^{+0.14}_{-0.05}M_{\odot}$ and a massive giant planet with a mass $M_{\rm p} = 7.70^{+11.34}_{-3.90}M_{\rm Jup}$. However, the interpretation of the secondary as a planet (i.e., as having $M_{\rm p} < 13 M_{\rm Jup}$) rests entirely on the Bayesian analysis. Motivated by this event, we conduct an investigation to determine which constraints meaningfully affect Bayesian analyses for microlensing events. We find that the masses inferred from such a Bayesian analysis are determined almost entirely by the measured value of $\theta_{\rm E}$ and are relatively insensitive to other factors such as the direction of the event $(\ell, b)$, the lens-source relative proper motion $\mu_{\rm rel}$, or the specific Galactic model prior.

We compare the star forming main sequence (SFMS) -- both integrated and resolved on 1kpc scales -- between the high-resolution TNG50 simulation of IllustrisTNG and observations from the 3D-HST slitless spectroscopic survey at z~1. Contrasting integrated star formation rates (SFRs), we find that the slope and normalization of the star-forming main sequence in TNG50 are quantitatively consistent with values derived by fitting observations from 3D-HST with the Prospector Bayesian inference framework. The previous offsets of 0.2-1dex between observed and simulated main sequence normalizations are resolved when using the updated masses and SFRs from Prospector. The scatter is generically smaller in TNG50 than in 3D-HST for more massive galaxies with M_*>10^10Msun, even after accounting for observational uncertainties. When comparing resolved star formation, we also find good agreement between TNG50 and 3D-HST: average specific star formation rate (sSFR) radial profiles of galaxies at all masses and radii below, on, and above the SFMS are similar in both normalization and shape. Most noteworthy, massive galaxies with M_*>10^10.5Msun, which have fallen below the SFMS due to ongoing quenching, exhibit a clear central SFR suppression, in both TNG50 and 3D-HST. In TNG this inside-out quenching is due to the supermassive black hole (SMBH) feedback model operating at low accretion rates. In contrast, the original Illustris simulation, without this same physical SMBH mechanism, does not reproduce the central SFR profile suppression seen in data. The observed sSFR profiles provide support for the TNG quenching mechanism and how it affects gas on kiloparsec scales in the centers of galaxies.

Indirect detection experiments typically measure the flux of annihilating dark matter (DM) particles propagating freely through galactic halos. We consider a new scenario where celestial bodies "focus" DM annihilation events, increasing the efficiency of halo annihilation. In this setup, DM is first captured by celestial bodies, such as neutron stars or brown dwarfs, and then annihilates within them. If DM annihilates to sufficiently long-lived particles, they can escape and subsequently decay into detectable radiation. This produces a distinctive annihilation morphology, which scales as the product of the DM and celestial body densities, rather than as DM density squared. We show that this signal can dominate over the halo annihilation rate in $\gamma$-ray observations in both the Milky Way Galactic center and globular clusters. We use \textit{Fermi} and H.E.S.S. data to constrain the DM-nucleon scattering cross section, setting powerful new limits down to $\sim10^{-39}~$cm$^2$ for sub-GeV DM using brown dwarfs, which is up to nine orders of magnitude stronger than existing limits. We demonstrate that neutron stars can set limits for TeV-scale DM down to about $10^{-47}~$cm$^2$.

We study the multimessenger signals from the merger of a black hole with a magnetized neutron star using resistive magneto-hydrodynamics simulations coupled to full general relativity. We focus on a case with a 5:1 mass-ratio, where only a small amount of the neutron star matter remains post-merger, but we nevertheless find that significant electromagnetic radiation can be powered by the interaction of the neutron star's magnetosphere with the black hole. In the lead-up to merger, strong twisting of magnetic field lines from the inspiral leads to plasmoid emission and results in a luminosity in excess of that expected from unipolar induction. We find that the strongest emission occurs shortly after merger during a transitory period in which magnetic loops form and escape the central region. The remaining magnetic field collimates around the spin axis of the remnant black hole before dissipating, an indication that, in more favorable scenarios (higher black hole spin/lower mass ratio) with larger accretion disks, a jet would form.

We study the abundance of satellite galaxies around 198 Milky Way- (MW) and M31-like hosts in TNG50, the final instalment in the IllustrisTNG suite of cosmological magnetohydrodynamical simulations. MW/M31-like analogues are defined as disky galaxies with stellar masses of Mstar = 10^10.5-11.2 Msun in relative isolation at z = 0. By defining satellites as galaxies with Mstar > 5*10^6 Msun within 300 kpc (3D) of their host, we find a remarkable level of diversity and host-to-host scatter across individual host galaxies. The median (16th - 84th percentiles) TNG50 MW/M31-like galaxy hosts a total of 5 (2-11) satellites with Mstar > 5*10^6 Msun, reaching up to Mstar ~ 10^8.5 Msun (10^7.4-9.4 Msun). The abundance of subhaloes with Mdyn > 5*10^7 Msun is larger by a factor of more than 10. The number of all satellites (subhaloes) ever accreted is larger by a factor of 4-5 (3-5) than those surviving to z = 0. Hosts with larger galaxy stellar mass, brighter K-band luminosity, larger total halo mass, and more recent halo assembly typically have a larger number of surviving satellites. The satellite abundances around TNG50 MW/M31-like galaxies are consistent with similar hosts from observational surveys (e.g. SAGA) and previous simulations (e.g. Latte). While the observed MW satellite system falls within the TNG50 scatter across all stellar masses considered, M31 is slightly more satellite-rich than our 1 sigma scatter, possibly due to volume and mass limitations. We find a handful of systems with both a Large and a Small Magellanic Cloud-like satellite. There is no missing satellites problem with TNG50.

In the two-phase scenario of galaxy formation, a galaxy's stellar mass growth is first dominated by in-situ star formation, and subsequently by accretion. We analyse the radial distribution of the accreted stellar mass in ~500 galaxies from the hydrodynamical cosmological simulation Magneticum. Generally, we find good agreement with other simulations in that higher mass galaxies have larger accreted fractions, but we predict higher accretion fractions for low-mass galaxies. Based on the radial distribution of the accreted and in-situ components, we define 6 galaxy classes, from completely accretion dominated to completely in-situ dominated, and measure the transition radii between in-situ and accretion-dominated regions for galaxies that have such a transition. About 70% of our galaxies have one transition radius. However, we also find about 10% of the galaxies to be accretion dominated everywhere, and about 13% to have two transition radii, with the centre and the outskirts both being accretion dominated. We show that these classes are strongly correlated with the galaxy merger histories, especially with the mergers' cold gas fractions. We find high total in-situ (low accretion) fractions to be associated with smaller, lower mass galaxies, lower central dark matter fractions, and larger transition radii. Finally, we show that the dips in observed surface brightness profiles seen in many early-type galaxies do not correspond to the transition from in-situ to accretion-dominated regions, and any inferred mass fractions are not indicative of the true accreted mass. Instead, these dips contain information about the galaxies' dry minor merger assembly history.

We present predictions for the extent of the dust-continuum emission of thousands of main-sequence galaxies drawn from the TNG50 simulation between $z=1-5$. To this aim, we couple the radiative transfer code SKIRT to the output of the TNG50 simulation and measure the dust-continuum half-light radius of the modeled galaxies, assuming a Milky Way dust type and a metallicity dependent dust-to-metal ratio. The dust-continuum half-light radius at observed-frame 850 $\mu$m is up to $\sim$75 per cent larger than the stellar half-mass radius, but significantly more compact than the observed-frame 1.6 $\mu$m (roughly corresponding to H-band) half-light radius, particularly towards high redshifts: the compactness compared to the 1.6 $\mu$m emission increases with redshift. This is driven by obscuration of stellar light from the galaxy centres, which increases the apparent extent of 1.6 $\mu$m disk sizes relative to that at 850 $\mu$m. The difference in relative extents increases with redshift because the observed-frame 1.6 $\mu$m emission stems from ever shorter wavelength stellar emission. These results suggest that the compact dust-continuum emission observed in $z>1$ galaxies is not (necessarily) evidence of the buildup of a dense central stellar component. We also find that the dust-continuum half-light radius very closely follows the radius containing half the star formation in galaxies, indicating that single band dust-continuum emission is a good tracer of the location of (obscured) star formation. The dust-continuum emission is more compact than the H2 mass (for galaxies at $z\geq 2$) and the underlying dust mass. The dust emission strongly correlates with locations with the highest dust temperatures, which do not need to be the locations where most H$_2$ and/or dust is located. The presented results are a common feature of main-sequence galaxies.

We explore in general relativity the survival time of neutron stars that host an endoparasitic, possibly primordial, black hole at their center. Corresponding to the minimum steady-state Bondi accretion rate for adiabatic flow that we found earlier for stiff nuclear equations of state (EOSs), we derive analytically the maximum survival time after which the entire star will be consumed by the black hole. We also show that this maximum survival time depends only weakly on the stiffness for polytropic EOSs with $\Gamma \geq 5/3$, so that this survival time assumes a nearly universal value that depends on the initial black hole mass alone. Establishing such a value is important for constraining the contribution of primordial black holes in the mass range $10^{-16} M_\odot \lesssim M \lesssim 10^{-10} M_\odot$ to the dark-matter content of the Universe.

We discuss the possible interpretation of the recently observed transient, GN-z11-flash as originating from a shock-breakout in a Population III supernova occurring in the GN-z11 galaxy at $z \sim 11$. We find that the parameters of the explosion are fully consistent with those expected from the shock breakout associated with a Type II supernova of a progenitor star of $\sim 300$ solar masses in this galaxy, with of order unity such events expected over an observing timescale of a few years. We forecast the expected number of such transients from $z > 10$ galaxies as a function of their host stellar mass and star formation rate.

We have entered a new era where integral-field spectroscopic surveys of galaxies are sufficiently large to adequately sample large-scale structure over a cosmologically significant volume. This was the primary design goal of the SAMI Galaxy Survey. Here, in Data Release 3 (DR3), we release data for the full sample of 3068 unique galaxies observed. This includes the SAMI cluster sample of 888 unique galaxies for the first time. For each galaxy, there are two primary spectral cubes covering the blue (370-570nm) and red (630-740nm) optical wavelength ranges at spectral resolving power of R=1808 and 4304 respectively. For each primary cube, we also provide three spatially binned spectral cubes and a set of standardized aperture spectra. For each galaxy, we include complete 2D maps from parameterized fitting to the emission-line and absorption-line spectral data. These maps provide information on the gas ionization and kinematics, stellar kinematics and populations, and more. All data are available online through Australian Astronomical Optics (AAO) Data Central.

We present the cumulative star-formation histories (SFHs) of >15000 dwarf galaxies ($M_{*}=10^{7-10}M_{\odot}$) from the TNG50 run of the IllustrisTNG suite across a vast range of environments. The key factors determining the dwarfs' SFHs are their status as central or satellite and their stellar mass, with centrals and more massive dwarfs assembling their stellar mass at later times on average compared to satellites and lower mass dwarfs. The satellites (in hosts of total mass $M_{200c,\,host}=10^{12-14.3}M_{\odot}$) assembled 90% of their z=0 stellar mass ~$7.0_{-5.5}^{+3.3}$ Gyr ago, while the centrals did so only ~$1.0_{-0.5}^{+4.0}$ Gyr ago. TNG50 predicts a large diversity in SFHs for both centrals and satellites, so that the stacked cumulative SFHs are representative of the TNG50 dwarf populations only in an average sense and individual dwarfs can have significantly different cumulative SFHs. Satellite dwarfs with the highest stellar mass to host mass ratios have the latest stellar mass assembly. Satellites at fixed stellar and host halo mass, found closer to the cluster centre, or accreted at earlier times, show significantly earlier stellar mass assembly. These trends, as well as the shapes of the SFHs themselves, are a manifestation of the varying proportions within a given subsample of quenched vs. star-forming galaxies, which exhibit markedly distinct SFH shapes. We also find a subtle effect whereby satellite dwarfs in the most massive hosts at z=0 have higher SFRs at early times, well before final infall into their z=0 host, compared to a control sample of centrals mass-matched at the time of accretion. This suggests that the large-scale environment can have a mild effect even on future satellites by providing the conditions for enhanced SF at early epochs. Our results are useful theoretical predictions for comparison to future resolved-stellar-population observations.

Many non-minimal dark matter scenarios lead to oscillatory features in the matter power spectrum induced by interactions either within the dark sector or with particles from the standard model. Observing such dark acoustic oscillations would therefore be a major step towards understanding dark matter. We investigate what happens to oscillatory features during the process of nonlinear structure formation. We show that at the level of the power spectrum, oscillations are smoothed out by nonlinear mode coupling, gradually disappearing towards lower redshifts. In the halo mass function, however, the same oscillations remain visible until the present epoch. As a consequence, dark acoustic oscillations could be detectable in observations that are either based on the halo mass function or on the high-redshift power spectrum. We investigate the effect of such oscillations on different observables, namely, the cluster mass function, the stellar-to-halo mass relation, and the Lyman-$\alpha$ flux power spectrum. We find that dark acoustic oscillations remain visible in all of these observables, but they are very extended and of low amplitude, making it challenging to detect them as distinct features in the data.

We present GLEAM (Galaxy Line Emission & Absorption Modeling), a Python tool for fitting Gaussian models to emission and absorption lines in large samples of 1D extragalactic spectra. GLEAM is tailored to work well in batch mode without much human interaction. With GLEAM, users can uniformly process a variety of spectra, including galaxies and active galactic nuclei, in a wide range of instrument setups and signal-to-noise regimes. GLEAM also takes advantage of multiprocessing capabilities to process spectra in parallel. With the goal of enabling reproducible workflows for its users, GLEAM employs a small number of input files, including a central, user-friendly configuration in which fitting constraints can be defined for groups of spectra and overrides can be specified for edge cases. For each spectrum, GLEAM produces a table containing measurements and error bars for the detected spectral lines and continuum, and upper limits for non-detections. For visual inspection and publishing, GLEAM can also produce plots of the data with fitted lines overlaid. In the present paper, we describe GLEAM's main features, the necessary inputs, expected outputs, and some example applications, including thorough tests on a large sample of optical/infra-red multi-object spectroscopic observations and integral field spectroscopic data. gleam is developed as an open-source project hosted at https://github.com/multiwavelength/gleam and welcomes community contributions.

The unprecedented sky coverage and observing cadence of the All-Sky Automated Survey for SuperNovae (ASAS-SN) has resulted in the discovery and continued monitoring of a large sample of Galactic transients. The vast majority of these are accretion-powered dwarf nova outbursts in cataclysmic variable systems, but a small subset are thermonuclear-powered classical novae. Despite improved monitoring of the Galaxy for novae from ASAS-SN and other surveys, the observed Galactic nova rate is still lower than predictions. One way classical novae could be missed is if they are confused with the much larger population of dwarf novae. Here, we examine the properties of 1617 dwarf nova outbursts detected by ASAS-SN and compare them to classical novae. We find that the mean classical nova brightens by ~11 magnitudes during outburst, while the mean dwarf nova brightens by only ~5 magnitudes, with the outburst amplitude distributions overlapping by roughly 15%. For the first time, we show that the amplitude of an outburst and the time it takes to decline by two magnitudes from maximum are positively correlated for dwarf nova outbursts. For classical novae, we find that these quantities are negatively correlated, but only weakly, compared to the strong anti-correlation of these quantities found in some previous work. We show that, even if located at large distances, only a small number of putative dwarf novae could be mis-classified classical novae suggesting that there is minimal confusion between these populations. Future spectroscopic follow-up of these candidates can show whether any are indeed classical novae.

We investigate the specific angular momentum (sAM) $ j(<r)$ profiles of intermediate redshift ($0.4<z<1.4$) star-forming galaxies (SFGs) in the relatively unexplored regime of low masses (down to $M_\star\sim 10^8$M$_{\odot}$) and small sizes (down to $R_{\rm e}\sim 1.5$ kpc) and characterize the sAM scaling relation and its redshift evolution. We have developed a 3D methodology to constrain sAM profiles of the star-forming gas using a forward modeling approach with \galpak{} that incorporates the effects of beam smearing, yielding the intrinsic morpho-kinematic properties even with limited spatial resolution data. Using mock observations from the TNG50 simulation, we find that our 3D methodology robustly recovers the SFR-weighted $j(<r)$ profiles down to low effective signal-to-noise ratio (SNR) of $\gtrapprox3$. We apply our methodology blindly to a sample of 494 \OII{}-selected SFGs in the MUSE Ultra Deep Field (UDF) 9~arcmin$^2$ mosaic data, covering the unexplored $8<\log M_*/$M$_{\odot}<9$ mass range. We find that the (SFR-weighted) sAM relation follows $j\propto M_\star^{\alpha}$ with an index $\alpha$ varying from $\alpha=0.3$ to $\alpha=0.5$, from $\log M_\star/$M$_{\odot}=8$ to $\log M_*/$M$_{\odot}=10.5$. The UDF sample supports a redshift evolution consistent with the $(1+z)^{-0.5}$ expectation from a Universe in expansion. The scatter of the sAM sequence is a strong function of the dynamical state with $\log j|_{M_*}\propto 0.65 \times \log(V_{\rm max}/\sigma)$ where $\sigma$ is the velocity dispersion at $2 R_{\rm e}$. In TNG50, SFGs also form a $j-M_{\star}-(V/\sigma)$ plane but correlates more with galaxy size than with morphological parameters. Our results suggest that SFGs might experience a dynamical transformation before their morphological transformation to becoming passive via either merging or secular evolution.

We derive the mass-radius relation and mass function of molecular clumps in the Large Magellanic Cloud (LMC) and interpret them in terms of the simple feedback model proposed by Fall, Krumholz, and Matzner (FKM). Our work utilizes the dendrogram-based catalog of clumps compiled by Wong et al. from $^{12}$CO and $^{13}$CO maps of six giant molecular clouds in the LMC observed with the Atacama Large Millimeter Array (ALMA). The Magellanic Clouds are the only external galaxies for which this type of analysis is possible at the necessary spatial resolution ($\sim1$ pc). We find that the mass-radius relation and mass function of LMC clumps have power-law forms, $R \propto M^{\alpha}$ and $dN/dM \propto M^{\beta}$, with indices $\alpha = 0.36 \pm 0.03$ and $\beta= -1.8 \pm 0.1 $ over the mass ranges $10^2 M_\odot \lesssim M \lesssim 10^5 M_\odot$ and $10^2 M_\odot \lesssim M \lesssim 10^4 M_\odot$, respectively. With these values of $\alpha$ and $\beta$ for the clumps (i.e., protoclusters), the predicted index for the mass function of young LMC clusters from the FKM model is $\beta \approx 1.7$, in good agreement with the observed index. The situation portrayed here for clumps and clusters in the LMC replicates that in the Milky Way.

We use a simulation-based modelling approach to analyse the anisotropic clustering of the BOSS LOWZ sample over the radial range $0.4 \, h^{-1} \, \mathrm{Mpc}$ to $63 \, h^{-1} \, \mathrm{Mpc}$, significantly extending what is possible with a purely analytic modelling framework. Our full-scale analysis yields constraints on the growth of structure that are a factor of two more stringent than any other study on large scales at similar redshifts. We infer $f \sigma_8 = 0.471 \pm 0.024$ at $z \approx 0.25$, and $f \sigma_8 = 0.431 \pm 0.025$ at $z \approx 0.40$; the corresponding $\Lambda$CDM predictions of the Planck CMB analysis are $0.470 \pm 0.006$ and $0.476 \pm 0.005$, respectively. Our results are thus consistent with Planck, but also follow the trend seen in previous low-redshift measurements of $f \sigma_8$ falling slightly below the $\Lambda$CDM+CMB prediction. We find that small and large radial scales yield mutually consistent values of $f \sigma_8$, but there are $1-2.5 \sigma$ hints of small scales ($< 10 \, h^{-1} \, \mathrm{Mpc}$) preferring lower values for $f \sigma_8$ relative to larger scales. We analyse the constraining power of the full range of radial scales, finding that most of the multipole information about $f\sigma_8$ is contained in the scales $2 \, h^{-1} \, \mathrm{Mpc} \lesssim s \lesssim 20 \, h^{-1} \, \mathrm{Mpc}$. Evidently, once the cosmological information of the quasi-to-nonlinear regime has been harvested, large-scale modes contain only modest additional information about structure growth. Finally, we compare predictions for the galaxy-galaxy lensing amplitude of the two samples against measurements from SDSS and assess the lensing-is-low effect in light of our findings.

Sub-solar mass black hole binaries, due to their light mass, would have to be primordial in origin instead of the result of stellar evolution. Soon after formation in the early Universe, primordial black holes can form binaries after decoupling from the cosmic expansion. Alternatively, primordial black holes as dark matter could also form binaries in the late Universe due to dynamical encounters and gravitational-wave braking. A significant feature for this channel is the possibility that some sources retain non-zero eccentricity in the LIGO/Virgo band. Assuming all dark matter is primordial black holes with a delta function mass distribution, $1M_\odot-1M_\odot$ binaries formed in this late Universe channel can be detected by Advanced LIGO and Virgo with at their design sensitivities with a rate of $\mathcal{O}(1)$/year, where $12\%(3\%)$ events have eccentricity at gravitational-wave frequency 10 Hz, $e^\mathrm{10Hz}\geq0.01(0.1)$, and non-detection can constrain the binary formation rate within this model. Third generation detectors would be expected to detect sub-solar mass eccentric binaries as light as $0.01 M_\odot$ within this channel, if they account for the majority of the dark matter. Furthermore, we use simulated gravitational-wave data to study the ability to search for eccentric gravitational-wave signals using quasi-circular waveform template bank with Advanced LIGO design sensitivity. Assuming binaries with a delta function mass of $0.1(1)M_\odot$ and the eccentricity distribution derived from this late Universe formation channel, for a match-filtering targeted search, $41\%(6\%)$ of the signals would be missed compared to ideal detection rate due to the mismatch in the gravitational-wave signal from eccentricity.

I searched for the ground state 6.8 and 9.2 GHz hyperfine transitions of rubidium and cesium toward M- and L-dwarfs that show Rb and Cs optical resonance lines. The optical lines can pump the hyperfine transitions, potentially forming masers. These spin-flip transitions of Rb and Cs are the principal transitions used in atomic clocks (the $^{133}$Cs hyperfine transition defines the second). If they are detected in stellar atmospheres, these transitions would provide exceptionally precise clocks that can be used as accelerometers, as exoplanet detectors, as probes of the predictions of general relativity, as probes of light propagation effects, and as a means to do fundamental physics with telescopes. Observations of 21 M- and L-dwarfs, however, show no evidence for Rb or Cs maser action, and a previous survey of giant stars made no Rb maser detections.

Despite the importance of Type Ia supernovae (SNe Ia) throughout astronomy, the precise progenitor systems and explosion mechanisms that drive SNe Ia are still unknown. An explosion scenario that has gained traction recently is the double detonation in which an accreted shell of helium detonates and triggers a secondary detonation in the underlying white dwarf. Our research presents a number of high resolution, multi-dimensional, full star simulations of thin-helium-shell, sub-Chandrasekhar-mass white dwarf progenitors that undergo a double detonation. We confirm the viability of the double detonation across a range of helium shell parameter space as well as present bulk yields and ejecta profiles for each progenitor. The yields obtained are generally consistent with previous works and indicate the likelihood of producing observables that resemble SNe Ia. The dimensionality of our simulations allow us to examine features of the double detonation more closely, including the details of the off-center secondary ignition and asymmetric ejecta. We find considerable differences in the high-velocity extent of post-detonation products across different lines of sight. The data from this work will be used to generate predicted observables and may further support the viability of the double detonation scenario as a SNe Ia channel as well as show how properties of the progenitor or viewing angle may influence trends in observable characteristics.

Based on the photometric redshift catalog of Zou H. et al. (2019), we apply a fast clustering algorithm to identify 540,432 galaxy clusters at $z\lesssim1$ in the DESI legacy imaging surveys, which cover a sky area of about 20,000 deg$^2$. Monte-Carlo simulations indicate that the false detection rate of our detecting method is about 3.1\%. The total masses of galaxy clusters are derived using a calibrated richness--mass relation that are based on the observations of X-ray emission and Sunyaev \& Zel'dovich effect. The median redshift and mass of our detected clusters are about 0.53 and $1.23\times10^{14} M_\odot$, respectively. Comparing with previous clusters identified using the data of the Sloan Digital Sky Survey (SDSS), we can recognize most of them, especially those with high richness. Our catalog will be used for further statistical studies on galaxy clusters and environmental effects on the galaxy evolution, etc.

Recent Juno observations have greatly extended the temporal and spatial coverage of lightning detection on Jupiter. We use these data to constrain a model of moist convection and lightning generation in Jupiter's atmosphere, and derive a roughly solar abundance of water at the base of the water cloud. Shallow lightning, observed by Juno (Becker et al., 2020, Nature, 584, 55-58) and defined as flashes originating at altitudes corresponding to pressure less than 2 bars, is reproduced, as is lightning at a deeper range of pressures, including those below the water cloud base. It is found that the generation of lightning requires ammonia to stabilize liquid water at altitudes corresponding to sub-freezing temperatures. We find a range of local water abundances in which lightning is possible, including subsolar values of water--consistent with other determinations of deep oxygen abundance.

Binary neutron star associated with short Gamma ray burst has drawn wide attention ever since the observation of GW170817, due to its potential application in cosmology. While, further application of this sort of event suffers from the problem of degeneracy between luminosity distance and inclination angle, especially in the face-on limit. In this paper, we aim to address this issue by taking into account a Gaussian prior on the inclination angle. To test the property of resulting posterior distribution, we generate four catalogues of 1000 events by varying the number of third-generation detectors and the scale of prior. It turns out that a network of detectors tends to recognize more and farther events than a single detector. Besides, adopting tighter prior and employing multiple detectors both lead to lower error at given redshifts. Also considered is the validity of a widely adopted formula $\sigma_{\Delta d_{\rm L}0} = 2d_{\rm L} / \rho$, which undergoes the change from overestimating to underestimating with the increase of redshift. In the case of $\Lambda$CDM cosmology, 800, 300, 500 and 200 events are required for the four configurations to achieve $1\%$ $H_0$ accuracy. With all 1000 events in each catalogue, $H_0$ can be constrained to $0.9\%$, $0.4\%$, $0.7\%$ and $0.4\%$, while the errors of $\Omega_m$ are 0.014, 0.008, 0.016 and 0.012 respectively. Besides, adopting $\sigma_{\Delta d_{\rm L}0}$ leads to underestimation on the errors of cosmological parameters for tighter prior and overestimation in the opposite case. Results of Gaussian process also show that gravitational wave standard siren can reach higher redshift than traditional standard candles, especially when a network of detector is available, while alteration of the prior only has insignificant impact.

We have recently developed a method to kinematically decompose simulated galaxies that helps to break the degeneracy among galactic stellar structures. For example, the concept of stellar halos is generalized to weakly-rotating structures that are composed of loosely bound stars, which can hence be associated to both disk and elliptical type morphologies. By applying this method to about 500 central galaxies with stellar mass $10^{10-11.5}\ M_\odot$ from the TNG50 simulation at $z=0$, we identify three broadly-defined types of galaxies: ones dominated by disk, by bulge, or by stellar halo structures. We then use the simulation to infer the underlying connection between the growth of structures and physical processes over cosmic time. Tracing galaxies back in time, we recognize three fundamental regimes: an early phase of evolution ($z\gtrsim2$), and internal and external (mainly mergers) processes that act at later times. We find that disk- and bulge-dominated galaxies are not significantly affected by mergers since $z\sim2$; the difference in their present-day structures originates from two distinct evolutionary pathways, extended vs. compact, that are likely determined by their parent dark matter halos in the early phase; i.e., nature. On the other hand, normal elliptical galaxies are typically halo-dominated, forming by external processes (e.g. major mergers) in the later phase, i.e., nurture. This picture challenges the general idea that elliptical galaxies are the same objects as classical bulges in disk galaxies. In observations, both bulge- and halo-dominated galaxies are likely to be classified as early-type galaxies with compact morphology and quiescent star formation. However, here we find them to have very different evolutionary histories.

Quasi-periodic fast propagating (QFP) waves are often excited by solar flares, and could be trapped in the coronal structure with low Alfv\'en speed, so they could be used as a diagnosing tool for both the flaring core and magnetic waveguide. As the periodicity of a QFP wave could originate from a periodic source or be dispersively waveguided, it is a key parameter for diagnosing the flaring core and waveguide. In this paper, we study two QFP waves excited by a GOES-class C1.3 solar flare occurring at active region NOAA 12734 on 8 March 2019. Two QFP waves were guided by two oppositely oriented coronal funnel. The periods of two QFP waves were identical and were roughly equal to the period of the oscillatory signal in the X-ray and 17 GHz radio emission released by the flaring core. It is very likely that the two QFP waves could be periodically excited by the flaring core. Many features of this QFP wave event is consistent with the magnetic tuning fork model. We also investigated the seismological application with QFP waves, and found that the magnetic field inferred with magnetohydrodynamic seismology was consistent with that obtained in magnetic extrapolation model. Our study suggest that the QFP wave is a good tool for diagnosing both the flaring core and the magnetic waveguide.

Solar flares and plasma eruptions are sudden releases of magnetic energy stored in the plasma atmosphere. To understand the physical mechanisms governing their occurrences, three-dimensional magnetic fields from the photosphere up to the corona must be studied. The solar photospheric magnetic fields are observable, whereas the coronal magnetic fields cannot be measured. One method for inferring coronal magnetic fields is performing data-driven simulations, which involves time-series observational data of the photospheric magnetic fields with the bottom boundary of magnetohydrodynamic simulations. We developed a data-driven method in which temporal evolutions of the observational vector magnetic field can be reproduced at the bottom boundary in the simulation by introducing an inverted velocity field. This velocity field is obtained by inversely solving the induction equation and applying an appropriate gauge transformation. Using this method, we performed a data-driven simulation of successive small eruptions observed by the Solar Dynamics Observatory and the Solar Magnetic Activity Telescope in November 2017. The simulation well reproduced the converging motion between opposite-polarity magnetic patches, demonstrating successive formation and eruptions of helical flux ropes.

We present the updated photometric calibration of the twelve optical passbands for the Javalambre Photometric Local Universe Survey (J-PLUS) second data release (DR2), comprising 1088 pointings of two square degrees, and study the systematic impact of metallicity in the stellar locus technique. The [Fe/H] metallicity from LAMOST DR5 for 146184 high-quality calibration stars, defined with S/N > 10 in J-PLUS passbands and S/N > 3 in Gaia parallax, was used to compute the metallicity-dependent stellar locus (ZSL). The initial homogenization of J-PLUS photometry, performed with a unique stellar locus, was refined by including the metallicity effect in colours via the ZSL. The variation of the average metallicity along the Milky Way produces a systematic offset in J-PLUS calibration. This effect is well above 1% for the bluer passbands and amounts 0.07, 0.07, 0.05, 0.03, and 0.02 mag in u, J0378, J0395, J0410, and J0430, respectively. We modelled this effect with the Milky Way location of the J-PLUS pointing, providing also an updated calibration for those observations without LAMOST information. The estimated accuracy in the calibration after including the metallicity effect is at 1% level for the bluer J-PLUS passbands and below for the rest. We conclude that photometric calibration with the stellar locus technique is prone to significant systematic bias along the Milky Way location for passbands bluer than lambda = 4500 A. The updated calibration method for J-PLUS DR2 reaches 1-2% precision and 1% accuracy for twelve optical filters within an area of 2176 square degrees.

It is an open challenge to estimate systematically the physical parameters of neutron star interiors from pulsar timing data while separating spin wandering intrinsic to the pulsar (achromatic timing noise) from measurement noise and chromatic timing noise (due to propagation effects). In this paper we formulate the classic two-component, crust-superfluid model of neutron star interiors as a noise-driven, linear dynamical system and use a state-space-based expectation-maximization method to estimate the system parameters using gravitational-wave and electromagnetic timing data. Monte Carlo simulations show that we can accurately estimate all six parameters of the two-component model provided that electromagnetic measurements of the crust angular velocity, and gravitational-wave measurements of the core angular velocity, are both available. When only electromagnetic data are available we can recover the overall relaxation time-scale, the ensemble-averaged spin-down rate, and the strength of the white-noise torque on the crust. However, the estimates of the secular torques on the two components and white noise torque on the superfluid are biased significantly.

Magnetic braking (MB) likely plays a vital role in the evolution of low-mass X-ray binaries (LMXBs). However, it is still uncertain about the physics of MB, and there are various proposed scenarios for MB in the literature. To examine and discriminate the efficiency of MB, we investigate the LMXB evolution with five proposed MB laws. Combining detailed binary evolution calculation with binary population synthesis, we obtain the expected properties of LMXBs and their descendants binary millisecond pulsars. We then discuss the strength and weakness of each MB law by comparing the calculated results with observations. We conclude that the $\tau$-boosted MB law seems to best match the observational characteristics.

Binary millisecond pulsars (MSPs) are believed to have descended from low-mass X-ray binaries (LMXBs), which have experienced substantial mass transfer and tidal circularization. Therefore, they should have very circular orbits. However, the discovery of several eccentric binary MSPs (with eccentricity $e\sim 0.01-0.1$) challenges this standard picture. Three models have been proposed thus far based on accretion-induced collapse of massive white dwarfs (WDs), neutron star-strange star transition, and formation of circumbinary disks. All of them are subject to various uncertainties, and are not entirely consistent with observations. Here we propose an alternative model taking into account the influence of thermonuclear flashes on proto-WDs. We assume that the flashes lead to asymmetrical mass ejection, which imparts a mild kick on the proto-WDs. By simulating orbital changes of binary MSPs with multiple shell flashes, we show that it is possible to reproduce the observed eccentricities, provided that the kick velocities are around a few kms$^{-1}$.

This paper is the first paper of a series that will present the derivation of the modal mineralogy of Mars (M3 project) at a global scale from the near-infrared dataset acquired by the imaging spectrometer OMEGA (Observatoire pour la Min\'eralogie, l'Eau, les Glaces et l'Activit\'e) on board ESA/Mars Express. The objective is to create and provide a global 3-D image-cube of Mars at 32px/{\deg} covering most of Mars surface. This product has several advantages. First, it can be used to instantaneously extract atmospheric- and aerosol-corrected near-infrared (NIR) spectra from any location on Mars. Second, several new data maps can be built as discussed here. That includes new global mineral distributions, quantitative mineral abundance distributions and maps of Martian surface chemistry (wt % oxide) detailed in a companion paper (Riu et al., submitted). Here we present the method to derive the global hyperspectral cube from several hundred millions of spectra. Global maps of some mafic minerals are then shown, and compared to previous works.

LiteBIRD, the Lite (Light) satellite for the study of B-mode polarization and Inflation from cosmic background Radiation Detection, is a space mission for primordial cosmology and fundamental physics. JAXA selected LiteBIRD in May 2019 as a strategic large-class (L-class) mission, with its expected launch in the late 2020s using JAXA's H3 rocket. LiteBIRD plans to map the cosmic microwave background (CMB) polarization over the full sky with unprecedented precision. Its main scientific objective is to carry out a definitive search for the signal from cosmic inflation, either making a discovery or ruling out well-motivated inflationary models. The measurements of LiteBIRD will also provide us with an insight into the quantum nature of gravity and other new physics beyond the standard models of particle physics and cosmology. To this end, LiteBIRD will perform full-sky surveys for three years at the Sun-Earth Lagrangian point L2 for 15 frequency bands between 34 and 448 GHz with three telescopes, to achieve a total sensitivity of 2.16 micro K-arcmin with a typical angular resolution of 0.5 deg. at 100GHz. We provide an overview of the LiteBIRD project, including scientific objectives, mission requirements, top-level system requirements, operation concept, and expected scientific outcomes.

A radiative transfer model was used to reproduce several millions of OMEGA (Observatoire pour la Min\'eralogie, l'Eau, les Glaces et l'Activit\'e) spectra representative of igneous terrains of Mars. This task provided the modal composition and grain sizes at a planetary scale. The lithology can be summarized in five mineral maps at km-scale. We found that the low albedo equatorial regions of the Martian surface (from 60{\deg}S to 30{\deg}N) are globally dominated by plagioclase with average abundance ~50 vol% and pyroxenes with total averaged abundance close to 40 vol%. An evolution of the LCP/(LCP+HCP) ratio is observed with time at the global scale, suggesting an evolution of the degree of partial melting throughout the Martian eras. Olivine and Martian dust are minor components of the modelled terrains. The olivine distribution is quite different from the other minerals because it is found on localized areas with abundance reaching 20 vol%. A statistical approach, to classify the pixels of the abundances maps, using k-means clustering, highlighted seven distinct mineral assemblages on the surface. This classification illustrates that diverse mineralogical units are found in the Noachian and Hesperian terrains, which suggests the presence of various and complex magmatic processes at a global scale during the two oldest eras. The chemical composition was derived from the modal composition maps. The OMEGA-derived chemical composition is quite consistent with several distinctive geochemical characteristics previously considered as fingerprints of the Martian surface. A major discrepancy is in regards to the Fe content that is significantly smaller than soil and rock analyses from GRS and in situ measurements. The discrepancy could be partly explained by the assumptions used for the spectral modelling or could also indicate surface alteration rinds.

Coronal magnetic flux ropes are generally considered to be the core structure of large-scale solar eruptions. Recent observations found that solar eruptions could be initiated by a sequence of "flux feeding," during which chromospheric fibrils rise upward from below, and merge with a pre-existing prominence. Further theoretical study has confirmed that the flux feeding mechanism is efficient in causing the eruption of flux ropes that are wrapped by bald patch separatrix surfaces. But it is unclear how flux feeding influences coronal flux ropes that are wrapped by hyperbolic flux tubes (HFT), and whether it is able to cause the flux-rope eruption. In this paper, we use a 2.5-dimensional magnetohydrodynamic model to simulate the flux feeding processes in HFT configurations. It is found that flux feeding injects axial magnetic flux into the flux rope, whereas the poloidal flux of the rope is reduced after flux feeding. Flux feeding is able to cause the flux rope to erupt, provided that the injected axial flux is large enough so that the critical axial flux of the rope is reached. Otherwise, the flux rope system evolves to a stable equilibrium state after flux feeding, which might be even farther away from the onset of the eruption, indicating that flux feeding could stabilize the rope system with the HFT configuration in this circumstance.

After the successful landing of the Mars Science Laboratory rover, both NASA and ESA initiated a selection process for potential landing sites for the Mars2020 and ExoMars missions, respectively. Two ellipses located in the Mawrth Vallis region were proposed and evaluated during a series of meetings (3 for Mars2020 mission and 5 for ExoMars). We describe here the regional context of the two proposed ellipses as well as the framework of the objectives of these two missions. Key science targets of the ellipses and their astrobiological interests are reported. This work confirms the proposed ellipses contain multiple past Martian wet environments of subaerial, subsurface and/or subaqueous character, in which to probe the past climate of Mars, build a broad picture of possible past habitable environments, evaluate their exobiological potentials and search for biosignatures in well-preserved rocks. A mission scenario covering several key investigations during the nominal mission of each rover is also presented, as well as descriptions of how the site fulfills the science requirements and expectations of in situ martian exploration. These serve as a basis for potential future exploration of the Mawrth Vallis region with new missions and describe opportunities for human exploration of Mars in terms of resources and science discoveries.

In the framework of the ESA Athena mission, the X-ray Integral Field Unit (X-IFU) instrument to be on board the X-ray Athena Observatory is a cryogenic micro-calorimeter array of Transition Edge Sensor (TES) detectors aimed at providing spatially resolved high-resolution spectroscopy. As a part of the on-board Event Processor (EP), the reconstruction software will provide the energy, spatial location and arrival time of the incoming X-ray photons hitting the detector and inducing current pulses on it. Being the standard optimal filtering technique the chosen baseline reconstruction algorithm, different modifications have been analyzed to process pulses shorter than those considered of high resolution (those where the full length is not available due to a close pulse after them) in order to select the best option based on energy resolution and computing performance results. It can be concluded that the best approach to optimize the energy resolution for short filters is the 0-padding filtering technique, benefiting also from a reduction in the computational resources. However, its high sensitivity to offset fluctuations currently prevents its use as the baseline treatment for the X-IFU application for lack of consolidated information on the actual stability it will get in flight.

Most neutron stars are expected to be born in supernovae, but only about half of supernova remnants (SNRs) are associated with a compact object. In many cases, a supernova progenitor may have resulted in a black hole. However, there are several possible reasons why true pulsar-SNR associations may have been missed in previous surveys: The pulsar's radio beam may not be oriented towards us; the pulsar may be too faint to be detectable; or there may be an offset in the pulsar position caused by a kick. Our goal is to find new pulsars in SNRs and explore their possible association with the remnant. The search and selection of the remnants presented in this paper was inspired by the non-detection of any X-ray bright compact objects in these remnants when previously studied. Five SNRs were searched for radio pulsars with the Green Bank Telescope at 820 MHz with multiple pointings to cover the full spatial extent of the remnants. A periodicity search plus an acceleration search up to 500 m/s^2 and a single pulse search were performed for each pointing in order to detect potential isolated binary pulsars and single pulses, respectively. No new pulsars were detected in the survey. However, we were able to re-detect a known pulsar, PSR J2047+5029, near SNR G89.0+4.7. We were unable to detect the radio-quiet gamma-ray pulsar PSR J2021+4026, but we do find a flux density limit of 0.08 mJy. Our flux density limits make our survey two to 16 times more sensitive than previous surveys, while also covering the whole spatial extent of the same remnants. We discuss potential explanations for the non-detection of a pulsar in the studied SNRs and conclude that sensitivity is still the most likely factor responsible for the lack of pulsars in some remnants.

Here we aim to explore the origin of the strong C2H lines to reimagine the chemistry of protoplanetary disks. There are a few key aspects that drive our analysis. First, C2H is detected in young and old systems, hinting at a long-lived chemistry. Second, as a radical, C2H is rapidly destroyed, within <1000 yr. These two statements hint that the chemistry responsible for C2H emission must be predominantly in the gas-phase and must be in equilibrium. Combining new and published chemical models we find that elevating the total volatile (gas and ice) C/O ratio is the only natural way to create a long lived, high C2H abundance. Most of the \ce{C2H} resides in gas with a Fuv/n-gas ~ 10^-7 G0 cm^3. To elevate the volatile C/O ratio, additional carbon has to be released into the gas to enable an equilibrium chemistry under oxygen-poor conditions. Photo-ablation of carbon-rich grains seems the most straightforward way to elevate the C/O ratio above 1.5, powering a long-lived equilibrium cycle. The regions at which the conditions are optimal for the presence of high C/O ratio and elevated C2H abundances in the gas disk set by the Fuv/n-gas ~ 10^-7 G0 cm^3 condition lie just outside the pebble disk as well as possibly in disk gaps. This process can thus also explain the (hints of) structure seen in C2H observations.

Transmission spectroscopy of WASP-107b revealed 7-8% absorption at the position of metastable HeI triplet at 10830 {\AA} in Doppler velocity range of [-20; 10] km/s, which is stronger than that measured in other exoplanets. With a dedicated 3D self-consistent hydrodynamic multi-fluid model we calculated the expanding upper atmosphere of WASP-107b and reproduced within the observations accuracy the measured HeI absorption profiles, constraining the stellar XUV flux to (6-10) erg cm-2 s-1 at 1 a.u., and the upper atmosphere helium abundance He/H to 0.075-0.15. The radiation pressure acting on the metastable HeI atoms was shown to be an important factor affecting the shape of the absorption profiles. Its effect is counterbalanced by the processes of collisional depopulation of the HeI metastable state. Altogether, the observed HeI absorption in WASP-107b can be interpreted with the expected reasonable parameters of the stellar-planetary system and appropriate account of the electron and atom impact processes.

Constraining Jupiter's internal structure is crucial for understanding its formation and evolution history. Recent interior models of Jupiter that fit Juno's measured gravitational field suggest an inhomogeneous interior and potentially the existence of a diluted core. These models, however, strongly depend on the model assumptions and the equations of state used. A complementary modelling approach is to use empirical structure models. These can later be used to reveal new insights on the planetary interior and be compared to standard models. Here we present empirical structure models of Jupiter where the density profile is constructed by piecewise-polytropic equations. With these models we investigate the relation between the normalized moment of inertia (MoI) and the gravitational moments $J_2$ and $J_4$. Given that only the first few gravitational moments of Jupiter are measured with high precision, we show that an accurate and independent measurement of the MoI value could be used to further constrain Jupiter's interior. An independent measurement of the MoI with an accuracy better than $\sim 0.1\%$ could constrain Jupiter's core region and density discontinuities in its envelope. We find that models with a density discontinuity at $\sim$ 1 Mbar, as would produce a presumed hydrogen-helium separation, correspond to a fuzzy core in Jupiter. We next test the appropriateness of using polytropes, by comparing them with empirical models based on polynomials. We conclude that both representations result in similar density profiles and ranges of values for quantities like core mass and MoI.

We consider the possibility that LIGO events GW190521, GW190425 and GW190814 may have emerged from the mirror world binaries. Theories of star evolution predict so called upper and lower mass gaps and masses of these merger components lie in that gaps. In order to explain these challenging events very specific assumptions are required and we argue that such scenarios are order of magnitude more probable in mirror world, where star formation begins earlier and matter density can exceed 5 times the ordinary matter density.

Space weather phenomena such as solar flares, have massive destructive power when reaches certain amount of magnitude. Such high magnitude solar flare event can interfere space-earth radio communications and neutralize space-earth electronics equipment. In the current study, we explorer the deep learning approach to build a solar flare forecasting model and examine its limitations along with the ability of features extraction, based on the available time-series data. For that purpose, we present a multi-layer 1D Convolutional Neural Network (CNN) to forecast solar flare events probability occurrence of M and X classes at 1,3,6,12,24,48,72,96 hours time frame. In order to train and evaluate the performance of the model, we utilised the available Geostationary Operational Environmental Satellite (GOES) X-ray time series data, ranged between July 1998 and January 2019, covering almost entirely the solar cycles 23 and 24. The forecasting model were trained and evaluated in two different scenarios (1) random selection and (2) chronological selection, which were compare afterward. Moreover we compare our results to those considered as state-of-the-art flare forecasting models, both with similar approaches and different ones.The majority of the results indicates that (1) chronological selection obtain a degradation factor of 3\% versus the random selection for the M class model and elevation factor of 2\% for the X class model. (2) When consider utilizing only X-ray time-series data, the suggested model achieve high score results compare to other studies. (3) The suggested model combined with solely X-ray time-series fails to distinguish between M class magnitude and X class magnitude solar flare events. All source code are available at https://github.com/vladlanda/Low-Dimensional-Convolutional-Neural-Network-For-Solar-Flares-GOES-Time-Series-Classification

TAUKAM stands for "TAUtenburg KAMera", which will become the new prime-focus imager for the Tautenburg Schmidt telescope. It employs an e2v 6kx6k CCD and is under manufacture by Spectral Instruments Inc. We describe the design of the instrument and the auxiliary components, its specifications as well as the concept for integrating the device into the telescope infrastructure. First light is foreseen in 2017. TAUKAM will boost the observational capabilities of the telescope for what concerns optical wide-field surveys.

This study investigate the effectiveness of using Deep Learning (DL) for the classification of planetary nebulae (PNe). It focusses on distinguishing PNe from other types of objects, as well as their morphological classification. We adopted the deep transfer learning approach using three ImageNet pre-trained algorithms. This study was conducted using images from the Hong Kong/Australian Astronomical Observatory/Strasbourg Observatory H-alpha Planetary Nebula research platform database (HASH DB) and the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS). We found that the algorithm has high success in distinguishing True PNe from other types of objects even without any parameter tuning. The Matthews correlation coefficient is 0.9. Our analysis shows that DenseNet201 is the most effective DL algorithm. For the morphological classification, we found for three classes, Bipolar, Elliptical and Round, half of objects are correctly classified. Further improvement may require more data and/or training. We discuss the trade-offs and potential avenues for future work and conclude that deep transfer learning can be utilized to classify wide-field astronomical images.

Galaxy clusters exhibit a rich morphology during the early and intermediate stages of mass assembly, especially beyond their boundary. A classification scheme based on shapefinders deduced from the Minkowski functionals is examined to fully account for the morphological diversity of galaxy clusters, including relaxed and merging clusters, clusters fed by filamentary structures, and cluster-pair bridges. These configurations are conveniently treated with idealised geometric models and analytical formulae, some of which are novel. Examples from CLASH and LC$^2$ clusters and observed cluster-pair bridges are discussed.

We present an investigation of the low-frequency radio and ultraviolet properties of a sample of $\simeq$10,500 quasars from the Sloan Digital Sky Survey Data Release 14, observed as part of the first data release of the Low-Frequency-Array (LOFAR) Two-metre Sky Survey. The quasars have redshifts $1.5 < z < 3.5$ and luminosities $44.6 < \log(L_{\text{bol}}/\text{erg s}^{-1}) < 47.2$. We employ ultraviolet spectral reconstructions based on an independent component analysis to parametrize the CIV$\lambda$1549-emission line which is used to infer the strength of accretion disc winds, and the HeII$\lambda$1640 line, an indicator of the soft X-ray flux. We find that radio-detected quasars are found in the same region of CIV blueshift versus equivalent-width space as radio-undetected quasars, but that the loudest, most luminous and largest radio sources exist preferentially at low CIV blueshifts. Additionally, the radio-detection fraction increases with blueshift whereas the radio-loud fraction decreases. In the radio-quiet population, we observe a range of HeII equivalent widths as well as a Baldwin effect with bolometric luminosity, whilst the radio-loud population has mostly strong HeII, consistent with a stronger soft X-ray flux. The presence of strong HeII is a necessary but not sufficient condition to detect radio-loud emission suggesting some degree of stochasticity in jet formation. Using energetic arguments and Monte Carlo simulations, we explore the plausibility of winds, compact jets and star formation as sources of the radio quiet emission, ruling out none. The existence of quasars with similar ultraviolet properties but differing radio properties suggests, perhaps, that the radio and ultraviolet emission is tracing activity occurring on different timescales.

Martian araneiform terrain, located in the Southern polar regions, consists of features with central pits and radial troughs which are thought to be associated with the solid state greenhouse effect under a CO$_{2}$ ice sheet. Sublimation at the base of this ice leads to gas buildup, fracturing of the ice and the flow of gas and entrained regolith out of vents and onto the surface. There are two possible pathways for the gas: through the gap between the ice slab and the underlying regolith, as proposed by Kieffer et al (2007), or through the pores of a permeable regolith layer, which would imply that regolith properties can control the spacing between adjacent spiders, as suggested by Hao et al. We test this hypothesis quantitatively in order to place constraints on the regolith properties. Based on previously estimated flow rates and thermophysical arguments, we suggest that there is insufficient depth of porous regolith to support the full gas flow through the regolith. By contrast, free gas flow through a regolith--ice gap is capable of supplying the likely flow rates for gap sizes on the order of a centimetre. This size of gap can be opened in the centre of a spider feature by gas pressure bending the overlying ice slab upwards, or by levitating it entirely as suggested in the original Kieffer et al (2007) model. Our calculations therefore support at least some of the gas flowing through a gap opened between the regolith and ice. Regolith properties most likely still play a role in the evolution of spider morphology, by regolith cohesion controlling the erosion of the central pit and troughs, for example.

We present a detailed study of the methodology for correlating `dark sirens' (compact binaries coalescences without electromagnetic counterpart) with galaxy catalogs. We propose several improvements on the current state of the art, and we apply them to the published LIGO/Virgo gravitational wave (GW) detections, performing a detailed study of several sources of systematic errors that, with the expected increase in statistics, will eventually become the dominant limitation. We provide a measurement of $H_0$ from dark sirens alone, finding as the best result $H_0=75^{+25}_{-22}\,\,{\rm km}\, {\rm s}^{-1}\, {\rm Mpc}^{-1}$ ($68\%$ c.l., for a flat prior in the range $[20,140] \,\,{\rm km}\, {\rm s}^{-1}\, {\rm Mpc}^{-1}$) which is, currently, the most stringent constraint obtained using only dark sirens. Combining dark sirens with the counterpart for GW170817 we find $H_0=70^{+11}_{-7} \,{\rm km}\, {\rm s}^{-1}\, {\rm Mpc}^{-1}$. We also study modified GW propagation, which is a smoking gun of dark energy and modifications of gravity at cosmological scales, and we show that current observations of dark sirens already start to provide interesting limits. From dark sirens alone, our best result for the parameter $\Xi_0$ that measures deviations from GR (with $\Xi_0=1$ in GR) is $\Xi_0=1.88^{+3.83}_{-1.10}$. We finally discuss limits on modified GW propagation under the tentative identification of the flare ZTF19abanrhr as the electromagnetic counterpart of the binary black hole coalescence GW190521, in which case our most stringent result is $\Xi_0=1.6^{+1.0}_{-0.6}$. We release the publicly available code $\tt{DarkSirensStat}$, which is available under open source license at https://github.com/CosmoStatGW/DarkSirensStat.

We present an extensive search in the literature and Gaia DR2 for visual co-moving binary companions to stars hosting exoplanets and brown dwarfs within 200 pc. We found 218 planet hosts out of 938 to be part of multiple-star systems, with 10 newly discovered binaries and 2 new tertiary stellar components. This represents an overall raw multiplicity rate of 23.2$\pm$1.6% for hosts to exoplanets across all spectral types, with multi-planet systems found to have a lower duplicity frequency at the 2.2$\sigma$ level. We found that more massive hosts are more often in binary configurations, and that planet-bearing stars in multiple systems are predominantly the most massive component of stellar binaries. Investigations of multiplicity as a function of planet mass and separation revealed that giant planets with masses >0.1 MJup are more frequently seen in stellar binaries than small sub-Jovian planets with a 3.6$\sigma$ difference, a trend enhanced for the most massive (>7 MJup) short-period (<0.5 AU) planets and brown dwarf companions. Binarity was found to have no significant effect on the demographics of low-mass planets (<0.1 MJup) or warm and cool gas giants (>0.5 AU). While stellar companion mass appears to have no impact on planet properties, binary separation seems to be an important factor in the resulting structure of planetary systems. Stellar companions on separations <1000 AU can play a role in the formation or evolution of massive close-in planets, while planets in wider binaries show similar properties to planets orbiting single stars. Finally, numerous stellar companions on separations <1-3 arcsec likely remain undiscovered to this date. Continuous efforts to complete our knowledge of stellar multiplicity on separations of tens to hundreds of AU are essential to confirm the reported trends and further our understanding of the roles played by multiplicity on exoplanets.

We study the effect of different Type Ia SN nucleosynthesis prescriptions on the Milky Way chemical evolution. To this aim, we run detailed one-infall and two-infall chemical evolution models, adopting a large compilation of yield sets corresponding to different white dwarf progenitors (near-Chandrasekar and sub-Chandrasekar) taken from the literature. We adopt a fixed delay time distribution function for Type Ia SNe , in order to avoid degeneracies in the analysis of the different nucleosynthesis channels. We also combine yields for different Type Ia SN progenitors in order to test the contribution to chemical evolution of different Type Ia SN channels. The results of the models are compared with recent LTE and NLTE observational data. We find that "classical" W7 and WDD2 models produce Fe masses and [$\alpha$/Fe] abundance patterns similar to more recent and physical near-Chandrasekar and sub- Chandrasekar models. For Fe-peak elements, we find that the results strongly depend either on the white dwarf explosion mechanism (deflagration-to-detonation, pure deflagration, double detonation) or on the initial white dwarf conditions (central density, explosion pattern). The comparison of chemical evolution model results with observations suggests that a combination of near-Chandrasekar and sub-Chandrasekar yields is necessary to reproduce the data of V, Cr, Mn and Ni, with different fractions depending on the adopted massive stars stellar yields. This comparison also suggests that NLTE and singly ionised abundances should be definitely preferred when dealing with most of Fe-peak elements at low metallicity.

We present the second data release for the HI-MaNGA programme of HI follow-up observations for the SDSS-IV MaNGA survey. This release contains HI measurements from 2135 Green Bank Telescope (GBT) observations and an updated crossmatch of the SDSS-IV MaNGA sample with the ALFALFA survey. We combine these data with MaNGA spectroscopic measurements to examine relationships between HI-to-stellar mass ratio (M_HI/M_*) and average ISM/star formation properties probed by optical emission lines. M_HI/M_* is very weakly correlated with the equivalent width of Halpha, implying a loose connection between the instantaneous star formation rate and the HI reservoir, although the link between M_HI/M_* and star formation strengthens when averaged even over only moderate timescales (~30 Myrs). Galaxies with elevated HI depletion times have enhanced [OI]/Halpha and depressed Halpha surface brightness, consistent with more HI residing in a diffuse and/or shock heated phase which is less capable of condensing into molecular clouds. Of all optical lines, M_HI/M_* correlates most strongly with oxygen equivalent width, EW(O), which is likely a result of the existing correlation between M_HI/M_* and gas-phase metallicity. Residuals in the M_HI/M_*-EW(O) relation are again correlated with [OI]/Halpha and Halpha surface brightness, suggesting they are also driven by variations in the fraction of diffuse and/or shock-heated gas. We recover the strong anti-correlation between M_HI/M_* and gas-phase metallicity seen in previous studies. We also find a positive relationship between M_HI/M_* and [OI]/Halpha, independent of metallicity, which may suggest that higher fractions of diffuse and/or shock heated gas are more prevalent in gas-rich galaxies.

We report the relationship between the $Gaia$--VLBI position differences and the magnitudes of source structure effects in VLBI observations. Because the $Gaia$--VLBI position differences are statistically significant for a considerable number of common sources, we attempt to discuss and explain these position differences based on VLBI observations and available source images at cm-wavelengths. Based on the derived closure amplitude root-mean-square (CARMS), which quantifies the magnitudes of source structure effects in the VLBI observations used for building the third realization of the International Celestial Reference Frame, the arc lengths and normalized arc lengths of the position differences are examined in detail. The radio jet directions and the directions of the $Gaia$--VLBI position differences are investigated for a small sample of sources. Both the arc lengths and normalized arc lengths of the $Gaia$ and VLBI positions are found to increase with the CARMS values. The majority of the sources with statistically significant position differences are associated with the sources having extended structure. Radio source structure is the one of the major factors of these position differences, and it can be the dominate factor for a number of sources. The vectors of the $Gaia$ and VLBI position differences are parallel to the radio-jet directions, which is confirmed with stronger evidence.

Recent developments in astronomical observations enable direct imaging of circumstellar disks. Precise characterization of such extended structure is essential to our understanding of stellar systems. However, the faint intensity of the circumstellar disks compared to the brightness of the host star compels astronomers to use tailored observation strategies, in addition to state-of-the-art optical devices. Even then, extracting the signal of circumstellar disks heavily relies on post-processing techniques. In this work, we propose a morphological component analysis (MCA) approach that leverages low-complexity models of both the disks and the stellar light corrupting the data. In addition to disks, our method allows to image exoplanets. Our approach is tested through numerical experiments.

In this work, we explore the dynamical impacts and observable signatures of two-fluid effects in the parameter regimes when ion-neutral collisions do not fully couple the neutral and charged fluids. The purpose of this study is to deepen our understanding of the RTI and the effects of the partial ionization on the development of RTI using non-linear two-fluid numerical simulations. Our two-fluid model takes into account neutral viscosity, thermal conductivity, and collisional interaction between neutrals and charges: ionization/recombination, energy and momentum transfer, and frictional heating. In this paper II, the sensitivity of the RTI dynamics to collisional effects for different magnetic field configurations supporting the prominence thread is explored. This is done by artificially varying, or eliminating, effects of both elastic and inelastic collisions by modifying the model equations. We find that ionization and recombination reactions between ionized and neutral fluids, if in equilibrium prior to the onset of the instability, do not substantially impact the development of the primary RTI. However, such reactions can impact development of secondary structures during mixing of the cold prominence and hotter surrounding coronal material. We find that collisionality within and between ionized and neutral particle populations play an important role in both linear and non-linear development of RTI, with ion-neutral collision frequency as the primary determining factor in development or damping of small scale structures. We also observe that degree and signatures of flow decoupling between ion and neutral fluids can depend both on the inter-particle collisionality and the magnetic field configuration of the prominence thread.

In this work we study 35 stellar clusters in the Small Magellanic Cloud (SMC) in order to provide their mean metallicities and ages. We also provide mean metallicities of the fields surrounding the clusters. We used Str\"omgren photometry obtained with the 4.1 m SOAR telescope and take advantage of $(b - y)$ and $m1$ colors for which there is a metallicity calibration presented in the literature. The spatial metallicity and age distributions of clusters across the SMC are investigated using the results obtained by Str\"omgren photometry. We confirm earlier observations that younger, more metal-rich star clusters are concentrated in the central regions of the galaxy, while older, more metal-poor clusters are located farther from the SMC center. We construct the age-metallicity relation for the studied clusters and find good agreement with theoretical models of chemical enrichment, and with other literature age and metallicity values for those clusters. We also provide the mean metallicities for old and young populations of the field stars surrounding the clusters, and find the latter to be in good agreement with recent studies of the SMC Cepheid population. Finally, the Str\"omgren photometry obtained for this study is made publicly available.

Since its launch, the Alpha Magnetic Spectrometer - 02 (AMS-02) has delivered outstanding quality measurements of the spectra of cosmic-ray (CR) species, $\bar{p}$, $e^{\pm}$, and nuclei, $_1$H-$_8$O, $_{10}$Ne, $_{12}$Mg, $_{14}$Si, which resulted in a number of breakthroughs. One of the latest long awaited surprises is the spectrum of $_{26}$Fe just published by AMS-02. Because of the large fragmentation cross section and large ionization energy losses, most of CR iron at low energies is local, and may harbor some features associated with relatively recent supernova (SN) activity in the solar neighborhood. Our analysis of new iron spectrum together with Voyager 1 and ACE-CRIS data reveals an unexpected bump in the iron spectrum and in the Fe/He, Fe/O, and Fe/Si ratios at 1-2 GV, while a similar feature in the spectra of He, O, Si, and in their ratios is absent, hinting at a local source of low-energy CRs. The found excess fits well with recent discoveries of radioactive $^{60}$Fe deposits in terrestrial and lunar samples, and in CRs. We provide an updated local interstellar spectrum (LIS) of iron in the energy range from 1 MeV nucleon$^{-1}$ to $\sim$10 TeV nucleon$^{-1}$. Our calculations employ the GalProp-HelMod framework that is proved to be a reliable tool in deriving the LIS of CR $\bar{p}$, $e^{-}$, and nuclei $Z\le28$.

Recent work reported the discovery of a gamma-ray burst (GRB) associated with the galaxy GN-z11 at $z\sim 11$. The extreme improbability of the transient source being a GRB in the very early Universe requires robust elimination of all plausible alternative hypotheses. We identify numerous examples of similar transient signals in separate archival MOSFIRE observations and argue that Solar system objects -- natural or artificial -- are a far more probable explanation for these phenomena.

We develop new strategies to build numerical relativity surrogate models for eccentric binary black hole systems, which are expected to play an increasingly important role in current and future gravitational-wave detectors. We introduce a new surrogate waveform model, \texttt{NRSur2dq1Ecc}, using 47 nonspinning, equal-mass waveforms with eccentricities up to $0.2$ when measured at a reference time of $5500M$ before merger. This is the first waveform model that is directly trained on eccentric numerical relativity simulations and does not require that the binary circularizes before merger. The model includes the $(2,2)$, $(3,2)$, and $(4,4)$ spin-weighted spherical harmonic modes. We also build a final black hole model, \texttt{NRSur2dq1EccRemnant}, which models the mass, and spin of the remnant black hole. We show that our waveform model can accurately predict numerical relativity waveforms with mismatches $\approx 10^{-3}$, while the remnant model can recover the final mass and dimensionless spin with errors smaller than $\approx 5 \times 10^{-4}M$ and $\approx 2 \times10^{-3}$ respectively. We demonstrate that the waveform model can also recover subtle effects like mode-mixing in the ringdown signal without any special ad-hoc modeling steps. Finally, we show that despite being trained only on equal-mass binaries, \texttt{NRSur2dq1Ecc} can be reasonably extended up to mass ratio $q\approx3$ with mismatches $\simeq 10^{-2}$ for eccentricities smaller than $\sim 0.05$ as measured at a reference time of $2000M$ before merger. The methods developed here should prove useful in the building of future eccentric surrogate models over larger regions of the parameter space.

We search for gravitational-wave signals produced by cosmic strings in the Advanced LIGO and Virgo full O3 data set. Search results are presented for gravitational waves produced by cosmic string loop features such as cusps, kinks and, for the first time, kink-kink collisions.cA template-based search for short-duration transient signals does not yield a detection. We also use the stochastic gravitational-wave background energy density upper limits derived from the O3 data to constrain the cosmic string tension, $G\mu$, as a function of the number of kinks, or the number of cusps, for two cosmic string loop distribution models.cAdditionally, we develop and test a third model which interpolates between these two models. Our results improve upon the previous LIGO-Virgo constraints on $G\mu$ by one to two orders of magnitude depending on the model which is tested. In particular, for one loop distribution model, we set the most competitive constraints to date, $G\mu\lesssim 4\times 10^{-15}$.

In this paper, we continue studying the thermodynamics of Hayward black hole, which has been recently approached by Molina \& Villanueva regarding the laws of black hole thermodynamics, by introducing the Hayward's parameter as being responsible for a possible regularization of the Schwarzschild black hole. Here, we show that the adiabatic foliations of the thermodynamic manifold are confined by the extremal subspace, and therefore, the latter cannot be reached adiabatically. A direct consequence of this features, is the impossibility of the merger of two extremal Hayward black holes.

In the same base setup as Sakharov's induced gravity, we investigate emergence of gravity in effective quantum field theories (QFT), with particular emphasis on the gauge sector in which gauge bosons acquire anomalous masses in proportion to the ultraviolet cutoff $\Lambda_\wp$. Drawing on the fact that $\Lambda_\wp^2$ corrections explicitly break the gauge and Poincare symmetries, we find that it is possible to map $\Lambda_\wp^2$ to spacetime curvature as a covariance relation and we find also that this map erases the anomalous gauge boson masses. The resulting framework describes gravity by the general relativity (GR) and matter by the QFT itself with $\log\Lambda_\wp$ corrections (dimensional regularization). This QFT-GR concord predicts existence of new physics beyond the Standard Model such that the new physics can be a weakly-interacting or even a non-interacting sector comprising the dark matter, dark energy and possibly more. The concord has consequential implications for collider, astrophysical and cosmological phenomena.

We determine the complete space-time metric from the bootstrapped Newtonian potential generated by a static spherically symmetric source in the surrounding vacuum. This metric contains post-Newtonian parameters which can be further used to constrain the complete underlying dynamical theory. For values of the post-Newtonian parameters within experimental bounds, the reconstructed metric appears very close to the Schwarzschild solution of General Relativity in the whole region outside the event horizon. The latter is however larger in size for the same value of the mass compared to the Schwarzschild case.

Recent developments in the gravitational waves interferometry require more pertinent theoretical models of gravitational waves generation and propagation. Untouched possible mechanisms of spin-2 spacetime perturbations production, we will consider their subsequent scattering on other black holes (BHs). Specifically, we consider a generalization of the Regge-Wheeler-Zerilli equations for the case of distorted BHs (BHs surrounded with matter) in Minkowski and Anti-de Sitter spacetimes, the metric potential of which obeys the Liouville equation. We establish significant differences in scattering characteristics of waves of different spins and angular momenta, including the gravitational waves, caused by losing the spherical symmetry of their propagation background. In particular, we demonstrate the strong impact of the background geometry deformation on the grey-body factors, hence on the absorption cross-sections of scattering waves, and explore the issue of stability of the background geometry upon changing the deformation degree parameters.

Using the framework of higher-form global symmetries, we examine the regime of validity of force-free electrodynamics by evaluating the lifetime of the electric field operator, which is non-conserved due to screening effects. We focus on a holographic model which has the same global symmetry as that of low energy plasma and obtain the lifetime of (non-conserved) electric flux in a strong magnetic field regime. The lifetime is inversely correlated to the magnetic field strength and thus suppressed in the strong field regime.