Papers for Tuesday, Jun 28 2022

We combine the isothermal Jeans model and the model of adiabatic halo contraction into a simple semi-analytic procedure for computing the density profile of self-interacting dark-matter (SIDM) haloes with the gravitational influence from the inhabitant galaxies. We show that the model agrees well with cosmological SIDM simulations over the entire core-forming stage and up to the onset of gravothermal core-collapse. Using this model, we show that the halo response to baryons is more diverse in SIDM than in CDM and depends sensitively on galaxy size, a desirable link in the context of the structural diversity of bright dwarf galaxies. The fast speed of the method facilitates analyses that would be challenging for numerical simulations -- notably, 1) we quantify the SIDM halo response as functions of the baryonic properties, on a fine mesh grid spanned by the baryon-to-total-mass ratio, $M_{\rm b}/M_{\rm vir}$, and galaxy compactness, $r_{1/2}/R_{\rm vir}$; 2) we show with high statistical precision that for typical Milky-Way-like systems, the SIDM profiles are similar to their CDM counterparts; and 3) we delineate the regime of gravothermal core-collapse in the $M_{\rm b}/M_{\rm vir}$-$r_{1/2}/R_{\rm vir}$ space, for a given cross section and a given halo concentration. Finally, we compare the isothermal Jeans model with the more sophisticated gravothermal fluid model, and show that the former yields faster core formation and agrees better with cosmological simulations. We attribute the difference to whether the target CDM halo is used as a boundary condition or as the initial condition for the gravothermal evolution, and thus comment on possible future improvements of the fluid model. We have made our programs for the model publicly available at https://github.com/JiangFangzhou/SIDM.

The ``Spectroscopy and H-band Imaging of Virgo cluster galaxies'' (SHIVir) survey is an optical and near-infrared survey which combines SDSS photometry, deep H-band photometry, and long-slit optical spectroscopy for 190 Virgo cluster galaxies (VCGs) covering all morphological types over the stellar mass range log (M_*/M_Sun) = 7.8-11.5$. We present the spectroscopic sample selection, data reduction, and analysis for this SHIVir sample. We have used and optimised the \texttt{pPXF} routine to extract stellar kinematics from our data. Ultimately, resolved kinematic profiles (rotation curves and velocity dispersion profiles) are available for 133 SHIVir galaxies. A comprehensive database of photometric and kinematic parameters for the SHIVir sample is presented with: grizH magnitudes, effective surface brightnesses, effective and isophotal radii, rotational velocities, velocity dispersions, and stellar and dynamical masses. Parameter distributions highlight some bimodal distributions and possible sample biases. A qualitative study of resolved extended velocity dispersion profiles suggests a link between the so-called ``sigma-drop'' kinematic profile and the presence of rings in lenticular S0 galaxies. Rising dispersion profiles are linked to early-type spirals or dwarf ellipticals for which a rotational component is significant, whereas peaked profiles are tied to featureless giant ellipticals.

The recent Gaia Data Release 3 (DR3) represents an unparalleled revolution in Galactic Archaeology, providing us with numerous radial velocities chemical abundances for million of stars, overcoming the spatial biases that ground-based spectroscopic surveys suffer from ground-based surveys. We present a new chemical evolution model for the Galactic disc components (high- and low- $\alpha$ sequence stars) designed to reproduce the new abundance ratios provided by the GSP-spec module for the Gaia DR3 and also constrained by i) stellar ages and ii) detailed star formation histories for both the thick and thin disc stars inferred from previous Gaia releases. Gaia DR3 highlighted the presence of young (massive) low-$\alpha$ disc stars which show evidence of a recent chemical impoverishment in several elements. Additionally, previous Gaia releases have found evidence for narrow episodes of enhanced SF inferred in recent time. In order to reproduce these observables, we propose a new chemical evolution model in which the low-$\alpha$ sequence is generated by two distinct infall episodes. Hence, in this study we compare Gaia DR3 chemical abundances with the predictions of a three-infall chemical evolution model for the high- and low-$\alpha$ components. The proposed three-infall chemical evolution model nicely reproduces the main features of the abundance ratio [X/Fe] versus [M/H] (X=Mg, Si, Ca, Ti, $\alpha$) of Gaia DR3 stars in different age bins for the considered $\alpha$ elements. Moreover, the most recent gas infall - which started $\sim$ 2.7 Gyr ago - allows us to predict well the Gaia DR3 young population which has experienced a recent chemical impoverishment.

Historically, FU Orionis-type stars are low-mass, pre-main sequence stars. The members of this class experience powerful accretion outbursts and remain in an enhanced accretion state for decades or centuries. V1515 Cyg, a classical FUor, started brightening in the 1940s and reached its peak brightness in the late 1970s. Following a sudden decrease in brightness it stayed in a minimum state for a few months, then started a brightening for several years. We present results of our ground-based photometric monitoring complemented with optical/NIR spectroscopic monitoring. Our light curves show a long-term fading with strong variability on weekly and monthly time scales. The optical spectra show P Cygni profiles and broad blue-shifted absorption lines, common properties of FUors. However, V1515 Cyg lacks the P Cygni profile in the Ca II 8498 \r{A} line, a part of the Ca infrared triplet (IRT), formed by an outflowing wind, suggesting that the absorbing gas in the wind is optically thin. The newly obtained near-infrared spectrum shows the strengthening of the CO bandhead and the FeH molecular band, indicating that the disk has become cooler since the last spectroscopic observation in 2015. The current luminosity of the accretion disk dropped from the peak value of 138 $L_{\odot}$ to about 45 $L_{\odot}$, suggesting that the long-term fading is also partly caused by the dropping of the accretion rate.

CzeV343 was previously identified as a candidate double eclipsing binary (2+2 quadruple), where the orbital periods of the two eclipsing binaries ($P_A \approx 1.2$ days and $P_B \approx 0.8$ days) lie very close to 3:2 resonance. Here, we analyze 11 years of ground-based photometry, 4 sectors of TESS 2-minute and full-frame photometry, and two optical spectra. We construct a global model of our photometry, including apsidal motion of binary A and light-travel time effect (LTTE) of the mutual outer orbit, and explore the parameter space with Markov Chain Monte Carlo. We estimate component masses for binary A ($1.8+1.3 M_\odot$) and binary B ($1.4+1.2 M_\odot$). We identify pseudo-synchronous rotation signal of binary A in TESS photometry. We detect apsidal motion in binary A with a period of about 33 years, which is fully explained by tidal and rotational contributions of stars aligned with the orbit. The mutual orbit has a period of about 1450 days and eccentricity of about 0.7. The LTTE amplitude is small, which points to low inclination of the outer orbit and a high degree of misalignment with the inner orbits. We find that when apsidal motion and mutual orbit are taken into account the orbital period resonance is exact to within $10^{-5}$ cycles/day. Many properties of CzeV343 are not compatible with requirements of the 3:2 resonance capture theory for coplanar orbits. Future evolution of CzeV343 can lead to mergers, triple common envelope, double white dwarf binaries, or a Type Ia supernova. More complex evolutionary pathways will likely arise from dynamical instability caused by orbital expansion when either of the binaries undergoes mass transfer. This instability has not been so far explored in 2+2 quadruples.

We present a deep neural network Real/Bogus classifier that improves classification performance in the Tomo-e Gozen transient survey by handling label errors in the training data. In the wide-field, high-frequency transient survey with Tomo-e Gozen, the performance of conventional convolutional neural network classifier is not sufficient as about $10^6$ bogus detections appear every night. In need of a better classifier, we have developed a new two-stage training method. In this training method, label errors in the training data are first detected by normal supervised learning classification, and then they are unlabeled and used for training of semi-supervised learning. For actual observed data, the classifier with this method achieves an area under the curve (AUC) of 0.9998 and a false positive rate (FPR) of 0.0002 at true positive rate (TPR) of 0.9. This training method saves relabeling effort by humans and works better on training data with a high fraction of label errors. By implementing the developed classifier in the Tomo-e Gozen pipeline, the number of transient candidates was reduced to $\sim$40 objects per night, which is $\sim$1/130 of the previous version, while maintaining the recovery rate of real transients. This enables more efficient selection of targets for follow-up observations.

As part of our program to identify host galaxies of known z=2-3 MgII absorbers with the Keck Cosmic Web Imager (KCWI), we discovered a compact group giving rise to a z=2.431 DLA with ultra-strong MgII absorption in quasar field J234628+124859. The group consists of four star-forming galaxies within 8-28 kpc and $v\sim50-300$ km s$^{-1}$ of each other, where tidal streams are weakly visible in deep HST imaging. The group geometric centre is D=25 kpc from the quasar (D=20-40 kpc for each galaxy). Galaxy G1 dominates the group ($1.66L_{\ast}$, ${\rm SFR}_{\rm FUV}=11.6$ M$_{\odot}$ yr$^{-1}$) while G2, G3, and G4 are less massive ($0.1-0.3L_{\ast}$, ${\rm SFR}_{\rm FUV}=1.4-2.0$ M$_{\odot}$ yr$^{-1}$). Using a VLT/UVES quasar spectrum covering the HI Lyman series and metal lines such as MgII, SiIII, and CIV, we characterised the kinematic structure and physical conditions along the line-of-sight with cloud-by-cloud multiphase Bayesian modelling. The absorption system has a total $\log(N(HI)/{\rm cm}^{-2})=20.53$ and an $N(HI)$-weighted mean metallicity of $\log(Z/Z_{\odot})=-0.68$, with a very large MgII linewidth of $\Delta v\sim700$ km s$^{-1}$. The highly kinematically complex profile is well-modelled with 30 clouds across low and intermediate ionisation phases with values ${13\lesssim\log(N(HI)/{\rm cm}^{-2})\lesssim20}$ and $-3\lesssim\log(Z/Z_{\odot})\lesssim1$. Comparing these properties to the galaxy properties, we infer a wide range of gaseous environments, including metal-rich outflows, metal-poor IGM accretion, and tidal streams from galaxy--galaxy interactions. This diversity of structures forms the intragroup medium around a complex compact group environment at the epoch of peak star formation activity. Surveys of low redshift compact groups would benefit from obtaining a more complete census of this medium for characterising evolutionary pathways.

One of the most important questions in blazar physics is the origin of broadband emission and fast-flux variation. In this work, we studied the broadband temporal and spectral properties of a TeV blazar 1ES 1727+502 and explore the one-zone synchrotron-self Compton (SSC) model to fit the broadband spectral energy distribution (SED). We collected the long-term (2014-2021) multiband data which includes both the low and high flux states of the source. The entire light curve is divided into three segments of different flux states and the best-fit parameters obtained by broadband SED modeling corresponding to three flux states were then compared. The TeV blazar 1ES 1727+502 has been observed to show the brightest flaring episode in X-ray followed by optical-UV and gamma-ray. The fractional variability estimated during various segments behaves differently in multiple wavebands, suggesting a complex nature of emission in this source. This source has shown a range of variability time from days scale to month scale during this long period of observations between 2014-2021. A "harder-when-brighter" trend is not prominent in X-ray but seen in optical-UV and an opposite trend is observed in gamma-ray. The complex nature of correlation among various bands is observed. The SED modeling suggests that the one-zone SSC emission model can reproduce the broadband spectrum in the energy range from optical-UV to very high energy gamma-ray.

