Papers for Thursday, Jul 01 2021

Context. The [CII] 158micron far-infrared fine-structure line is one of the dominant cooling lines of the star-forming interstellar medium (ISM). Hence [CII] emission originates in and thus can be used to trace a range of ISM processes. Velocity-resolved large-scale mapping of [CII] in star-forming regions provides a unique perspective of the kinematics of these regions and their interactions with the exciting source of radiation. Aims. We explore the scientific applications of large-scale mapping of velocity-resolved [CII] observations. With the [CII] observations, we investigate the effect of stellar feedback on the ISM. We present the details of observation, calibration, and data reduction using a heterodyne array receiver mounted on an airborne observatory. Results. A square-degree [CII] map with a spectral resolution of 0.3 km/s is presented. The scientific potential of this data is summarized with discussion of mechanical and radiative stellar feedback, filament tracing using [CII], [CII] opacity effects, [CII] and carbon recombination lines, and [CII] interaction with the large molecular cloud. The data quality and calibration is discussed in detail, and new techniques are presented to mitigate the effects of unavoidable instrument deficiencies (e.g. baseline stability) and thus to improve the data quality. A comparison with a smaller [CII] map taken with the Herschel/Heterodyne Instrument for the Far-Infrared (HIFI) spectrometer is presented.

We investigate the attenuation law in $z\sim 6$ quasars by combining cosmological zoom-in hydrodynamical simulations of quasar host galaxies, with multi-frequency radiative transfer calculations. We consider several dust models differing in terms of grain size distributions, dust mass and chemical composition, and compare the resulting synthetic Spectral Energy Distributions (SEDs) with data from bright, early quasars. We show that only dust models with grain size distributions in which small grains ($a < 0.1~\mu$m, corresponding to $\approx 60\%$ of the total dust mass) are selectively removed from the dusty medium provide a good fit to the data. Removal can occur if small grains are efficiently destroyed in quasar environments and/or early dust production preferentially results in large grains. Attenuation curves for these models are close to flat, and consistent with recent data; they correspond to an effective dust-to-metal ratio $f_d \simeq 0.38$, i.e. close to the Milky Way value.

A growing number of eclipsing binary systems of the "HW Vir" kind (i. e., composed by a subdwarf-B/O primary star and an M dwarf secondary) show variations in their orbital period, also called Eclipse Time Variations (ETVs). Their physical origin is not yet known with certainty: while some ETVs have been claimed to arise from dynamical perturbations due to the presence of circumbinary planetary companions, other authors suggest that the Applegate effect or other unknown stellar mechanisms could be responsible for them. In this work, we present twenty-eight unpublished high-precision light curves of one of the most controversial of these systems, the prototype HW Virginis. We homogeneously analysed the new eclipse timings together with historical data obtained between 1983 and 2012, demonstrating that the planetary models previously claimed do not fit the new photometric data, besides being dynamically unstable. In an effort to find a new model able to fit all the available data, we developed a new approach based on a global-search genetic algorithm and eventually found two new distinct families of solutions that fit the observed timings very well, yet dynamically unstable at the 10^5-year time scale. This serves as a cautionary tale on the existence of formal solutions that apparently explain ETVs but are not physically meaningful, and on the need of carefully testing their stability. On the other hand, our data confirm the presence of an ETV on HW Vir that known stellar mechanisms are unable to explain, pushing towards further observing and modelling efforts.

We present H-band interferometric observations of the red supergiant (RSG) AZ Cyg made with the Michigan Infra-Red Combiner (MIRC) at the six-telescope Center for High Angular Resolution Astronomy (CHARA) Array. The observations span 5 years (2011-2016), offering insight into the short and long-term evolution of surface features on RSGs. Using a spectrum of AZ Cyg obtained with SpeX on the NASA InfraRed Telescope Facility (IRTF) and synthetic spectra calculated from spherical MARCS, spherical PHOENIX, and SAtlas model atmospheres, we derive $T_{\text{eff}}$ is between $3972 K$ and $4000 K$ and $\log~g$ between $-0.50$ and $0.00$, depending on the stellar model used. Using fits to the squared visibility and Gaia parallaxes we measure its average radius $R=911^{+57}_{-50}~R_{\odot}$. Reconstructions of the stellar surface using our model-independent imaging codes SQUEEZE and OITOOLS.jl show a complex surface with small bright features that appear to vary on a timescale of less than one year and larger features that persist for more than one year. 1D power spectra of these images suggest a characteristic size of $0.52-0.69~R_{\star}$ for the larger, long lived features. This is close to the values of $0.51-0.53~R_{\star}$ derived from 3D RHD models of stellar surfaces. We conclude that interferometric imaging of this star is in line with predictions of 3D RHD models but that short-term imaging is needed to more stringently test predictions of convection in RSGs.

Stellar RV jitter due to surface activity may bias the RV semi-amplitude and mass of rocky planets. The amplitude of the jitter may be estimated from the uncertainty in the rotation period, allowing the mass to be more accurately obtained. We find candidate rotation periods for 17 out of 35 TESS Objects of Interest (TOI) hosting <3 R_Earth planets as part of the Magellan-TESS Survey, which is the first-ever statistically robust study of exoplanet masses and radii across the photo-evaporation gap. Seven periods are 3+ sigma detections, two are 1.5+ sigma, and 8 show plausible variability but the periods remain unconfirmed. The other 18 TOIs are non-detections. Candidate rotators include the host stars of the confirmed planets L 168-9 b, the HD 21749 system, LTT 1445 A b, TOI 1062 b, and the L 98-59 system. 13 candidates have no counterpart in the 1000 TOI rotation catalog of Canto Martins et al. (2020). We find periods for G3-M3 dwarfs using combined light curves from TESS and the Evryscope all-sky array of small telescopes, sometimes with longer periods than would be possible with TESS alone. Secure periods range from 1.4 to 26 d with Evryscope-measured photometric amplitudes as small as 2.1 mmag in g'. We also apply Monte Carlo sampling and a Gaussian Process stellar activity model from the code exoplanet to the TESS light curves of 6 TOIs to confirm the Evryscope periods.

When exoplanets pass in front of their stars, they imprint a transit signature on the stellar light curve which to date has been assumed to be symmetric in time, owing to the planet being modelled as a circular area occulting the stellar surface. However, this signature might be asymmetric due to different temperature/pressure and/or chemical compositions in the different terminator regions of the transiting planet. catwoman is a Python package that allows to model these asymmetric transit lightcurves, calculating lightcurves for any radially symmetric stellar limb darkening law, and where planets are modelled as two semi-circles, of different radii, using the integration algorithm developed in arXiv:1507.08285 and implemented in the batman library, from which catwoman builds upon.

