Papers for Wednesday, Sep 08 2021

Context: The current sample of known brown dwarfs (BDs) around FGK-stars is only of the order of a hundred. The ongoing ESA mission Gaia has already collected its nominal 5 yr of mission data and might operate up to 10 yr. Aims: Using detailed simulations, we estimate the number of BDs that could be discovered by Gaia astrometry, radial velocity, and photometric transits around main sequence (V) and subgiants (IV) FGK host stars for a 5 and 10-yr mission. Methods: Using a robust $\Delta \chi^2$ statistic we analyse the BD companion detectability from the Besan\c{c}on Galaxy population synthesis model complemented by Gaia DR2 data for the bright end, using the latest Gaia performance and scanning law, and literature-based BD-parameter distributions. Results: We report here reliable detection numbers ($\Delta \chi^2$>50) for a 5-yr [10-yr] mission. Astrometry alone yields 28,000-42,000 [45,000-55,000] detections out to several hundred pc [>kiloparsec], with the majority around G magnitude 14-15 [14-16] and P>200 d. Gaia radial velocity time series allow detection of 830-1100 [1500-1900] mainly massive BDs (55-80 M_J), most having P <10 d. Systems with at least 3 photometric transits (S/N>3) are expected for 720-1100 [1400-2300] BDs, averaging at 4-5 [5-6] transits per source. Overlap of astrometric and radial velocity detection yield 370-410 [870-950] candidates, transit and radial velocity 17-27 [35-56], and transit and astrometric detection 1-3 [4-6]. Conclusions: Though above numbers have +/- 50% uncertainty due to the uncertain occurrence rate and period distribution of BDs around FGK host stars, Gaia detections will number in the tens-of-thousands, enlarging the current sample by at least two orders of magnitude, allowing to investigate the BD fraction and orbital architectures as a function of host stellar parameters in greater detail than every before.

The tidal disruption event, AT2019dsg, was observed across a broad range of electromagnetic-wavelengths from radio to X-rays, and it was possibly accompanied by a high-energy neutrino. We study the nature of the radio-emitting outflow by analyzing the synchrotron self-absorbed spectra in terms of the equipartition model. We find that the time evolution of the outflow radius can be interpreted as either free expansion or deceleration. If the former, the outflow was launched $\simeq$ 40 days before the optical peak; if the latter, the launch was $\simeq$ 10 days after the optical peak. In addition, the energy in the radio-emitting region increases over time. This second conclusion is most naturally interpreted by a scenario resembling the earliest stage of a supernova remnant: as more and more material is swept up, it is heated by the forward shock at the expense of the outflow's kinetic energy. Energy injection from an accreting BH cannot be completely excluded, but the injection rate is very different from the fallback luminosity, requiring further physical explanation. If the neutrino association is real, the scale of energy injection needed is much greater than for the radio emission, suggesting that the detected neutrino did not arise from the radio emitting region.

The new generation of wide-field time-domain surveys has made it feasible to study the clustering of supernova (SN) host galaxies in the large-scale structure (LSS) for the first time. We investigate the LSS environment of SN populations, using 106 dark matter density realisations with a resolution of $\sim$ 3.8 Mpc, constrained by the 2M++ galaxy survey. We limit our analysis to redshift $z<0.036$, using samples of 498 thermonuclear and 782 core-collapse SNe from the Zwicky Transient Facility's Bright Transient Survey and Census of the Local Universe catalogues. We detect clustering of SNe with high significance; the observed clustering of the two SNe populations is consistent with each other. Further, the clustering of SN hosts is consistent with that of the Sloan Digital Sky Survey (SDSS) Baryon Oscillation Spectroscopic Survey (BOSS) DR12 spectroscopic galaxy sample in the same redshift range. Using a tidal shear classifier, we classify the LSS into voids, sheets, filaments and knots. We find that both SNe and SDSS galaxies are predominantly found in sheets and filaments. SNe are significantly under-represented in voids and over-represented in knots compared to the volume fraction in these structures. This work opens the potential for using forthcoming wide-field deep SN surveys as a complementary LSS probe.

WL 17 is a young transition disk in the Ophiuchus L1688 molecular cloud complex. Even though WL 17 is among the brightest disks in L1688 and massive enough to expect dust self-scattering, it was undetected in polarization down to ALMA's instrument sensitivity limit. Such low polarization fractions could indicate unresolved polarization within the beam or optically thin dust emission. We test the latter case by combining the high sensitivity 233 GHz Stokes I data from the polarization observations with previous ALMA data at 345 GHz and 100 GHz. We use simple geometric shapes to fit the observed disk visibilities in each band. Using our simple models and assumed dust temperature profiles, we estimate the optical depth in all three bands. The optical depth at 233 GHz peaks at $\tau_{233} \sim 0.3$, which suggests the dust emission may not be optically thick enough for dust self-scattering to be efficient. We also find the higher sensitivity 233 GHz data show substructure in the disk for the first time. The substructure appears as brighter lobes along the major axis, on either side of the star. We attempt to fit the lobes with a simple geometric model, but they are unresolved in the 233 GHz data. We propose that the disk may be flared at 1 mm such that there is a higher column of dust along the major axis than the minor axis when viewed at an inclination. These observations highlight the strength of high sensitivity continuum data from dust polarization observations to study disk structures.

Driven by increasingly complete observational knowledge of systems of satellite galaxies, mutual spatial alignments and relations in velocities among satellites belonging to a common host have become a productive field of research. Numerous studies have investigated different types of such phase-space correlations, and were met with varying degrees of attention by the community. The Planes of Satellite Galaxies issue is maybe the best-known example, with a rich field of research literature and an ongoing, controversial debate on how much of a challenge it poses to the LCDM model of cosmology. Another type of correlation, the apparent excess of close pairs of dwarf galaxies, has received considerably less attention despite its reported tension with LCDM expectations. With the expansion of proper motion measurements, largely driven by Gaia, other peculiar phase-space correlations were uncovered among the Milky Way satellites. Examples are the tangential velocity excess of satellites compared to cosmological expectations, and the unexpected preference of satellites to be close to pericenter. At the same time, other kinds of correlations have been found to be more in line with cosmological expectations, specifically lopsided satellite galaxy systems and the accretion of groups of satellite galaxies. The latter has mostly been studied in cosmological simulations thus far, but offers the potential to address some of the other issues by providing a way to produce correlations among the orbits of a group's satellite galaxy members. This review is the first to provide an introduction to the highly active field of phase-space correlations among satellite galaxy systems. It summarizes existing, recent research and highlights interdependencies between the different, currently almost exclusively individually considered types of correlations. Future prospects are also briefly discussed. (abridged)

The formation of stars and planetary systems is a complex phenomenon, which relies on the interplay of multiple physical processes. Nonetheless, it represents a crucial stage for our understanding of the Universe, and in particular of the conditions leading to the formation of key molecules (e.g. water) on comets and planets. {\it Herschel} observations demonstrated that stars form out of gaseous filamentary structures in which the main constituent is molecular hydrogen (H$_2$). Depending on its nuclear spin H$_2$ can be found in two forms: `ortho' with parallel spins and `para' where the spins are anti-parallel. The relative ratio among these isomers, i.e.\,the ortho-to-para ratio (OPR), plays a crucial role in a variety of processes related to the thermodynamics of star-forming gas and to the fundamental chemistry affecting the formation of water in molecular clouds and its subsequent deuteration, commonly used to determine the origin of water in Solar System's bodies. Here, for the first time, we assess the evolution of the OPR starting from the warm neutral medium, by means of state-of-the-art three-dimensional magneto-hydrodynamic simulations of turbulent molecular clouds. Our results show that star-forming clouds exhibit a low OPR ($\ll 0.1$) already at moderate densities ($\sim$1000 cm$^{-3}$). We also constrain the cosmic rays ionisation rate, finding that $10^{-16}\,\rm s^{-1}$ is the lower limit required to explain the observations of diffuse clouds. Our results represent a step forward in the understanding of the star and planet formation process providing a robust determination of the chemical initial conditions for both theoretical and observational studies.

Making use of both MUSE observations of 85 galaxies from the survey GASP (GAs Stripping Phenomena in galaxies with MUSE) and a large sample from MaNGA (Mapping Nearby Galaxies at Apache Point Observatory survey) we investigate the distribution of gas metallicity gradients as a function of stellar mass, for local cluster and field galaxies. Overall, metallicity profiles steepen with increasing stellar mass up to $10^{10.3}\,{\rm M_\odot}$ and flatten out at higher masses. Combining the results from the metallicity profiles and the stellar mass surface density gradients, we propose that the observed steepening is a consequence of local metal enrichment due to in-situ star formation during the inside-out formation of disk galaxies. The metallicity gradient-stellar mass relation is characterized by a rather large scatter, especially for $10^{9.8}<{\rm M_\star/M_\odot}<10^{10.5}$, and we demonstrate that metallicity gradients anti-correlate with the galaxy gas fraction. Focusing on the galaxy environment, at any given stellar mass, cluster galaxies have systematically flatter metallicity profiles than their field counterparts. Many sub-populations coexist in clusters: galaxies with shallower metallicity profiles appear to have fallen into their present host halo sooner and have experienced the environmental effects for a longer time than cluster galaxies with steeper metallicity profiles. Recent galaxy infallers, like galaxies currently undergoing ram-pressure stripping, show metallicity gradients more similar to those of field galaxies, suggesting they have not felt the effect of the cluster yet.