Rotation causes an increase in a neutron star's mass and equatorial radius. The mass and radius depend sensitively on the unknown equation of state (EOS) of cold, dense matter. However, the increases in mass and radius due to rotation are almost independent of the EOS. The EOS independence leads to the idea of neutron star universality. In this paper, we compute sequences of rotating neutron stars with constant central density. We use a collection of randomly generated EOS to construct simple correction factors to the mass and radius computed from the equations of hydrostatic equilibrium for non-rotating neutron stars. The correction factors depend only on the non-rotating star's mass and radius and are almost independent of the EOS. This makes it computationally inexpensive to include observations of rotating neutron stars in EOS inference codes. We also construct a mapping from the measured mass and radius of a rotating neutron star to a corresponding non-rotating star. The mapping makes it possible to construct a zero-spin mass-radius curve if the masses and radii of many neutron stars with different spins are measured. We show that the changes in polar and equatorial radii are symmetric, in that the polar radius shrinks at the same rate that the equatorial radius grows. This symmetry is related to the observation that the equatorial compactness (the ratio of mass to radius) is almost constant on one of the constant-density sequences.

Massive Early-Type Galaxies (ETGs) show several strong CO absorption features in their H- and K-band spectra that cannot be explained by state-of-the-art stellar population models. For many years, the disagreement has been attributed to the presence of intermediate-age stellar components that are dominated by stars in the Asymptotic Giant Branch (AGB) phase. However, no robust evidence of this scenario has been provided so far. One way to test this claim is by comparison of CO indices for ETGs and for relic galaxies. Lacking the intermediate-age stellar populations, relic galaxies provide us with a unique opportunity to address the origin of strong CO absorptions in ETGs. Here, we utilize the prototype relic galaxy NGC 1277 and compare the CO absorption features of this galaxy with the ones of a representative sample of massive ETGs. We show that the CO lines in both systems have similar strengths, significantly stronger than the predictions of stellar population synthesis models. We conclude that intermediate-age stellar populations in massive ETGs are not the culprit of the strong CO absorptions.

Forecasting the solar cycle amplitude is important for a better understanding of the solar dynamo as well as for many space weather applications. We demonstrated a steady relationship between the maximal growth rate of sunspot activity in the ascending phase of a cycle and the subsequent cycle amplitude on the basis of four data sets of solar activity indices: total sunspot numbers, hemispheric sunspot numbers from the new catalogue from 1874 onwards, total sunspot areas, and hemispheric sunspot areas. For all the data sets, a linear regression based on the maximal growth rate precursor shows a significant correlation. Validation of predictions for cycles 1-24 shows high correlations between the true and predicted cycle amplitudes reaching r = 0.93 for the total sunspot numbers. The lead time of the predictions varies from 2 to 49 months, with a mean value of 21 months. Furthermore, we demonstrated that the sum of maximal growth rate indicators determined separately for the north and the south hemispheric sunspot numbers provides more accurate predictions than that using total sunspot numbers. The advantages reach 27% and 11% on average in terms of rms and correlation coefficient, respectively. The superior performance is also confirmed with hemispheric sunspot areas with respect to total sunspot areas. The maximal growth rate of sunspot activity in the ascending phase of a solar cycle serves as a reliable precursor of the subsequent cycle amplitude. Furthermore, our findings provide a strong foundation for supporting regular monitoring, recording, and predictions of solar activity with hemispheric sunspot data, which capture the asymmetric behaviour of the solar activity and solar magnetic field and enhance solar cycle prediction methods.

We have analysed two cosmological zoom simulations with $M_{\rm vir}\sim 10^{12}{\rm\,M}_\odot$ from the NIHAO series, both with and without feedback. We show that an entropy criterion based on the equation of state of the intergalactic medium can successfully separate cold- and hot-mode accretion. The shock-heated gas has non-negligible turbulent support and cools inefficiently. In the simulations without feedback, only a small fraction ($\sim 20$ per cent) of the stellar mass comes from baryons that have been in the hot circumgalactic medium, although quantitative conclusions should be taken with caution due to our small-number statistics. With feedback, the fraction is larger because of the reaccretion of gas heated by supernovae, which has lower entropies and shorter cooling times than the gas heated by accretion shocks. We have compared the results of NIHAO to predictions of the GalICS 2.1 semianalytic model of galaxy formation. The shock-stability criterion implemented in GalICS 2.1 successfully reproduces the transition from cold- to hot-mode accretion.

We present the fifth incarnation of the Mid-Infrared Array Camera (MIRAC-5) instrument which will use a new GeoSnap (3 - 13 microns) detector. Advances in adaptive optics (AO) systems and detectors are enabling ground-based mid-infrared systems capable of high spatial resolution and deep contrast. As one of the only 3 - 13 micron cameras used in tandem with AO, MIRAC-5 will be complementary to the James Webb Space Telescope (JWST) and capable of characterizing gas giant exoplanets and imaging forming protoplanets (helping to characterize their circumplanetary disks). We describe key features of the MIRAC-5 GeoSnap detector, a long-wave Mercury-Cadmium-Telluride (MCT) array produced by Teledyne Imaging Sensors (TIS), including its high quantum efficiency (> 65%), large well-depth, and low noise. We summarize MIRAC-5's important capabilities, including prospects for obtaining the first continuum mid-infrared measurements for several gas giants and the first 10.2-10.8 micron NH3 detection in the atmosphere of the warm companion GJ 504b (Teff ~550 K) within 8 hours of observing time. Finally, we describe plans for future upgrades to MIRAC-5 such as adding a coronagraph. MIRAC-5 will be commissioned on the MMT utilizing the new MAPS AO system in late 2022 with plans to move to Magellan with the MagAO system in the future.

Phosphorus is a necessary element for life on Earth, but at present we have limited constraints on its chemistry in star- and planet-forming regions: to date, phosphorus carriers have only been detected towards a few low-mass protostars. Motivated by an apparent association between phosphorus molecule emission and outflow shocking, we used the IRAM 30m telescope to target PN and PO lines towards seven Solar-type protostars with well-characterized outflows, and firmly detected phosphorus molecules in three new sources. This sample, combined with archival observations of three additional sources, enables the first exploration of the demographics of phosphorus chemistry in low-mass protostars. The sources with PN detections show evidence for strong outflow shocks based on their H$_2$O 1$_{10}$-1$_{01}$ fluxes. On the other hand, no protostellar properties or bulk outflow mechanical properties are found to correlate with the detection of PN. This implies that gas-phase phosphorus is specifically linked to shocked gas within the outflows. Still, the PN and PO line kinematics suggest an emission origin in post-shocked gas rather than directly shocked material. Despite sampling a wide range of protostellar properties and outflow characteristics, we find a fairly narrow range of source-averaged PO/PN ratios (0.6-2.2) and volatile P abundances as traced by (PN+PO)/CH$_3$OH ($\sim$1-3%). Spatially resolved observations are needed to further constrain the emission origins and environmental drivers of the phosphorus chemistry in these sources.

A statistical model of discrete finite length random processes with negative power law spectral densities is presented. The definition of terms is followed by a description of the spectral density trend. An algorithmic construction of random process, and a short block of computer code is given to implement the construction of the random process. The relationship between the second order properties of the random processes and the parameters of the construction is developed and demonstrated. The paper ends with a demonstration of the connection between the frequency of the random process sign changes and the power law exponent.

We present the results of high-resolution spectropolarimetric observations of the optically dominant component in the rare hydrogen-deficient binary system upsilon Sgr. Only a small number of such systems in a very late phase of helium shell burning are currently known. The mass transfer from the donor star in binary systems usually leads to the stripping of its hydrogen envelope. Consequently, since the mass of the secondary increases, it appears rejuvenated. Using a few ESO FORS1 low-resolution spectropolarimetric observations of this system, Hubrig et al. announced in 2009 the presence of a magnetic field of the order of -70 - -80G. Here we report on more recent high-resolution ESO HARPS spectropolarimetric observations showing that the primary in upsilon Sgr is a spectrum variable star and possesses a weak magnetic field of the order of a few tens of Gauss. The detection of a magnetic field in this rare hydrogen-deficient binary is of particular interest, as such systems are frequently discussed as probable progenitors of core-collapse supernovae and gravitational-wave sources. Future magnetic studies of such systems will be worthwhile to gain deeper insights into the role of magnetic fields in the evolution of massive stars in binary systems.

Previously, we proposed a new model of the origin of life, named Nebula-Relay Nebula-Relay, which assumed that the life on Earth originated in the planetary system of the sun's predecessor star and then filled in the pre-solar nebula after its death. What is the source of life's energy in the molecular clouds? This draft discussed two possible mechanisms: methanogenesis and cosmic ray-driven bioenergetics. We found that enough free energy is released from the chemical reaction for methanogens. But the scarcity of carbon compounds is a possible limiting factor. The second one is driven by cosmic ray ionization, which means that radiation hazards become the source of life energy. Protons are naturally produced in this scenario, which may be chemiosmosis's origin.

The Blanco DECam Bulge Survey (BDBS) has imaged more than 200 square degrees of the southern Galactic bulge, providing photometry in the ugrizy filters for $\sim 250$ million unique stars. The presence of a strong foreground disk population, along with complex reddening and extreme image crowding, has made it difficult to constrain the presence of young and intermediate age stars in the bulge population. We employed an accurate cross-match of BDBS with the latest data release (EDR3) from the Gaia mission, matching more than 140 million sources with BDBS photometry and Gaia EDR3 photometry and astrometry. We relied on Gaia EDR3 astrometry, without any photometric selection, to produce clean BDBS bulge colour-magnitude diagrams (CMDs). Gaia parallaxes were used to filter out bright foreground sources, and a Gaussian mixture model fit to Galactic proper motions could identify stars kinematically consistent with bulge membership. We applied this method to 127 different bulge fields of $1$ deg$^2$ each, with $|\ell| \leq 9.5^\circ$ and $-9.5^\circ \leq b \leq -2.5^\circ$. The astrometric cleaning procedure removes the majority of blue stars in each field, especially near the Galactic plane, where the ratio of blue to red stars is $\lesssim 10\%$, increasing to values $\sim 20\%$ at higher Galactic latitudes. We rule out the presence of a widespread population of stars younger than 2 Gyr. The vast majority of blue stars brighter than the turnoff belong to the foreground population, according to their measured astrometry. We introduce the distance between the observed red giant branch bump and the red clump as a simple age proxy for the dominant population in the field, and we confirm the picture of a predominantly old bulge. Further work is needed to apply the method to estimate ages to fields at higher latitudes, and to model the complex morphology of the Galactic bulge.

Recent asteroseismic observations by the $\it Kepler$ space mission have revealed the dip fine structure in the period-spacing versus period diagram of rapidly rotating $\gamma$ Doradus stars. Following the successful reproduction of the dip structure by numerical calculations in previous studies, we present in this paper the physical mechanism of how the dip is formed as a result of the interaction between the gravito-inertial waves in the radiative envelope and the pure inertial waves in the convective core. We analytically describe the wave solutions in both of the radiative envelope and the convective core, and match them at the interface to construct an eigenmode. We have found from the analysis the following points: the dip structure is mainly controlled by a parameter that has an inverse correlation with Brunt-V\"ais\"al\"a frequency at the interface; the depth and the width of the dip is shallower and larger, respectively, as the parameter gets large; the shape of the dip can be approximated by the Lorentzian function; the period at the central position of the dip is equal to or slightly smaller than that of the involved pure inertial mode in the convective core. We have also understood based on the evolutionary models of main-sequence stars that the parameter is inversely correlated with the chemical composition gradient at the convective-core boundary. The dip structure thus would provide information about the poorly-understood physical processes, such as diffusion, convective overshooting and rotational mixing, around the boundary between the convective core and the radiative envelope.