We present a search for transient radio sources on timescales of 2-9 years at 150 MHz. This search is conducted by comparing the first Alternative Data Release of the TIFR GMRT Sky Survey (TGSS ADR1) and the second data release of the LOFAR Two-metre Sky Survey (LoTSS DR2). The overlapping survey area covers 5570 $\rm{deg}^2$ on the sky, or 14% of the entire hemisphere. We introduce a method to compare the source catalogues that involves a pair match of sources, a flux density cutoff to meet the survey completeness limit and a newly developed compactness criterion. This method is used to identify both transient candidates in the TGSS source catalogue that have no counterpart in the LoTSS catalogue and transient candidates in LoTSS without a counterpart in TGSS. We find that imaging artefacts and uncertainties and variations in the flux density scales complicate the transient search. Our method to search for transients by comparing two different surveys, while taking into account imaging artefacts around bright sources and misaligned flux scales between surveys, is universally applicable to future radio transient searches. No transient sources were identified, but we are able to place an upper limit on the transient surface density of $<5.4 \cdot 10^{-4}\ \text{deg}^{-2}$ at 150 MHz for compact sources with an integrated flux density over 100 mJy. Here we define a transient as a compact source with flux greater than 100 mJy that appears in the catalogue of one survey without a counterpart in the other survey.

Measurement of the distances to nearby galaxies have improved rapidly in recent decades. The ever-present challenge is to reduce systematic effects, especially as greater distances are probed, and the uncertainties become larger. In this paper, we combine several recent calibrations of the Tip of the Red Giant Branch (TRGB) method. These calibrations are internally self-consistent at the 1% level. New Gaia Early Data Release 3 (EDR3) data provide an additional consistency check, at a (lower) 5% level of accuracy, a result of the well-documented Gaia angular covariance bias. The updated TRGB calibration applied to a distant sample of Type Ia supernovae from the Carnegie Supernova Project results in a value of the Hubble constant of Ho = 69.8 $\pm$ 0.6 (stat) $\pm$ 1.6 (sys) km/s/Mpc. No statistically significant difference is found between the value of Ho based on the TRGB and that determined from measurements of the cosmic microwave background. The TRGB results are also consistent to within 2$\sigma$ with the SHoES and Spitzer plus HST Key Project Cepheid calibrations. The TRGB results alone do not demand additional new physics beyond the standard Lambda-CDM cosmological model. They have the advantage of simplicity of the underlying physics (the core He flash) and small systematic uncertainties (from extinction, metallicity and crowding). Finally, the strengths and weaknesses of both the TRGB and Cepheids are reviewed, and prospects for addressing the current discrepancy with future Gaia, HST and JWST observations are discussed. Resolving this discrepancy is essential for ascertaining if the claimed tension in Ho between the locally-measured and the CMB-inferred value is physically motivated.

Asteroseismic studies of red giants generally assume that the oscillation modes can be treated as linear perturbations to the background star. However, observations by the Kepler mission show that the oscillation amplitudes increase dramatically as stars ascend the red giant branch. The importance of nonlinear effects should therefore be assessed. In previous work, we found that mixed modes in red giants are unstable to nonlinear three-wave interactions over a broad range of stellar mass and evolutionary state. Here we solve the amplitude equations that describe the mode dynamics for large networks of nonlinearly coupled modes. The networks consist of stochastically driven parent modes coupled to resonant secondary modes (daughters, granddaughters, etc.). We find that nonlinear interactions can lower the energy of gravity-dominated mixed modes by $\gtrsim 80\%$ compared to linear theory. However, they have only a mild influence on the energy of pressure-dominated mixed modes. Expressed in terms of the dipole mode visibility $V^2$, i.e., the summed amplitudes of dipole modes relative to radial modes, we find that $V^2$ can be suppressed by $50-80\%$ relative to the linear value for highly-evolved red giants whose frequency of maximum power $\nu_{\rm max} \lesssim 100\,\mu\textrm{Hz}$. However, for less evolved red giants with $150\lesssim \nu_{\rm max} \lesssim 200\,\mu\textrm{Hz}$, $V^2$ is suppressed by only $10-20\%$. We conclude that resonant mode coupling can have a potentially detectable effect on oscillations at $\nu_{\rm max} \lesssim 100\,\mu\textrm{Hz}$ but it cannot account for the population of red giants that exhibit dipole modes with unusually small amplitudes at high $\nu_{\rm max}$.

The technique of transmission spectroscopy -- the variation of the planetary radius with wavelength due to opacity sources in the planet's terminator region -- has been to date one of the most successful in the characterization of exoplanet atmospheres, providing key insights into the composition and structure of these distant worlds. A common assumption made when using this technique, however, is that the variations are the same in the entire terminator region. In reality, the morning and evening terminators might have distinct temperature, pressure and thus compositional profiles due to the inherent 3-D nature of the planet which would, in turn, give rise to different spectra on each side of it. Constraining those might be fundamental for our understanding of not only the weather patterns in these distant worlds, but also of planetary formation signatures which might only be possible to extract once these features are well understood. Motivated by this physical picture, in this work we perform a detailed study on the observational prospects of detecting this effect. We present an open-source semi-analytical framework with which this information can be extracted directly from transit lightcurves, and perform a detailed study on the prospects of detecting the effect with current missions such as TESS and upcoming ones such as JWST. Our results show that these missions show great promise for the detection of this effect. Transmission spectroscopy studies with JWST, in particular, could provide spectra of each of the limbs allowing us to convey 3-D information previously accessible only via phase-curves.

A novel equatorial platform is described. While the invention lacks the compactness some desire in an equatorial platform, it holds promise of high precision, large payload weight and inexpensive construction.

We present the results of the asteroseismic analysis of the hydrogen-deficient white dwarf PG 0112+104 from the $Kepler$-2 field. Our seismic procedure using the forward method based on physically sound, static models, includes the new core parameterization leading us to reproduce the periods of this star near the precision of the observations. This new fit outperforms current state-of-the-art standards by order of magnitudes. We precisely establish the internal structure and unravel the inner C/O stratification of its core. This opens up interesting perspectives on better constraining key processes in stellar physics such as nuclear burning, convection, and mixing, that shape this stratification over time.

We have performed a series of numerical experiments aimed at studying the activation of Kerr black holes (BHs) by advection of small scale magnetic fields. Such configurations may potentially give rise to the formation of quasi-striped Blandford-Znajek jets. It can also lead to enhanced dissipation and generation of plasmoids in current sheets formed in the vicinity of the BH horizon, which may constitute a mechanism to power the hard X-ray emission seen in many accreting BH systems (a la lamppost models). Our analysis suggests that formation of quasi-striped jets with significant power may be possible provided loops with alternating polarity having sizes larger than $\sim 10 r_g$ or so can be maintained (either form sporadically or advected from outside) at a radius $\lesssim 10^2 r_g$. This conclusion is consistent with recent results of general relativistic force-free simulations. We also find that the accretion dynamics exhibits cyclic behaviour in MAD states, alternating between high accretion phases and quenched accretion phases during which the magnetosphere becomes force-free out to radii $\gtrsim 10r_g$. We suggest that such a behaviour should lead to notable variations of the observed luminosity and image of the inner disc (BH shadow image). Finally, we find that the transition between accreted loops on the BH gives rise to the formation of current sheets and energetic plasmoids on the jet boundary during intermittent periods when the jet becomes inactive, in addition to an equatorial current sheet that forms during peaks in the jet activity.