We present a follow-up analysis examining the dynamics and structures of 41 massive, large star-forming galaxies at z~0.67-2.45 using both ionized and molecular gas kinematics. We fit the galaxy dynamics with models consisting of a bulge, a thick, turbulent disk, and a NFW dark matter halo, using code that fully forward models the kinematics, including all observational and instrumental effects. We explore the parameter space using Markov Chain Monte Carlo (MCMC) sampling, including priors based on stellar and gas masses and disk sizes. We fit the full sample using extracted 1D kinematic profiles. For a subset of 14 well-resolved galaxies, we also fit the 2D kinematics. The MCMC approach robustly confirms the results from least-squares fitting presented in Paper I (Genzel et al. 2020): the sample galaxies tend to be baryon-rich on galactic scales (within one effective radius). The 1D and 2D MCMC results are also in good agreement for the subset, demonstrating that much of the galaxy dynamical information is captured along the major axis. The 2D kinematics are more affected by the presence of non-circular motions, which we illustrate by constructing a toy model with constant inflow for one galaxy that exhibits residual signatures consistent with radial motions. This analysis, together with results from Paper I and other studies, strengthens the finding that massive, star-forming galaxies at z~1-2 are baryon-dominated on galactic scales, with lower dark matter fractions towards higher baryonic surface densities. Finally, we present details of the kinematic fitting code used in this analysis.

We derive abundance ratios for nine stars in the relatively high-metallicity bulge globular cluster NGC 6380. We find a mean cluster metallicity between [Fe/H]$= -0.80$ and $-0.73$, with no clear evidence for a variation in iron abundances beyond the observational errors. Stars with strongly enhanced in [N/Fe] abundance ratios populate the cluster, and are anti-correlated with [C/Fe], trends that are considered a signal of the multiple-population phenomenon in this cluster. We detect an apparent intrinsic star-to-star spread ($\gtrsim 0.27$ dex) in the slow neutron-capture process element (s-element) Ce II. Moreover, the [Ce/Fe] abundance ratio exhibits a likely correlation with [N/Fe], and a somewhat weaker correlation with [Al/Fe]. If confirmed, NGC 6380 could be the first high-metallicity globular cluster where a N-Ce correlation is detected. Furthermore, this correlation suggests that Ce may also be an element involved in the multiple-population phenomenon. Currently, a consensus interpretation for the origin of the this apparent N-Ce correlation in high-metallicity clusters is lacking. We tentatively suggest that it could be reproduced by different channels - low-mass asymptotic giant-branch stars in the high-metallicity regime or fast-rotating massive stars ("spinstars"), due to the rotational mixing. It may also be the cumulative effect of several pollution events including the occurrence of peculiar stars. Our findings should guide stellar nucleosynthesis models, in order to understand the reasons for its apparent exclusivity in relatively high-metallicity globular clusters.

Fast radio bursts (FRBs) are enigmatic transients with very short-duration radio emission. Their nature is still unknown and is widely debated. I provide the first analysis of atomic gas properties of FRB hosts to provide constraints on their nature. HI observations exist for NGC3252, the host of FRB 181030A, M81, the host of FRB 200120E, and the Milky Way, the host of FRB 200428. I report three observables: i) all three FRB hosts are interacting galaxies; ii) the HI spectra of both FRB hosts with such data available are highly asymmetric, several standard deviations above the general population of galaxies; iii) two FRB hosts have normal atomic gas properties and one is strongly deficient in atomic gas. This indicates that FRBs are connected with a recent enhancement of star formation due to interaction. This supports fast FRB channels, for example a massive star with a short delay time so that interaction signatures giving rise to the birth of the progenitor are still visible. Long gamma-ray burst (GRB) and broad-lined type Ic supernova (SN) hosts exhibit much more symmetric spectra, even though they were claimed to experience gas inflow from the intergalactic medium. The difference can be explained by the interactions experienced by FRB hosts being more disruptive than these gas inflow, or by the mass effect, with GRB/SN host at lower masses having less organized gas motions, so with HI lines closer to a symmetrical Gaussian. This also suggests that the emission mechanisms of FRBs and GRBs are likely different.

The co-evolution between central supermassive black holes (BH), their host galaxies, and dark matter halos is still a matter of intense debate. Present theoretical models suffer from large uncertainties and degeneracies, for example, between the fraction of accreting sources and their characteristic accretion rate. In recent work we showed that Active Galactic Nuclei (AGN) clustering represents a powerful tool to break degeneracies when analysed in terms of mean BH mass, and that AGN bias at fixed stellar mass is largely independent of most of the input parameters, such as the AGN duty cycle and the mean scaling between BH mass and host galaxy stellar mass. In this paper we take advantage of our improved semi-empirical methodology and recent clustering data derived from large AGN samples at $z \sim 1.2$, demonstrate that the AGN bias as a function of host galaxy stellar mass is a crucial diagnostic of the BH--galaxy connection, and is highly dependent on the scatter around the BH mass--galaxy mass scaling relation and on the relative fraction of satellite and central active BHs. Current data at $z \sim 1.2$ favour relatively high values of AGN in satellites, pointing to a major role of disc instabilities in triggering AGN, unless a high minimum host halo mass is assumed. The data are not decisive on the magnitude/covariance of the BH-galaxy scatter at $z \sim 1.2$ and intermediate host masses $M_\mathrm{star} \lesssim 10^{11} \,\mathrm{M}_\odot$. However, future surveys like Euclid/LSST will be pivotal in shedding light on the BH--galaxy co-evolution.

We study the $f(R,T)$ cosmological models under the self-similarity hypothesis. We determine the exact form that each physical and geometrical quantity may take in order that the Field Equations (FE) admit exact self-similar solutions through the matter collineation approach. We study two models: the case$\ f(R,T)=f_{1}(R)+f_{2}(T)$ and the case $f(R,T)=f_{1}(R)+f_{2} (R)f_{3}(T)$. In each case, we state general theorems which determine completely the form of the unknown functions $f_{i}$ such that the field equations admit self-similar solutions. We also state some corollaries as limiting cases. These results are quite general and valid for any homogeneous self-similar metric$.$ In this way, we are able to generate new cosmological scenarios. As examples, we study two cases by finding exact solutions to these particular models.

We report the discovery of a new, chemically distinct population of relatively high-metallicity ([Fe/H] $> -0.7$) red giant stars with super-solar [N/Fe] ($\gtrsim +0.75$) identified within the bulge, disk, and halo of the Milky Way. This sample of stars was observed during the second phase of the Apache Point Observatory Galactic Evolution Experiment (APOGEE-2); the spectra of these stars are part of the seventeenth Data Release (DR 17) of the Sloan Digital Sky Survey. We hypothesize that this newly identified population was formed in a variety of progenitors, and are likely made up of either fully or partially destroyed metal-rich globular clusters, which we refer to as Globular Cluster Debris (GCD), identified by their unusual photospheric nitrogen abundances. It is likely that some of the GCD stars were probable members of the Gaia-Enceladus-Sausage accretion event, along with clusters formed in situ.

We study the impact of AGN ionised outflows on star formation in high-redshift AGN hosts, by combining NIR IFS observations, mapping the H$\alpha$ emission and [OIII] outflows, with matched-resolution observations of the rest-frame FIR emission. We present high-resolution ALMA Band 7 observations of eight X-ray selected AGN at z~2 from the SUPER sample, targeting the rest-frame ~260 um continuum at ~2 kpc (0.2'') resolution. We detected 6 out of 8 targets with S/N>10 in the ALMA maps, with continuum flux densities F = 0.27-2.58 mJy and FIR half-light radii Re = 0.8-2.1 kpc. The FIR Re of our sample are comparable to other AGN and star-forming galaxies at a similar redshift from the literature. However, we find that the mean FIR size in X-ray AGN (Re = 1.16+/- 0.11 kpc) is slightly smaller than in non-AGN (Re = 1.69+/-0.13 kpc). From SED fitting, we find that the main contribution to the 260 um flux density is dust heated by star formation, with < 4% contribution from AGN-heated dust and < 1% from synchrotron emission. The majority of our sample show different morphologies for the FIR (mostly due to reprocessed stellar emission) and the ionised gas emission (H$\alpha$ and [OIII], mostly due to AGN emission). This could be due to the different locations of dust and ionised gas, the different sources of the emission (stars and AGN), or the effect of dust obscuration. We are unable to identify any residual H$\alpha$ emission, above that dominated by AGN, that could be attributed to star formation. Under the assumption that the FIR emission is a reliable tracer of obscured star formation, we find that the obscured star formation activity in these AGN host galaxies is not clearly affected by the ionised outflows. However, we cannot rule out that star formation suppression is happening on smaller spatial scales than the ones we probe with our observations (< 2 kpc) or on different timescales.

We provide the first full K2 transiting exoplanet sample, using photometry from Campaigns 1-8 and 10-18, derived through an entirely automated procedure. This homogeneous planet candidate catalog is a crucial to perform a robust demographic analysis of transiting exoplanets with K2. We identify 747 unique planet candidates and 57 multi-planet systems. Of these candidates, 366 have not been previously identified, including one resonant multi-planet system and one system with two short-period gas giants. By automating the construction of this list, measurements of sample biases (completeness and reliability) can be quantified. We carried out a light curve-level injection/recovery test of artificial transit signals and found a maximum completeness of 61%, a consequence of the significant detrending required for K2 data analysis. Through this operation we attained measurements of the detection efficiency as a function of signal strength, enabling future population analysis using this sample. We assessed the reliability of our planet sample by testing our vetting software EDI-Vetter against inverted transit-free light curves. We estimate 91% of our planet candidates are real astrophysical signals, increasing up to 94% when limited to the FGKM dwarf stellar population. We also constrain the contamination rate from background eclipsing binaries to less than 5%. The presented catalog, along with the completeness and reliability measurements, enable robust exoplanet demographic studies to be carried out across the fields observed by the K2 mission for the first time.

The numerical solution of relativistic hydrodynamics equations in conservative form requires root-finding algorithms that invert the conservative-to-primitive variables map. These algorithms employ the equation of state of the fluid and can be computationally demanding for applications involving sophisticated microphysics models. This work explores the use of machine learning methods to speed up the recovery of primitives in relativistic hydrodynamics. Artificial neural networks are trained to replace either the interpolations of a tabulated equation of state or directly the conservative-to-primitive map. The application of these neural networks to simple benchmark problems show that both approaches improve over traditional root finders with tabular equation-of-state and multi-dimensional interpolations. In particular, the neural networks for the conservative-to-primitive map accelerate the variable recovery by more than an order of magnitude over standard methods while maintaining accuracy. Neural networks are thus an interesting option to improve the speed and robustness of relativistic hydrodynamics algorithms.