We aim at establishing safe membership and evolutionary status of 11 chemically peculiar (CP) stars that are residing in the domain of the open cluster NGC2516 and are frequently referred to as cluster members. We queried the Gaia EDR3 catalogue in an area with a radius of 1deg and selected 37508 stars brighter than G=19mag. The cluster membership was determined in parallax-proper motion-space and 719 probable and 764 possible members were found. The obtained average astrometric and photometric parameters of the cluster are in good agreement with the most recent literature data. The evolutionary status of the target stars was determined with respect to Padova isochrones. After minor adjustments including the metallicity, the reddening, and the transformation scale variation, a perfect fit of the model to the observations over the whole observed magnitude range was achieved. Only 5 of the 11 considered CP stars could be classified as highly probable cluster members. Among the Ap/Bp stars with previously detected magnetic fields HD65987 and HD65712 have a high membership probability and the magnetic star CPD-60 944B is a possible cluster member. Further we discuss the blue straggler nature of HD66194 and the magnetic star HD65987. To our knowledge, HD65987 is currently the only known blue straggler, with a field of the order of a few hundred Gauss. The most striking result of our study is that the strongly magnetic A0p star HD66318 with previously reported very low fractional age does not belong to the NGC2516 cluster at a high level of confidence.

Several space missions and instruments for UV spectropolarimetry are in preparation, such as the proposed NASA MIDEX Polstar project, the proposed ESA M mission Arago, and the Pollux instrument on the future LUVOIR-like NASA flagship mission. In the frame of Polstar, we have studied the capabilities these observatories would offer to gain information on the magnetic and plasma properties of the magnetospheres of hot stars, helping us test the fundamental hypothesis that magnetospheres should act to rapidly drain angular momentum, thereby spinning the star down, whilst simultaneously reducing the net mass-loss rate. Both effects are expected to lead to dramatic differences in the evolution of magnetic vs. non-magnetic stars.

We report a one-second-cadence wide-field survey for M-dwarf flares using the Tomo-e Gozen camera mounted on the Kiso Schmidt telescope. We detect 22 flares from M3-M5 dwarfs with rise times and amplitudes ranging from $5\, \mathrm{sec} \lesssim t_\mathrm{rise} \lesssim 100\,\mathrm{sec}$ and $0.5 \lesssim \Delta F/F_{\star} \lesssim 20$, respectively. The flare light curves mostly show steeper rises and shallower decays than those obtained from the Kepler one-minute cadence data and tend to have flat peak structures. Assuming a blackbody spectrum with temperatures of $9,000-15,000\,\mathrm{K}$, the peak luminosities and bolometric energies are estimated to be $10^{29}\,\mathrm{erg\,sec^{-1}} \lesssim L_\mathrm{peak} \lesssim 10^{31}\,\mathrm{erg\,sec^{-1}}$ and $10^{31}\,\mathrm{erg} \lesssim E_{\rm bol} \lesssim 10^{34}\,\mathrm{erg}$, which constitutes the bright end of fast optical flares for M dwarfs. We confirm that more than 90\% of the host stars of the detected flares are magnetically active based on their H$\alpha$ emission line intensities obtained by LAMOST. The estimated occurrence rate of the detected flares is $\sim 0.7$ per day per an active star, indicating they are common in magnetically active M dwarfs. We argue that the flare light curves can be explained by the chromospheric compression model; the rise time is broadly consistent with the Alfv\'en transit time of a magnetic loop with a length scale of $l_\mathrm{loop} \sim 10^4\,\mathrm{km}$ and a field strength of $1,000\,\mathrm{G}$, while the decay time is likely determined by the radiative cooling of the compressed chromosphere down to near the photosphere with a temperature of $\gtrsim 10,000\,\mathrm{K}$. These flares from M dwarfs could be a major contamination source for a future search of fast optical transients of unknown types.

The additional signals observed in the frequency spectra of the first-overtone RR Lyrae stars, that form a period ratio around 0.61 with the period of the first overtone, are a common phenomenon for RRc and RRd stars, as well as for first-overtone classical Cepheids. The recently proposed model explains these signals as harmonics of non-radial modes of degrees 8 or 9 in the case of RR Lyrae stars and 7, 8, or 9 in the case of classical Cepheids. We selected at least triple-mode RR Lyrae stars pulsating in radial and non-radial modes for asteroseismic modeling. We assume the identification of the non-radial modes as predicted by the model. We calculated a dense grid of models for RR Lyrae stars using envelope pulsation code. By matching first overtone period and period ratios we obtained physical parameters for the selected sample of triple-mode stars. It is the very first attempt of modeling RR Lyrae stars with non-radial modes. We compared our results with predictions of stellar evolution theory, which resulted in a mass discrepancy more noticeable for long-period stars: pulsation masses seem higher than evolutionary masses. We compared metallicity estimates for RRc stars from modeled sample with metallicities determined spectroscopically for a sample of RRc stars in the solar neighbourhood: both distributions are consistent.

It is reasonable to assume that the structure of a planet and the interior distribution of its components are determined by its formation history. We thus follow the growth of a planet from a small embryo through its subsequent evolution. We estimate the accretion rate range based on a protoplanetary disk model at a large enough distance from the central star, for water ice to be a major component. We assume the accreted material to be a mixture of silicate rock and ice, with no H-He envelope, as the accretion timescale is much longer than the time required for the nebular gas to dissipate. We adopt a thermal evolution model that includes accretional heating, radioactive energy release, and separation of ice and rock. Taking the Safronov parameter and the ice-to-rock ratio as free parameters, we compute growth and evolutionary sequences for different parameter combinations, for 4.6 Gyr. We find the final structure to depend significantly on both parameters. Low initial ice to rock ratios and high accretion rates, each resulting in increased heating rate, lead to the formation of extended rocky cores, while the opposite conditions leave the composition almost unchanged and result in relatively low internal temperatures. When rocky cores form, the ice-rich outer mantles still contain rock mixed with the ice. We find that a considerable fraction of the ice evaporates upon accretion, depending on parameters, and assume it is lost, thus the final surface composition and bulk density of the planet do not necessarily reflect the protoplanetary disk composition.

The X-ray plateau emission observed in many Long gamma-ray bursts (LGRBs) has been usually interpreted as the spin-down luminosity of a rapidly spinning, highly magnetized neutron star (millisecond magnetar). If this is true, then the magnetar may emit extended gravitational wave (GW) emission associated with the X-ray plateau due to non-axisymmetric deformation or various stellar oscillations. The advanced LIGO and Virgo detectors have searched for long-duration GW transients for several years, no evidence of GWs from any magnetar has been found until now. In this work, we attempt to search for signature of GW radiation in the electromagnetic observation of 30 LGRBs under the assumption of the magnetar model. We utilize the observations of the LGRB plateau to constrain the properties of the new-born magnetar, including the initial spin period $P_0$, diploe magnetic field strength $B_p$ and the ellipticity $\epsilon$. We find that there are some tight relations between magnetar parameters, e.g., $\epsilon \propto B_p^{1.29}$ and $B_p\propto P_0^{1.14}$. In addition, we derive the GW strain for magnetar sample via their spin-down processes, and find that the GWs from these objects may not be detectable by the aLIGO and ET detectors. For a rapidly spinning magnetar ($P\sim 1\mbox{ ms}$, $ B \sim10^{15}\mbox{ G}$), the detection horizon for advanced LIGO O5 detector is $\sim 180\mbox{ Mpc}$. The detection of such GW signal associated with the X-ray plateau would be a smoking gun that the central engine of GRB is a magnetar.

Accurate calculation of the phase lags of quasi-periodic oscillations (QPOs) will provide insight into their origin. In this paper we investigate the phase lag correction method which has been applied to calculate the intrinsic phase lags of the QPOs in MAXI J1820+070. We find that the traditional additive model between BBN and QPOs in the time domain is rejected, but the convolution model is accepted. By introducing a convolution mechanism in the time domain, the Fourier cross-spectrum analysis shows that the phase lags between QPOs components in different energy bands will have a simple linear relationship with the phase lags between the total signals, so that the intrinsic phase lags of the QPOs can be obtained by linear correction. The power density spectrum (PDS) thus requires a multiplicative model to interpret the data. We briefly discuss a physical scenario for interpreting the convolution. In this scenario, the corona acts as a low-pass filter, the Green's function containing the noise is convolved with the QPOs to form the low-frequency part of the PDS, while the high-frequency part requires an additive component. We use a multiplicative PDS model to fit the data observed by Insight-HXMT. The overall fitting results are similar compared to the traditional additive PDS model. Neither the width nor the centroid frequency of the QPOs obtained from each of the two PDS models were significantly different, except for the r.m.s. of the QPOs. Our work thus provides a new perspective on the coupling of noise and QPOs.

A major outstanding problem in solar physics is the confinement of the solar tachocline, the thin shear layer that separates nearly solid-body rotation in the radiative interior from strong differential rotation in the convection zone. Here, we present the first 3-D, global solar simulation in which a tachocline is confined by a self-excited dynamo. The non-axisymmetric magnetism is initially built in the convection zone and then diffusively imprints downward. Additionally, the field is locally amplified throughout the radiative interior by vigorous horizontal motions that arise from equatorial Rossby waves and possibly shear instabilities. Our work thus challenges the long-held notion that the Sun's dynamo magnetic field is amplified only as deep as the tachocline and stored in a quiescent radiative interior.

The distribution of young stars into OB associations has long been in need of updating. High-precision {\sl Gaia} early Data Release 3 astrometry, coupled with modern machine-learning methods, allows this to be done. We have compiled a well-defined sample which includes OB stars and young open clusters, in total comprising about 47,700 objects. To break the sample down into groupings resembling associations, we applied the HDBSCAN$^{*}$ clustering algorithm. We used a Monte Carlo method to estimate the kinematic ages of the resulting clusters and the Student's $t$-test to assess the significance of the linear correlations between proper motions and coordinates, indicating the presence of possible cluster expansion signatures. The ages of the majority of clusters demonstrating a general expansion at a 1$\sigma$ confidence level are several tens of Myr, which is in agreement with the expected ages of OB associations. We found 32 open clusters which turned out to be members of the resulting groupings; their ages are consistent with one another within the uncertainties. Comparison of the clusters thus obtained with the historical composition of OB associations in the literature shows a correspondence between their positions in the Galaxy but an apparent absence of good one-to-one stellar matches. Therefore, we suggest that the historical composition of OB associations needs to be revised.

I present progenitor luminosities ($L$) for a sample of 112 Type II supernovae (SNe II), computed directly from progenitor photometry and the bolometric correction technique, or indirectly from empirical correlations between progenitor luminosity and [OI] $\lambda\lambda$6300, 6364 line luminosity at 350 d since explosion, $^{56}$Ni mass, or absolute $V$-band magnitude at 50 d since explosion. To calibrate these correlations, I use twelve SNe II with progenitor luminosities measured from progenitor photometry. I find that the correlations mentioned above are strong and statistically significant, and allow to estimate progenitor luminosities to a precision between 20 and 24 per cent. I correct the SN sample for selection bias and define a subsample of 112 SNe II with progenitor luminosities between $\log(L/L_{\odot})=4.6$ dex, corresponding to the completeness limit of the corrected sample, and the maximum observed progenitor luminosity of $\log(L/L_{\odot})=5.091$ dex. The luminosity distribution for this subsample is statistically consistent with those for red supergiants (RSGs) in LMC, SMC, M31, and M33 with $4.6\leq\log(L/L_{\odot})\leq5.091$. This supports that SN II progenitors correspond to RSGs. The conspicuous absence of SN II progenitors with $\log(L/L_{\odot})>5.1$ dex with respect to what is observed in RSG luminosity distributions, known as the RSG problem, is significant at a $5.2\pm0.5\,\sigma$ level.