Gravitational wave observations of binary neutron star mergers provide valuable information about neutron star structure and the equation of state of dense nuclear matter. Numerous methods have been proposed to analyze the population of observed neutron stars and previous work has demonstrated the necessity of jointly fitting the astrophysical distribution and the equation of state in order to accurately constrain the equation of state. In this work, we introduce a new framework to simultaneously infer the distribution of binary neutron star masses and the nuclear equation of state using Gaussian mixture model density estimates which mitigates some of the limitations previously-used methods suffer from. Using our method, we reproduce previous projections for the expected precision of our joint mass distribution and equation of state inference with tens of observations. We also show that mismodeling the equation of state can bias our inference of the neutron star mass distribution. While we focus on neutron star masses and matter effects, our method is widely applicable to population inference problems.

Magnetosheath jets are currently an important topic in the field of magnetosheath physics. It is thought that 97~\% of the jets are produced by the shock rippling at quasi-parallel shocks. Recently, large statistical studies of magnetosheath jets have been performed, however it is not clear whether rippling also produces jets found downstream of quasi-perpendicular shocks. We analyze four types of events in the quasi-perpendicular magnetosheath with signatures characteristic of magnetosheath jets, namely increased density and/or dynamic pressure, that were not produced by the shock rippling: 1) magnetic flux tubes connected to the quasi-parallel bow-shock, 2) non-reconnecting current sheets, 3) reconnection exhausts and 4) mirror mode waves. The flux tubes are downstream equivalents of the upstream traveling foreshocks. Magnetosheath jets can impact the magnetopause, so knowing the conditions under which they form may enable us to understand their signatures in the magnetosphere.

SS Cyg has long been recognized as the prototype of a group of dwarf novae that show only outbursts. However, this object has entered a quite anomalous event in 2021, which at first appeared to be standstill, i.e., an almost constant luminosity state, observed in Z Cam-type dwarf novae. This unexpected event gives us a great opportunity to reconsider the nature of standstill in cataclysmic variables. We have observed this anomalous event and its forerunner, a gradual and simultaneous increase in the optical and X-ray flux during quiescence, through many optical telescopes and the X-ray telescopes NICER and NuSTAR. We have not found any amplification of the orbital hump during quiescence before the anomalous event, which suggests that the mass transfer rate did not significantly fluctuate on average. The estimated X-ray flux was not enough to explain the increment of the optical flux during quiescence via X-ray irradiation of the disk and the secondary star. It would be natural to consider that viscosity in the quiescent disk was enhanced before the anomalous event, which increased mass accretion rates in the disk and raised not only the optical flux but also the X-ray flux. We suggest that enhanced viscosity also triggered the standstill-like phenomenon in SS Cyg, which is considered to be a series of small outbursts. The inner part of the disk would always stay in the outburst state and only its outer part would be unstable against the thermal-viscous instability during this phenomenon, which is consistent with the observed optical color variations. This scenario is in line with our X-ray spectral analyses which imply that the X-ray emitting inner accretion flow became hotter than usual and vertically expanded and that it became denser and was cooled down after the onset of the standstill-like state.

We observe the density wave angular pattern speed OMEGA-p to be near 12 to 17 km / s / kpc, by the separation between a typical optical HII region (from the spiral arm dust lane) and using a HII evolution time model to yield its relative speed, and independently by the separation between a typical radio maser (from the spiral arm dust lane) with a maser model.

We analyse the tidal disruption probability of potential neutron star--black hole (NSBH) merger gravitational wave (GW) events, including GW190426_152155, GW190814, GW200105_162426 and GW200115_042309, detected during the third observing run of the LIGO/Virgo Collaboration, and the detectability of kilonova emission in connection with these events. The posterior distributions of GW190814 and GW200105_162426 show that they must be plunging events and hence no kilonova signal is expected from these events. With the stiffest NS equation of state allowed by the constraint of GW170817 taken into account, the probability that GW190426_152155 and GW200115_042309 can make tidal disruption is $\sim24\%$ and $\sim3\%$, respectively. However, the predicted kilonova brightness is too faint to be detected for present follow-up search campaigns, which explains the lack of electromagnetic (EM) counterpart detection after triggers of these GW events. Based on the best constrained population synthesis simulation results, we find that disrupted events account for only $\lesssim20\%$ of cosmological NSBH mergers since most of the primary BHs could have low spins. The associated kilonovae for those disrupted events are still difficult to be discovered by LSST after GW triggers in the future, because of their low brightness and larger distances. For future GW-triggered multi-messenger observations, potential short-duration gamma-ray bursts and afterglows are more probable EM counterparts of NSBH GW events.

The discovery that many classical novae produce detectable GeV $\gamma$-ray emission has raised the question of the role of shocks in nova eruptions. Here we use radio observations of nova V809 Cep (Nova Cep 2013) with the Jansky Very Large Array to show that it produced non-thermal emission indicative of particle acceleration in strong shocks for more than a month starting about six weeks into the eruption, quasi-simultaneous with the production of dust. Broadly speaking, the radio emission at late times -- more than a six months or so into the eruption -- is consistent with thermal emission from $10^{-4} M_\odot$ of freely expanding, $10^4$~K ejecta. At 4.6 and 7.4 GHz, however, the radio light-curves display an initial early-time peak 76 days after the discovery of the eruption in the optical ($t_0$). The brightness temperature at 4.6 GHz on day 76 was greater than $10^5 K$, an order of magnitude above what is expected for thermal emission. We argue that the brightness temperature is the result of synchrotron emission due to internal shocks within the ejecta. The evolution of the radio spectrum was consistent with synchrotron emission that peaked at high frequencies before low frequencies, suggesting that the synchrotron from the shock was initially subject to free-free absorption by optically thick ionized material in front of the shock. Dust formation began around day 37, and we suggest that internal shocks in the ejecta were established prior to dust formation and caused the nucleation of dust.

We present SoFiA 2, the fully automated 3D source finding pipeline for the WALLABY extragalactic HI survey with the Australian SKA Pathfinder (ASKAP). SoFiA 2 is a reimplementation of parts of the original SoFiA pipeline in the C programming language and makes use of OpenMP for multi-threading of the most time-critical algorithms. In addition, we have developed a parallel framework called SoFiA-X that allows the processing of large data cubes to be split across multiple computing nodes. As a result of these efforts, SoFiA 2 is substantially faster and comes with a much reduced memory footprint compared to its predecessor, thus allowing the large WALLABY data volumes of hundreds of gigabytes of imaging data per epoch to be processed in real-time. The source code has been made publicly available to the entire community under an open-source licence. Performance tests using mock galaxies injected into genuine ASKAP data suggest that in the absence of significant imaging artefacts SoFiA 2 is capable of achieving near-100% completeness and reliability above an integrated signal-to-noise ratio of about 5-6. We also demonstrate that SoFiA 2 generally recovers the location, integrated flux and w20 line width of galaxies with high accuracy. Other parameters, including the peak flux density and w50 line width, are more strongly biased due to the influence of the noise on the measurement. In addition, very faint galaxies below an integrated signal-to-noise ratio of about 10 may get broken up into multiple components, thus requiring a strategy to identify fragmented sources and ensure that they do not affect the integrity of any scientific analysis based on the SoFiA 2 output.