We examine the consequences of, and apply, the formalism developed in Terquem (2021) for calculating the rate $D_R$ at which energy is exchanged between fast tides and convection. In this previous work, $D_R$ (which is proportional to the gradient of the convective velocity) was assumed to be positive in order to dissipate the tidal energy. Here we argue that, even if energy is intermittently transferred from convection to the tides, it must ultimately return to the convective flow and transported efficiently to the stellar surface on the convective timescale. This is consistent with, but much less restrictive than, enforcing $D_R>0$. Our principle result is a calculation of the circularization timescale of late-type binaries, taking into account the full time evolution of the stellar structure. We find that circularization is very efficient during the PMS phase, inefficient during the MS, and once again efficient when the star approaches the RGB. These results are in much better agreement with observations than earlier theories. We also apply our formalism to hot Jupiters, and find that tidal dissipation in a Jupiter mass planet yields a circularization timescale of 1 Gyr for an orbital period of 3 d, also in good overall agreement with observations. The approach here is novel, and the apparent success of the theory in resolving longstanding timescale puzzles is compelling.

The reconstruction of the solar spectral irradiance (SSI) on various time scales is essential for the understanding of the Earth's climate response to the SSI variability. The driver of the SSI variability is understood to be the intensity contrast of magnetic features present on the Sun with respect to the largely non-magnetic quiet Sun. However, different spectral synthesis codes lead to diverging projections of SSI variability. In this study we compare three different radiative transfer codes and carry out a detailed analysis of their performance. We perform the spectral synthesis at the continuum wavelength of 665 nm with the Code for Solar Irradiance (COSI), and the Rybicki-Hummer (RH), and Max Planck University of Chicago Radiative MHD (MURaM) codes for three 3D MHD simulations snapshots, a non-magnetic case, and MHD simulations with 100 G, and 200 G magnetic field strength. We determine the intensity distributions, the intensity differences and ratios for the spectral synthesis codes. We identify that the largest discrepancies originate in the intergranular lanes where the most field concentration occurs. Overall, the applied radiative transfer codes give consistent intensity distributions. Also, the intensity variation as a function of magnetic field strength for the particular 100 G and 200 G snapshots agree within the 2-3% range.

The Diffuse Ionized Gas (DIG) contributes to the nebular emission of galaxies, resulting in emission line flux ratios that can be significantly different from those produced by HII regions. Comparing the emission of [SII]6717,31 between pointed observations of HII regions in nearby galaxies and integrated spectra of more distant galaxies, it has been recently claimed that the DIG can also deeply affect the emission of bright, star-forming galaxies, and that a large correction must be applied to observed line ratios to recover the genuine contribution from HII regions. Here we show instead that the effect of DIG on the integrated spectra of star-forming galaxies is lower than assumed in previous work. Indeed, aperture effects on the spectroscopy of nearby HII regions are largely responsible for the observed difference: when spectra of local HII regions are extracted using large enough apertures while still avoiding the DIG, the observed line ratios are the same as in more distant galaxies. This result is highly relevant for the use of strong-line methods to measure metallicity.

Millimeter continuum imaging of protoplanetary disks reveals the distribution of solid particles and the presence of substructures (gaps and rings) beyond 5-10 au, while infrared (IR) spectra provide access to abundances of gaseous species at smaller disk radii. Building on recent observational findings of an anti-correlation between the inner disk water luminosity and outer dust disk radius, we aim here at investigating the dynamics of icy solids that drift from the outer disk and sublimate their ice inside the snow line, enriching the water vapor that is observed in the IR. We use a volatile-inclusive disk evolution model to explore a range of conditions (gap location, particle size, disk mass, and alpha-viscosity) under which gaps in the outer disk efficiently block the inward drift of icy solids. We find that inner-disk vapor enrichment is highly sensitive to the location of a disk gap, yielding for each particle size a radial "sweet spot" that reduces the inner-disk vapor enrichment to a minimum. For pebbles of 1-10 mm in size, which carry the most mass, this sweet spot is at 7-15 au, suggesting that inner gaps may have a key role in reducing ice delivery to the inner disk and may not allow the formation of Earths and super-Earths. This highlights the importance of observationally determining the presence and properties of inner gaps in disks. Finally, we argue that the inner water vapor abundance can be used as a proxy for estimating the pebble drift efficiency and mass-flux entering the inner disk.

We generalize the Laplace resonance among three satellites, S1, S2 , and S3, by considering different ratios of the mean-longitude variations. These resonances, which we call Laplace-like, are classified as first order in the cases of the 2:1&2:1, 3:2&3:2, and 2:1&3:2 resonances, second order in the case of the 3:1&3:1 resonance, and mixed order in the case of the 2:1&3:1 resonance. We consider a model that includes the gravitational interaction with the central body together with the effect due to its oblateness, the mutual gravitational influence of the satellites S1, S2, and S3 and the secular gravitational effect of a fourth satellite S 4 , which plays the role of Callisto in the Galilean system. In addition, we consider the dissipative effect due to the tidal torque between the inner satellite and the central body. We investigate these Laplace-like resonances by studying different aspects: (i) we study the survival of the resonances when the dissipation is included, taking two different expressions for the dissipative effect in the case of a fast- or a slowly rotating central body, (ii) we investigate the behavior of the Laplace-like resonances when some parameters are varied, specifically, the oblateness coefficient, the semimajor axes, and the eccentricities of the satellites, (iii) we analyze the linear stability of first-order resonances for different values of the parameters, and (iv) we also include the full gravitational interaction with S 4 to analyze its possible capture into resonance. The results show a marked difference between first-, second-, and mixed-order resonances, which might find applications when the evolutionary history of the satellites in the Solar System are studied, and also in possible actual configurations of extrasolar planetary systems.

In a purely cold dark matter universe, the initial matter power spectrum and its subsequent gravitational growth contain no special mass- or time-scales, and so neither do the emergent population statistics of internal dark matter (DM) halo properties. Using 1.5 million halos from three IllustrisTNG realizations of a LambdaCDM universe, we show that galaxy formation physics drives non-monotonic features ("wiggles") into DM property statistics across six decades in halo mass, from dwarf galaxies to galaxy clusters. We characterize these features by extracting the halo mass-dependent statistics of five DM halo properties -- velocity dispersion, NFW concentration, density- and velocity-space shapes, and formation time -- using kernel-localized linear regression (KLLR). Comparing precise estimates of normalizations, slopes, and covariances between realizations with and without galaxy formation, we find systematic deviations across all mass-scales, with maximum deviations of 25% at the Milky-Way mass of 1e12 Msun. The mass-dependence of the wiggles is set by the interplay between different cooling and feedback mechanisms, and we discuss its observational implications. The property covariances depend strongly on halo mass and physics treatment, but the correlations are mostly robust. Using multivariate KLLR and interpretable machine learning, we show the halo concentration and velocity-space shape are principal contributors, at different mass, to the velocity dispersion variance. Statistics of mass accretion rate and DM surface pressure energy are provided in an appendix. We publicly release halo property catalogs and KLLR parameters for the TNG runs at twenty epochs up to z = 12.

The majority of exoplanets found to date have been discovered via the transit method, and transmission spectroscopy represents the primary method of studying these distant worlds. Currently, in-depth atmospheric characterization of transiting exoplanets entails the use of spectrographs on large telescopes, requiring significant observing time to study each planet. Previous studies have demonstrated trends for solar system worlds using color-color photometry of reflectance spectra, as well as trends within transmission spectra for hot Jupiters. Building on these concepts, we have investigated the use of transmission color photometric analysis for efficient, coarse categorization of exoplanets and for assessing the nature of these worlds, with a focus on resolving the bulk composition degeneracy to aid in discriminating super-Earths and sub-Neptunes. We present our methodology and first results, including spectrum models, model comparison frameworks, and wave band selection criteria. We present our results for different transmission "color" metrics, filter selection methods, and numbers of filters. Assuming noise-free spectra of isothermal atmospheres in chemical equilibrium, with our pipeline, we are able to constrain atmospheric mean molecular weight in order to distinguish between super-Earth and sub-Neptune atmospheres with >90$\%$ overall accuracy using as few as two specific low-resolution filter combinations. We also found that increasing the number of filters does not substantially impact this performance. This method could allow for broad characterization of large numbers of planets much more efficiently than current methods permit, enabling population and system-level studies. Additionally, data collected via this method could inform follow-up observing time by large telescopes for more detailed studies of worlds of interest.

Star formation induced by a spiral shock wave, which in turn is generated by a spiral density wave, produces an azimuthal age gradient across the spiral arm, which has opposite signs on either side of the corotational resonance. An analysis of the spatial separation between young star clusters and nearby HII regions made it possible to determine the position of the corotation radius in the studied galaxies. Fourier analysis of the gas velocity field in the same galaxies independently confirmed the corotation radius estimates obtained by the morphological method presented here.

We used archival data from the \emph{Chandra X-ray Observatory} (\emph{Chandra}) and the \emph{Hubble Space Telescope}, to identify 334 candidate X-ray binary (XRB) systems and their potential optical counterparts in the interacting galaxy pair NGC 5194/5195 (M51). We present the catalog and data analysis of X-ray and optical properties for those sources, from the deep $892$ ks \emph{Chandra} observations, along with the magnitudes of candidate optical sources as measured in the $8.16$ ks \emph{HST} observations. The X-ray luminosity function of the X-ray sources above a few times $10^{36}\, {\rm erg\,s^{-1}}$ follows a power law $N(>L_{X,b})\propto L_{X,b}^{1-\alpha}$ with $\alpha=1.65\pm0.03$. Aproximately 80\% of sources are variable over a 30 day window. Nearly half of the X-ray sources (173/334) have an optical counterparts within $0{\mbox{$.\!\!^{\prime\prime}$}}5$.