We discuss the secondary cosmic microwave background (CMB) anisotropy due to kinetic Sunyaev-Zel'dovich (kSZ) effect from ionized bubbles around individual quasars prior to the reionization of the Universe. The bubbles create local ionization modulations which move with the large scale structure linear bulk flow and acts as source kSZ. We improve upon previous calculations of this effect, using a halo model based description of quasar abundance, and find that the kSZ distortion power spectrum, $C_\ell$, from the bubbles to be sub-dominant to kSZ from patchy reionization driven by galaxies. However, the shape of the two $C_\ell$'s are very different with the quasar bubble $C_\ell$ having a peak at $\ell \approx 500-700$ whereas the $C_\ell$ due to patchy reionization flattening out at $\ell > 1000$ thus making it plausible to separate the two using $C_\ell$ template-fitting in a future survey like CMB-HD. We also show that the amplitude of the quasar bubble induced power spectrum has a strong dependence on the epoch of reionization, and can be negligible for early reionization models. Next, we look at the imprint of a single quasar bubble on the CMB and show that it can be detected in a high-resolution, ambitious effort like CMB-HD. The amplitude and the oscillatory features of a single bubble power spectrum can be used to constrain the apparent asymmetric shape of the relativistically expanding ionized bubble, and depends on intrinsic quasar physics parameters such as the quasar photon emission rate and the quasar lifetime and cosmological probes like high redshift large scale linear velocity and the neutral hydrogen fraction. A detection of a high redshift quasar bubble in the CMB would carry complimentary information to its detection in HI or Lyman-$\alpha$ and a joint analysis can be used to break parameter degeneracies.

The General Antiparticle Spectrometer (GAPS) is an upcoming balloon mission to measure low-energy cosmic-ray antinuclei during at least three ~35-day Antarctic flights. With its large geometric acceptance and novel exotic atom-based particle identification, GAPS will detect ~500 cosmic antiprotons per flight and produce a precision cosmic antiproton spectrum in the kinetic energy range of ~0.07-0.21 GeV/n at the top of the atmosphere. With these high statistics extending to lower energies than any previous experiment, and with complementary sources of experimental uncertainty compared to traditional magnetic spectrometers, the GAPS antiproton measurement will be sensitive to dark matter, primordial black holes, and cosmic ray propagation. The antiproton measurement will also validate the GAPS antinucleus identification technique for the antideuteron and antihelium rare-event searches. This analysis demonstrates the GAPS sensitivity to cosmic-ray antiprotons using a full instrument simulation and event reconstruction, and including solar and atmospheric effects.

Using proper motions from Gaia Early Data Release 3 (Gaia EDR 3) and radial velocities from several surveys, we identify 60 candidate high-velocity stars with total velocity greater than 75\% escape velocity that probably origin from Sagittarius dwarf spheroidal galaxy (Sgr) by orbital analysis. Sgr's gravity has little effect on the results and the Large Magellanic Cloud's gravity has non-negligible effect on only a few stars. The closest approach of these stars to the Sgr occurs when the Sgr passed its pericenter ($\sim$ 38.2 Myr ago), which suggest they were tidally stripped from the Sgr. The positions of these stars in the HR diagram and the chemical properties of 19 of them with available [Fe/H] are similar with the Sgr stream member stars. This is consistent with the assumption of their accretion origin. Two of the 60 are hypervelocity stars, which may also be produced by Hills mechanism.

The Square Kilometre Array (SKA) project consists of delivering two largest radio telescope arrays being built by the SKA Observatory (SKAO), which is an intergovernmental organization bringing together nations from around the world with China being one of the major member countries. The computing resources needed to process, distribute, curate and use the vast amount of data that will be generated by the SKA telescopes are too large for the SKAO to manage on its own. To address this challenge, the SKAO is working with the international community to create a shared, distributed data, computing and networking capability called the SKA Regional Centre Alliance. In this model, the SKAO will be supported by a global network of SKA Regional Centres (SRCs) distributed around the world in its member countries to build an end-to-end science data system that will provide astronomers with high-quality science products. SRCs undertake deep processing, scientific analysis, and long-term storage of the SKA data, as well as user support. China has been actively participating in and promoting the construction of SRCs. This paper introduces the international cooperation and ongoing prototyping of the global SRC network, the construction plan of the China SRC and describes in detail the China SRC prototype. The paper presents examples of scientific applications of SKA precursor and pathfinder telescopes completed using resources from the China SRC prototype. Finally, the future prospects of the China SRC are presented.

Ultraviolet (UV) bursts and Ellerman bombs (EBs) are small-scale magnetic reconnection events taking place in the highly stratified, low solar atmosphere. It is still not clear whether UV bursts have to be generated at a higher atmospheric layer than EBs or whether both UV bursts and EBs can occur in the low chromosphere. We numerically studied the low $\beta$ magnetic reconnection process around the solar temperature minimum region (TMR). The time-dependent ionization degrees of hydrogen and helium are included in the MHD code, which lead to a more realistic magnetic diffusion caused by electron-neutral collision and ambipolar diffusion. A more realistic radiative cooling model from Carlsson & Leenaarts 2012 is included in the simulations. Our results in high resolution indicate that the plasmas in the reconnection region are heated up to more than $20,000$ K if the reconnecting magnetic field is as strong as $500$ G, which suggests that UV bursts can be generated in the dense low chromosphere. The dominant mechanism for producing the UV burst in the low chromosphere is heating, as a result of the local compression in the reconnection process. The thermal energy occurring in the reconnection region rapidly increases after the turbulent reconnection mediated by plasmoids is invoked. The average power density of the generated thermal energy in the reconnection region can reach over $1000$ erg cm$^{-3}$ s$^{-1}$, which is comparable to the average power density accounting for a UV burst. With the strength of the reconnecting magnetic field exceeding $900$ G, the width of the synthesized Si IV 1394 A line profile with multiple peaks can reach up to $100$ km s$^{-1}$, which is consistent with observations.

Our goal is to characterise the physical properties of the metal-poor brown dwarf population. In particular, we focus on the recently discovered peculiar dwarf WISE J1810055$-$1010023. We collected optical iz and near-infrared J-band imaging on multiple occasions over 1.5 years to derive accurate trigonometric parallax and proper motion of the metal-depleted ultra-cool dwarf candidate WISE1810. We also acquired low-resolution optical spectroscopy (0.6$-$1.0 $\mu$m) and new infrared (0.9$-$1.3 $\mu$m) spectra of WISE1810 that were combined with our photometry, other existing data from the literature and our trigonometric distance to determine the object's luminosity from the integration of the observed spectral energy distribution covering from 0.6 through 16$\mu$m. We compared the full optical and infrared spectrum with state-of-the-art atmosphere models to further constrain its effective temperature, surface gravity and metallicity. WISE1810 is detected in the $iz$ bands with AB magnitudes of $i$=23.871$\pm$0.104 and $z$=20.147$\pm$0.083 mag in the PanSTARRS system. It does not show any obvious photometric variability beyond 0.1$-$0.2 mag in any of the $z$- and $J$-band filters. The very red $z-J \approx 2.9$ mag colour is compatible with an ultra-cool dwarf nature. Fitting for parallax and proper motion, we measure a trigonometric parallax of 112.5$^{+8.1}_{-8.0}$ mas for WISE1810, placing the object at only 8.9$^{+0.7}_{-0.6}$ pc, about three times closer than previously thought. We employed Monte Carlo methods to estimate the error on the parallax and proper motion. The object's luminosity was determined at log$L/L_\odot$=$-$5.78$\pm$0.11 dex. From the comparison to atmospheric models, we infer a likely metallicity of [Fe/H] $\approx -1.5$ and an effective temperature cooler than 1000K. Abridged

We consider the acceleration of charged particles in relativistic shearing flows, with Lorentz factor up to $\Gamma_0 \sim 20$. We present numerical solutions to the particle transport equation and compare these with results from analytical calculations. We show that in the highly relativistic limit the particle energy spectrum that results from acceleration approaches a power law, $N(E)\propto E^{-\tilde{q}}$, with a universal value $\tilde{q}=(1+\alpha)$ for the slope of this power law, where $\alpha$ parameterizes the power-law momentum dependence of the particle mean free path. At mildly relativistic flow speeds, the energy spectrum becomes softer and sensitive to the underlying flow profile. We explore different flow examples, including Gaussian and power-law-type velocity profiles, showing that the latter yield comparatively harder spectra, producing $\tilde{q}\simeq 2$ for $\Gamma_0 \simeq 3$ and Kolmogorov turbulence. We provide a comparison with a simplified leaky-box approach and derive an approximate relation for estimating the spectral index as a function of the maximum shear flow speed. These results are of relevance for jetted, high-energy astrophysical sources such as active galactic nuclei, since shear acceleration is a promising mechanism for the acceleration of charged particles to relativistic energies and is likely to contribute to the high-energy radiation observed.

Double-lined eclipsing binaries allow the direct determination of masses and radii, which are key to test stellar models. With the launch of the TESS mission, many well-known eclipsing binaries have been observed at higher photometric precision, permitting the improvement of the absolute dimensions determinations. Using TESS data and newly-obtained spectroscopic observations, we aim at determining the masses and radii of the eccentric eclipsing binary systems V889 Aql and V402 Lac, together with their apsidal motion parameters. We modelled simultaneously radial velocity curves and times of eclipse for each target to precisely determine the orbital parameters of the systems, which we used to analyse the light curves and then obtain their absolute dimensions. We compared the obtained values with those predicted by theoretical models. We determined masses and radii of the components of both systems with relative uncertainties lower than 2%. V889 Aql is composed of two stars with masses $2.17\pm0.02$ M$_{\odot}$ and $2.13\pm0.01$ M$_{\odot}$ and radii $1.87\pm0.04$ R$_{\odot}$ and $1.85\pm0.04$ R$_{\odot}$. We found conclusive evidence of the presence of a third body orbiting V889 Aql with a period of 67 years. Based on the detected third light and the absence of signal in the spectra, we suggest that this third body could in turn be a binary composed by two $\sim$1.4 M$_{\odot}$ stars. V402 Lac is composed by two stars with masses $2.80\pm0.05$ M$_{\odot}$ and $2.78\pm0.05$ M$_{\odot}$ and radii $2.38\pm0.03$ R$_{\odot}$ and $2.36\pm0.03$ R$_{\odot}$. The times of minimum light are compatible with the presence of a third body for this system too, although its period is not yet fully sampled. In both cases we have found a good agreement between the observed apsidal motion rates and the model predictions.

Since the first gravitational wave (GW) event from binary black hole (BBH) was detected by LIGO-Virgo, GWs have become a useful probe on astrophysics and cosmology. If primordial black holes (PBHs) contribute a significant fraction of dark matter at wide mass range, they will cause microlensing in the GW signals with long wavelengths that are distinct from the lensing effects of electromagnetic signals from astrophysical objects. In this paper, we apply the lensing effect of GW from BBH to derive constraints on the abundance of PBHs for the Cosmic Explorer, a third-generation ground-based GW detector. We firstly consider two channels of formation of BBH that contribute to low and high redshift GW sources, including the astrophysical origin BBH scenario, and the primordial origin BBH scenario. Secondly, comparing with the method of optical depth, we use the Bayesian analysis to derive constraints on the abundance of PBHs with different mass function of lens taken into consideration. For a null search with $1000$ detected GW events of BBH, we find that the abundance of PBHs could be constrained to $\leq0.1\%$ in the mass range $\geq500~M_{\odot}$ at $68\%$ confidence level. In addition, if a GW event lensed by $M_{\rm PBH}=100-300~M_{\odot}$ is detected in $100$ detected GW events of BBH, we can derive that the abundance of PBHs is from $2.3\%$ to $25.2\%$ in this mass range with the Bayesian analysis.