A tidal disruption event (TDE) involves the tidal shredding of a star in the vicinity of a dormant supermassive black hole. The nearby ($\approx$230 mega-parsec) radio-quiet (radio luminosity of $4 \times 10^{38}$ erg s$^{-1}$) AT2019dsg is the first TDE potentially associated with a neutrino event. The origin of the non-thermal emission in AT2019dsg remains inconclusive; possibilities include a relativistic jet or a sub-relativistic outflow. Distinguishing between them can address neutrino production mechanisms. High resolution very long baseline interferometry monitoring provides uniquely constraining flux densities and proper motion of the ejecta. A non-relativistic (outflow velocity of $\approx$0.1 $c$) decelerated expansion in a relatively dense environment is found to produce the radio emission. Neutrino production may be related to the acceleration of protons by the outflow. The present study thus helps exclude jet-related origins for the non-thermal emission and neutrino production, and constrains non-jetted scenarios.

Outbursts of young stellar objects occur when the mass accretion rate suddenly increases. However, such outbursts are difficult to detect for deeply embedded protostars due to their thick envelope and the rarity of outbursts. The near-IR spectroscopy is a useful tool to identify ongoing outburst candidates by the characteristic absorption features that indicate a disk origin. However, without high-resolution spectroscopy, the spectra of outburst candidates can be confused with the late-type stars since they have similar spectral features. For the protostar IRAS 16316-1540, the near-IR spectrum has line equivalent widths that are consistent with M-dwarf photospheres. However, our high-resolution IGRINS spectra reveal that the absorption lines have boxy and/or double-peaked profiles, as expected from a disk and not the star. The continuum emission source is likely the hot, optically thick disk, heated by viscous accretion. The projected disk rotation velocity of 41$\pm$5 km s$^{-1}$ corresponds to $\sim 0.1$ AU. Based on the result, we suggest IRAS 16316-1540 as an ongoing outburst candidate. Viscous heating of disks is usually interpreted as evidence for ongoing bursts, which may be more common than previously estimated from low-resolution near-IR spectra.

In this letter, we implement a test of the standard law for the dark matter density evolution. For this purpose, only a flat universe and the validity of the FRW metric are assumed. A deformed dark matter density evolution law is considered, given by $\rho_c(z) \propto (1+z)^{3+\epsilon}$, and constraints on $\epsilon$ are obtained by using galaxy cluster gas mass fractions, and cosmic chronometers measurements. We find that $\epsilon =0$ within 2$\sigma$ c.l., in full agreement with other recent analyses.

Galactic electron density distribution models are crucial tools for estimating the impact of the ionised interstellar medium on the impulsive signals from radio pulsars and fast radio bursts. The two prevailing Galactic electron density models are YMW16 (Yao et al., 2017) and NE2001 (Cordes & Lazio, 2002). Here, we introduce a software package PyGEDM which provides a unified application programming interface (API) for these models and the YT20 (Yamasaki & Totani, 2020) model of the Galactic halo. We use PyGEDM to compute all-sky maps of Galactic dispersion measure (DM) for YMW16 and NE2001, and compare the large-scale differences between the two. In general, YMW16 predicts higher DM values toward the Galactic anticentre. YMW16 predicts higher DMs at low Galactic latitudes, but NE2001 predicts higher DMs in most other directions. We identify lines of sight for which the models are most discrepant, using pulsars with independent distance measurements. YMW16 performs better on average than NE2001, but both models show significant outliers. We suggest that future campaigns to determine pulsar distances should focus on targets where the models show large discrepancies, so future models can use those measurements to better estimate distances along those line of sight. We also suggest that the Galactic halo should be considered as a component in future GEDMs, to avoid overestimating the Galactic DM contribution for extragalactic sources such as FRBs.

It is generally accepted that CMEs result from eruptions of magnetic flux ropes, which are dubbed as magnetic clouds in interplanetary space. The composition (including the ionic charge states and elemental abundances) is determined prior to and/or during CME eruptions in the solar atmosphere, and does not alter during magnetic cloud propagation to 1 AU and beyond. It has been known that the composition is not uniform within a cross section perpendicular to magnetic cloud axis, and the distribution of ionic charge states within a cross section provides us an important clue to investigate the formation and eruption processes of flux ropes due to the freeze-in effect. The flux rope is a three dimensional magnetic structure intrinsically, and it remains unclear whether the composition is uniform along the flux rope axis as most magnetic clouds are only detected by one spacecraft. In this paper we report a magnetic cloud that was observed by ACE near 1 AU on 1998 March 4--6 and Ulysses near 5.4 AU on March 24--28 sequentially. At these times, both spacecraft were located around the ecliptic plane, and the latitudinal and longitudinal separations between them were $\sim$2.2$^{\circ}$ and $\sim$5.5$^{\circ}$, respectively. It provides us an excellent opportunity to explore the axial inhomogeneity of flux rope composition, as both spacecraft almost intersected the cloud center at different sites along its axis. Our study shows that the average values of ionic charge states exhibit significant difference along the axis for carbon, and the differences are relatively slight but still obvious for charge states of oxygen and iron, as well as the elemental abundances of iron and helium. Besides the means, the composition profiles within the cloud measured by both spacecraft also exhibit some discrepancies. We conclude that the inhomogeneity of composition exists along the cloud axis.

Long gamma-ray bursts are associated with the core-collapse of massive, rapidly spinning stars. However, the believed efficient angular momentum transport in stellar interiors leads to predominantly slowly-spinning stellar cores. Here, we report on binary stellar evolution and population synthesis calculations, showing that tidal interactions in close binaries not only can explain the observed sub-population of spinning, merging binary black holes, but also lead to long gamma-ray bursts at the time of black-hole formation, with rates matching the empirical ones. We find that $\approx$10% of the GWTC-2 reported binary black holes had a long gamma-ray burst associated with their formation, with GW190517 and GW190719 having a probability of $\approx$85% and $\approx$60%, respectively, being among them.

Optical, near-infrared (NIR) photometric and spectroscopic studies along with the optical imaging polarimetric results for SN 2012au to constrain the nature of the progenitor and other properties are presented in this paper. Well-calibrated multi-band optical photometric data (from $-$0.2 to +413 d since $B$-band maximum) was used to compute the bolometric light curve and to perform semi$-$analytical light-curve modelling using the $\texttt{MINIM}$ code. A spin-down millisecond magnetar-powered model reasonably explains the observed photometric evolution of SN 2012au. Early time imaging polarimetric follow-up observations ($-$2 to +31 d) and their comparison with other similar cases indicate signatures of asphericity in the ejecta. Good spectral coverage of SN 2012au (from $-$5 to +391 d) allows us to trace the evolution of layers of SN ejecta in detail. SN 2012au exhibits higher line velocities in comparison to the other SNe~Ib. Late nebular phase spectra of SN 2012au indicate for a Wolf-Rayet star as the possible progenitor for SN 2012au, with oxygen, He-core, and main sequence masses of $\sim$1.62 $\pm$ 0.15 M$_\odot$, $\sim$4$-$8 M$_\odot$, and $\sim$17$-$25 M$_\odot$, respectively. There is a clear absence of first-overtone of the carbon monoxide (CO) features up to +319 d in the $K$-band region of the NIR spectra. Overall analysis suggests that SN 2012au is one of the most luminous slow-decaying Type Ib SN having comparatively higher ejecta mass ($\sim$4.7$-$8.3 M$_\odot$) and kinetic energy ($\sim$[4.8 $-$ 5.4] $\times$ 10$^{51}$ erg). Detailed modelling using the $\texttt{MESA}$ and the results obtained through $\texttt{STELLA}$ and $\texttt{SNEC}$ explosions also strongly support spin-down of a magnetar as a possible powering source of SN 2012au having mass around 20 M$_\odot$ and metallicity of Z = 0.04.