We report results obtained from the study of 12 thermonuclear X-ray bursts in 6 AstroSat observations of a neutron star X-ray binary and well-known X-ray burster, 4U 1636$-$536. Burst oscillations at $\sim$581 Hz are observed with 4$-$5$\sigma$ confidence in three of these X-ray bursts. The rising phase burst oscillations show a decreasing trend of the fractional rms amplitude at 3$\sigma$ confidence,by far the strongest evidence of thermonuclear flame spreading observed with AstroSat. During the initial 0.25 second of the rise a very high value (34.0$\pm$6.7%) is observed. The concave shape of the fractional amplitude profile provides a strong evidence of latitude-dependent flame speeds, possibly due to the effects of the Coriolis force. We observe decay phase oscillations with amplitudes comparable to that observed during the rising phase, plausibly due to the combined effect of both surface modes as well as the cooling wake. The Doppler shifts due to the rapid rotation of the neutron star might cause hard pulses to precede the soft pulses, resulting in a soft lag. The distance to the source estimated using the PRE bursts is consistent with the known value of $\sim$6 kpc.

Based on a 14-moment closure for non-resistive (general-) relativistic viscous plasmas, we describe a new numerical scheme that is able to handle all first-order dissipative effects (heat conduction, bulk and shear viscosities), as well the anisotropies induced by the presence of magnetic fields. The latter is parameterized in terms of a thermal gyrofrequency or, equivalently, a thermal Larmor radius and allows to correctly capture the thermal Hall effect. By solving an extended Israel-Stewart-like system for the dissipative quantities that enforces algebraic constraints via stiff-relaxation, we are able to cast all first-order dissipative terms in flux-divergence form. This allows us to apply traditional high-resolution shock capturing methods to the equations, making the system suitable for the numerical study of highly turbulent flows. We present several numerical tests to assess the robustness of our numerical scheme in flat spacetime. The 14-moment closure can seamlessly interpolate between the highly collisional limit found in neutron star mergers, and the highly anisotropic limit of relativistic Braginskii magnetohydrodynamics appropriate for weakly collisional plasmas in black-hole accretion problems. We believe that this new formulation and numerical scheme will be useful for a broad class of relativistic magnetized flows.

We use the idea of repeat measurement to determine the mean value and error of light curve solution parameters of contact binaries. Our simulation is realized by the Monte Carlo algorithm and Wilson-Devinney code. This method can obtain the systematic and random error simultaneously. Within our 48 models, the systematic errors are smaller than the random errors in most case. According to the numerical calculations, it is found that the relative errors of photometric mass ratios are less than 1% for totally eclipsing contact binaries, while they are generally between 10% and 20% for partly eclipsing ones. The effect of third light on the errors of photometric solution is also investigated. With a third light, these errors are close to 10% for totally eclipsing contact binaries. Specially, it is better to set the third light to zero in flux if that light is very faint (e.g., less than 1% contribution in luminosity), because such faint third light will bring big errors to the light curve solutions.

HD 98800 is a young ($\sim10$ Myr old) and nearby ($\sim45$ pc) quadruple system, composed of two spectroscopic binaries orbiting around each other (AaAb and BaBb), with a gas-rich disk in polar configuration around BaBb. While the orbital parameters of BaBb and AB are relatively well constrained, this is not the case for AaAb. A full characterisation of this quadruple system can provide insights on the formation of such a complex system. The goal of this work is to determine the orbit of the AaAb subsystem and refine the orbital solution of BaBb using multi-epoch interferometric observations with the VLTI/PIONIER and radial velocities. The PIONIER observations provide relative astrometric positions and flux ratios for both AaAa and BaBb subsystems. Combining the astrometric points with radial velocity measurements, we determine the orbital parameters of both subsystems. We refined the orbital solution of BaBb and derived, for the first time, the full orbital solution of AaAb. We confirmed the polar configuration of the circumbinary disk around BaBb. From our solutions, we also inferred the dynamical masses of AaAb ($M_{Aa} = 0.93 \pm 0.09$ and $M_{Ab} = 0.29 \pm 0.02 M$_{\odot}$). We also revisited the parameters of the AB outer orbit. Using the N-body simulation, we show that the system should be dynamically stable over thousands of orbital periods and that it made preliminary predictions for the transit of the disk in front of AaAb which is estimated to start around 2026. We discuss the lack of a disk around AaAb, which can be explained by the larger X-ray luminosity of AaAb, promoting faster photo-evaporation of the disk. High-resolution infrared spectroscopic observations would provide radial velocities of Aa and Ab (blended lines in contemporary observations), which would allow us to calculate the dynamical masses of Aa and Ab independently of the parallax of BaBb.

We present SDO/AIA observation of three types of fast-mode propagating magnetosonic waves in a GOES C3.0 flare on 2013 April 23, which was accompanied by a prominence eruption and a broad coronal mass ejection (CME). During the fast rising phase of the prominence, a large-scale dome-shaped extreme ultraviolet (EUV) wave firstly formed ahead of the CME bubble and propagated at a speed of about 430 km/s in the CME's lateral direction. One can identify the separation process of the EUV wave from the CME bubble. The reflection effect of the on-disk counterpart of this EUV wave was also observed when it interacted with a remote active region. Six minutes after the first appearance of the EUV wave, a large-scale quasi-periodic EUV train with a period of about 120 seconds appeared inside the CME bubble, which emanated from the flare epicenter and propagated outward at an average speed up to 1100 km/s. In addition, another narrow quasi-periodic EUV wave train was observed along a closed-loop system connecting two adjacent active regions, which also emanated from the flare epicenter, propagated at a speed of about475 km/s and with a period of about 110 seconds. We propose that all the observed waves are fast-mode magnetosonic waves, in which the large-scale dome-shaped EUV wave ahead of the CME bubble was driven by the expansion of the CME bubble, while the large-scale quasi-periodic EUV train within the CME bubble and the narrow quasi-periodic EUV wave train along the closed-loop system were excited by the intermittent energy-releasing process in the flare. Coronal seismology application and energy carried by the waves are also estimated based on the measured wave parameters.

We present a host morphological study of 1,265 far-infrared galaxies (FIRGs) and sub-millimeter galaxies (SMGs) in the Cosmic Evolution Survey field using the F160W and F814W images obtained by the Hubble Space Telescope. The FIRGs and the SMGs are selected from the Herschel Multi-tiered Extragalactic Survey and the SCUBA-2 Cosmology Legacy Survey, respectively. Their precise locations are based on the interferometry data from the Atacama Large Millimeter/submillimeter Array and the Very Large Array. The vast majority of these objects are at $0.1\lesssim z\lesssim 3$. While the SMGs do not constitute a subset of the FIRGs in our selection due to the signal-to-noise ratio thresholds, SMGs can be regarded as the population at the high-redshift tail of FIRGs. Most of our FIRGs/SMGs have total infrared luminosity ($L_{\rm IR}$) in the regimes of luminous and ultra-luminous infrared galaxies (LIRGs, $L_{\rm IR} = 10^{11-12}L_\odot$; ULIRGs, $L_{\rm IR}>10^{12}L_\odot$). The hosts of the SMG ULIRGs, FIRG ULIRGs and FIRG LIRGs are of sufficient numbers to allow for detailed analysis, and they are only modestly different in their stellar masses. Their morphological types are predominantly disk galaxies (type "D") and irregular/interacting systems (type "Irr/Int"). There is a morphological transition at $z\approx 1.25$ for the FIRG ULIRG hosts, above which the Irr/Int galaxies dominate and below which the D and the Irr/Int galaxies have nearly the same contributions. The SMG ULIRG hosts seem to experience a similar transition. This suggests a shift in the relative importance of galaxy mergers/interactions versus secular gas accretions in "normal" disk galaxies as the possible triggering mechanisms of ULIRGs. The FIRG LIRG hosts are predominantly D galaxies over $z=0.25-1.25$ where they are of sufficient statistics.

Based on observations with the $Chandra$ X-ray Observatory, we present the latest spectral evolution of the X-ray remnant of SN 1987A (SNR 1987A). We present a high-resolution spectroscopic analysis using our new deep ($\sim$312 ks) $Chandra$ HETG observation taken in March 2018, as well as archival $Chandra$ gratings spectroscopic data taken in 2004, 2007, and 2011 with similarly deep exposures ($\sim$170 - 350 ks). We perform detailed spectral model fits to quantify changing plasma conditions over the last 14 years. Recent changes in electron temperatures and volume emission measures suggest that the shocks moving through the inner ring have started interacting with less dense circumstellar material, probably beyond the inner ring. We find significant changes in the X-ray line flux ratios (among H- and He-like Si and Mg ions) in 2018, consistent with changes in the thermal conditions of the X-ray emitting plasma that we infer based on the broadband spectral analysis. Post-shock electron temperatures suggested by line flux ratios are in the range $\sim$0.8 - 2.5 keV as of 2018. We do not yet observe any evidence of substantial abundance enhancement, suggesting that the X-ray emission component from the reverse-shocked metal-rich ejecta is not yet significant in the observed X-ray spectrum.

In this work, a mass-effective temperature-surface gravity relation (MTGR) is developed for main sequence stars in the range of 6400 K < $T_{\rm eff}$ < 20000 K with log$g$ > 3.44. The MTGR allows the simple estimation of the masses of stars from their effective temperatures and surface gravities. It can be used for solar metallicity and can be rescaled for any metallicity within -1.00 < [Fe/H] < 0.7. The effect of alpha-enhanced compositions can also be considered with the help of correction terms. It is aimed to develop an MTGR that can estimate the masses of main-sequence stars from their atmospheric parameters. One advantage of an MTGR over the classical mass-luminosity relations is that its mass estimation is based on parameters that can be obtained by purely spectroscopic methods and, therefore, the interstellar extinction or reddening do not have to be known. The use of surface gravity ($g$) also relates an MTGR with stellar evolution and provides a more reliable mass estimation. A synthetical MTGR is obtained from theoretical isochrones using a Levenberg-Marquardt chi-square minimization algorithm. The validity of the MTGR is then checked by testing over 278 binary components with precise absolute masses. Very good agreement has been obtained between the absolute masses of 278 binary star components and their masses estimated from the MTGR. A mathematical expression is also given to calculate the propagated uncertainties of the MTGR masses. For the typical uncertainties in atmospheric parameters and metallicity, the typical uncertainties in the masses estimated from the MTGR mostly remain around 5-9%. The fact that this uncertainty level is only on average about three times as large as that of the absolute masses indicates that the MTGR is a very powerful tool for stellar mass estimation. A computer code, mtgr.pro, written in GDL/IDL is also provided for the relation.