We present a study of photochemical hazes of exoplanet atmospheres based on a self-consistent model including haze microphysics, disequilibrium chemistry, and radiative feedbacks. We derive the haze properties required to match HST observations of ten hot-Jupiters. HAT-P-12b, HD-189733b, HD-209458b and WASP-6b require haze mass fluxes between 5x10$^{-15}$ and 9x10$^{-12} g.cm^{-2}.s^{-1}$ to match the observations. WASP-12b and WASP-19b with equilibrium temperatures above 2000 K are incompatible with the presence of haze and are better fitted by heavy metals. HAT-P-1b and WASP-31b do not show clear evidence for the presence of hazes with upper mass fluxes of 10$^{-15}$ and 10$^{-16}g.cm^{-2}.s^{-1}$, respectively, while WASP-17b and WASP-39b present an upper mass flux limit of 10$^{-16} g.cm^{-2}.s^{-1}$. We discuss the implications of the self-consistent model and we derive upper limits for the haze abundances based on photochemistry results. Our results suggest HCN as the main haze precursor up to 1300 K effective temperatures and CO above. Our derived haze mass fluxes based on the fit to the observations are consistent with the photochemistry with formation yields up to $\sim$6.4\%. Disequilibrium chemistry has negligible impact on the spectra considering the low resolution observations used but impacts the chemical composition and temperature profiles. We find that hazes produce hotter upper atmosphere temperatures with a detectable impact on the spectra. Clouds may have implications for interpreting the transit spectra of HD-209458b, WASP-31b and WASP-39b. Nevertheless, the presence of silicate and iron clouds is expected in all studied atmospheres except WASP-12b and WASP-19b.

Simplified assumptions about the thermodynamics of the electrons are normally employed in general-relativistic magnetohydrodynamic (GRMHD) simulations of accretion onto black holes. To counter this, we have developed a self-consistent approach to study magnetised and radiatively cooled two-temperature accretion flows around a Kerr black hole in two spatial dimensions. The approach includes several heating processes, radiative cooling, and a coupling between the electrons and the ions via Coulomb interaction. We test our approach by performing axisymmetric GRMHD simulations of magnetised tori accreting onto a Kerr black hole under various astrophysical scenarios. In this way, we find that the inclusion of the Coulomb interaction and the radiative cooling impacts the thermodynamical properties of both the ions and electrons, changing significantly the temperature distribution of the latter. Since the electrons are responsible for the electromagnetic emission from these accretion flows, our results underline the importance of a two-temperature approach when imaging these flows. In addition, we find that the accretion rate influences the bulk properties of the flow and the thermodynamics of the electrons and ions. Interestingly, we observe qualitatively distinct temperature properties for SANE and MAD accretion modes while maintaining the same accretion rates, which could help distinguishing MAD and SANE accretion flows via observations. Finally, we propose two new relations for the temperature ratios of the electrons, ions, and of the gas in terms of the plasma-$\beta$ parameter. The new relations not only provide a more accurate description of the thermodynamics of the accretion flow in all relevant regimes, but they also represent a simple and effective approach to treat two-temperature accretion flows on supermassive black holes such as Sgr A* and M\,87*.

We analyze Multi Unit Spectroscopic Explorer observations of nine low-redshift (z < 0.1) Palomar-Green quasar host galaxies to investigate the spatial distribution and kinematics of the warm, ionized interstellar medium, with the goal of searching for and constraining the efficiency of active galactic nucleus (AGN) feedback. After separating the bright AGN from the starlight and nebular emission, we use pixel-wise, kpc-scale diagnostics to determine the underlying excitation mechanism of the line emission, and we measure the kinematics of the narrow-line region (NLR) to estimate the physical properties of the ionized outflows. The radial size of the NLR correlates with the AGN luminosity, reaching scales of $\sim 5\,$kpc and beyond. The geometry of the NLR is well-represented by a projected biconical structure, suggesting that the AGN radiation preferably escapes through the ionization cone. We find enhanced velocity dispersions ($\sim 100\,$km$\,$s$^{-1}$) traced by the H$\alpha$ emission line in localized zones within the ionization cones. Interpreting these kinematic features as signatures of interaction between an AGN-driven ionized gas outflow and the host galaxy interstellar medium, we derive mass outflow rates of $\sim 0.008-1.6\, M_\odot \,$yr$^{-1}$ and kinetic injection rates of $\sim 10^{39}-10^{42} \,$erg$\,$s$^{-1}$, which yield extremely low coupling efficiencies of $\lesssim 10^{-3}$. These findings add to the growing body of recent observational evidence that AGN feedback is highly ineffective in the host galaxies of nearby AGNs.

A cosmological model with a time-varying mass of electrons seems a promising solution for the so-called Hubble tension. We examine the big bang nucleosynthesis (BBN) constraints on the time-varying electron mass model, because a larger electron mass gives rise to the smaller neutron decay rate which could affect the light element abundance. Additionally, different inferred cosmological parameters, primarily baryon asymmetry, to keep the cosmic background power spectrum unchanged could affect the abundance of light element. We find that the predicted helium fraction becomes larger and the deuterium abundance becomes smaller as the electron mass at the BBN time becomes larger. Thus, we conclude that an acceptable electron mass at the BBN time would be only approximately 1% greater than the current electron mass.

Using Multi-Unit Spectroscopic Explorer (MUSE) spectroscopy, we analyse the stellar kinematics of 18 brightest group early-type (BGEs) galaxies, selected from the Complete Local-Volume Groups Sample (CLoGS). We analyse the kinematic maps for distinct features, and measure specific stellar angular momentum within one effective radius ($\lambda_{e}$). We classify the BGEs as fast (10/18) or slow (8/18) rotators, suggesting at least two different evolution paths. We quantify the anti-correlation between higher-order kinematic moment $h_{3}$ and V/$\sigma$ (using the $\xi_{3}$ parameter), and the kinematic misalignment angle between the photometric and kinematic position angles (using the $\Psi$ parameter), and note clear differences between these parameter distributions of the fast and slow rotating BGEs. We find that all 10 of our fast rotators are aligned between the morphological and kinematical axis, consistent with an oblate galaxy shape, whereas the slow rotators are spread over all three classes: oblate (1/8), triaxial (4/8), and prolate (3/8). We place the results into context using known radio properties, X-ray properties, and observations of molecular gas. We find consistent merger histories inferred from observations for the fast-rotating BGEs, indicating that they experienced gas-rich mergers or interactions, and these are very likely the origin of the cold gas. Observational evidence for the slow rotators are consistent with gas-poor mergers. For the slow rotators with cold gas, all evidence point to cold gas cooling from the intragroup medium.

We present the first dynamical simulation that recreates the major properties of the archetypal nearby spiral galaxy M101. Our model describes a grazing but relatively close (14 kpc) passage of the companion galaxy NGC 5474 through M101's outer disk approximately 200 Myr ago. The passage is retrograde for both disks, which yields a relatively strong gravitational response while suppressing the formation of long tidal tails. The simulation reproduces M101's overall lopsidedness, as well as the extended NE Plume and sharp western edge of the galaxy's disk. The post-starburst populations observed in M101's NE Plume are likely a result of star formation triggered at the point of contact where the galaxies collided. Over time, this material will mix azimuthally, leaving behind diffuse, kinematically coherent stellar streams in M101's outer disk. At late times after the encounter, the density profile of M101's disk shows a broken "upbending" profile similar to those seen in spiral galaxies in denser environments, further demonstrating the connection between interactions and long-term structural changes in galaxy disks.

Planetary systems such as our own are formed after a long process where matter condenses from diffuse clouds to stars, planets, asteroids, comets and residual dust, undergoing dramatic changes in physical and chemical state in less than a few million years. Several studies have shown that the chemical composition during the early formation of a Solar-type planetary system is a powerful diagnostic to track the history of the system itself. Among the approximately 270 molecules so far detected in the ISM, the so-called interstellar complex organic molecules (iCOMs) are of particular interest both because of their evolutionary diagnostic power and because they might be potential precursors of biomolecules, which are at the basis of terrestrial life. This Chapter focuses on the evolution of organic molecules during the early stages of a Solar-type planetary system, represented by the prestellar, Class 0/I and protoplanetary disk phases, and compares them with what is observed presently in Solar System comets. Our twofold goal is to review the processes at the base of organic chemistry during Solar-type star formation and, in addition, to possibly provide constraints on the early history of our own planetary system.

Photometric data-driven classification of supernovae becomes a challenge due to the appearance of real-time processing of big data in astronomy. Recent studies have demonstrated the superior quality of solutions based on various machine learning models. These models learn to classify supernova types using their light curves as inputs. Preprocessing these curves is a crucial step that significantly affects the final quality. In this talk, we study the application of multilayer perceptron (MLP), bayesian neural network (BNN), and normalizing flows (NF) to approximate observations for a single light curve. We use these approximations as inputs for supernovae classification models and demonstrate that the proposed methods outperform the state-of-the-art based on Gaussian processes applying to the Zwicky Transient Facility Bright Transient Survey light curves. MLP demonstrates similar quality as Gaussian processes and speed increase. Normalizing Flows exceeds Gaussian processes in terms of approximation quality as well.

Cataloguing and classifying celestial objects is one of the fundamental activities of observational astrophysics. In this work, we compare the contents of two comprehensive databases, the NASA Extragalactic Database (NED) and Set of Identifications, Measurements and Bibliography for Astronomical Data (SIMBAD) in the vicinity of nearby galaxies. These two databases employ different classification schemes -- one flat and one hierarchical -- and our goal was to determine the compatibility of classifications for objects in common. Searching both databases for objects within the respective isophotal radius of each of the ~1300 individual galaxies in the Local Volume Galaxy sample, we found that on average, NED contains about ten times as many entries as SIMBAD and about two thirds of SIMBAD objects are matched by position to a NED object, at 5 arcsecond tolerance. These quantities do not depend strongly on the properties of the parent galaxies. We developed an algorithm to compare individual object classifications between the two databases and found that 88% of the classifications agree; we conclude that NED and SIMBAD contain consistent information for sources in common in the vicinity of nearby galaxies. Because many galaxies have numerous sources contained only in one of NED or SIMBAD, researchers seeking the most complete picture of an individual galaxy's contents are best served by using both databases.

We study the observational implications of a class of inflationary models wherein the inflaton is coupled to the Einstein tensor through a generalized non-minimal derivative coupling (GNMDC). Such a coupling can be realized in the framework of Horndeski theories or generalized Galileon theories and leads to novel and distinguishable inflationary predictions. In particular, we explore whether such models can provide a possible explanation to the large scale anomalies such as the power suppression and other localized features associated with the CMB temperature and polarisation anisotropies at low multipoles or large angular scales. For a specific choice of the GNMDC coupling function, we find that these models can lead to suitable localized features in the power spectrum on large scales. To our knowledge, such models have not been used earlier to explain the origin of these modulations in the CMB power spectrum. We work in the regime of parameter space such that we avoid the gradient instability and the superluminal propagation of scalar perturbations. A very interesting aspect of our analysis is that a class of inflationary models such as the hilltop-quartic model results in a better agreement with the Planck data in the presence of an additional GNMDC term. Further, we compare the GNMDC model with the data using CosmoMC and find that these models provide a considerable improvement over the best fit reference $\Lambda$CDM model with a featureless, power law, primordial spectrum. Future CMB experiments such as CMB-S4 will certainly impose better constraints on various parameters of the GNMDC scenario. Finally, we discuss the wider implications of our results and argue that these models may also be useful in providing an explanation to other anomalies associated with the CMB observations.