Expanding on the comment by Lyne et al (2017), that intermittent pulsars tend to congregate near a stripe in the logarithmic period versus period-derivative diagram, representing a small range of polar cap electric potential, as well as the fact (already apparent in their Fig.~7, but not explicitly stated there) that high-fraction nulling pulsars also tend to reside within this and an additional stripe, we make the observation that the two stripes further match the "death lines" for double and single-pole interpulses, associated with nearly orthogonal and aligned rotators respectively. These extreme inclinations are known to suffer from pair production deficiencies, so we propose to explain intermittency and high-fraction nulling by reinvigorating some older quiescent (no pulsar wind or radio emission) "electrosphere" solutions. Specifically, as the polar potential drops below the two threshold bands (i.e., the two stripes), corresponding to the aligned and orthogonal rotators, their respective magnetospheres transition from being of the active pair-production-sustained type into becoming the electrospheres, in which charges are only lifted from the star. The borderline cases sitting in the gap outside of the stable regime of either case manifest as high-fraction nullers. Hall evolution of the magnetic field inside orthogonally rotating neutron stars can furthermore drive secular regime changes, resulting in intermittent pulsars.

Dark energy has been introduced to explain the present accelerating expansion of the universe. In the LambdaCDM model, the present standard model of cosmology, dark energy is described as a cosmological constant which is time independent. The possibility of the time dependence of dark energy has been investigated to obtain deeper understanding of it by looking redshit-magnitude relation of type Ia supernova, for example. We investigate the constraints to the time dependence of dark energy with several phenomenological models by fitting their parameters to the redshit-magnitude relation of the Pantheon supernova catalog, and confirm that little time dependence can not be excluded by the goodness-of-fit criterion only. Such a non-trivial time dependence causes different expansion rate in the period of reionization, which affects the low-l polarization power spectrum of CMB. We investigate the possibility to detect the effect in future probe like LiteBIRD, for example. We find that the low-l EE polarization power spectrum is enhanced in general, but it will be difficult to be detected as far as looking the EE polarization power spectrum only because of the limitation by cosmic variance.

Direct measurements of the masses of supermassive black holes (SMBHs) are key to understanding their growth and constrain their symbiotic relationship to their host galaxies. However, current methods used to directly measure black hole masses in active quasars become challenging or impossible beyond $z\gtrsim0.2$. Spectroastrometry (SA) measures the spatial centroid of an object's spectrum as a function of wavelength, delivering angular resolution far better than PSF for high signal-to-noise ratio observations. We observed the luminous quasar SDSS J212329.47--005052.9 at $z=2.279$ with the aim of resolving its $\sim100\mu\mathrm{as}$ H$\alpha$ broad emission line region (BLR), and present the first SA constraints on the size and kinematic structure of the BLR. Using a novel pipeline to extract the SA signal and reliable uncertainties, we achieved a centroiding precision of $\simeq100\mu\mathrm{as}$, or $>2000\times$ smaller than the $K$-band AO-corrected PSF, yielding a tentative $3.2\sigma$ detection of an SA signal from the BLR. Modeling the BLR emission as arising from an inclined rotating disk with a mixture of coherent and random motions we constrain $r_\mathrm{BLR}=454^{+565}_{-162}\,\mu\mathrm{as}$ ($3.71^{+4.65}_{-1.28}\,\mathrm{pc}$), providing a $95\%$ confidence upper limit on the black hole mass $M_\mathrm{BH}\,\sin^2\,i \leq1.8 \times10^9\,\mathrm{M}_\odot$. Our results agree with the $r_\mathrm{BLR}-L$ relation measured for lower-$z$ quasars, but expands its dynamic range by an order of magnitude in luminosity. We did not detect the potentially stronger SA signal from the narrow line region, but discuss in detail why it may be absent. Already with existing instrumentation, SA can deliver $\sim6\times$ smaller uncertainties ($\sim15\,\mu\mathrm{as}$) than achieved here, enabling $\sim10\%$ measurements of SMBH masses in high-$z$ quasars.

In 2007, the young star 1SWASP J140747.93-394542.6 (V1400 Cen) underwent a complex series of deep eclipses over 56 days. This was attributed to the transit of a ring system filling a large fraction of the Hill sphere of an unseen substellar companion. Subsequent photometric monitoring has not found any other deep transits from this candidate ring system, but if there are more substellar companions and they are coplanar with the potential ring system, there is a chance that they will transit the star as well. This young star is active and the light curves show a 5% modulation in amplitude with a dominant rotation period of 3.2 days due to star spots rotating in and out of view. We model and remove the rotational modulation of the J1407 light curve and search for additional transit signatures of substellar companions orbiting around J1407. We combine the photometry of J1407 from several observatories, spanning a 19 year baseline. We remove the rotational modulation by modeling the variability as a periodic signal, whose periodicity changes slowly with time over several years due to the activity cycle of the star. A Transit Least Squares (TLS) analysis searches for any periodic transiting signals within the cleaned light curve. We identify an activity cycle of J1407 with a period of 5.4 years. A Transit Least Squares search does not find any plausible periodic eclipses in the light curve, from 1.2% amplitude at 5 days up to 1.9% at 20 days. This sensitivity is confirmed by injecting artificial transits into the light curve and determining the recovery fraction as a function of transit depth and orbital period. J1407 is confirmed as a young active star with an activity cycle consistent with a rapidly rotating solar mass star. With the rotational modulation removed, the TLS analysis rules out transiting companions with radii larger than about 1 Jupiter.

The role of solar wind expansion in generating whistler waves is investigated using the EB-iPic3D code, which models solar wind expansion self-consistently within a fully kinetic semi-implicit approach. The simulation is initialized with an electron velocity distribution function modeled after Parker Solar Probe observations during its first perihelion at 0.166 au, consisting of a dense core and an anti-sunward strahl. This distribution function is initially stable with respect to kinetic instabilities. Expansion drives the solar wind into successive regimes where whistler heat flux instabilities are triggered. These instabilities produce sunward whistler waves initially characterized by predominantly oblique propagation with respect to the interplanetary magnetic field. The excited waves interact with the electrons via resonant scattering processes. As a consequence, the strahl pitch angle distribution broadens and its drift velocity reduces. Strahl electrons are scattered in the direction perpendicular to the magnetic field, and an electron halo is formed. At a later stage, resonant electron firehose instability is triggered and further affects the electron temperature anisotropy as the solar wind expands. Wave-particle interaction processes are accompanied by a substantial reduction of the solar wind heat flux. The simulated whistler waves are in qualitative agreement with observations in terms of wave frequencies, amplitudes and propagation angles. Our work proposes an explanation for the observations of oblique and parallel whistler waves in the solar wind. We conclude that solar wind expansion has to be factored in when trying to explain kinetic processes at different heliocentric distances.