Understanding the mechanisms of accretion-ejection during X-ray binaries outbursts has been a problem for several decades. It is still not clear yet what controls the spectral evolution of these objects from the hard to the soft states and then back to the hard states at the end of the outburst, tracing the well-known hysteresis cycle in the hardness-intensity diagram. Moreover, the link between the spectral states and the presence/absence of radio emission is still highly debated. In a series of papers, we developed a model composed of a truncated outer standard accretion disk (SAD, from the solution of Shakura and Sunyaev) and an inner jet emitting disk (JED). In this paradigm, the JED plays the role of the hot corona while simultaneously explaining the presence of a radio jet. Our goal is to apply for the first time direct fitting procedures of the JED-SAD model to the hard states of four outbursts of GX 339-4 observed during the 2000-2010 decade by RXTE, combined with simultaneous or quasi simultaneous ATCA observations. We built JED-SAD model tables usable in Xspec as well as a reflection model table based on the Xillver model. We apply our model to the 452 hard state observations obtained with RXTE/PCA. We were able to correctly fit the X-ray spectra and simultaneously reproduce the radio flux. We show that the functional dependency of the radio emission on the model parameters (mainly the accretion rate and the transition radius between the JED and the SAD) is similar between all the rising phases of the different outbursts of GX 339-4. But it is significantly different from the functional dependency obtained in the decaying phases. This result strongly suggests a change in the radiative and/or dynamical properties of the ejection between the beginning and the end of the outburst. We discuss possible scenarios that could explain these differences.

Generating pre-initial conditions (or particle loads) is the very first step to set up a cosmological N-body simulation. In this work, we revisit the numerical convergence of pre-initial conditions on dark matter halo properties using a set of simulations which only differs in initial particle loads, i.e. grid, glass, and the newly introduced capacity constrained Voronoi tessellation (CCVT). We find that the median halo properties agree fairly well (i.e. within a convergence level of a few per cent) among simulations running from different initial loads. We also notice that for some individual haloes cross-matched among different simulations, the relative difference of their properties sometimes can be several tens of per cent. By looking at the evolution history of these poorly converged haloes, we find that they are usually merging haloes or haloes have experienced recent merger events, and their merging processes in different simulations are out-of-sync, making the convergence of halo properties become poor temporarily. We show that, comparing to the simulation starting with an anisotropic grid load, the simulation with an isotropic CCVT load converges slightly better to the simulation with a glass load, which is also isotropic. Among simulations with different pre-initial conditions, haloes in higher density environments tend to have their properties converged slightly better. Our results confirm that CCVT loads behave as well as the widely used grid and glass loads at small scales, and for the first time we quantify the convergence of two independent isotropic particle loads (i.e. glass and CCVT) on halo properties.

Prominence eruption is closely related to coronal mass ejections and is an important topic in solar physics. Spectroscopic observation is an effective way to explore the plasma properties, but the spectral observations of eruptive prominences are rare. In this paper we will introduce an eruptive polar crown prominence with spectral observations from the Interface Region Imaging Spectrograph (IRIS), and try to explain some phenomena that are rarely reported in previous works. The eruptive prominence experiences a slow-rise and fast-rise phase, while the line-of-sight motions of the prominence plasma could be divided into three periods: two hours before the fast-rise phase, opposite Doppler shifts are found at the two sides of the prominence axis;then, red shifts dominate the prominence gradually; in the fast-rise phase, the prominence gets to be blue-shifted. During the second period, a faint component appears in Mg II k window with a narrow line width and a large red shift. A faint region is also found in AIA 304-angstrom images along the prominence spine, and the faint region gets darker during the expansion of the spine. We propose that the opposite Doppler shifts in the first period is a feature of the polar crown prominence that we studied. The red shifts in the second period is possibly due to mass drainage during the elevation of the prominence spine, which could accelerate the eruption in return. The blue shifts in the third period is due to that the prominence erupts toward the observer. We suggest that the faint component appears due to the decreasing of the plasma density, and the latter results from the expansion of the prominence spine.

The cut-off mass ratio is under debate for contact binaries. In this paper, we present the investigation of two contact binaries with mass ratios close to the low mass ratio limit. It is found that the mass ratios of VSX J082700.8+462850 (hereafter J082700) and 1SWASP J132829.37+555246.1 (hereafter J132829) are both less than 0.1 ($q\sim0.055$ for J082700, and $q\sim0.089$ for J132829). J082700 is a shallow contact binary with a contact degree of $\sim$19\%, and J132829 is a deep contact system with a fillout factor of $\sim$70\%. The $O-C$ diagram analysis indicated that both the two systems manifest long-term period decrease. In addition, J082700 exhibits a cyclic modulation which is more likely resulted from Applegate mechanism. In order to explore the properties of extremely low mass ratio contact binaries (ELMRCBs), we carried out a statistical analysis on contact binaries with mass ratios of $q\lesssim0.1$ and discovered that the values of $J_{spin}/J_{orb}$ of three systems are greater than 1/3. Two possible explanations can interpret this phenomenon. One is that some physical processes, unknown to date, are not considered when Hut presented the dynamically instability criterion. The other is that the dimensionless gyration radius ($k$) should be smaller than the value we used ($k^2=0.06$). We also found that the formation of ELMRCBs possibly has two channels. The study of evolutionary states of ELMRCBs reveals that their evolutionary states are similar with those of normal W UMa contact binaries.

Due to the increasing number of observations Type Ia supernovae are nowadays regarded as a heterogeneous class of objects consisting of several subclasses. One of the largest of these is the class of Type Iax supernovae (SNe Iax) which have been suggested to originate from pure deflagrations in CO Chandrasekhar-mass white dwarfs (WDs). Although a few deflagration studies have been carried out, the full diversity of the class is not captured yet. We therefore present a parameter study of single-spot ignited deflagrations with varying ignition locations, central densities, metallicities and compositions. We also explore a rigidly rotating progenitor and carry out 3D hydrodynamic simulations, nuclear network calculations and radiative transfer. The new models extend the range in brightness covered by previous studies to the lower end. Our explosions produce $^{56}$Ni masses from $5.8 \times 10^{-3}$ to $9.2 \times 10^{-2}\,M_\odot$. In spite of the wide exploration of the parameter space the main characteristics of the models are primarily driven by the mass of $^{56}$Ni. Secondary parameters have too little impact to explain the observed trend among faint SNe~Iax. We report kick velocities of the bound explosion remnants from $6.9$ to $369.8\,$km$\,s^{-1}$. The wide exploration of the parameter space and viewing-angle effects in the radiative transfer lead to a significant spread in the synthetic observables. The trends towards the faint end of the class are, however, not reproduced. This motivates a quantification of the systematic uncertainties in the modeling procedure and the influence of the $^{56}$Ni-rich bound remnant. While the pure deflagration scenario remains a favorable explanation for bright and intermediate luminosity SNe~Iax, the possibility that SNe~Iax do not consist of a single explosion scenario needs to be considered.

We use a semi-analytic galaxy formation model to study the co-evolution of supermassive black holes (SMBHs) with their host galaxies. Although the coalescence of SMBHs is not important, the quasar-mode accretion induced by mergers plays a dominant role in the growth of SMBHs. Mergers play a more important role in the growth of SMBH host galaxies than in the SMBH growth. It is the combined contribution from quasar mode accretion and mergers to the SMBH growth and the combined contribution from starburst and mergers to their host galaxy growth that determine the observed scaling relation between the SMBH masses and their host galaxy masses. We also find that mergers are more important in the growth of SMBH host galaxies compared to normal galaxies which share the same stellar mass range as the SMBH host galaxies.

Context. A proper estimate of the chromospheric magnetic fields is believed to improve modelling of both active region and coronal mass ejection evolution. Aims. We investigate the similarity between the chromospheric magnetic field inferred from observations and the field obtained from a magnetohydrostatic (MHS) extrapolation. Methods. Based Fe i 6173 {\AA} and Ca ii 8542 {\AA} observations of NOAA active region 12723, we employed the spatially-regularised weak-field approximation (WFA) to derive the vector magnetic field in the chromosphere from Ca ii, as well as non-LTE inversions of Fe i and Ca ii to infer a model atmosphere for selected regions. Milne-Eddington inversions of Fe i serve as photospheric boundary for the MHS model that delivers the three-dimensional field, gas pressure and density. Results. For the line-of-sight component, the MHS chromospheric field generally agrees with the non-LTE inversions and WFA, but tends to be weaker than those when larger in magnitude than 300 G. The observationally inferred transverse component is stronger, especially in magnetically weaker regions, yet the qualitative distribution with height is similar to the MHS results. For either field component the MHS chromospheric field lacks the fine structure derived from the inversions. Furthermore, the MHS model does not recover the magnetic imprint from a set of high fibrils connecting the main polarities. Conclusions. The MHS extrapolation and WFA provide a qualitatively similar chromospheric field, where the azimuth of the former is better aligned with Ca ii 8542 {\AA} fibrils than that of the WFA, especially outside strong-field concentrations. The amount of structure as well as the transverse field strengths are underestimated by the MHS extrapolation. This underscores the importance of considering a chromospheric magnetic field constraint in data-driven modelling of active regions.