If primordial black holes (PBH) within the asteroid/moon-mass range indeed reside in galactic dark matter halos, they must necessarily collide with galactic neutron stars (NS). These collisions must, again necessarily, form light black holes (LBHs) with masses of typical NS, $M_{\rm LBH} \approx \,1.5\,M_{\odot}$. LBHs may be behind events already detected by ground-based gravitational-wave detectors (GW170817, GW190425, and others such as a mixed stellar black hole-neutron star mass event GW191219\_163120), and most recently by microlensing (OGLE-BLG-2011-0462). Although the status of these observations as containing LBHs is not confirmed, there is no question that gravitational-wave detectors and microlensing are in principle and in practice capable of detecting LBHs. We have calculated the creation rate of LBHs resulting from these light primordial black hole collisions with neutron stars. On this basis, we claim that if improved gravitational-wave detectors and microlensing statistics of the LBH events would indicate that the number of LBHs is significantly lower that what follows from the calculated creation rate, then this would be an unambiguous proof that there is no significant light PBH contribution to the galactic dark matter halos. Otherwise, if observed and calculated numbers of LBHs roughly agree, then the hypothesis of primordial black hole existence gets strong observational support, and in addition their collisions with neutron stars may be considered a natural creation channel for the LBHs, solving the problem of their origin, as it is known that they cannot be a product of standard stellar evolution.

The ALMA interferometer, with its unprecedented combination of high-sensitivity and high-angular resolution, allows for (sub-)mm wavelength mapping of protostellar systems at Solar System scales. Astrochemistry has benefited from imaging interstellar complex organic molecules in these jet-disk systems. Here we report the first detection of methanol (CH3OH) and methyl formate (HCOOCH3) emission towards the triple protostellar system VLA1623-2417 A1+A2+B, obtained in the context of the ALMA Large Program FAUST. Compact methanol emission is detected in lines from Eu = 45 K up to 61 K and 537 K towards components A1 and B, respectively. LVG analysis of the CH3OH lines towards VLA1623-2417 B indicates a size of 0.11-0.34 arcsec (14-45 au), a column density N(CH3OH) = 10^16-10^17 cm-2, kinetic temperature > 170 K, and volume density > 10^8 cm-3. An LTE approach is used for VLA1623-2417 A1, given the limited Eu range, and yields Trot < 135 K. The methanol emission around both VLA1623-2417 A1 and B shows velocity gradients along the main axis of each disk. Although the axial geometry of the two disks is similar, the observed velocity gradients are reversed. The CH3OH spectra from B shows two broad (4-5 km s-1) peaks, which are red- and blue-shifted by about 6-7 km s-1 from the systemic velocity. Assuming a chemically enriched ring within the accretion disk, close to the centrifugal barrier, its radius is calculated to be 33 au. The methanol spectra towards A1 are somewhat narrower (about 4 km s-1), implying a radius of 12-24 au.

The atmospheres of exoplanets harbor critical information about their habitability. However, extracting and interpreting that information requires both high-quality spectroscopic data and a comparative analysis to characterize the findings. Looking forward to data availability, we propose a novel, assumption-free approach adapting the Jensen-Shannon divergence ($\mathcal{D}_{JS}$) information measure to identify Earth-like planets through their transmission spectra. We apply this method to simulated Earth-like and Jupiter-like planets, including high-interest observed exoplanets such as Trappist-1e and GJ 667 Cc, and demonstrate that $\mathcal{D}_{JS}$ can discriminate between different planet types. We argue that this method can be used to identify habitable and even inhabited planets as more precise transit spectroscopy data becomes available in the coming years.

We present a comprehensive study of the ionization structure and kinematics of the planetary nebula (PN) NGC40 (a.k.a. the Bow-tie Nebula). A set of narrow-band images obtained with the ALhambra Faint Object Spectrograph and Camera (ALFOSC) at the Nordic Optical Telescope (NOT) are used to study the turbulent distribution of gas in the main cavity, the ionization stratification and the density of this PN. High-resolution Manchester Echelle Spectrograph (MES) observations obtained at 2.1m telescope of the San Pedro M\'{a}rtir (SPM) Observatory are used to unveil in great detail the kinematic signatures of all morphological features in NGC40. The images and spectra suggest that NGC40 had multiple mass ejections in its recent formation history. We found 4 jet-like ejections not aligned with the main axis of NGC40 (PA=20$^{\circ}$), some of them having pierced the main cavity along the SW-NE direction as well as the southern lobe. Using a tailor-made morpho-kinematic model of NGC40 produced with SHAPE we found that the main cavity has a kinematic age of 6,500 yr and the two pairs of lobes expanding towards the N and S directions have an averaged age of 4,100$\pm$550 yr. NGC40 thus adds to the group of PNe with multiple ejections along different axes that challenge the models of PN formation.

Context: The far-infrared (FIR) and sub-millimeter (submm) emissivity of the Milky Way (MW) cirrus is an important benchmark for dust grain models. Dust masses in other galaxies are generally derived from the FIR/submm using the emission properties of these MW-calibrated models. Aims: We seek to derive the FIR/submm emissivity in nine nearby spiral galaxies to check its compatibility with MW cirrus measurements. Methods: We obtained values of the emissivity at 70 to 500 um, using maps of dust emission from the Herschel satellite and of gas surface density from the THINGS and HERACLES surveys on a scale generally corresponding to 440 pc. We studied the variation of the emissivity with the surface brightness ratio I(250um)/I(500um), a proxy for the intensity of the interstellar radiation field heating the dust. Results: We find that the average value of the emissivity agrees with MW estimates for pixels sharing the same color as the cirrus, namely, for I(250um)/I(500um) = 4.5. For I(250um)/I(500um) > 5, the measured emissivity is instead up to a factor ~2 lower than predicted from MW dust models heated by stronger radiation fields. Regions with higher I(250um)/I(500um) are preferentially closer to the galactic center and have a higher overall (stellar+gas) surface density and molecular fraction. The results do not depend strongly on the adopted CO-to-molecular conversion factor and do not appear to be affected by the mixing of heating conditions. Conclusions: Our results confirm the validity of MW dust models at low density, but are at odds with predictions for grain evolution in higher density environments. If the lower-than-expected emissivity at high I(250um)/I(500um) is the result of intrinsic variations in the dust properties, it would imply an underestimation of the dust mass surface density of up to a factor ~2 when using current dust models.

We examine the observed properties of the Nili Fossae olivine-clay-carbonate lithology from orbital data and in situ by the Mars 2020 rover at the S\'e\'itah unit in Jezero crater, including: 1) composition (Liu, 2022) 2) grain size (Tice, 2022) 3) inferred viscosity (calculated based on geochemistry collected by SuperCam (Wiens, 2022)). Based on the low viscosity and distribution of the unit we postulate a flood lava origin for the olivine-clay-carbonate at S\'e\'itah. We include a new CRISM map of the clay 2.38 {\mu}m band and use in situ data to show that the clay in the olivine cumulate in the S\'e\'itah formation is consistent with talc or serpentine from Mars 2020 SuperCam LIBS and VISIR and MastCam-Z observations. We discuss two intertwining aspects of the history of the lithology: 1) the emplacement and properties of the cumulate layer within a lava lake, based on terrestrial analogs in the Pilbara, Western Australia, and using previously published models of flood lavas and lava lakes, and 2) the limited extent of post emplacement alteration, including clay and carbonate alteration (Clave, 2022; Mandon, 2022).

We present a first theoretical modeling of joint Parker Solar Probe (PSP) - Metis/Solar Orbiter (SolO) quadrature observations Telloni et al 2022c. The combined observations describe the evolution of a slow solar wind plasma parcel from the extended solar corona ($3.5-6.3$ R$_\odot$) to the very inner heliosphere (23.2 R$_\odot$). The Metis/SolO instrument remotely measures the solar wind speed finding a range from $96-201$ kms$^{-1}$, and PSP measures the solar wind plasma in situ, observing a radial speed of 219.34 kms$^{-1}$. We find theoretically and observationally that the solar wind speed accelerates rapidly within 3.3 -- 4 R$_\odot$, and then increases more gradually with distance. Similarly, we find that the theoretical solar wind density is consistent with the remotely and in situ observed solar wind density. The normalized cross-helicity and normalized residual energy observed by PSP are 0.96 and -0.07, respectively, indicating that the slow solar wind is very Alfv\'enic. The theoretical NI/slab results are very similar to PSP measurements, which is a consequence of the highly magnetic field-aligned radial flow ensuring that PSP can measure slab fluctuations and not 2D. Finally, we calculate the theoretical 2D and slab turbulence pressure, finding that the theoretical slab pressure is very similar to that observed by PSP.

It is speculated that the high star-formation efficiency observed in spiral-arm molecular clouds is linked to the prevalence of compressive (curl-free) turbulent modes, while the shear-driven solenoidal (divergence-free) modes appear to be the main cause of the low star-formation efficiency that characterises clouds in the Central Molecular Zone. Similarly, analysis of the Orion B molecular cloud has confirmed that, although turbulent modes vary locally and at different scales within the cloud, the dominant solenoidal turbulence is compatible with its low star formation rate. This evidence points to inter-and intra-cloud fluctuations of the solenoidal modes being an agent for the variability of star formation efficiency. We present a quantitative estimation of the relative fractions of momentum density in the solenoidal modes of turbulence in a large sample of plane molecular clouds in the \ce{^{13}CO}/\ce{C^{18}O} ($J=3\rightarrow 2$) Heterodyne Inner Milky Way Plane Survey (CHIMPS). We find a negative correlation between the solenoidal fraction and star-formation efficiency. This feature is consistent with the hypothesis that solenoidal modes prevent or slow down the collapse of dense cores. In addition, the relative power in the solenoidal modes of turbulence (solenoidal fraction) appears to be higher in the Inner Galaxy declining with a shallow gradient with increasing Galactocentric distance. Outside the Inner Galaxy, the slowly, monotonically declining values suggest that the solenoidal fraction is unaffected by the spiral arms.

Recent high-spatial-resolution observations have revealed dust substructures in protoplanetary disks such as rings and gaps, which do not always correlate with gas. Because radial gas flow induced by low-mass, non-gas-gap-opening planets could affect the radial drift of dust, it potentially forms these dust substructures in disks. We investigate the potential of gas flow induced by low-mass planets to sculpt the rings and gaps in the dust profiles. We first perform three-dimensional hydrodynamical simulations, which resolve the local gas flow past a planet. We then calculate the trajectories of dust influenced by the planet-induced gas flow. Finally, we compute the steady-state dust surface density by incorporating the influences of the planet-induced gas flow into a one-dimensional dust advection-diffusion model. The outflow of the gas toward the outside of the planetary orbit inhibits the radial drift of dust, leading to dust accumulation (the dust ring). The outflow toward the inside of the planetary orbit enhances the inward drift of dust, causing dust depletion around the planetary orbit (the dust gap). Under weak turbulence ($\alpha_{\rm diff}\lesssim10^{-4}$, where $\alpha_{\rm diff}$ is the turbulence strength parameter), the gas flow induced by the planet with $\gtrsim1\,M_{\oplus}$ (Earth mass) generates the dust ring and gap in the distribution of small dust grains ($\lesssim1$ cm) with the radial extent of $\sim1\text{--}10$ times gas scale height around the planetary orbit without creating a gas gap and pressure bump. The gas flow induced by low-mass, non-gas-gap-opening planets can be considered a possible origin of the observed dust substructures in disks. Our results may be helpful to explain the disks whose dust substructures were found not to correlate with those of the gas.