Recent work has revealed two classes of Globular Clusters (GCs), dubbed Type-I and Type-II. Type-II GCs are characterized by a blue- and a red- red giant branch composed of stars with different metallicities, often coupled with distinct abundances in the slow-neutron capture elements (s-elements). Here we continue the chemical tagging of Type-II GCs by adding the two least-massive clusters of this class, NGC1261 and NGC6934. Based on both spectroscopy and photometry, we find that red stars in NGC1261 are slightly enhanced in [Fe/H] by ~0.1 dex and confirm that red stars of NGC 6934 are enhanced in iron by ~0.2 dex. Neither NGC1261 nor NGC6934 show internal variations in the s-elements, which suggests a GC mass threshold for the occurrence of s-process enrichment. We found a significant correlation between the additional Fe locked in the red stars of Type-II GCs and the present-day mass of the cluster. Nevertheless, most Type II GCs retained a small fraction of Fe produced by SNe II, lower than the 2%; NGC6273, M54 and omega Centauri are remarkable exceptions. In the appendix, we infer for the first time chemical abundances of Lanthanum, assumed as representative of the s-elements, in M54, the GC located in the nucleus of the Sagittarius dwarf galaxy. Red-sequence stars are marginally enhanced in [La/Fe] by 0.10\pm0.06 dex, in contrast with the large [La/Fe] spread of most Type II GCs. We suggest that different processes are responsible for the enrichment in iron and s-elements in Type-II GCs.

The Galactic B[e] supergiant MWC 137 is surrounded by a large-scale optical nebula. To shed light on the physical conditions and kinematics of the nebula, we analyze the optical forbidden emission lines [NII] 6548,6583 and [SII] 6716,6731 in long-slit spectra taken with ALFOSC at the Nordic Optical Telescope. The radial velocities display a complex behavior but, in general, the northern nebular features are predominantly approaching while the southern ones are mostly receding. The electron density shows strong variations across the nebula with values spreading from about zero to ~800 cm$^{-3}$. Higher densities are found closer to MWC~137 and in regions of intense emission, whereas in regions with high radial velocities the density decreases significantly. We also observe the entire nebula in the two [SII] lines with the scanning Fabry-Perot interferometer attached to the 6-m telescope of the Special Astrophysical Observatory. These data reveal a new bow-shaped feature at PA = 225-245 and a distance 80" from MWC 137. A new H$\alpha$ image has been taken with the Danish 1.54-m telescope on La Silla. No expansion or changes in the nebular morphology appear within 18.1 years. We derive a mass of 37 (+9/-5) solar masses and an age of $4.7\pm0.8$ Myr for MWC 137. Furthermore, we detect a period of 1.93 d in the time series photometry collected with the TESS satellite, which could suggest stellar pulsations. Other, low-frequency variability is seen as well. Whether these signals are caused by internal gravity waves in the early-type star or by variability in the wind and circumstellar matter currently cannot be distinguished.

The origin of magnetic fields in the Universe is an open problem. Seed magnetic fields possibly produced in early times may have survived up to the present day close to their original form, providing an untapped window to the primeval Universe. The recent observations of high-energy neutrinos from the blazar TXS 0506+056 in association with an electromagnetic counterpart in a broad range of wavelengths can be used to probe intrinsic properties of this object and the traversed medium. Here we show that intergalactic magnetic fields (IGMFs) can affect the intrinsic spectral properties of this object reconstructed from observations. In particular, we point out that the reconstructed maximum gamma-ray energy of TXS 0506+056 can be significantly higher if IGMFs are strong. Finally, we use this flare to constrain both the magnetic-field strength and the coherence length of IGMFs.

We present a deep learning model to predict the r-band bulge-to-total light ratio (B/T) of nearby galaxies using their multi-band JPEG images alone. Our Convolutional Neural Network (CNN) based regression model is trained on a large sample of galaxies with reliable decomposition into the bulge and disk components. The existing approaches to estimate the B/T use galaxy light-profile modelling to find the best fit. This method is computationally expensive, prohibitively so for large samples of galaxies, and requires a significant amount of human intervention. Machine learning models have the potential to overcome these shortcomings. In our CNN model, for a test set of 20000 galaxies, 85.7 per cent of the predicted B/T values have absolute error (AE) less than 0.1. We see further improvement to 87.5 per cent if, while testing, we only consider brighter galaxies (with r-band apparent magnitude < 17) with no bright neighbours. Our model estimates B/T for the 20000 test galaxies in less than a minute. This is a significant improvement in inference time from the conventional fitting pipelines, which manage around 2-3 estimates per minute. Thus, the proposed machine learning approach could potentially save a tremendous amount of time, effort and computational resources while predicting B/T reliably, particularly in the era of next-generation sky surveys such as the Legacy Survey of Space and Time (LSST) and the Euclid sky survey which will produce extremely large samples of galaxies.

We present here the discovery of a new class of super-slow rotating asteroids (P>1000 hours) in data extracted from the Asteroid Terrestrial-impact Last Alert System (ATLAS) and Zwicky Transient Facility (ZTF) all-sky surveys. Of the 39 rotation periods we report here, 32 have periods longer than any previously reported unambiguous rotation periods currently in the Asteroid Light Curve Database. In our sample, 7 objects have a rotation period > 4000 hours and the longest period we report here is 4812 hours (~200 days). We do not observe any correlation between taxonomy, albedo, or orbital properties with super-slow rotating status. The most plausible mechanism for the creation of these very slow rotators is if their rotations were slowed by YORP spin-down. Super-slow rotating asteroids may be common, with at least 0.4% of the main-belt asteroid population with a size range between 2 and 20 km in diameter rotating with periods longer than 1000 hours.

For the first time, a direct NLTE analysis of carbon and nitrogen lines in the spectra of nine RR Lyrae stars was carried out. We have determined the abundances of these elements together with oxygen, and have shown that the nitrogen content is increased in metallicity deficient program stars. We conclude that this is a sign of the first dredge up, which occurred at the previous stage of the red giant branch, and brought material processed in an incomplete CNO cycle to the surface of the star. This effect is significantly enhanced by thermohaline (extra-) mixing, which is more effective for metal-poor RR Lyrae stars. This is clearly seen in the plot showing that C/N ration in our sample of stars gradually decreases as metallicity decreases from about --0.2 to --2. Oxygen abundance depends on metallicity in a similar way to what we see in the Population II stars.