We study the statistics of various large-scale structure tracers in gravity-only cosmological simulations including baryons and cold dark matter (CDM) initialized with two different transfer functions, and simulated as two distinct fluids. This allows us to study the impact of baryon-CDM relative perturbations on these statistics. In particular, we focus on the statistics of cosmic voids, as well as on the matter and halo real-space 2-point correlation function and baryon acoustic oscillations (BAO) peak. We find that the void size function is affected at the 1-2% level at maximum, and that the impact is more important at higher redshift, while the void density profile and void bias are roughly unaffected. We do not detect a sizeable impact of relative baryon-CDM perturbations on the real-space correlation functions of matter and halos or the BAO peak, which is in line with results from previous works. Our results imply that it would be hard to use voids or real-space correlation functions to constrain baryon-CDM relative perturbations, but also that we might not have to include them in models for the analysis of future cosmological surveys data.

Acoustic waves excited in the photosphere and below might play an integral part in the heating of the solar chromosphere and corona. However, it is yet not fully clear how much of the initially acoustic wave flux reaches the corona and in what form. We investigate the wave propagation, damping, transmission, and conversion in the lower layers of the solar atmosphere using 3D numerical MHD simulations. A model of a gravitationally stratified expanding straight coronal loop, stretching from photosphere to photosphere, is perturbed at one footpoint by an acoustic driver with a period of 370 seconds. For this period acoustic cutoff regions are present below the transition region (TR). About 2% of the initial energy from the driver reach the corona. The shape of the cutoff regions and the height of the TR show a highly dynamic behavior. Taking only the driven waves into account, the waves have a propagating nature below and above the cutoff region, but are standing and evanescent within the cutoff region. Studying the driven waves together with the background motions in the model reveals standing waves between the cutoff region and the TR. These standing waves cause an oscillation of the TR height. In addition, fast or leaky sausage body-like waves might have been excited close to the base of the loop. These waves then possibly convert to fast or leaky sausage surface-like waves at the top of the main cutoff region, followed by a conversion to slow sausage body-like waves around the TR.

Dust has been detected in high-redshift ($z>5$) galaxies but its origin is still being debated. Dust production in high-redshift galaxies could be dominated by stellar production or by accretion (dust growth) in the interstellar medium. Previous studies have shown that these two dust sources predict different grain size distributions, which lead to significantly different extinction curves. In this paper, we investigate how the difference in the extinction curves affects the dust attenuation properties of galaxies by performing radiative transfer calculations. To examine the major effects of the dust--stars distribution geometry, we adopt two representative cases in spherical symmetry: the well-mixed geometry (stars and dust are homogeneously mixed) and the two-layer geometry (young stars are more concentrated in the centre). In both cases, we confirm that the attenuation curve can be drastically steepened by scattering and by different optical depths between young and old stellar populations, and can be flattened by the existence of unobscured stellar populations. We can reproduce similar attenuation curves even with very different extinction curves. Thus, we conclude that it is difficult to distinguish the dust sources only with attenuation curves. However, if we include information on dust emission and plot the IRX (infrared excess)--$\beta$ (ultraviolet spectral slope) relation, different dust sources predict different positions in the IRX--$\beta$ diagram. A larger $\beta$ is preferred under a similar IRX if dust growth is the dominant dust source.

In this work we analyse in detail the possibility of using small and intermediate-scale gravitational wave anisotropies to constrain the inflationary particle content. First, we develop a phenomenological approach focusing on anisotropies generated by primordial tensor-tensor-scalar and purely gravitational non-Gaussianities. We highlight the quantities that play a key role in determining the detectability of the signal. To amplify the power of anisotropies as a probe of early universe physics, we consider cross-correlations with CMB temperature anisotropies. We assess the size of the signal from inflationary interactions against so-called induced anisotropies. In order to arrive at realistic estimates, we obtain the projected constraints on the non-linear primordial parameter $F_{\rm NL}$ for several upcoming gravitational wave probes in the presence of the astrophysical gravitational wave background. We further illustrate our findings by considering a concrete inflationary realisation and use it to underscore a few subtleties in the phenomenological analysis.

We measure the rate of environmentally-driven star formation quenching in galaxies at $z\sim 1$, using eleven massive ($M\approx 2\times10^{14}\,\mathrm{M}_\odot$) galaxy clusters spanning a redshift range $1.0<z<1.4$ from the GOGREEN sample. We identify three different types of transition galaxies: "green valley" (GV) galaxies identified from their rest-frame $(NUV-V)$ and $(V-J)$ colours; "blue quiescent" (BQ) galaxies, found at the blue end of the quiescent sequence in $(U-V)$ and $(V-J)$ colour; and spectroscopic post-starburst (PSB) galaxies. We measure the abundance of these galaxies as a function of stellar mass and environment. For high stellar mass galaxies ($\log{M/\mathrm{M}_\odot}>10.5$) we do not find any significant excess of transition galaxies in clusters, relative to a comparison field sample at the same redshift. It is likely that such galaxies were quenched prior to their accretion in the cluster, in group, filament or protocluster environments. For lower stellar mass galaxies ($9.5<\log{M/\mathrm{M}_\odot}<10.5$) there is a small but significant excess of transition galaxies in clusters, accounting for an additional $\sim 5-10$ per cent of the population compared with the field. We show that our data are consistent with a scenario in which 20--30 per cent of low-mass, star-forming galaxies in clusters are environmentally quenched every Gyr, and that this rate slowly declines from $z=1$ to $z=0$. While environmental quenching of these galaxies may include a long delay time during which star formation declines slowly, in most cases this must end with a rapid ($\tau<1$ Gyr) decline in star formation rate.

Earth-skimming neutrinos are those which travel through the Earth's crust at a shallow angle. For Ultra-High-Energy (E > 1 PeV; UHE) earth-skimming tau neutrinos, there is a high-probability that the tau lepton created by a neutrino-Earth interaction will emerge from the ground before it decays. When this happens, the decaying tau particle initiates an air shower of relativistic sub-atomic particles which emit Cherenkov radiation. To observe this Cherenkov radiation, we propose the Trinity Observatory. Using a novel optical structure design, pointing at the horizon, Trinity will observe the Cherenkov radiation from upward-going neutrino-induced air showers. Being sensitive to neutrinos in the 1-10,000 PeV energy range, Trinity's expected sensitivity will have a unique role to play filling the gap between the observed astrophysical neutrinos observed by IceCube and the expected sensitivity of radio UHE neutrino detectors.

We measured the equivalent width of a large set of near-infrared (NIR, 0.8--2.4$ \ \mu$m) line-strength indices in the XShooter medium-resolution spectra of the central regions of 14 galaxies. We found that two aluminum indices Al at 1.31 $\mu$m and Al1 at 1.67 $\mu$m and the two CO indices CO1 at 1.56 $\mu$m and CO4 at 1.64 $\mu$m are tightly correlated with the velocity dispersion. Moreover, the NIR Al and CO1 indices show strong correlations with the optical Mg2 and Mgb indices, which are usually adopted as $\alpha$/Fe-enhancement diagnostics. The molecular FeH1 index at 1.58 $\mu$m tightly correlates with the optical <Fe> and [MgFe]' indices, which are used as total metallicity diagnostics. The NIR Pa$\beta$ index at 1.28 $\mu$m has a behaviour similar to the optical H$\beta$ index, which is a diagnostic of mean age. We defined two new composite indices, <Al> and [AlFeH], as possible candidates to be used as NIR diagnostics of total metallicity and $\alpha$/Fe enhancement. The NIR <Al> index has a strong correlation with the optical Mg2 and Mgb indices, while the [AlFeH] index is tightly correlated with the optical <Fe> and [MgFe]' indices. The distribution of the data points in the NIR Pa$\beta$-<Al> and Pa$\beta$-[AlFeH] diagrams mimic that in the optical [MgFe]'-H$\beta$ and the Mgb-<Fe> diagrams, which are widely used to constraint the properties of the unresolved stellar populations. We concluded that some NIR line-strength indices could be useful in studying stellar populations as well as in fine-tuning stellar population models.

Measuring the Cosmic X-ray Background (CXB) is a key to understand the Active Galactic Nuclei population, their absorption distribution and their average spectra. However, hard X-ray instruments suffer from time-dependent backgrounds and cross-calibration issues. The uncertainty of the CXB normalization remain of the order of 20%. To obtain a more accurate measurement, the Monitor Vsego Neba (MVN) instrument was built in Russia but not yet launched to the ISS (arXiv:1410.3284). We follow the same ideas to develop a CXB detector made of four collimated spectrometers with a rotating obturator on top. The collimators block off-axis photons below 100 keV and the obturator modulates on-axis photons allowing to separate the CXB from the instrumental background. Our spectrometers are made of 20 mm thick CeBr$_{3}$ crystals on top of a SiPM array. One tube features a $\sim$20 cm$^2$ effective area and more energy coverage than MVN, leading to a CXB count rate improved by a factor of $\sim$10 and a statistical uncertainty $\sim$0.5% on the CXB flux. A prototype is being built and we are seeking for a launch opportunity.

POLAR is a dedicated Gamma-Ray Burst polarimeter making use of Compton-scattering which took data from the second Chinese spacelab, the Tiangong-2 from September 2016 to April 2017. It has a wide Field of View of $\sim6$ steradians and an effective area of $\sim400\ cm^2$ at 300 keV. These features make it one of the most sensitive instruments in its energy range (15-500 keV), and therefore capable of almost continuously monitoring persistent sources such as pulsars. Significant folded pulsation from both PSR B0531+21 (the Crab Pulsar) and PSR B1509-58 has been observed. Observations of the Crab Pulsar with POLAR have previously been used for phase-resolved spectroscopy of the Crab Pulsar to calibrate the instrumental responses of POLAR. In this work, we investigate a polarimetric joint-fitting method for observations of the Crab Pulsar with POLAR. Unlike a GRB observation with POLAR, the observations of the Crab Pulsar are complicated by multiple observational datasets during which the polarization plane rotates as well. So before fitting, we have to correct the modulation curves under different datasets, by taking into account the rotations of the Crab Pulsar's relative position in the detctor's local coordinate, and the changes of detector response in different datasets. Despite these difficulties and the low signal to background for such sources constraining, polarization measurements were possible with the POLAR data. We will present the methodology briefly, which could be applied to any wide FoV polarimeter, and polarization results of the Crab pulsar with POLAR. Finally, the inferred ability of pulsar detection with POLAR-2 (the successor of POLAR) will also be discussed.