We present an analytical proof assisted by computer calculations for the dynamical stability of the eight main planets and Pluto for the next 100,000 years. It means that the semi-major axes of the planets will not change significantly during this period. Also the eccentricities and inclinations of the orbits will remain sufficiently small. A standard linear four-step numerical method is used to integrate approximately the orbits of Mercury, Venus, Earth, Mars, Jupiter, Saturn, Uranus, Neptune and Pluto. Written in orbital elements, the dynamics of the nine planets manifests a system of 54 first-order ordinary differential equations. The step-size of the numerical method -- about six days, has been performed 6,290,000 times. We estimate the total accumulation of rounding-off errors, deviations related to possible uncertainty in the astronomical data and the accuracy of the computer calculations.

Nano-diamonds remain an intriguing component of the dust in the few sources where they have been observed in emission. This work focusses on the nano-diamonds observed in circumstellar discs and is an attempt to derive critical information about their possible sizes, compositions, and evolution using a recently-derived set of optical constants. The complex indices of refraction of nano-diamonds and their optical properties (the efficiency factors Qext, Qsca, Qabs, and Qpr) were used to determine their temperatures, lifetimes, and drift velocities as a function of their radii (0.5-100 nm), composition (surface hydrogenation and irradiated states), and distance from the central stars in circumstellar regions. The nano-diamond temperature profiles were determined for the stars HR 4049, Elias 1, and HD 97048 in the optically-thin limit. The results indicate that large nano-diamonds (a = 30 - 100 nm) are the hottest and therefore the least resistant in the inner disc regions (~ 10-50 AU), while small (a < 10 nm) fully-hydrogenated nano-diamonds remain significantly cooler in these same regions. We discuss these results within the context of nano-diamond formation in circumstellar discs. Large nano-diamonds, being the hottest, are most affected by the stellar radiation field, however, the effects of radiation pressure appear to be insufficient to move them out of harm's way. The nano-diamonds that best survive and therefore shine in the inner regions of proto-planetary discs are then seemingly small (a < 10 nm), hydrogenated, and close in size to pre-solar nano-diamonds (a ~ 1.4 nm). Nevertheless, it does not yet appear possible to reconcile their existence there with their seemingly short lifetimes in such regions.

Red Clump stars have been found to be enhanced in lithium relative to stars at the tip of the Red Giant Branch (TRGB), which is unexpected in current stellar models. At the TRGB, stars undergo the helium flash, during which helium burning briefly generates roughly $10^9 \, L_\odot$ of power and drives vigorous convection within the star's core. The helium-burning shell excites large fluxes of internal gravity waves. Here we investigate whether or not these waves can deposit enough heat to destabilize the hydrogen-burning shell, generate a convection zone there, and thereby drive the Cameron-Fowler process to enhance surface $^{7}\mathrm{Li}$. We study this with detailed stellar evolution models, and find that while the waves deposit $\sim 10^6 \, L_\odot$ near the hydrogen-burning shell, and while this generally does produce a convection zone, the resulting convection does not reach high enough to merge with the envelope, and so cannot explain enhancements to surface $^{7}\mathrm{Li}$.

A promising energy range to look for angular correlation between cosmic rays of extragalactic origin and their sources is at the highest energies, above few tens of EeV ($1\:{\rm EeV}\equiv 10^{18}\:$eV). Despite the flux of these particles being extremely low, the area of ${\sim}\:3{,}000 \: \text{km}^2$ covered at the Pierre Auger Observatory, and the 17-year data-taking period of the Phase 1 of its operations, have enabled us to measure the arrival directions of more than 2,600 ultra-high energy cosmic rays above $32\:\text{EeV}$. We publish this data set, the largest available at such energies from an integrated exposure of $122{,}000 \: \text{km}^2\:\text{sr}\:\text{yr}$, and search it for anisotropies over the $3.4\pi$ steradians covered with the Observatory. Evidence for a deviation in excess of isotropy at intermediate angular scale, with ${\sim}\:15^\circ$ Gaussian spread or ${\sim}\:25^\circ$ top-hat radius, is obtained at the $4\:\sigma$ significance level for cosmic-ray energies above ${\sim}\:40\:\text{EeV}$.

Stars grazing supermassive black holes (SMBHs) on bound orbits may survive tidal disruption, causing periodic flares. Inspired by the recent discovery of the periodic nuclear transient ASASSN-14ko, a promising candidate for a repeating tidal disruption event (TDE), we study the tidal deformation of stars approaching SMBHs on eccentric orbits. With both analytical and hydrodynamics methods, we show the overall tidal deformation of a star is similar to that in a parabolic orbit provided that the eccentricity is above a critical value. This allows one to make use of existing simulation libraries from parabolic encounters to calculate the mass fallback rate in eccentric TDEs. We find the flare structures of eccentric TDEs show a complicated dependence on both the SMBH mass and the orbital period. For stars orbiting SMBHs with relatively short periods, we predict significantly shorter-lived duration flares than those in parabolic TDEs, which can be used to predict repeating events if the mass of the SMBH can be independently measured. Using an adiabatic mass loss model, we study the flare evolution over multiple passages, and show the evolved stars can survive many more passages than main sequence stars. We apply this theoretical framework to the repeating TDE candidate ASASSN-14ko and suggest that its recurrent flares originate from a moderately massive ($M\gtrsim 1\,\mathrm{M_\odot}$), extended ($\approx$ a few $\mathrm{R_\odot}$), evolved star on a grazing, bound orbit around the SMBH. Future hydrodynamics simulations of multiple tidal interactions will enable realistic models on the individual flare structure and the evolution over multiple flares.

Measurement of spin-precession in black hole binary mergers observed with gravitational waves is an exciting milestone as it relates to both general relativistic dynamics and astrophysical binary formation scenarios. In this study, we revisit the evidence for spin-precession in GW200129 and localize its origin to data in LIGO Livingston in the 20--50\,Hz frequency range where the signal amplitude is lower than what expected from a non-precessing binary given all the other data. These data are subject to known data quality issues as a glitch was subtracted from the detector's strain data. The lack of evidence for spin-precession in LIGO Hanford leads to a noticeable inconsistency between the inferred binary mass ratio and precessing spin in the two LIGO detectors, something not expected from solely different Gaussian noise realizations. We revisit the LIGO Livingston glitch mitigation and show that the difference between a spin-precessing and a non-precessing interpretation for GW200129 is smaller than the statistical and systematic uncertainty of the glitch subtraction, finding that the support for spin-precession depends sensitively on the glitch modeling. We also investigate the signal-to-noise ratio $\sim7$ trigger in the less sensitive Virgo detector. Though not influencing the spin-precession studies, the Virgo trigger is grossly inconsistent with the ones in LIGO Hanford and LIGO Livingston as it points to a much heavier system. We interpret the Virgo data in the context of further data quality issues. While our results do not disprove the presence of spin-precession in GW200129, we argue that any such inference is contingent upon the statistical and systematic uncertainty of the glitch mitigation. Our study highlights the role of data quality investigations when inferring subtle effects such as spin-precession for short signals such as the ones produced by high-mass systems.

We study the formation and evolution of vortices in $U(1)$ dark photon dark matter and dark photon clouds that arise through black hole superradiance. We show how the production of both longitudinal mode and transverse mode dark photon dark matter can lead to the formation of vortices. After vortex formation, the energy stored in the dark photon dark matter will be transformed into a large number of vortex strings. In the case where a dark photon magnetic field is produced, bundles of vortex strings are formed in a superheated phase transition, and evolve towards a configuration consisting of many string loops that are uncorrelated on large scales, analogous to a melting phase transition in condensed matter. In the process, they dissipate via dark photon and gravitational wave emission, offering a target for experimental searches. Vortex strings were also recently shown to form in dark photon superradiance clouds around black holes, and we discuss the dynamics and observational consequences of this phenomenon with phenomenologically motivated parameters. In that case, the string loops ejected from the superradiance cloud, apart from producing gravitational waves, are also quantised magnetic flux lines and can be looked for with magnetometers. We discuss the connection between the dynamics in these scenarios and similar vortex dynamics found in type II superconductors.

The dressed metric and the hybrid approach to perturbations are the two main approaches to capture the effects of quantum geometry in the primordial power spectrum in loop quantum cosmology. Both consider Fock quantized perturbations over a loop quantized background and result in very similar predictions except for the modes which exit the horizon in the effective spacetime in the Planck regime. Understanding precise relationship between both approaches has so far remained obscured due to differences in construction and technical assumptions. We explore this issue at the classical and effective spacetime level for linear perturbations, ignoring backreaction, which is the level at which practical computations of the power spectrum in both of the approaches have so far been performed. We first show that at the classical level both the approaches lead to the same Hamiltonian up to the second order in perturbations and result in the same classical mass functions in the Mukhanov-Sasaki equation on the physical solutions. At the effective spacetime level, the difference in phenomenological predictions between the two approaches in the Planck regime can be traced to whether one uses the Mukhanov-Sasaki variable $Q_{\vec k}$ (the dressed metric approach) or its rescaled version $\nu_{\vec k}=aQ_{\vec k}$ (the hybrid approach) to write the Hamiltonian of the perturbations, and associated polymerization ambiguities. It turns out that if in the dressed metric approach one chooses to work with $\nu_{\vec{k}}$, the effective mass function can be written exactly as in the hybrid approach, thus leading to identical phenomenological predictions in all regimes. Our results explicitly show that the dressed metric and the hybrid approaches for linear perturbations, at a practical computational level, can be seen as two sides of the same coin.

Multiband observations of coalescing stellar-mass black holes binaries could deliver valuable information on the formation of those sources and potential deviations from General Relativity. Some of these binaries might be first detected by the space-based detector LISA and, then, several years later, observed with ground-based detectors. Due to large uncertainties in astrophysical models, it is hard to predict the population of such binaries that LISA could observe. In this work, we assess the ability of LISA to detect the events of the third catalogue of gravitational wave sources released by the LIGO/Virgo/KAGRA collaboration. We consider the possibility of directly detecting the source with LISA and performing archival searches in the LISA data stream, after the event has been observed with ground-based detectors. We also assess how much could LISA improve the determination of source parameters. We find that it is not guaranteed that any event other than GW150914 would have been detected. Nevertheless, if any event is detected by LISA, even with a very low signal-to-noise ratio, the measurement of source parameters would improve by combining observations of LISA and ground based detectors, in particular for the chirp mass.

The IKKT matrix model yields an emergent space-time. We further develop these ideas and give a proposal for an emergent metric. Based on previous numerical studies of this model, we provide evidence that the emergent space-time is continuous and infinite in extent, both in space and in time, and that the metric is spatially flat. The time evolution describes the transition from a string-theoretic emergent phase to a phase in which the $SO(9)$ symmetry of the model is spontaneously broken to $SO(6) \times SO(3)$, with three dimensions of space expanding, becoming classical and at later times evolving like in a radiation-dominated universe, and the remaining six dimensions of space stabilized at the string scale. We speculate on how this analysis can be extended to yield an early universe cosmology which, in addition to the above-mentioned properties, also leads to a roughly scale-invariant spectrum of cosmological fluctuations and gravitational waves.

We study the cosmological inflation and dark matter (DM) in a unified way within a $Z_3$ complex scalar model. The real and imaginary parts of the complex scalar act as the inflaton and DM respectively. The slow-rolling inflation with non-minimal coupling in both the metric and Palatini formalisms can be realized. We examine the whole parameters space by fully considering the theoretical and experimental constraints. We find that in the low-energy scale, the DM relic density and the DM-nucleon direct scattering experiments favor the mixing angle $|\theta| \lesssim 0.25$, the DM mass $m_\chi \gtrsim 80\rm{GeV}$, and the mass of Higgs-like scalar $m_{h_2} \gtrsim 300\rm{GeV}$. In the high-energy scale, after further considering the cosmological constraints of the scalar spectral index and the tensor-to-scalar ratio for the two forms of inflation, the scalar spectral indices are both $\sim 0.965$, the non-minimum coupling coefficients are $\sim 10^4$ and $\sim 10^9$, and the tensor-to-scalar ratios are $\sim 10^{-3}$ and $\lesssim 10^{-11}$ respectively, which suggests that the inflation under the two formalisms can be distinguished by measuring the tensor-to-scalar ratio with higher precision.