Recent observations have detected excess H$\alpha$ emission from young stellar systems with an age of several Myr such as PDS 70. One-dimensional radiation-hydrodynamic models of shock-heated flows that we developed previously demonstrate that planetary accretion flows of $>$ a few ten km s$^{-1}$ can produce H$\alpha$ emission. It is, however, a challenge to understand the accretion process of proto-giant planets from observations of such shock-originated emission because of a huge gap in scale between the circumplanetary disk (CPD) and the microscopic accretion shock. To overcome the scale gap problem, we combine two-dimensional, high-spatial-resolution global hydrodynamic simulations and the one-dimensional local radiation hydrodynamic model of the shock-heated flow. From such combined simulations for the protoplanet-CPD system, we find that the H$\alpha$ emission is mainly produced in localized areas on the protoplanetary surface. The accretion shocks above CPD produce much weaker H$\alpha$ emission (approximately 1-2 orders of magnitude smaller in luminosity). Nevertheless, the accretion shocks above CPD significantly affect the accretion process onto the protoplanet. The accretion occurs at a quasi-steady rate, if averaged on a 10-day timescale, but its rate shows variability on shorter timescales. The disk surface accretion layers including the CPD-shocks largely fluctuate, which results in the time-variable accretion rate and H$\alpha$ luminosity of the protoplanet. We also model the spectral emission profile of the H$\alpha$ line and find that the line profile is less time-variable, despite the large variability in luminosity. High-spectral resolution spectroscopic observation and monitoring will be key to reveal the property of the accretion process.

The age-metallicity relation is fundamental to study the formation and evolution of the disk. Observations have shown that this relation has a large scatter which can not be explained by observational errors only. That scatter is hence attributed to the effects of radial migration in which stars tracing different chemical evolution histories in the disk get mixed. However, the recent study of Nissen et al. 2020, using high precision observational data of solar type stars, found two relatively tight age-metallicity relations. One sequence of older and metal-richer stars probably traces the chemical enrichment history of the inner disk while the other sequence of younger and metal-poorer stars the chemical enrichment history of the outer disk. If uncertainties in age measurements increase, these sequences mix explaining the scatter of the one relation observed in other studies. This work follows up on these results, by analysing an independent sample of red clump giants observed by APOGEE. Since ages for red giants are significantly more uncertain, the [C/N] ratios are considered as a proxy for age. This larger dataset is used to investigate these relations at different Galactic radii, finding that these distinct sequences exist only in the solar neighbourhood. The APOGEE dataset is further used to explore different abundance and kinematical planes to shed light on the nature of these populations.

One of the major discoveries of Hinode's Extreme-ultraviolet Imaging Spectrometer (EIS) is the presence of upflows at the edges of active regions. As active regions are magnetically connected to the large-scale field of the corona, these upflows are a likely contributor to the global mass cycle in the corona. Here we examine the driving mechanism(s) of the very strong upflows with velocities in excess of 70 km/s, known as blue-wing asymmetries, observed during the eruption of a flux rope in AR 10977 (eruptive flare SOL2007-12-07T04:50). We use Hinode/EIS spectroscopic observations combined with magnetic-field modeling to investigate the possible link between the magnetic topology of the active region and the strong upflows. A Potential Field Source Surface (PFSS) extrapolation of the large-scale field shows a quadrupolar configuration with a separator lying above the flux rope. Field lines formed by induced reconnection along the separator before and during the flux-rope eruption are spatially linked to the strongest blue-wing asymmetries in the upflow regions. The flows are driven by the pressure gradient created when the dense and hot arcade loops of the active region reconnect with the extended and tenuous loops overlying it. In view of the fact that separator reconnection is a specific form of the more general quasi-separatrix (QSL) reconnection, we conclude that the mechanism driving the strongest upflows is, in fact, the same as the one driving the persistent upflows of approx. 10 - 20 km/s observed in all active regions.

A kilonova signal is generally expected after a Black Hole - Neutron Star merger. The strength of the signal is related to the equation of state of neutron star matter and it increases with the stiffness of the latter. The recent results obtained by NICER suggest a rather stiff equation of state and the expected kilonova signal is therefore strong, at least if the mass of the Black Hole does not exceed $\sim 10 M_\odot$. We compare the predictions obtained by considering equations of state of neutron star matter satisfying the most recent observations and assuming that only one family of compact stars exists with the results predicted in the two-families scenario. In the latter a soft hadronic equation of state produces very compact stellar objects while a rather stiff quark matter equation of state produces massive strange quark stars, satisfying NICER results. The expected kilonova signal in the two-families scenario is very weak: the Strange Quark Star - Black Hole merger does not produce a kilonova signal because, according to simulations, the amount of mass ejected is negligible and the Hadronic Star - Black Hole merger produces a much weaker signal than in the one-family scenario because the hadronic equation of state is very soft. This prediction will be easily tested with the new generation of detectors.

We present analysis of XMM-Newton Optical Monitor observations in the near-ultraviolet of HD 189733, covering twenty primary transits of its hot Jupiter planet. The transit is clearly detected with both the UVW2 and UVM2 filters, and our fits to the data reveal transit depths in agreement with that observed optically. The measured depths correspond to radii of $1.059^{+0.046}_{-0.050}$ and $0.94^{+0.15}_{-0.17}$ times the optically-measured radius (1.187 R$_{\rm J}$ at 4950 \r{A}) in the UVW2 and UVM2 bandpasses, respectively. We also find no statistically significant variation in the transit depth across the 8 year baseline of the observations. We rule out extended broadband absorption towards or beyond the Roche lobe at the wavelengths investigated, although observations with higher spectral resolution are required to determine if absorption out to those distances from the planet is present in individual near-UV lines.

With the aim of understanding the role of outflows in star formation, we performed a statistical study of the physical parameters of outflows in eleven massive protoclusters associated with ultra-compact HII regions. A total of 106 outflow lobes are identified in these protoclusters using the ALMA CO (3-2), HCN (4-3) and HCO+ (4-3) line observations. Although the position angles of outflow lobes do not differ in these three tracers, HCN and HCO+ tend to detect lower terminal velocity of the identified outflows compared to CO. The majority of the outflows in our targets are young with typical dynamical time-scales of 10^2-10^4 years, and are mostly composed of low-mass outflows along with at least one high-mass outflow in each target. An anti-correlation of outflow rate with dynamical time-scale indicates that the outflow rate possibly decreases with time. Also, a rising trend of dynamical time-scale with the mass of the associated core hints that the massive cores might have longer accretion histories than the low mass cores. Estimation of different energies in these protoclusters shows that outflows studied here cannot account for the generation of the observed turbulence, but can sustain the turbulence at the current epoch as the energy injection rate from the outflows is similar to the estimated dissipation rate.

We study the energy-dependent time lags and rms fractional amplitude of the kilohertz quasi-periodic oscillations (kHz QPOs) of a group of neutron-star low mass X-ray binaries (LMXBs). We find that for the lower kHz QPO the slope of the best-fitting linear model to the time-lag spectrum and the total rms amplitude integrated over the 2 to 25 keV energy band both decrease exponentially with the luminosity of the source. For the upper kHz QPO the slope of the time-lag spectrum is consistent with zero, while the total rms amplitude decreases exponentially with the luminosity of the source. We show that both the slope of the time-lag spectrum and the total rms amplitude of the lower kHz QPO are linearly correlated with a slope of ~1. Finally, we discuss the mechanism that could be responsible for the radiative properties of the kHz QPOs, with the variability originating in a Comptonising cloud or corona that is coupled to the innermost regions of the accretion disc, close to the neutron star.