Context. Wide-field spectrometers are needed to deal with current astrophysical challenges that require multiband observations at millimeter wavelengths. An example of these is the KIDs Interferometer Spectrum Survey (KISS), which uses two arrays of kinetic inductance detectors (KIDs) coupled to a Martin-Puplett interferometer (MPI). KISS has a wide instantaneous field of view (1 deg in diameter) and a spectral resolution up to 1.45 GHz in the 120-180 GHz electromagnetic band. The instrument is installed on the 2.25 m Q-U-I JOint TEnerife telescope in Teide Observatory (Tenerife, Canary Islands), at an altitude of 2395 m above sea level. Aims. This work presents an original readout modulation method developed to improve the sky signal reconstruction accuracy for types of instruments for which a fast sampling frequency is required both to remove atmospheric fluctuations and to perform full spectroscopic measurements on each sampled sky position. Methods. We first demonstrate the feasibility of this technique using simulations. Then, we apply such a scheme to on-sky calibration. Results. We show that the sky signal can be reconstructed to better than 0.5% for astrophysical sources, and to better than 2% for large background variations such as in "skydip", in an ideal noiseless scenario. The readout modulation method is validated by observations on-sky during the KISS commissioning campaign. Conclusions. We conclude that accurate photometry can be obtained for future KID-based MPI.

From Oct. 2019 to Apr. 2020, LAMOST performs a time-domain spectroscopic survey of four $K$2 plates with both low- and med-resolution observations. The low-resolution spectroscopic survey gains 282 exposures ($\approx$46.6 hours) over 25 nights, yielding a total of about 767,000 spectra, and the med-resolution survey takes 177 exposures ($\approx$49.1 hours) over 27 nights, collecting about 478,000 spectra. More than 70%/50% of low-resolution/med-resolution spectra have signal-to-noise ratio higher than 10. We determine stellar parameters (e.g., $T_{\rm eff}$, log$g$, [Fe/H]) and radial velocity (RV) with different methods, including LASP, DD-Payne, and SLAM. In general, these parameter estimations from different methods show good agreement, and the stellar parameter values are consistent with those of APOGEE. We use the $Gaia$ DR2 RV data to calculate a median RV zero point (RVZP) for each spectrograph exposure by exposure, and the RVZP-corrected RVs agree well with the APOGEE data. The stellar evolutionary and spectroscopic masses are estimated based on the stellar parameters, multi-band magnitudes, distances and extinction values. Finally, we construct a binary catalog including about 2700 candidates by analyzing their light curves, fitting the RV data, calculating the binarity parameters from med-resolution spectra, and cross-matching the spatially resolved binary catalog from $Gaia$ EDR3. The LAMOST TD survey is expected to get breakthrough in various scientific topics, such as binary system, stellar activity, and stellar pulsation, etc.

The recent high spatial/spectral resolution observations have enabled constraining formation mechanisms of giant planets, especially at the final stages. The current interpretation of such observations is that these planets undergo magnetospheric accretion, suggesting the importance of planetary magnetic fields. We explore the properties of accreting, magnetized giant planets surrounded by their circumplanetary disks, using the physical parameters inferred for PDS 70 b/c. We compute the magnetic field strength and the resulting spin rate of giant planets, and find that these planets may possess dipole magnetic fields of either a few 10 G or a few 100 G; the former is the natural outcome of planetary growth and radius evolution, while the resulting spin rate cannot reproduce the observations. For the latter, a consistent picture can be drawn, where strong magnetic fields induced by hot planetary interiors lead both to magnetospheric accretion and to spin-down due to disk locking. We also compute the properties of circumplanetary disks in the vicinity of these planets, taking into account planetary magnetic fields. The resulting surface density becomes very low, compared with the canonical models, implying the importance of radial movement of satellite-forming materials. Our model predicts a positive gradient of the surface density, which invokes the traps for both satellite migration and radially drifting dust particles. This work thus concludes that the final formation stages of giant planets are similar to those of low-mass stars such as brown dwarfs, as suggested by recent studies.

Radio-to-TeV observations of the bright nearby (z=0.034) blazar Markarian 501 (Mrk 501), performed from December 2012 to April 2018, are used to study the emission mechanisms in its relativistic jet. We examined the multi-wavelength variability and the correlations of the light curves obtained by eight different instruments, including the First G-APD Cherenkov Telescope (FACT), observing Mrk 501 in very high-energy (VHE) gamma-rays at TeV energies. We identified individual TeV and X-ray flares and found a sub-day lag between variability in these two bands. Simultaneous TeV and X-ray variations with almost zero lag are consistent with synchrotron self-Compton (SSC) emission, where TeV photons are produced through inverse Compton scattering. The characteristic time interval of 5-25 days between TeV flares is consistent with them being driven by Lense-Thirring precession.

Galactic nuclei are promising sites for stellar origin black hole (BH) mergers, as part of merger hierarchies in deep potential wells. We show that binary black hole (BBH) merger rates in active galactic nuclei (AGN) should always exceed merger rates in quiescent galactic nuclei (nuclear star clusters, NSCs) around supermassive BHs (SMBHs) without accretion disks. This is primarily due to average binary lifetimes in AGN that are significantly shorter than in NSCs. The lifetime difference comes from rapid hardening of BBHs in AGN, such that their semi-major axes are smaller than the hard-soft boundary of their parent NSC; this contrasts with the large average lifetime to merger for BBHs in NSCs around SMBHs, due to binary ionization mechanisms. Secondarily, merger rates in AGNs are enhanced by gas-driven binary formation mechanisms. Formation of new BHs in AGN disks are a minor contributor to the rate differences. With the gravitational wave detection of several BBHs with at least one progenitor in the upper mass gap, and signatures of dynamical formation channels in the $\chi_{\rm eff}$ distribution, we argue that AGN could contribute $\sim 25\%-80\%$ of the LIGO-Virgo measured rate of $\sim 24 \rm{Gpc}^{-3} \rm{yr}^{-1}$.

We study gravitational lensing of gravitational waves from compact object binaries as a probe of compact dark matter (DM) objects such as primordial black holes. Assuming a point mass lens, we perform parameter estimation of lensed gravitational wave signals from compact object binaries to determine the detectability of the lens with ground based laser interferometers. Then, considering binary populations that LIGO-Virgo has been probing, we derive a constraint on the abundance of compact DM from non-observation of lensed events. We find that the LIGO-Virgo observations imply that compact objects heavier than $M_l = 50M_\odot$ can not constitute all DM and less than $15\%$ of DM can be in compact objects heavier than $M_l = 200M_\odot$. We also show that the DM fraction in compact objects can be probed by LIGO in its final sensitivity for $M_l > 20M_\odot$ reaching $0.7\%$ of the DM abundance at $M_l > 100M_\odot$, and by ET for $M_l > 0.4M_\odot$ reaching DM fraction as low as $3\times 10^{-5}$ at $M_l > 20M_\odot$.

Works submitted to the 37th International Cosmic Ray Conference. Contents: Operation, Calibration, and Analysis Galactic Gamma-Ray Physics Extra Galactic Gamma-Ray Physics Cosmic-Ray Physics

We analyse field fluctuations during an Ultra Slow-Roll phase in the stochastic picture of inflation and the resulting non-Gaussian curvature perturbation, fully including the gravitational backreaction of the field's velocity. By working to leading order in a gradient expansion, we first demonstrate that consistency with the momentum constraint of General Relativity prevents the field velocity from having a stochastic source, reflecting the existence of a single scalar dynamical degree of freedom on long wavelengths. We then focus on a completely level potential surface, $V=V_0$, extending from a specified exit point $\phi_{\rm e}$, where slow roll resumes or inflation ends, to $\phi\rightarrow +\infty$. We compute the probability distribution in the number of e-folds $\mathcal{N}$ required to reach $\phi_{\rm e}$ which allows for the computation of the curvature perturbation. We find that, if the field's initial velocity is high enough, all points eventually exit through $\phi_{\rm e}$ and a finite curvature perturbation is generated. On the contrary, if the initial velocity is low, some points enter an eternally inflating regime despite the existence of $\phi_{\rm e}$. In that case the probability distribution for $\mathcal{N}$, although normalizable, does not possess finite moments, leading to a divergent curvature perturbation.

In this work, we describe an updated version of single arm locking, and the noise amplification due to the nulls can be flexibly restricted with the help of optical frequency comb. We show that, the laser phase noise can be divided by a specific factor with optical frequency comb as the bridge. The analytical results indicate that, the peaks in the science band have been greatly reduced. The performance of the noise suppression shows that the total noise after arm locking can well satisfy the requirement of time delay interferometry, even with the free-running laser source. We also estimate the frequency pulling characteristics of the updated single arm locking, and the results suggest that the pulling rate can be tolerated, without the risk of mode hopping. Arm locking will be a valuable solution for the noise reduction in the space-borne GW detectors. We demonstrate that, with the precise control of the returned laser phase noise, the noise amplification in the science band can be efficiently suppressed based on the updated single arm locking. Not only our method allows the suppression of the peaks, the high gain, low pulling rate, it can also serve for full year, without the potential risk of locking failure due to the arm length mismatch. We finally discuss the unified demonstration of the updated single arm locking, where both the local and the returned laser phase noises can be tuned to generate the expected arm-locking sensor actually. Our work could provide a powerful method for the arm locking in the future space-borne GW detectors.