For the inflaton field we determine a new exact solution by using the Lie symmetry analysis. Specifically, we construct a second-order differential master equation for arbitrary scalar field potential by assuming that the spectral index for the density perturbations $n_{s}$ and the scalar to tensor ratio $r$ are related as $n_{s}-1=h\left( r\right) $. Function $h\left( r\right) $ is classified according to the admitted Lie symmetries for the master equation. The possible admitted Lie symmetries form the $A_{2}$, $A_{3,2}$, $A_{3,3}$ and $sl\left( 3,R\right) $ Lie algebras. The new inflationary solution is recovered by the Lie symmetries of the $A_{3,3}$ algebra. Scalar field potential is derived explicitly, while we compare the resulting spectral indices with the observations.

We apply singularity analysis to investigate the integrability properties of the gravitational field equations in Weyl Integrable Spacetime for a spatially flat Friedmann--Lema\^itre--Robertson--Walker background spacetime induced with an ideal gas. We find that the field equations possess the Painlev\'e property in the presence of the cosmological constant and the analytic solution is given by a left Laurent expansion.

The experimental investigation of matter under extreme densities and temperatures as they occur for example in astrophysical objects and nuclear fusion applications constitutes one of the most active frontiers at the interface of material science, plasma physics, and engineering. The central obstacle is given by the rigorous interpretation of the experimental results, as even the diagnosis of basic parameters like the temperature T is rendered highly difficult by the extreme conditions. In this work, we present a simple, approximation-free method to extract the temperature of arbitrarily complex materials from scattering experiments, without the need for any simulations or an explicit deconvolution. This new paradigm can be readily implemented at modern facilities and corresponding experiments will have a profound impact on our understanding of warm dense matter and beyond, and open up a gamut of appealing possibilities in the context of thermonuclear fusion, laboratory astrophysics, and related disciplines.

Timing observations from the Nanshan 26-m radio telescope for nine pulsars between 2000 and 2014 have been used to search for glitches. The data span for nine pulsars ranges from 11.6 to 14.2 years. From the total of 114 yr of pulsar rotational history, 16 new glitches were identified in 9 pulsars. Glitch parameters were measured by fitting the timing residuals data. All 16 glitches have a small fractional size. Six new glitches have been detected in PSR J1833-0827, making it another frequent glitching pulsar. Some of the 16 glitches may experience exponential or linear recovery, but it is unlikely for us to make further analyses with the large gap in the data set. All the glitch rates obtained from Nanshan are higher than that from Jodrell Bank Observatory. The small glitch size and high glitch rate could possibly attribute to the high observation cadence.

The generation of gravitational waves from a post-Newtonian source endowed with a quantum spin, modeled by the Weyssenhoff fluid, is investigated in the context of Einstein-Cartan theory at the first post-Newtonian level by resorting to the Blanchet-Damour formalism. After having worked out the basic principles of the hydrodynamics in Einstein-Cartan framework, we study the Weyssenhoff fluid within the post-Newtonian approximation scheme. The complexity of the underlying dynamical equations suggests to employ a discrete description via the point-particle limit, a procedure which permits the analysis of inspiralling spinning compact binaries. We then provide a first application of our results by considering binary neutron star systems.

In this paper we employ a recent proposal of C. Tsallis and formulate the first law of thermodynamics for gravitating systems in terms of the extensive but non-additive entropy. We pay a particular attention to an integrating factor for the heat one-form and show that in contrast to conventional thermodynamics it factorizes into thermal and entropic part. Ensuing first law of thermodynamics implies Tsallis cosmology, which is then subsequently used to address the observed discrepancy between current bound on the Dark Matter relic abundance and present IceCube data on high-energy neutrinos. To resolve this contradiction we keep the conventional minimal Yukawa-type interaction between standard model and Dark Matter particles but replace the usual Friedmann field equations with Tsallis-cosmology-based modified Friedmann equations. We show that when the Tsallis scaling exponent $\delta \sim 1.57$ (or equivalently, the holographic scaling exponent $\alpha \sim 3.13$) the aforementioned discrepancy disappears.

Black holes can accumulate a large amount of energy, responsible for highly energetic astrophysical phenomena Recently, fast magnetic reconnection (MR) of the magnetic field was proposed as a new way to extract energy and in this paper, we investigate this phenomena in a bumblebee Kerr-Sen BH. We find that the presence of the charge parameter strongly changes the simple Kerr case, making this extraction mechanism possible even for not extremely rotating black holes ($a \sim 0.7$). We also show that, under appropriate circumstances, MR is more efficient compared to the Blandford-Znajek mechanism. We finally compare these results with quintessence black-hole solutions not finding and enhancement respect to Kerr solution.

We report measurements of the angular distributions of low momentum atmospheric muons at 38 m above sea level for zenith angles $\theta$ between -60 and 60 degrees in the south-north direction. The muon detection was performed with two NaI(Tl) scintillation detectors mounted in coincidence. An adjustable lead thickness placed between the detectors allowed to select muons with a minimal momentum ranging from 0.3 to 0.9 GeV$/$c. The integrated and the differential muon flux were determined by analyzing the deposited energy spectra in the scintillators backed up by a Geant4 simulation of the experimental setup. The results are consistent with the $\cos^{n}(\theta)$ distribution in good agreement with the literature. These data contribute to fill the gap in the geomagnetic cutoff rigidity interval 8~GV$<P_c<$14~GV where no similar measurements were performed before. We found that $n=1.88-0.12~P_{\mu}^c$ in this domain of muon momenta cutoff $P_{\mu}^c$<1 GeV$/$c. The present measurements are useful for many muon studies requiring an accurate integrated flux.

We investigate the inflation driven by a non-linear electromagnetic field based on a NLED lagrangian density ${\cal L}_{nled} = - {\cal F} f \left( {\cal F} \right)$, where $f \left( {\cal F}\right)$ is a generalized functional depending on ${\cal F}$. We first formulate an $f$-NLED cosmological model with a more general functional $f \left( {\cal F}\right)$ and show that all NLED models can be expressed in this framework; then, we investigate in details two interesting examples of the functional $f \left( {\cal F}\right)$. We present our phenomenological model based on a new Lagrangian for NLED. Solutions to the field equations with the physical properties of the cosmological parameters are obtained. We show that the early Universe had no Big-Bang singularity, which accelerated in the past. We also investigate the qualitative implications of NLED by studying the inflationary parameters, like the slow-roll parameters, spectral index $n_s$, and tensor-to-scalar ratio $r$ and compare our results with observational data. Detailed phase-space analysis of our NLED cosmological model is performed with and without matter source. As a first approach, we consider the motion of a particle of unit mass in an effective potential. Our systems correspond to fast-slow systems for physical values of the electromagnetic field and the energy densities at the end of inflation. We analyze a complementary system using Hubble-normalized variables to investigate the cosmological evolution previous to the matter-dominated Universe.

We study self-consistent approximations such as the $\Phi$-derivable and virial approaches to dilute strongly interacting systems in equilibrium. We consider a system of non-relativistic fermions of one kind interacting via a pair potential Thermodynamical quantities are expressed in terms of various spectral functions. We review the $\Phi$ derivable approximation scheme demonstrating the exact conservation of the Noether and the Botermans-Malfliet fermion number densities, and then for $\Phi$ described by the tadpole and sandwich diagrams we show the coincidence of these two number densities. As examples of test pair potentials, we consider the Yukawa central nucleon-nucleon potentials within Walecka, CD Bonn, and Reid parameterizations, and the corresponding classical Lennard-Jones potentials. Expressions for the second and third virial coefficients are derived and analyzed for $\Phi$ described by the tadpole and sandwich diagrams. Next, we focus on the virial approach to the equation of state. Classical, semiclassical, and purely quantum approaches are studied in detail. Then, different extrapolations of the virial equation of state are considered including the van~der~Waals form and excluded-volume models. We derive the expression for the second virial coefficient using the effective range approximation for the scattering amplitude and compare the result with the purely quantum result using the experimental phase shifts. Attention is focused on the problem of anomalously large value of the nucleon-nucleon scattering length appearing due to the presence of the quasi-bound state in nucleon-nucleon scattering, which can be destroyed in the matter because of the action of the Pauli blocking. We present results for the second virial coefficient subtracting this term. We discuss the validity of such a procedure to describe the equation of state of the nuclear matter in the virial limit.

The kilo-Hertz gravitational waves radiated by the neutron star merger remnants carry rich information about the physics of high-density nuclear matter states, and many important astrophysical phenomena such as gamma-ray bursts and black hole formation. Current laser interferometer gravitational wave detectors, such as LIGO, VIRGO, and KAGRA have limited signal response at the kilo-Hertz band, thereby unable to capture these important physical phenomena. This work proposes an alternative protocol for boosting the sensitivity of the gravitational wave detectors at high frequency by implementing an optomechanical quantum amplifier. With the auxiliary quantum amplifier, this design has the feature of Parity-Time (PT) symmetry so that the detection band will be significantly broadened within the kilo-Hertz range. In this work, we carefully analyze the quantum-noise-limited sensitivity and the dynamical stability of this design. Based on our protocol, our result shows that the quantum-noise-limited sensitivity will be improved by one order of magnitude around 3kHz, which indicates the potential of our design for a future search of neutron star merger signals.

In this chapter we present basic concepts of neutrino physics. We start with a brief introduction to the standard model electroweak sector, followed by calculations of some relevant neutrino interaction cross sections. We obtain the oscillation formalism from the Dirac equation, alongside with a self-contained derivation of matter effects from the standard model electroweak lagrangian. Some key features of neutrino oscillation phenomenology are discussed as well. We then review the most precise oscillation measurements and finalize with some broad discussion of the current open questions in neutrino physics.

It is well-known that gravitino propagation in standard supergravities is free of any causality problems. However, two issues related to gravitino propagation were recently uncovered in specific supergravities with nonlinear supersymmetry. One of them concerns potential acausality/superluminality, whereas the second one arises from the vanishing of the sound speed at specific points during inflation. The former is famously related to positivity constraints on specific EFT operators, derived from dispersion relations on the energy-growing part of scattering amplitudes, and indeed we show that subluminality constraints for the gravitino are related via the equivalence theorem to positivity bounds in low-energy goldstino actions. However, the former are stronger, in the sense that they apply to functions of the scalar fields not only in the ground state, but for any field values such as those scanned by time-dependent solutions, unlike bounds derived from $2\to 2$ scattering amplitudes in the vacuum. We also argue that nontrivial causality constraints arise only in the case where nonlinear supersymmetry in the matter sector is encoded into superfield constraints which do not seem to arise from microscopic two-derivative lagrangians, in particular for the orthogonal constraint used to build minimal models of inflation in supergravity. This allows us to propose simple alternatives which maintain the minimality of the spectra and are causal in all points of the theory parameter space. We also discuss minimal supergravity models of inflation along these lines.

Within the extension of the $\Lambda$CDM model, allowing for the presence of neutrinos or warm dark matter, we develop the analytical cosmological perturbation theory. It covers all spatial scales where the weak gravitational field regime represents a valid approximation. Discrete particles - the sources of the inhomogeneous gravitational field - may be relativistic. Similarly to the previously investigated case of nonrelativistic matter, the Yukawa interaction range is naturally incorporated into the first-order scalar metric corrections.