We present a set of ultramassive white dwarf models, focused on masses above $1.3\,M_\odot$. Given the uncertainties about the formation and compositions of such objects, we construct parameterized model sequences, guided by evolutionary calculations including both single star and double white dwarf merger formation channels. We demonstrate that the cooling of objects with central densities in excess of $10^9\,\rm g\,cm^{-3}$ is dominated by neutrino cooling via the Urca process in the first $\approx 100$ Myr after formation. Our models indicate that the recently discovered ultramassive white dwarf ZTF J190132.9+145808.7 is likely to have experienced this Urca-dominated cooling regime. We also show that the high densities imply that diffusion is unlikely to significantly alter the core compositions of these objects before they crystallize.

Variable red- and blue-shifted absorption features observed in the Ca ii K line towards the A-type shell star $\phi$ Leo have been suggested by us in a previous work to be likely due to solid, comet-like bodies in the circumstellar (CS) environment. Our aim is to expand our observational study of this object to other characteristic spectral lines of A-type photospheres as well as to lines arising in their CS shells. We have obtained more than 500 high-resolution optical spectra collected at different telescopes from December 2015 to January 2019. We have analysed some photospheric lines, in particular Ca i 4226 \AA ~and Mg ii 4481 \AA, as well as the circumstellar shell lines Ca ii H\&K, Ca ii IR triplet, Fe ii, Ti ii, and the Balmer lines H$\alpha$ and H$\beta$. Our observational study reveals that $\phi$ Leo is a variable $\delta$ Scuti star whose spectra show remarkable dumps and bumps superimposed on the photospheric line profiles, which vary their strength and sharpness, propagate from blue- to more red-shifted radial velocities and persisting during a few hours, likely produced by non-radial pulsations. At the same time, all shell lines present an emission at $\sim$3 km/s centered at the core of the CS features, and two variable absorption minima at both sides of the emission. The variations observed in the Ca ii H\&K, Fe ii and Ti ii lines occur at any time scale from minutes to days and observing run, but without any clear correlation or recognizable temporal pattern among the different lines. In the case of H$\alpha$ the CS contribution is also variable in just one of the observing runs. We suggest that $\phi$ Leo is a rapidly rotating $\delta$ Scuti star surrounded by a variable, (nearly) edge-on CS disk possibly re-supplied by the $\delta$ Scuti pulsations.

The big rip, the little rip and the little sibling of the big rip are cosmological doomsdays predicted by some phantom dark energy models that could describe the future evolution of our Universe. When the universe evolves towards either of these future cosmic events all bounded structures and, ultimately, space-time itself are ripped apart. Nevertheless, it is commonly belief that quantum gravity effects may smooth or even avoid these classically predicted singularities. In this review, we discuss the classical and quantum occurrence of these riplike events in the scheme of metric $f(R)$ theories of gravity. The quantum analysis is performed in the framework of $f(R)$ quantum geometrodynamics. In this context, we analyse the fulfilment of the DeWitt criterion for the avoidance of these singular fates.

A surprising consequence of non-linear flavor evolution is the spontaneous breaking of the initial symmetries of the neutrino gas propagating in a dense astrophysical environment. We explore the flavor conversion physics within a three flavor framework, by taking into account the polar and azimuthal angular distributions of neutrinos. The very first example of spontaneous symmetry breaking in the context of fast flavor mixing is presented. Intriguingly, in the presence of spontaneous symmetry breaking, fast flavor mixing does not always develop in the proximity of crossings in the electron lepton number, as commonly considered, and large flavor mixing can affect all neutrino modes. Such findings are peculiar to non-azimuthally symmetric systems and are not predictable from the linear regime of the flavor evolution; they can have major consequences on the physics of compact astrophysical objects.

We construct a generally-covariant formulation of non-equilibrium thermodynamics in General Relativity. We find covariant entropic forces arising from gradients of the entropy density, and a corresponding non-conservation of the energy momentum tensor in terms of these forces. We also provide a Hamiltonian formulation of General Relativity in the context of non-equilibrium phenomena and write the Raychaudhuri equations for a congruence of geodesics. We find that a fluid satisfying the strong energy condition could avoid collapse for a positive and sufficiently large entropic-force contribution. We then study the forces arising from gradients of the bulk entropy of hydrodynamical matter, as well as the entropy of boundary terms in the action, like those of black hole horizons. Finally, we apply the covariant formulation of non-equilibrium thermodynamics to the expanding universe and obtain the modified Friedmann equations, with an extra term corresponding to an entropic force satisfying the second law of thermodynamics.

General relativistic entropic acceleration theory may explain the present cosmic acceleration from first principles without the need of introducing a cosmological constant. Following the covariant formulation of non-equilibrium phenomena in the context of a homogeneous and isotropic Friedmann-Lemaitre-Robertson-Walker (FLRW) metric, we find that the growth of entropy associated with the causal horizon of our universe (inside a finite bubble in eternal inflation) induces an acceleration that is essentially indistinguishable from that of $\Lambda$CDM, except for a slightly larger present rate of expansion compared to what would be expected from the CMB in $\Lambda$CDM, possibly solving the so-called $H_0$ tension. The matter content of the universe is unchanged and the coincidence problem is resolved since it is the growth of the causal horizon of matter that introduces this new relativistic entropic force. The cosmological constant is made unnecessary and the future hypersurface is Minkowsky rather than de Sitter.

The detection of gravitational radiation, emitted in the aftermath of the excitation of neutron star quasi-normal modes, has the potential to provide unprecedented access to the properties of matter in the star interior, and shed new light on the dynamics of nuclear interactions at microscopic level. Of great importance, in this context, will be the sensitivity to themodelling of three-nucleon interactions, which are known to play a critical role in the high-density regime. We report the results of a calculation of the frequencies and damping times of the fundamental mode, carried out using the equation of state of Akmal, Pandharipande and Ravenhall as a baseline, and varying the strength of the isoscalar repulsive term the Urbana IX potential within a range consistent with multimessenger astrophysical observations. The results of our analysis indicate that repulsive three-nucleon interactions strongly affect the stiffness of the equation of state, which in turn determines the pattern of the gravitational radiation frequencies, largely independent of the mass of the source. The observational implications are also discussed.

We describe the application of the lattice covering problem to the placement of templates in a search for continuous gravitational waves from the low-mass X-Ray binary Scorpius X-1. Efficient placement of templates to cover the parameter space at a given maximum mismatch is an application of the sphere covering problem, for which an implementation is available in the LatticeTiling software library. In the case of Sco X-1, potential correlations, in both the prior uncertainty and the mismatch metric, between the orbital period and orbital phase, lead to complications in the efficient construction of the lattice. We define a shearing coordinate transformation which simultaneously minimizes both of these sources of correlation, and allows us to take advantage of the small prior orbital period uncertainty. The resulting lattices have a factor of about 3 fewer templates than the corresponding parameter space grids constructed by the prior straightforward method, allowing a more sensitive search at the same computing cost and maximum mismatch.