Nature's most powerful high-energy sources are capable of accelerating particles to high energy and radiate it away on extremely short timescales, even shorter than the light crossing time of the system. It is yet unclear what physical processes can produce such an efficient acceleration, despite the copious radiative losses. By means of radiative particle-in-cell simulations, we show that magnetically dominated turbulence in pair plasmas subject to strong synchrotron cooling generates a nonthermal particle spectrum with a hard power-law range (slope $p \sim 1$) within a few eddy turnover times. Low pitch-angle particles can significantly exceed the nominal radiation-reaction limit, before abruptly cooling down. The particle spectrum becomes even harder ($p < 1$) over time owing to particle cooling with an energy-dependent pitch-angle anisotropy. The resulting synchrotron spectrum is hard ($\nu F_\nu \propto \nu^s$ with $s \sim 1$). Our findings have important implications for understanding the nonthermal emission from high-energy astrophysical sources, most notably the prompt phase of gamma-ray bursts and gamma-ray flares from the Crab nebula.

These lectures, presented at the 2021 Les Houches Summer School on Dark Matter, provide an introduction to key methods and tools of indirect dark matter searches, as well as a status report on the field circa summer 2021. Topics covered include the possible effects of energy injection from dark matter on the early universe, methods to calculate both the expected energy distribution and spatial distribution of particles produced by dark matter interactions, an outline of theoretical models that predict diverse signals in indirect detection, and a discussion of current constraints and some claimed anomalies. These notes are intended as an introduction to indirect dark matter searches for graduate students, focusing primarily on intuition-building estimates and useful concepts and tools.

We study an extension of the minimal gauged $L_{\mu}-L_{\tau}$ model in order to explain the anomalous magnetic moments of muon and electron simultaneously. Presence of an additional scalar doublet $\eta$ and an in-built $Z_2$ symmetry under which the right handed singlet fermions and $\eta$ are odd, leads to light neutrino mass in scotogenic fashion along with a stable dark matter candidate. In spite of the possibility of having positive and negative contributions to $(g-2)$ from vector boson and charged scalar loops respectively, the minimal scotogenic $L_{\mu}-L_{\tau}$ model can not explain muon and electron $(g-2)$ simultaneously while being consistent with other experimental bounds. We then extend the model with a vector like lepton doublet which not only leads to a chirally enhanced negative contribution to electron $(g-2)$ but also leads to the popular singlet-doublet fermion dark matter scenario. With this extension, the model can explain both electron and muon $(g-2)$ while being consistent with neutrino mass, dark matter and other direct search bounds. The model remains predictive at high energy experiments like collider as well as low energy experiments looking for charged lepton flavour violation, dark photon searches, in addition to future $(g-2)$ measurements.

The observation of Earth matter effects in the spectrum of neutrinos coming from a next galactic core-collapse supernova (CCSN) could, in principle, reveal if neutrino mass ordering is normal or inverted. One of the possible ways to identify the mass ordering is through the observation of the modulations that appear in the spectrum when neutrinos travel through the Earth before they arrive at the detector. These features in the neutrino spectrum depend on two factors, the average neutrino energies, and the difference between the primary neutrino fluxes of electron and other flavors produced inside the supernova. However, recent studies indicate that the Earth matter effect for CCSN neutrinos is expected to be rather small and difficult to be observed by currently operating or planned neutrino detectors mainly because of the similarity of average energies and fluxes between electron and other flavors of neutrinos, unless the distance to CCSN is significantly smaller than the typically expected one, $\sim 10$ kpc. Here, we are looking towards the possibility if the non-standard neutrino properties such as decay of neutrinos can enhance the Earth matter effect. In this work we show that invisible neutrino decay can potentially enhance significantly the Earth matter effect for both $\nu_e$ and $\bar{\nu}_e$ channels at the same time for both mass orderings, even if the neutrino spectra between electron and other flavors of neutrinos are very similar, which is a different feature not expected for CCSN neutrinos with standard oscillation without the decay effect.

This paper presents a relativistic version of Newtonian Fractional-Dimension Gravity (NFDG), an alternative gravitational model recently introduced and based on the theory of fractional-dimension spaces. This extended version - Relativistic Fractional-Dimension Gravity (RFDG) - is based on other existing theories in the literature and might be useful for astrophysical and cosmological applications. In particular, we review the mathematical theory for spaces with non-integer dimensions and its connections with the non-relativistic NFDG. The Euler-Lagrange equations for scalar fields can also be extended to spaces with fractional dimensions, by adding an appropriate weight factor, and then can be used to generalize the Laplacian operator for rectangular, spherical, and cylindrical coordinates. In addition, the same weight factor can be added to the standard Hilbert action in order to obtain the field equations, following methods used for scalar-tensor models of gravity, multi-scale spacetimes, and fractional gravity theories. We then apply the field equations to standard cosmology and to the Friedmann-Lema\^itre-Robertson-Walker metric. Using a suitable weight $v_{t}\left(t\right)$, depending on the synchronous time $t$ and on a single time-dimension parameter $\alpha _{t}$, we extend the Friedmann equations to the RFDG case. This allows for the computation of the scale factor $a\left(t\right)$ for different values of the fractional time-dimension $\alpha _{t}$ and the comparison with standard cosmology results. Future additional work on the subject, including studies of the cosmological late-time acceleration, type Ia supernovae data, and related dark energy theory will be needed to establish this model as a relativistic alternative theory of gravity.

We forecast constraints on neutrino decay via capture of the Cosmic Neutrino Background on tritium, with emphasis on the PTOLEMY-type experiment. In particular, in the case of invisible neutrino decay into lighter neutrinos in the Standard Model and invisible particles, we can constrain not only the neutrino lifetime but also the masses of the invisible particles. For this purpose, we also formulate the energy spectra of the lighter neutrinos produced by 2-body and 3-body decays, and those of the electrons emitted in the process of the detection of the lighter neutrinos.

We review recent work on the possibility of primordial synchronization of different Mixmaster regions in generic inhomogeneous spacetime. It is shown that inhomogeneous domains undergoing chaotic oscillations may synchronize themselves exponentially fast and proceed in perfect symphony asymptotically in the past direction. Implications of this result for the structure and evolution of an early phase of the universe are briefly discussed.

Chemical constants extracted from $^{124}$Xe+ $^{124}$Sn collisions at 32 AMeV are compared to the predictions of an extended Nuclear Statistical Equilibrium model including mean-field interactions and in-medium binding energy shifts for the light ($Z\leq 2$) clusters. The ion species and density dependence of the in-medium modification is directly extracted from the experimental data. We show that the shift increases with the mass of the cluster and the density of the medium, and we provide a simple linear fit for future use in astrophysical simulations in the framework of the CompOSE data base. The resulting mass fractions are computed in representative thermodynamic conditions relevant for supernova and neutron star mergers. A comparison to the results of a similar analysis of the same data performed in the framework of a relativistic mean-field model shows a good agreement at low density, but significant discrepancies close to the Mott dissolution of clusters in the dense medium.

The ringdown data of both GW150914 and GW190521\_074359 (GW190521r) show evidence supporting the presence of overtone. Previous studies all adopt a fundamental assumptions, which were motivated more by convenience than by first principles, that the first overtone and the fundamental mode share a same onset. In this work, for the first time we relax such assumption, and we aim to probe the possible chromatic onsets of these two components within the GW150914 and GW190521r ringdown data. For both events, we bound the onset lags to be $\Delta t_0\leq 10M_f$ at high probabilities (i.e., $\geq 96.5\%$), where $M_f$ is the mass of the remnant black hole formed in the merger. Moreover, for GW190521r (GW150914) we have $\Delta t_0\geq 3M_f$ at $98.3\%$ ($69.5\%$) credibility, indicating the non-simultaneous excitation between the fundamental mode and the first overtone in the ringdown.

The neutrino floor is a theoretical lower limit on WIMP-like dark matter models that are discoverable in direct detection experiments. It is commonly interpreted as the point at which dark matter signals become hidden underneath a remarkably similar-looking background from neutrinos. However, it has been known for some time that the neutrino floor is not a hard limit, but can be pushed past with sufficient statistics. As a consequence, some have recently advocated for calling it the "neutrino fog" instead. The downside of current methods of deriving the neutrino floor are that they rely on arbitrary choices of experimental exposure and energy threshold. Here we propose to define the neutrino floor as the boundary of the neutrino fog, and develop a calculation free from these assumptions. The technique is based on the derivative of a hypothetical experimental discovery limit as a function of exposure, and leads to a neutrino floor that is only influenced by the systematic uncertainties on the neutrino flux normalisations. Our floor is broadly similar to those found in the literature, but differs by almost an order of magnitude in the sub-GeV range, and above 20 GeV.

Using an exact law for incompressible Hall magnetohydrodynamics (HMHD) turbulence, the energy cascade rate is computed from three-dimensional HMHD-CGL (bi-adiabatic ions and isothermal electrons) and Landau fluid (LF) numerical simulations that feature different intensities of Landau damping over a broad range of wavenumbers, typically $0.05\lesssim k_\perp d_i \lesssim100$. Using three sets of cross-scale simulations where turbulence is initiated at large, medium and small scales, the ability of the fluid energy cascade to "sense" the kinetic Landau damping at different scales is tested. The cascade rate estimated from the exact law and the dissipation calculated directly from the simulation are shown to reflect the role of Landau damping in dissipating energy at all scales, with an emphasis on the kinetic ones. This result provides new prospects on using exact laws for simplified fluid models to analyze dissipation in kinetic simulations and spacecraft observations, and new insights into theoretical description of collisionless magnetized plasmas.

We present a modified commercial L-4C geophone with interferometric readout that demonstrated a resolution 60 times lower than the included coil-magnet readout at low frequencies. The intended application for the modified sensor is in vibration isolation platforms that require improved performance at frequencies lower than 1 Hz. A controls and noise-model of an Advanced LIGO 'HAM-ISI' vibration isolation system was developed, and it shows that our sensor can reduce the residual vibration by a factor of 70 at 0.1 Hz