Papers for Friday, May 07 2021

Context. Infrared spectroscopy of star and planet forming regions is at the dawn of a new age with the upcoming James Webb Space Telescope. In support of these observations, laboratory spectra are required to identify complex organic molecules in the ices that cover the dust grains in these regions. Aims. This study aims to provide reference spectra to firmly detect icy methyl formate in the different stages of star and planet forming regions. Methyl formate is mixed in astronomically relevant matrices, and the peak positions, FWHMs, and relative band intensities are characterized for different temperatures to provide an analytical tool for astronomers. Methods. Methyl formate is deposited at 15 K under high-vacuum conditions. Specifically, methyl formate is deposited pure and mixed with CO, H$_2$CO, CH$_3$OH, H$_2$O, and CO:H$_2$CO:CH$_3$OH combined. Throughout the experiment infrared spectra are acquired with a FTIR spectrometer in the range from 4000-500 cm$^{-1}$ (2.5-20 $\mu$m) at a spectral resolution of 0.5 cm$^{-1}$. Results. We present the characterization of five solid-state methyl formate vibrational modes in pure and astronomically relevant ice matrices. The five selected vibrational modes, namely the C=O stretch, C$-$O stretch, CH$_3$ rocking, O$-$CH$_3$ stretching, and OCO deformation, are best suited for a JWST identification of methyl formate. For each of these vibrational modes, and each of the mixtures the TvS heatmaps, peak position versus FWHM, and relative band intensities are given. Additionally, the acquired reference spectra of methyl formate are compared with Spitzer observations of HH 46. A tentative detection of methyl formate provides an upper limit to the column density of $1.7\times10^{17}$ cm$^{-2}$, corresponding to an upper limit relative to water of $\leq 2.2\%$ and $\leq 40\%$ with respect to methanol.

Common envelope evolution (CEE) is believed to be an important stage in the evolution of binary/multiple stellar systems. Following this stage, the CE is thought to be ejected, leaving behind a compact binary (or a merger product). Although extensively studied, the CEE process is still little understood, and although most binaries have non-negligible eccentricity, the effect of initial eccentricity on the CEE has been little explored. Moreover, most studies assume a complete circularization of the orbit by the CE onset, while observationally such eccentricities are detected in many post-CE binaries. Here we use smoothed particle hydro-dynamical simulations (SPH) to study the evolution of initially eccentric ($0\le e\le0.95$) CE-systems. We find that initially eccentric binaries only partially circularize. In addition, higher initial eccentricity leads to a higher post-CE eccentricity, with eccentricities of post-CE binaries as high as 0.2 in the most eccentric cases, and even higher if the initial peri-center of the orbit is located inside the star (e.g. following a kick into an eccentric orbit, rather than a smooth transition). CEE of more eccentric binaries leads to enhanced dynamical mass-loss of the CE compared with more circular binaries, and depends on the initial closest approach of the binary. We show that our results and the observed eccentricities of post-CE binaries suggest that the typical assumptions of circular orbits following CEE should potentially be revised. We expect post-CE eccentricities to affect the delay time distributions of various transients such as supernovae, gamma-ray bursts and gravitational-wave sources by up to tens of percents.

Merging compact black-hole (BH) binaries are likely to exist in the nuclear star clusters around supermassive BHs (SMBHs), such as Sgr A$^\ast$. They may also form in the accretion disks of active galactic nuclei. Such compact binaries can emit gravitational waves (GWs) in the low-frequency band (0.001-1 Hz) that are detectable by several planned space-borne GW observatories. We show that the orbital axis of the compact binary may experience significant variation due to the frame-dragging effect associated with the spin of the SMBH. The dynamical behavior of the orbital axis can be understood analytically as a resonance phenomenon. We show that rate of change of the binary orbital axis encodes the information on the spin of the SMBH. Therefore detecting GWs from compact binaries around SMBHs, particularly the modulation of the waveform associated with the variation of the binary orbital axis, can provide a new probe on the spins of SMBHs.

Using the IllustrisTNG cosmological galaxy formation simulations, we analyze the physical properties of young quiescent galaxies at z=2 with stellar masses above 10^10.5 solar masses. This key population provides an unaltered probe into the evolution of galaxies from star-forming to quiescent, and has been recently targeted by several observational studies. Young quiescent galaxies in the simulations do not appear unusually compact, in tension with observations, but they show unique age profiles that are qualitatively consistent with the observed color gradients. In particular, more than half of the simulated young quiescent galaxies show positive age profiles due to recent intense central starbursts, which are triggered by significant mergers. Yet, there is a sizable population of recently quenched galaxies without significant mergers and with flat age profiles. Our results suggest that mergers play a fundamental role in structural transformation, but are not the only available pathway to quench a z=2 galaxy.

We use the induced geometry on the two dimensional transverse cross section of a photon beam propagating on a perturbed Friedmann-Robertson-Walker (FRW) spacetime to find the Cosmic Microwave Background (CMB) photon distribution over a telescope's collecting area today. It turns out that at each line of sight the photons are diluted along a transverse direction due to gravitational shearing. The effect can be characterized by two spin-weight-two variables, which are reminiscent of the Stokes polarization parameters. Similar to that case, one can construct a scalar and a pseudo-scalar function where the latter only gets contributions from the tensor modes. We analytically determine the power spectrum of the pseudo-scalar at superhorizon scales in a simple inflationary model and briefly discuss possible observational consequences.

In the era of precision cosmology and ever-improving cosmological simulations, a better understanding of different galaxy components such as bulges and discs will give us new insight into galactic formation and evolution. Based on the fact that the stellar populations of the constituent components of galaxies differ by their dynamical properties, we develop two simple models for galaxy decomposition using the IllustrisTNG cosmological hydrodynamical simulation. The first model uses a single dynamical parameter and can distinguish 4 components: thin disc, thick disc, counter-rotating disc and bulge. The second model uses one more dynamical parameter, was defined in a probabilistic manner, and distinguishes two components: bulge and disc. The number fraction of disc-dominated galaxies at a given stellar mass obtained by our models agrees well with observations for masses exceeding $ \log_{10}(M_*/M_\odot)=10$. S\'ersic indices and half-mass radii for the bulge components agree well with those for real galaxies. The mode of the distribution of S\'ersic indices for the disc components is at the expected value of $n=1$. However, disc half-mass radii are smaller than those for real galaxies, in accordance with previous findings that the IllustrisTNG simulation produces undersized discs.

We explore variations of the dust extinction law of the Milky Way by selecting stars from the Swift/UVOT Serendipitous Source Catalogue, cross-matched with Gaia DR2 and 2MASS to produce a sample of 10,452 stars out to ~4kpc with photometry covering a wide spectral window. The near ultraviolet passbands optimally encompass the 2175A bump, so that we can simultaneously fit the net extinction, quoted in the V band (A$_V$), the steepness of the wavelength dependence ($\delta$) and the bump strength (E$_b$). The methodology compares the observed magnitudes with theoretical stellar atmospheres from the models of Coelho. Significant correlations are found between these parameters, related to variations in dust composition, that are complementary to similar scaling relations found in the more complex dust attenuation law of galaxies - that also depend on the distribution of dust among the stellar populations within the galaxy. We recover the strong anticorrelation between A$_V$ and Galactic latitude, as well as a weaker bump strength at higher extinction. $\delta$ is also found to correlate with latitude, with steeper laws towards the Galactic plane. Our results suggest that variations in the attenuation law of galaxies cannot be fully explained by dust geometry.

There has been much interest in novel models of dark matter that exhibit interesting behavior on galactic scales. A primary motivation is the observed Baryonic Tully-Fisher Relation in which the mass of galaxies increases as the quartic power of rotation speed. This scaling is not obviously accounted for by standard cold dark matter. This has prompted the development of dark matter models that exhibit some form of so-called MONDian phenomenology to account for this galactic scaling, while also recovering the success of cold dark matter on large scales. A beautiful example of this are the so-called superfluid dark matter models, in which a complex bosonic field undergoes spontaneous symmetry breaking on galactic scales, entering a superfluid phase with a 3/2 kinetic scaling in the low energy effective theory, that mediates a long-ranged MONDian force. In this work we examine the causality and locality properties of these and other related models. We show that the Lorentz invariant completions of the superfluid models exhibit high energy perturbations that violate global hyperbolicity of the equations of motion in the MOND regime and can be superluminal in other parts of phase space. We also examine a range of alternate models, finding that they also exhibit forms of non-locality.

The majority of non-merging stellar mass black holes are discovered by observing high energy emission from accretion processes. Here we pursue the large, but still mostly unstudied population of non-interacting black holes and neutron stars by searching for the tidally-induced ellipsoidal variability of their stellar companions. We start from a sample of about 200,000 rotational variables, semi-regular variables, and eclipsing binary stars from the All-Sky Automated Survey for Supernovae (ASAS-SN). We use a $\chi^2$ ratio test followed by visual inspection to identify 369 candidates for ellipsoidal variability. We also discuss how to combine the amplitude of the variability with mass and radius estimates for observed stars to calculate a minimum companion mass, identifying the most promising candidates for high mass companions.

The noteworthy four-planet HR 8799 system teeters on the brink of gravitational instability and contains an A-type host star which is characteristic of the progenitors of the majority of known white dwarf planetary system hosts. Gozdziewski and Migaszewski (2020) have demonstrated that the system can retain all four planets for at least 1 Gyr along the main sequence if the planets evolve within an externally unperturbed 8:4:2:1 mean motion resonance configuration. Here we propagate forward their most stable fit beyond the main sequence, and incorporate external effects from Galactic tides and stellar flybys. We find that (i) giant branch mass loss always breaks the resonance, and usually triggers the ejection of two of the planets, (ii) stellar flybys and Galactic tides rarely break the resonance during the main-sequence and giant branch phases, but play a crucial role in determining the final planetary configurations around the eventual white dwarf host star, and (iii) the meanderings of the surviving planets vary significantly, occupying regions from under 1 au to thousands of au. The ubiquitous survival of at least one planet and the presence of the debris discs in the system should allow for dynamical pathways for the white dwarf to be metal-polluted.

The atmospheres of synchronously rotating exoplanets are intrinsically three-dimensional, and fast vertical and horizontal winds are expected to mix the atmosphere, driving the chemical composition out of equilibrium. Due to the longer computation times associated with multi-dimensional forward models, horizontal mixing has only been investigated for a few case studies. In this paper, we aim to generalize the impact of horizontal and vertical mixing on the chemistry of exoplanet atmospheres over a large parameter space. We do this by applying a sequence of post-processed forward models for a large grid of synchronously rotating gaseous exoplanets, where we vary the effective temperature (between 400 K and 2600 K), surface gravity, and rotation rate. We find that there is a dichotomy in the horizontal homogeneity of the chemical abundances. Planets with effective temperatures below 1400 K tend to have horizontally homogeneous, vertically quenched chemical compositions, while planets hotter than 1400 K exhibit large compositional day-night differences for molecules such as methane. Furthermore, we find that the planet's rotation rate impacts the planetary climate, and thus also the molecular abundances and transmission spectrum. By employing a hierarchical modelling approach, we assess the relative importance of disequilibrium chemistry on the exoplanet transmission spectrum, and conclude that the temperature has the most profound impact. Temperature differences are also the main cause of limb asymmetries, which we estimate could be observable with the James Webb Space Telescope. This work highlights the value of applying a consistent modelling setup to a broad parameter space in exploratory theoretical research.

We present a detailed 3D kinematic analysis of the central regions ($R<30''$) of the low-mass and dynamically evolved galactic globular cluster NGC 6362. The study is based on data obtained with ESO-VLT/MUSE used in combination with the adaptive optics module and providing $\sim3000$ line-of-sight radial velocities, which have been complemented with Hubble Space Telescope proper motions. The quality of the data and the number of available radial velocities allowed us to detect for the first time a significant rotation signal along the line of sight in the cluster core with amplitude of $\sim 1$ km/s and with a peak located at only $\sim20''$ from the cluster center, corresponding to only $\sim10\%$ of the cluster half-light radius. This result is further supported by the detection of a central and significant tangential anisotropy in the cluster innermost regions. This is one of the most central rotation signals ever observed in a globular cluster to date. We also explore the rotational properties of the multiple populations hosted by this cluster and find that Na-rich stars rotate about two times more rapidly than the Na-poor sub-population thus suggesting that the interpretation of the present-day globular cluster properties require a multi-component chemo-dynamical approach. Both the rotation amplitude and peak position would fit qualitatively the theoretical expectations for a system that lost a significant fraction of its original mass because of the long-term dynamical evolution and interaction with the Galaxy. However, to match the observations more quantitatively further theoretical studies to explore the initial dynamical properties of the cluster are needed.

We present a detailed analysis of time variability of two distinct C IV broad absorption line (BAL) components seen in the spectrum of J162122.54+075808.4 ($z_{em}$ = 2.1394) using observations from SDSS, NTT and SALT taken at seven different epochs spanning about 15 years. The blue-BAL component (with an ejection velocity, $v_{\rm e}\sim37,500$ kms$^{-1}$) is an emerging absorption that shows equivalent width variations and kinematic shifts consistent with acceleration. The red-BAL component ($v_{\rm e} \sim 15,400$ kms$^{-1}$) is a three component absorption. One of the components is emerging and subsequently disappearing. The two other components show kinematic shifts consistent with acceleration coupled with equivalent width variability. Interestingly, we find the kinematic shifts and equivalent width variability of the blue- and red-BAL components to be correlated. While the C IV emission line flux varies by more than 17% during our monitoring period, the available light-curves (covering rest frame 1300-2300 angstrom do not show more than a 0.1 mag variability in the continuum. This suggests that the variations in the ionizing flux are larger than that of the near-UV flux. However, the correlated variability seen between different BAL components cannot be explained solely by photoionization models without structural changes. In the framework of disk wind models, any changes in the radial profiles of density and/or velocity triggered either by disk instabilities or by changes in the ionizing radiation can explain our observations. High resolution spectroscopic monitoring of J1621+0758 is important to understand the physical conditions of the absorbing gas and thereby to constrain the parameters of disk-wind models.

A late end to reionisation at redshift $z\simeq 5.3$ is consistent with observed spatial variations in the Ly$\alpha$ forest transmission and the deficit of Ly$\alpha$ emitting galaxies around extended Ly$\alpha$ absorption troughs at $z=5.5$. In this model, large islands of neutral hydrogen should persist in the diffuse intergalactic medium (IGM) until $z\simeq 6$. We use a novel, hybrid approach that combines high resolution cosmological hydrodynamical simulations with radiative transfer to predict the incidence of strong 21 cm forest absorbers with optical depths $\tau_{21}>10^{-2}$ from the diffuse IGM in these late reionisation models. We include the effect of redshift space distortions on the simulated 21 cm forest spectra, and treat the highly uncertain heating of the pre-reionisation IGM by soft X-rays as a free parameter. For a model with only modest IGM pre-heating, such that average gas kinetic temperatures in the diffuse IGM remain below $T_{\rm K}\simeq 10^{2} \rm\, K$, we find that strong 21 cm forest absorption lines should persist until $z=6$. For a sample of $\sim 10$ sufficiently radio loud background sources, a null-detection of 21 cm forest absorbers at $z\simeq 6$ with SKA1-low or possibly LOFAR should provide an informative lower limit on the still largely unconstrained soft X-ray background at high redshift and the temperature of the pre-reionisation IGM.

Sun-like stars form from the contraction of cold and dense interstellar clouds. How the collapse proceeds and the main physical processes driving it, however, are still under debate and a final consensus on the timescale of the process has not been reached. Does this contraction proceed slowly, sustained by strong magnetic fields and ambipolar diffusion, or is it driven by fast collapse with gravity dominating the entire process? One way to answer this question is to measure the age of prestellar cores through statistical methods based on observations or via reliable chemical chronometers, which should better reflect the physical conditions of the cores. Here we report APEX observations of ortho-H$_2$D$^+$ and para-D$_2$H$^+$ for six cores in the Ophiuchus complex and combine them with detailed three-dimensional magneto-hydrodynamical simulations including chemistry, providing a range of ages for the observed cores of 100-200 kyr. The outcome of our simulations and subsequent analysis provides a good matching with the observational results in terms of physical (core masses and volume densities) and dynamical parameters such as the Mach number and the virial parameter. We show that models of fast collapse successfully reproduce the observed range of chemical abundance ratios as the timescales to reach the observed stages is shorter than the free-fall time of the cores and much shorter than the ambipolar diffusion time, measured from the electron fraction in the simulations. Our work establishes the ortho-H$_2$D$^+$/para-D$_2$H$^+$ ratio as a reliable chemical clock and opens up to the possibility of exploring the star formation process in a statistically relevant sample through observations of these tracers.

The 3D matter power spectrum, $P_{\delta}(k,z)$ is a fundamental quantity in the analysis of cosmological data such as large-scale structure, 21cm observations, and weak lensing. Existing computer models (Boltzmann codes) such as CLASS can provide it at the expense of immoderate computational cost. In this paper, we propose a fast Bayesian method to generate the 3D matter power spectrum, for a given set of wavenumbers, $k$ and redshifts, $z$. Our code allows one to calculate the following quantities: the linear matter power spectrum at a given redshift (the default is set to 0); the non-linear 3D matter power spectrum with/without baryon feedback; the weak lensing power spectrum. The gradient of the 3D matter power spectrum with respect to the input cosmological parameters is also returned and this is useful for Hamiltonian Monte Carlo samplers. The derivatives are also useful for Fisher matrix calculations. In our application, the emulator is accurate when evaluated at a set of cosmological parameters, drawn from the prior, with the fractional uncertainty, $\Delta P_{\delta}/P_{\delta}$ centered on 0. It is also $\sim 300$ times faster compared to CLASS, hence making the emulator amenable to sampling cosmological and nuisance parameters in a Monte Carlo routine. In addition, once the 3D matter power spectrum is calculated, it can be used with a specific redshift distribution, $n(z)$ to calculate the weak lensing and intrinsic alignment power spectra, which can then be used to derive constraints on cosmological parameters in a weak lensing data analysis problem. The software ($\texttt{emuPK}$) can be trained with any set of points and is distributed on Github, and comes with with a pre-trained set of Gaussian Process (GP) models, based on 1000 Latin Hypercube (LH) samples, which follow roughly the current priors for current weak lensing analyses.

AM CVn systems are a rare type of accreting binary that consists of a white dwarf and a helium-rich, degenerate donor star. Using the Zwicky Transient Facility (ZTF), we searched for new AM CVn systems by focusing on blue, outbursting stars. We first selected outbursting stars using the ZTF alerts. We cross-matched the candidates with $Gaia$ and Pan-STARRS catalogs. The initial selection of candidates based on the $Gaia$ $BP$-$RP$ contains 1751 unknown objects. We used the Pan-STARRS $g$-$r$ and $r$-$i$ color in combination with the $Gaia$ color to identify 59 high-priority candidates. We obtained identification spectra of 35 sources, of which 18 are high priority candidates, and discovered 9 new AM CVn systems and one magnetic CV which shows only He-II lines. Using the outburst recurrence time, we estimate the orbital periods which are in the range of 29 to 50 minutes. We conclude that targeted followup of blue, outbursting sources is an efficient method to find new AM CVn systems, and we plan to followup all candidates we identified to systematically study the population of outbursting AM CVn systems.

The role of large amplitude whistler waves in the energization and scattering of solar wind electrons has long been an interesting problem in Space Physics. To study this wave-particle interaction, we developed a vectorized test particle simulation with a variational calculation of the Lyapunov exponents. From using secular perturbation theory on this Hamiltonian system of wave and particle, we confirmed that the pitch angle diffusion of the particle was along the constant Hamiltonian surface and that it was driven by the interaction with the resonance surfaces. We also showed that oblique whistlers could efficiently scatter field-aligned strahl electrons into the halo population in the solar wind. We demonstrated through simulation that these waves were capable of generating horn-like features in the velocity distribution function, similar to recent PIC simulation results in the literature.

Near-Sun Comet C/2019 Y4 (ATLAS) is the first member of a long-period comet group observed to disintegrate well before perihelion. Here we present our investigation into this disintegration event using images obtained in a 3-day {\it Hubble Space Telescope} (\hst) campaign. We identify two fragment clusters produced by the initial disintegration event, corresponding to fragments C/2019 Y4-A and C/2019 Y4-B identified in ground-based data. These two clusters started with similar integrated brightness, but exhibit different evolutionary behavior. C/2019 Y4-A was much shorter-lived compared to C/2019 Y4-B, and showed signs of significant mass-loss and changes in size distribution throughout the 3-day campaign. The cause of the initial fragmentation is undetermined by the limited evidence but crudely compatible with either the spin-up disruption of the nucleus or runaway sublimation of sub-surface supervolatile ices, either of which would lead to the release of a large amount of gas as inferred from the significant bluing of the comet observed shortly before the disintegration. Gas can only be produced by the sublimation of volatile ices, which must have survived at least one perihelion passage at a perihelion distance of $q=0.25$~au. We speculate that Comet ATLAS is derived from the ice-rich interior of a non-uniform, kilometer-wide progenitor that split during its previous perihelion. This suggests that comets down to a few kilometers in diameter can still possess complex, non-uniform interiors that can protect ices against intense solar heating.

We present the results from a full polarization study carried out with ALMA during the first VLBI campaign, which was conducted in Apr 2017 in the $\lambda$3mm and $\lambda$1.3mm bands, in concert with the Global mm-VLBI Array (GMVA) and the Event Horizon Telescope (EHT), respectively. We determine the polarization and Faraday properties of all VLBI targets, including Sgr A*, M87, and a dozen radio-loud AGN. We detect high linear polarization fractions (2-15%) and large rotation measures (RM $>10^{3.3}-10^{5.5}$ rad m$^{-2}$). For Sgr A* we report a mean RM of $(-4.2\pm0.3) \times10^5$ rad m$^{-2}$ at 1.3 mm, consistent with measurements over the past decade, and, for the first time, an RM of $(-2.1\pm0.1) \times10^5$ rad m$^{-2}$ at 3 mm, suggesting that about half of the Faraday rotation at 1.3 mm may occur between the 3 mm photosphere and the 1.3 mm source. We also report the first unambiguous measurement of RM toward the M87 nucleus at mm wavelengths, which undergoes significant changes in magnitude and sign reversals on a one year time-scale, spanning the range from -1.2 to 0.3 $\times\,10^5$ rad m$^{-2}$ at 3 mm and -4.1 to 1.5 $\times\,10^5$ rad m$^{-2}$ at 1.3 mm. Given this time variability, we argue that, unlike the case of Sgr A*, the RM in M87 does not provide an accurate estimate of the mass accretion rate onto the black hole. We put forward a two-component model, comprised of a variable compact region and a static extended region, that can simultaneously explain the polarimetric properties observed by both the EHT and ALMA. These measurements provide critical constraints for the calibration, analysis, and interpretation of simultaneously obtained VLBI data with the EHT and GMVA.

We forecast the reionization history constraints, inferred from Lyman-alpha damping wing absorption features, for a future sample of $\sim 20$ $z \geq 6$ gamma-ray burst (GRB) afterglows. We describe each afterglow spectrum by a three-parameter model. First, L characterizes the size of the ionized region (the "bubble size") around a GRB host halo. Second, $\langle{x_{\rm HI}\rangle}$ is the volume-averaged neutral fraction outside of the ionized bubble around the GRB, which is approximated as spatially uniform. Finally, $N_{\mathrm{HI}}$ denotes the column-density of a local damped Lyman-alpha absorber (DLA) associated with the GRB host galaxy. The size distribution of ionized regions is extracted from a numerical simulation of reionization, and evolves strongly across the Epoch of Reionization (EoR). The model DLA column densities follow the empirical distribution determined from current GRB afterglow spectra. We use a Fisher matrix formalism to forecast the $\langle{x_{\rm HI}(z)\rangle}$ constraints that can be obtained from follow-up spectroscopy of afterglows with SNR = 20 per R=3,000 resolution element at the continuum. We find that the neutral fraction may be determined to better than 10-15\% (1-$\sigma$) accuracy from this data across multiple independent redshift bins at $z \sim 6-10$, spanning much of the EoR, although the precision degrades somewhat near the end of reionization. A more futuristic survey with $80$ GRB afterglows at $z \geq 6$ can improve the precision here by a factor of $2$ and extend measurements out to $z \sim 14$. We further discuss how these constraints may be combined with estimates of the escape fraction of ionizing photons, derived from the DLA column density distribution towards GRBs extracted at slightly lower redshift. This combination will help in testing whether we have an accurate census of the sources that reionized the universe.

We present the results of a search for optical candidates of Ultraluminous X-ray sources (ULXs) in two dwarf galaxies: NGC 4861 and NGC 4449 using Hubble Space Telescope {\it HST} archival data. With a precise astrometry, we confirm that NGC 4861 X1 associated with an HII complex as reported by \cite{2014MNRAS.441.1841T} and we conclude that NGC 4861 X2 resides in a young star group with 400$\pm$80 M$\odot$. We also find that NGC 4449 X7 is associated with three optical candidates within an error radius of 0$\farcs$2 at the 90$\%$ confidence level. Absolute magnitudes (M$_{v}$) of these candidates are determined as $-$5.0 and $-$4.1. The age and mass values for the three candidates by using stellar evolutionary tracks are estimated as 40-50 Myr and $\sim$8 M$\odot$, respectively. The locations of optical candidates suggest possible association with a nearby group of stars. In addition, we analyzed the previously unused archival data of {\it XMM-Newton}, {\it Chandra} and {\it Swift} where the sources detected. Although the X-ray spectral data do not allow us to discriminate between physical models, long-term data at hand are consistent with the sources being in luminous hard states.

Planar acoustically dominated magneto\-hydro\-dynamic waves are initiated at the high-$\beta$ base of a simulated 2D isothermal stratified atmosphere with potential magnetic field exhibiting both open and closed field regions as well as neutral points. They shock on their way upward toward the Alfv\'en-acoustic equipartition surface $a=c$, where $a$ and $c$ are the Alfv\'en and sound speeds respectively. Expanding on recent 1.5D findings that such shocks mode-convert to fast shocks and slow smoothed waves on passing through $a=c$, we explore the implications for these more complex magnetic geometries. It is found that the 1.5D behaviour carries over to the more complex case, with the fast shocks strongly attracted to neutral points, which are disrupted producing extensive fine structure. It is also observed that shocks moving in the opposite direction, from $a>c$ to $a<c$, split into fast and slow components too, and that again it is the slow component that is smoothed.

Wide field planetary camera 2 (WFPC2) exposures are already some 20 years older than Gaia epoch observations, or future JWST observations. As such, they offer an unprecedented time baseline for high-precision proper-motion studies, provided the full astrometric potential of these exposures is reached. We have started such a project with the work presented here being its first step. We explore geometric distortions beyond the well-known ones published in the early 2000s. This task is accomplished by using the entire database of WFPC2 exposures in filters F555W, F606W and F814W and three standard astrometric catalogs: Gaia EDR3, 47 Tuc and $\omega$Cen. The latter two were constructed using HST observations made with cameras other than WFPC2. We explore a suite of centering algorithms, and various distortion maps in order to understand and quantify their performance. We find no high-frequency systematics beyond the 34th-row correction, down to a resolution of 10 pixels. Low-frequency systematics starting at a resolution of 50-pixels are present at a level of 30-50 millipix (1.4-2.3 mas) for the PC and 20-30 millipix (2-3 mas) for the WF chips. We characterize these low-frequency systematics by providing correction maps and updated cubic-distortion coefficients for each filter.

We present a tidal model for treating the rotational evolution in the general three-body problem with arbitrary viscosities, in which all the masses are considered to be extended and all the tidal interactions between pairs are taken into account. Based on the creep tide theory, we present the set of differential equations that describes the rotational evolution of each body, in a formalism that is easily extensible to the N tidally-interacting body problem. We apply our model to the case of a circumbinary planet and use a Kepler-38 like binary system as a working example. We find that, in this low planetary eccentricity case, the most likely final stationary rotation state is the 1:1 spin-orbit resonance, considering an arbitrary planetary viscosity inside the estimated range for the solar system planets. We derive analytical expressions for the mean rotational stationary state, based on high-order power series of the semimajor axes ratio a1 /a2 and low-order expansions of the eccentricities. These are found to reproduce very accurately the mean behaviour of the low-eccentric numerical integrations for arbitrary planetary relaxation factors, and up to a1/a2 \sim 0.4. Our analytical model is used to predict the stationary rotation of the Kepler circumbinary planets and find that most of them are probably rotating in a sub-synchronous state, although the synchrony shift is much less important than the one estimated in Zoppetti et al. (2019, 2020). We present a comparison of our results with those obtained with the Constant Time Lag and find that, unlike what we assumed in our previous works, the cross torques have a non-negligible net secular contribution, and must be taken into account when computing the tides over each body in an N-extended-body system from an arbitrary reference frame. These torques are naturally taken into account in the creep theory.

Active galactic nuclei (AGNs) can act as black hole assembly lines, funneling some of the stellar-mass black holes from the vicinity of the galactic center into the inner plane of the AGN disk where the black holes can merge through dynamical friction and gravitational wave emission. Here, we show that stars near the galactic center are also brought into the AGN disk, where they can be tidally disrupted by the stellar-mass black holes in the disk. Such micro-tidal disruption events (micro-TDEs) could be useful probe of stellar interaction with the AGN disk. We find that micro-TDEs in AGNs occur at a rate of $\sim170$ Gpc$^{-3}$yr$^{-1}$. Their cleanest observational probe may be the detection of tidal disruption in AGNs by heavy supermassive black holes ($M_{\bullet}\gtrsim10^{8}$ M$_{\odot}$) so that cannot tidally disrupt solar-type stars. We discuss two such TDE candidates observed to date (ASASSN-15lh and ZTF19aailpwl).

We present the results of a 3D global magnetohydrodynamic (MHD) simulation of an AM CVn system that was aimed at exploring eccentricity growth in the accretion disc self-consistently from a first principles treatment of the MHD turbulence. No significant eccentricity growth occurs in the simulation. In order to investigate the reasons why, we ran 2D alpha disc simulations with alpha values of 0.01, 0.1, and 0.2, and found that only the latter two exhibit significant eccentricity growth. We present an equation expressing global eccentricity evolution in terms of contributing forces and use it to analyze the simulations. As expected, we find that the dominant term contributing to the growth of eccentricity is the tidal gravity of the companion star. In the 2D simulations, the alpha viscosity directly contributes to eccentricity growth. In contrast, the overall magnetic forces in the 3D simulation damp eccentricity. We also analyzed the mode-coupling mechanism of Lubow, and confirmed that the spiral wave excited by the 3:1 resonance was the dominant contributor to eccentricity growth in the 2D $\alpha=0.1$ simulations, but other waves also contribute significantly. We found that the $\alpha=0.1$ and 0.2 simulations had more relative mass at larger radii compared to the $\alpha=0.01$ and 3D MHD simulation, which also had an effective $\alpha$ of 0.01. This suggests that in 3D MHD simulations without sufficient poloidal magnetic flux, MRI turbulence does not saturate at a high enough $\alpha$ to spread the disc to large enough radii to reproduce the superhumps observed in real systems.

The vertical shear instability (VSI) is a hydrodynamical instability that requires rapid gas cooling and has been suggested to operate in outer regions of protoplanetary disks. The VSI drives turbulence with strong vertical motions, which could regulate the dust growth and settling. However, dust growth and settling can regulate the VSI because dust depletion makes gas cooling inefficient in outer disk regions that are optically thin to their own thermal emission. In this study, we quantify this potentially stabilizing effects of dust evolution on the VSI based on the linear analysis. We construct a model for calculating the cooling timescale, taking into account dust growth beyond micron sizes and size-dependent settling. Combining the model with the linear stability analysis, we map the region where the VSI operates, which we call the VSI zone, and estimate the maximum growth rate at each radial position. We find that dust growth as well as settling makes the VSI zone more confined around the midplane. This causes a decrease in the growth rate because the vertical shear of the rotation velocity, which is the source of the instability, is weaker at lower altitude. In our default disk model with 0.01 solar masses, dust growth from 10 micron to 1 mm causes a decrease in the growth rate by a factor of more than 10. The suppression of VSI-driven turbulence by dust evolution may promote further dust evolution in the outer regions and also explain a high degree of dust settling observed in the disk around HL Tau.

We perform a stacking analysis of the HI spectra from the Arecibo Legacy Fast ALFA (ALFALFA) survey for optically-selected local galaxies from the Sloan Digital Sky Survey (SDSS) to study the average gas fraction of galaxies at fixed stellar mass ($M_*$) and star formation rate (SFR). We first confirm that the average gas fraction strongly depends on the stellar mass and SFR of host galaxies; massive galaxies tend to have a lower gas fraction, and actively star-forming galaxies show higher gas fraction, which is consistent with many previous studies. Then we investigate the morphological dependence of the HI gas mass fraction at fixed $M_*$ and SFR to minimize the effects of these parameters. We use three morphological classifications based on parametric indicator (S\'{e}rsic index), non-parametric indicator (C-index), and visual inspection (smoothness from the Galaxy Zoo 2 project) on the optical image. We find that there is no significant morphological dependence of the HI gas mass fraction at fixed $M_*$ and SFR when we use C-index. In comparison, there exists a hint of diminishment in the HI gas mass fraction for "smooth" galaxies compared with "non-smooth" galaxies. We find that the visual smoothness is sensitive to the existence of small-scale structures in a galaxy. Our result suggests that even at fixed $M_*$ and SFR, the presence of such small-scale structures (seen in the optical image) is linked to their total HI gas content.

We present the results of multi-epoch, multi-frequency monitoring of a blazar 4C +29.45, which was regularly monitored as part of the Interferometric Monitoring of GAmma-ray Bright Active Galactic Nuclei (iMOGABA) program - a key science program of the Korean Very long baseline interferometry Network (KVN). Observations were conducted simultaneously at 22, 43, 86, and 129 GHz during the 4 years from 5 December 2012 to 28 December 2016. We also used additional data from the 15 GHz Owens Valley Radio Observatory (OVRO) monitoring program. From the 15 GHz light curve, we estimated the variability time scales of the source during several radio flux enhancements. We found that the source experienced 6 radio flux enhancements with variability time scales of 9-187 days during the observing period, yielding corresponding variability Doppler factors of 9-27. From the multi-frequency simultaneous KVN observations, we were able to obtain accurate radio spectra of the source and hence to more precisely measure the turnover frequencies $\nu_{\rm r}$ of synchrotron self-absorbed (SSA) emission with a mean value of $\nu_{rm r}$ = 28.9 GHz. Using jet geometry assumptions, we estimated the size of the emitting region at the turnover frequency. We found that the equipartition magnetic field strength is up to two orders of magnitudes higher than the SSA magnetic field strength (0.001-0.1 G). This is consistent with the source being particle dominated. We performed a careful analysis of the systematic errors related to making these estimations. From the results, we concluded that the equipartition region locates upstream the SSA region.

We report the result of a search for neutrinos in coincidence with solar flares from the GOES flare database. The search was performed on a 10.8 kton-year exposure of KamLAND collected from 2002 to 2019. We found no statistical excess of neutrinos and established 90\% confidence level upper limits of $8.4 \times 10^7$\,cm$^{-2}$ ($3.0 \times 10^{9}$\,cm$^{-2}$) on electron anti-neutrino (electron neutrino) fluence at 20\,MeV normalized to the X12 flare, assuming that the neutrino fluence is proportional to the X-ray intensity. The 90\% C.L. upper limits from this work exclude the entire region of parameter space associated with the Homestake event excess for the large solar flare in 1991.

The goal of this work is to study magnetic fields of the cool, eclipsing binary star UV Piscium (UV Psc). This system contains two active late-type stars, UV Psc A (G5V) and B (K3V). To obtain a complete picture, the properties of both global and local magnetic field structures are studied for both components. High-resolution intensity and circular polarisation spectra, collected in 2016 with the ESPaDOnS spectropolarimeter at the CFHT, were used to analyse the magnetic field of UV Psc. To increase the signal-to-noise ratio, the multi-line technique of least-squares deconvolution (LSD) was used to obtain average Stokes IV profiles. Then, a Zeeman-Doppler imaging (ZDI) code was employed to obtain the large-scale magnetic field topology and brightness distribution for both components of UV Psc. In addition, the small-scale magnetic fields, not visible to ZDI, were studied using the Zeeman intensification of Fe I lines. Maps of the surface magnetic field for both components of UV Psc were obtained, the large-scale magnetic fields feature strong toroidal and non-axisymetric components. UV Psc A and B have average global field strengths of 137 G and 88 G, respectively. The small-scale fields are notably stronger, with average strengths of 2.5 and 2.2 kG, respectively. Our results are in agreement with previous studies of partly-convective stars. Overall, UV Psc A has a stronger magnetic field compared to UV Psc B. Due to the eclipsing binary geometry, certain magnetic field features are not detectable using circular polarisation only. An investigation into theoretical linear polarisation profiles shows that they could be used to reveal antisymmetric components of the magnetic field. This result also has implications for the study of exoplanetary transit hosts.

Both classical and relativistic weak-field and slow-motion perturbations to planetary orbits can be treated as perturbative corrections to the Keplerian model. In particular, tidal forces and General Relativity (GR) induce small precession rates of the apsidal line. Accurate measurements of these effects in transiting exoplanets could be used to test GR and to gain information about the planetary interiors. Unfortunately, models for transiting planets have a high degree of degeneracy in the orbital parameters that, combined to the uncertainties of photometric transit observations, results in large errors on the determinations of the argument of periastron and precludes a direct evaluation of the apsidal line precession. Moreover, tidal and GR precession time-scales are many order of magnitudes larger than orbital periods, so that on the observational time-spans required to cumulate a precession signal enough strong to be detected, even small systematic errors in transit ephemerides add up to cancel out the tiny variations due to precession. Here we present a more feasible solution to detect tidal and GR precession rates through the observation of variations of the time interval ($\Delta \tau$) between primary and secondary transits of hot Jupiters and propose the most promising target for such detection, WASP-14 b. For this planet we expect a cumulated $\Delta \tau$ $\approx$ -250 s, due to tidal and relativistic precession, since its first photometric observations.

Milky Way's hot gaseous halo extends up to the Galactic virial radius ($\sim 200$ kpc) and contains a significant component of baryon mass of the Galaxy. The halo properties can be constrained from X-ray spectroscopic observations and from satellite galaxies' ram-pressure stripping studies. Results of the former method crucially depend on the gas metallicity assumptions while the latter one's are insensitive to them. Here, a joint analysis of both kinds of data is presented to constrain electron density and metallicity of the gas. The power law is assumed for the electron density radial profile, while for the metallicity, a common-used constant-metallicity assumption is relaxed by introducing of a physically motivated spherical profile. The model is fitted to a sample of 433 (20) sight lines for O VII emission (absorption) measurements and 7 electron density constraints from ram-pressure stripping studies. The best-fitting electron density profile of $n_e\propto r^{-(0.9...1.1)}$ (where $r\gg1$ kpc is the galactocentric radius) is found. The metallicity is constrained as $Z \approx (0.1...0.7) Z_{\odot}$ (subscript $\odot$ represents the solar values) at $r>50$ kpc. These imply a total hot gas mass of $M \approx (2.4...8.7) \times 10^{10} M_{\odot}$, which accounts for $\sim(17...100)$ per cent of the Milky Way's missing baryon mass. The model uncertainties are discussed, and the results are examined in the context of previous studies.

In this work, we reanalyzed the 11 years of spectral data from Fermi-LAT of starburst galaxies (SBGs) and star-forming galaxies (SFGs) currently observed. We used a one-zone model provided by \textbf{NAIMA} and the hadronic origin to explain the GeV observation data of SBGs and SFGs. We found that a protonic distribution of a power-law form with an exponential cutoff can explain the spectra of most SBGs and SFGs. However, it cannot well explain the spectral hardening components from NGC 1068 and NGC 4945 in the GeV energy band. Therefore, we considered the two-zone model to well explain these phenomena. We summarized the features of various model parameters including the spectral index, cutoff energy, protonic total energy. Similar to the evolution of supernova remnants (SNRs) in the Milky Way, we estimated that the protonic acceleration limitation inside SBGs likely reach the order of 10$^{2}$ TeV by the one-zone model, close to those of SNRs in the Milky Way.

We perform high-resolution cosmological hydrodynamic simulations to study the formation of first galaxies that reach the masses of $10^{8-9}~h^{-1}~M_\odot$ at $z=9$. The resolution of the simulations is high enough to resolve minihaloes and allow us to successfully pursue the formation of multiple Population (Pop) III stars, their supernova (SN) explosions, resultant metal-enrichment of the inter-galactic medium (IGM) in the course of the build-up of the system. Metals are ejected into the IGM by multiple Pop III SNe, but some of the metal-enriched gas falls back onto the halo after $\gtrsim 100~\rm Myr$. The star formation history of the first galaxy depends sensitively on the initial mass function (IMF) of Pop III stars. The dominant stellar population transits from Pop III to Pop II at $z\sim 12-15$ in the case of power-law Pop III IMF, ${\rm d}n/{\rm d}M \propto M^{-2.35}$ with the mass range $10-500~M_\odot$. At $z\lesssim 12$, stars are stably formed in the first galaxies with a star formation rate of $\sim 10^{-3}$-$10^{-1}~M_\odot/{\rm yr}$. In contrast, for the case with a flat IMF, gas-deprived first galaxies form due to frequent Pop III pair-instability SNe, resulting in the suppression of subsequent Pop II star formation. In addition, we calculate UV continuum, Ly$\alpha$- and H$\alpha$-line fluxes from the first galaxies. We show that the James Webb Space Telescope will be able to detect both UV continuum, Ly$\alpha$ and H$\alpha$ line emission from first galaxies with halo mass $\gtrsim 10^{9}~M_\odot$ at $z \gtrsim 10$.

In order for off-Earth top surface structures built from regolith to protect astronauts from radiation, they need to be several meters thick. Technical University Delft (TUD) proposes to excavate into the ground to create subsurface habitats. By excavating not only natural protection from radiation can be achieved but also thermal insulation because the temperature is more stable underground. At the same time through excavation valuable resources can be mined for through in situ resource utilization (ISRU). The idea is that a swarm of autonomous mobile robots excavate the ground in a sloped downwards spiral movement. The excavated regolith will be mixed with cement, which can be reproduced on Mars through ISRU, in order to create concrete. The concrete is 3D printed/sprayed on the excavated tunnel to reinforce it. As soon as the tunnels are reinforced, the material in-between the tunnels can be removed in order to create a larger cavity that can be used for inhabitation. Proposed approach relies on Design-to-Robotic-Production (D2RP) technology developed at TUD1 for on-Earth applications. The rhizomatic 3D printed structure is a structurally optimized porous shell structure with increased insulation properties. In order to regulate the indoor pressurised environment an inflatable structure is placed in the 3D printed cavity. This inflatable structure is made of materials, which can also be at some point reproduced on Mars through ISRU. Depending on location the habitat and the production system are powered by a system combining solar and kite power. The ultimate goal is to develop an autarkic D2RP system for building subsurface autarkic habitats on Mars from locally obtained materials.

Aims: We present a detailed X-ray study of the recently discovered supernova remnant (SNR) G53.41+0.03 that follows up and further expands on the previous, limited analysis of archival data covering a small portion of the SNR. Methods: With the new dedicated 70ks XMM-Newton observation we investigate the morphological structure of the SNR in X-rays, search for a presence of a young neutron star and characterise the plasma conditions in the selected regions by means of spectral fitting. Results: The first full view of SNR G53.41+0.03 shows an X-ray emission region well aligned with the reported half-shell radio morphology. We find three distinct regions of the remnant that vary in brightness and hardness of the spectra, and are all best characterised by a hot plasma model in a non-equilibrium ionisation state. Of the three regions, the brightest one contains the most mature plasma, with ionisation age $\tau \approx 4\times10^{10}$s cm$^{-3}$ (where $\tau = n_e t$), a lower electron temperature of kT$_\mathrm{e} \approx 1$ keV and the highest estimated gas density of n$_\mathrm{H}\approx 0.87$ cm$^{-3}$. The second, fainter but spectrally harder, region reveals a younger plasma ($\tau \approx 1.7\times10^{10}$s cm$^{-3}$) with higher temperature (kT$_\mathrm{e} \approx 2$ keV) and two to three times lower density (n$_\mathrm{H}\approx 0.34$ cm$^{-3}$). The third region is very faint, but we identify spectroscopically the presence of a hot plasma.Employing several methods for age estimation, we find the remnant to be $t \approx 1000-5000$ yrs old, confirming the earlier reports of a relatively young age. The environment of the remnant also contains a number of point sources, of which most are expected to be positioned in the foreground. Of the two point sources in the geometrical centre of the remnant one is consistent with the characteristics of a young neutron star.

This paper is the second of a series that tackles the properties of molecular gas in galaxies residing in clusters and their related large-scale structures. We report on the observation of 22 galaxies in the CO(3-2) line gathered with the Atacama Large Millimeter Array (ALMA). These galaxies are either bona fide members of the CL1301.7$-$1139 cluster ($z=0.4828$, $\sigma_{cl}=681$ km s$^{-1}$), or located within $\sim 7 \times R_{200}$, its virial radius. They have been selected to sample the range of photometric local densities around CL1301.7$-$1139, with stellar masses above log($M_{\rm star}$) = 10, and to be located in the blue clump of star forming galaxies derived from the $u$, $g$, and $i$ photometric bands. Unlike previous works, our sample selection does not impose a minimum star formation rate or detection in the far infrared. As such and as much as possible, it delivers an unbiased view of the gas content of normal star forming galaxies at $z \sim 0.5$. Our study highlights the variety of paths to star formation quenching, and most likely the variety of physical properties (i.e. temperature, density) of the corresponding galaxy cold molecular gas. Just as in the case of CL1411.1$-$1148, although to a smaller extent, we identify a number of galaxies with lower gas fraction than classically found in other surveys. These galaxies can still be on the star forming main sequence. When these galaxies are not inside the cluster virialized region, we provide hints that they are linked to their infall regions within $\sim 4 \times R_{200}$.

Recently, the observed equation of state for dark energy appears to favor values below $-1$. The tendency implies that the nature of dark energy may be quite different from that of the cosmological constant. In view of the adjustment on the equation of state keeps decreasing, the introduction of the phantom energy seems inevitable. By employing observational constraints from supernovae and from the acoustic scale in which the accuracy of the data has become extraordinary, we apply a phenomenological scenario to be acquainted with the evolution of our universe. The demonstration on the constrained unfolding of the phantom energy shows the model has high consistency with the current observation.

We study the spectral diversity of Type Ia supernovae (SNe Ia) at maximum light using high signal-to-noise spectrophotometry of 173 SNe Ia from the Nearby Supernova Factory. We decompose the diversity of these spectra into different extrinsic and intrinsic components, and we construct a nonlinear parameterization of the intrinsic diversity of SNe Ia that preserves pairings of "twin" SNe Ia. We call this parameterization the "Twins Embedding". Our methodology naturally handles highly nonlinear variability in spectra, such as changes in the photosphere expansion velocity, and uses the full spectrum rather than being limited to specific spectral line strengths, ratios or velocities. We find that the time evolution of SNe Ia near maximum light is remarkably similar, with 84.6% of the variance in common to all SNe Ia. After correcting for brightness and color, the intrinsic variability of SNe Ia is mostly restricted to specific spectral lines, and we find intrinsic dispersions as low as ~0.02 mag between 6600 and 7200 A. With a nonlinear three-dimensional model plus one dimension for color, we can explain 89.2% of the intrinsic diversity in our sample of SNe Ia, which includes several different kinds of "peculiar" SNe Ia. A linear model requires seven dimensions to explain a comparable fraction of the intrinsic diversity. We show how a wide range of previously-established indicators of diversity in SNe Ia can be recovered from the Twins Embedding. In a companion article, we discuss how these results an be applied to standardization of SNe Ia for cosmology.

The temperature in most parts of a protoplanetary disk is determined by irradiation from the central star. Numerical experiments of Watanabe & Lin (2008) suggested that such disks, also called `passive disks', suffer from a thermal instability. Here, we use analytical and numerical tools to reveal the nature of this instability. We find that it is related to the flaring of the optical surface, the layer at which starlight is intercepted by the disk. Whenever a disk annulus is perturbed thermally and acquires a larger scale height, disk flaring becomes steeper in the inner part, and flatter in the outer part. Starlight now shines more overhead for the inner part and can penetrate into deeper layers; conversely, it is absorbed more shallowly in the outer part. These geometric changes allow the annulus to intercept more starlight, and the perturbation grows. We call this the irradiation instability. It requires only ingredients known to exist in realistic disks (stellar illumination, optically thick), and operates best in parts that are very optically thick (inside 20 AU, but can extend to further reaches when, e.g., dust settling is considered). An unstable disk develops travelling thermal waves that reach order-unity in amplitude. In thermal radiation, such a disk should appear as a series of bright rings interleaved with dark shadowed gaps, while in scattered light it resembles a moving staircase. Depending on the gas and dust responses, this instability could lead to a wide range of interesting consequences, such as dust traps, vertical circulation, vortices and turbulence.

The space and ground-based observations have shown a lot of activities and instabilities in the atmosphere of the giant ice planet Neptune. Using the archival data of high resolution Atacama Large Millimeter/Submillimeter Array (ALMA) with band 7 observation, we present the spectroscopic detection of the rotational emission line of sulfur dioxide (SO$_{2}$) at frequency $\nu$ = 343.476 GHz with transition J=57$_{15,43}$$-$58$_{14,44}$. We also re-detect the emission line of carbon monoxide (CO) at frequency $\nu$ = 345.795 GHz with transition J=3$-$2. The molecular lines of SO$_{2}$ and CO in the atmosphere of Nepure are detected with the $\geq$4$\sigma$ statistical significance. The statistical column density of SO$_{2}$ is N(SO$_{2}$) = 2.61$\times$10$^{15}$ cm$^{-2}$ with rotational temperature $T_{SO_{2}}$ = 50 K and the statistical column density CO is N(CO) = 1.86$\times$10$^{19}$ cm$^{-2}$ with $T_{CO}$ = 29 K. The typical mixing ratio in the atmosphere of Neptune for SO$_{2}$ is 1.24$\times$10$^{-10}$ and CO is 0.88$\times$10$^{-6}$. The SO$_{2}$ and CO gas in the atmosphere of Neptune may create due to Shoemaker-Levy 9 impacts in Jovian planets since 1994.

We analyzed the TeV gamma-ray image of a supernova remnant RX J1713.7$-$3946 (RX J1713) through a comparison with the interstellar medium (ISM) and the non-thermal X-rays. The gamma-ray datasets at two energy bands of $>$2 TeV and $>$250-300 GeV were obtained with H.E.S.S. (H.E.S.S. Collaboration 2018; Aharonian et al. 2007) and utilized in the analysis. We employed a new methodology which assumes that the gamma-ray counts are expressed by a linear combination of two terms; one is proportional to the ISM column density and the other proportional to the X-ray count. We then assume these represent the hadronic and leptonic components, respectively. By fitting the expression to the data pixels, we find that the gamma-ray counts are well represented by a flat plane in a 3D space of the gamma-ray counts, the ISM column density and the X-ray counts. The results using the latest H.E.S.S. data at 4.8 arcmin resolution show that the hadronic and leptonic components occupy $(67\pm8)$% and $(33\pm8)$% of the total gamma rays, respectively, where the two components have been quantified for the first time. The hadronic component is greater than the leptonic component, which reflects the massive ISM of $\sim$10$^4$ $M_{\odot}$ associated with the SNR, lending support for the acceleration of the cosmic-ray protons. There is a marginal hint that the gamma rays are suppressed at high gamma-ray counts which may be ascribed to the second order effects including the shock-cloud interaction and the penetration-depth effect.

We present the goals, strategy and first results of the high-cadence Galactic plane survey using the Zwicky Transient Facility (ZTF). The goal of the survey is to unveil the Galactic population of short-period variable stars, including short period binaries and stellar pulsators with periods less than a few hours. Between June 2018 and January 2019, we observed 64 ZTF fields resulting in 2990 deg$^2$ of high stellar density in ZTF-$r$ band along the Galactic Plane. Each field was observed continuously for 1.5 to 6 hrs with a cadence of 40 sec. Most fields have between 200 and 400 observations obtained over 2-3 continuous nights. As part of this survey we extract a total of $\approx$230 million individual objects with at least 80 epochs obtained during the high-cadence Galactic Plane survey reaching an average depth of ZTF-$r$ $\approx$20.5 mag. For four selected fields with 2 million to 10 million individual objects per field we calculate different variability statistics and find that $\approx$1-2% of the objects are astrophysically variable over the observed period. We present a progress report on recent discoveries, including a new class of compact pulsators, the first members of a new class of Roche Lobe filling hot subdwarf binaries as well as new ultracompact double white dwarfs and flaring stars. Finally we present a sample of 12 new single-mode hot subdwarf B-star pulsators with pulsation amplitudes between ZTF-$r$ = 20-76 mmag and pulsation periods between $P$ = 5.8-16 min with a strong cluster of systems with periods $\approx$ 6 min. All of the data have now been released in either ZTF Data Release 3 or data release 4.

The structure of the icy shells of ocean worlds is important for understanding the stability of their underlying oceans as it controls the rate at which heat can be transported outward and radiated to space. Future spacecraft exploration of the ocean worlds (e.g., by NASA's Europa Clipper mission) will allow for higher-resolution measurements of gravity and shape than currently available. In this paper, we study the sensitivity of gravity-topography admittance to the structure of icy shells in preparation for future data analysis. An analytical viscous relaxation model is used to predict admittance spectra given different shell structures determined by the temperature-dependent viscosity of a tidally heated, conductive shell. We apply these methods to the ocean worlds of Europa and Enceladus. We find that admittance is sensitive to the mechanisms of topography support at different wavelengths and estimate the required gravity performance to resolve transitions between these mechanisms. We find that the Airy isostatic model is unable to accurately describe admittance except at the longest wavelengths due to a rapid relaxation of the shell's basal topography resulting from a large viscosity contrast. Our models suggest that measurements of admittance at low spherical harmonic degrees are more sensitive to thick shells with high tidal dissipation, and may complement ice-penetrating radar measurements in constraining shell thickness. Finally, we find that admittance may be used to constrain the tidal dissipation within the icy shell, which would be complementary to a more demanding measurement of the tidal phase lag.

We have simulated the evolution of non-thermal cosmic ray electrons (CREs) in 3D relativistic magneto hydrodynamic (MHD) jets evolved up to a height of 9 kpc. The CREs have been evolved in space and in energy concurrently with the relativistic jet fluid, duly accounting for radiative losses and acceleration at shocks. We show that jets stable to MHD instabilities show expected trends of regular flow of CREs in the jet spine and acceleration at a hotspot followed by a settling backflow. However, unstable jets create complex shock structures at the jet-head (kink instability), the jet spine-cocoon interface and the cocoon itself (Kelvin-Helmholtz modes). CREs after exiting jet-head undergo further shock crossings in such scenarios and are re-accelerated in the cocoon. CREs with different trajectories in turbulent cocoons have different evolutionary history with different spectral parameters. Thus at the same spatial location, there is mixing of different CRE populations, resulting in a complex total CRE spectrum when averaged over a given area. Cocoons of unstable jets can have an excess build up of energetic electrons due to re-acceleration at turbulence driven shocks and slowed expansion of the decelerated jet. This will add to the non-thermal energy budget of the cocoon.

We use field-level forward models of galaxy clustering and the EFT likelihood formalism to study, for the first time for self-consistently simulated galaxies, the relations between the linear $b_1$ and second-order bias parameters $b_2$ and $b_{K^2}$. The forward models utilize all of the information available in the galaxy distribution up to a given order in perturbation theory, which allows us to infer these bias parameters with high signal-to-noise, even from relatively small volumes ($L_{\rm box} = 205{\rm Mpc}/h$). We consider galaxies from the IllustrisTNG simulations, and our main result is that the $b_2(b_1)$ and $b_{K^2}(b_1)$ relations obtained from gravity-only simulations for total mass selected objects are broadly preserved for simulated galaxies selected by stellar mass, star formation rate, color and black hole accretion rate. This consistency under different galaxy selection criteria suggests that theoretical priors on these bias relations may be used to improve cosmological constraints based on observed galaxy samples. We do identify some small differences between the bias relations in the hydrodynamical and gravity-only simulations, which we show can be linked to the environmental dependence of the relation between galaxy properties and mass. We also show that the EFT likelihood recovers the value of $\sigma_8$ to percent-level from various galaxy samples (including splits by color and star formation rate) and after marginalizing over 8 bias parameters. This demonstration using simulated galaxies adds to previous works based on halos as tracers, and strengthens further the potential of forward models to infer cosmology from galaxy data.

The space-borne gravitational wave (GW) detectors, LISA and TAIJI, are planned to be launched in the 2030s. The dual detectors with comparable sensitivities will form a network observing GW with significant advantages. In this work, we report the investigations on the three possible LISA-TAIJI networks for different location and orientation compositions of LISA orbit ($60^\circ$ inclination and tailing the Earth by $20^\circ$) and alternative TAIJI orbit configurations including TAIJIp ($60^\circ$ inclination and leading the Earth by $20^\circ$), TAIJIc ($60^\circ$ inclination and co-located with LISA), TAIJIm ($-60^\circ$ inclination and leading the Earth by $20^\circ$). In the three LISA-TAIJI configurations, the LISA-TAIJIm network shows a best performance on the sky localization and polarization determination for the massive binary system due to their complementary antenna pattern, and LISA-TAIJIc could achieve the best cross correlation and observe the stochastic GW background with an optimal sensitivity.

The idea of possible modification to gravity theory, whether it is in the Newtonian or general relativistic premises, is there for quite sometime. Based on it, astrophysical and cosmological problems are targeted to solve. But none of the Newtonian theories of modification has been performed from the first principle. Here, we modify Poisson's equation and propose two possible ways to modify the law gravitation which, however, reduces to Newton's law far away from the center of source. Based on these modified Newton's laws, we attempt to solve problems lying with white dwarfs. There are observational evidences for possible violation of the Chandrasekhar mass-limit significantly: it could be sub- as well as super-Chandrasekhar. We show that modified Newton's law, either by modifying LHS or RHS of Poisson's equation, can explain them.

We develop a formalism to describe the scattering of dark matter (DM) particles by electrons bound in crystals for a general form of the underlying DM-electron interaction. Such a description is relevant for direct-detection experiments of DM particles lighter than a nucleon, which might be observed in operating DM experiments via electron excitations in semiconductor crystal detectors. Our formalism is based on an effective theory approach to general non-relativistic DM-electron interactions, including the anapole, and magnetic and electric dipole couplings, combined with crystal response functions defined in terms of electron wave function overlap integrals. Our main finding is that, for the usual simplification of the velocity integral, the rate of DM-induced electronic transitions in a semiconductor material depends on at most five independent crystal response functions, four of which were not known previously. We identify these crystal responses, and evaluate them using density functional theory for crystalline silicon and germanium, which are used in operating DM direct detection experiments. Our calculations allow us to set 90% confidence level limits on the strength of DM-electron interactions from data reported by the SENSEI and EDELWEISS experiments. The novel crystal response functions discovered in this work encode properties of crystalline solids that do not interact with conventional experimental probes, suggesting the use of the DM wind as a probe to reveal new kinds of hidden order in materials.

Electron-tracking Compton camera, which is a complete Compton camera with tracking Compton scattering electron by a gas micro time projection chamber, is expected to open up MeV gamma-ray astronomy. The technical challenge for achieving several degrees of the point spread function is the precise determination of the electron-recoil direction and the scattering position from track images. We attempted to reconstruct these parameters using convolutional neural networks. Two network models were designed to predict the recoil direction and the scattering position. These models marked 41$~$degrees of the angular resolution and 2.1$~$mm of the position resolution for 75$~$keV electron simulation data in Argon-based gas at 2$~$atm pressure. In addition, the point spread function of ETCC was improved to 15$~$degrees from 22$~$degrees for experimental data of 662$~$keV gamma-ray source. These performances greatly surpassed that using the traditional analysis.

A new cosmological scenario is proposed in which a light scalaron of $f (R)$ gravity plays the role of dark matter. In this scenario, the scalaron initially resides at the minimum of its effective potential while the electroweak symmetry is unbroken. At the beginning of the electroweak crossover, the evolving expectation value of the Higgs field triggers the evolution of the scalaron due to interaction between these fields. After the electroweak crossover, the oscillating scalaron can represent cold dark matter. Its current energy density depends on a single free parameter, the scalaron mass $m$, and the value $m \simeq 4.4 \times 10^{-3}\, \text{eV}$ is required to explain the observed dark-matter abundance. Larger mass values would be required in scenarios where the scalaron is excited before the electroweak crossover.

We assume the anisotropic model of the Universe in the framework of varying speed of light $c$ and varying gravitational constant $G$ theories and study different types of singularities. For the singularity models, we write the scale factors in terms of cosmic time and found some conditions for possible singularities. For future singularities, we assume the forms of varying speed of light and varying gravitational constant. For regularizing big bang singularity, we assume two forms of scale factors: sine model and tangent model. For both the models, we examine the validity of null energy condition and strong energy condition. Start from the first law of thermodynamics, we study the thermodynamic behaviours of $n$ number of Universes (i.e., Multiverse) for (i) varying $c$, (ii) varying $G$ and (iii) both varying $c$ and $G$ models. We found the total entropies for all the cases in the anisotropic Multiverse model. We also found the nature of the Multiverse if total entropy is constant.

We consider Einstein-Gauss-Bonnet (EGB) inflationary models using the effective potential approach. We present evolution equations in the slow-roll regime using the effective potential and the tensor-to-scalar ratio. The choice of the effective potential is related to an expression of the spectral index in terms of e-folding number $N_e$. The satisfaction of the slow-roll regime is mostly related to the form of the tensor-to-scalar ratio $r$. The case of $r\sim1/N^2_e$ leads to a generalization of $\alpha$-attractors inflationary parameters to Einstein-Gauss-Bonnet gravity with exponential effective potential. Also the cosmological attractors include models with $r\sim1/N_e$. We are checking the satisfaction of the slow-roll regime during inflation for models with $r\sim1/N_e$.

Observations of Interplanetary Scintillation (IPS) are an efficient remote-sensing method to study the solar wind and inner heliosphere. From 2016 to 2018, some distinctive observations of IPS sources like 3C 286 and 3C 279 were accomplished with the Five-hundred-meter Aperture Spherical radio Telescope (FAST), the largest single-dish telescope in the world. Due to the 270-1620 MHz wide frequency coverage of the Ultra-Wideband (UWB) receiver, one can use both single-frequency and dual-frequency analyses to determine the projected velocity of the solar wind. Moreover, based on the extraordinary sensitivity owing to the large collecting surface area of FAST, we can observe weak IPS signals. With the advantages of both the wider frequency coverage and high sensitivity, also with our radio frequency interference (RFI) mitigation strategy and an optimized model-fitting method developed, in this paper, we analyze the fitting confidence intervals of the solar wind velocity, and present some preliminary results achieved using FAST, which points to the current FAST system being highly capable of carrying out observations of IPS

Dynamical systems methods are used to investigate a cosmological model with non-minimally coupled scalar field and asymptotically quadratic potential function. We found that for values of the non-minimal coupling constant parameter $\frac{3}{16}<\xi<\frac{1}{4}$ there exists an unstable asymptotic de Sitter state giving rise to non-singular beginning of the universe. The energy density associated with this state depends on value of the non-minimal coupling constant and can be much smaller than the Planck energy density. For $\xi=\frac{1}{4}$ we found that the initial state is in form of the static Einstein universe. Proposed evolutional model, contrary to the seminal Starobinsky's model, do not depend on the specific choice of initial conditions in phase space, moreover, a small change in the model parameters do not change the evolution thus the model is generic and structurally stable. The values of the non-minimal coupling constant can indicate for a new fundamental symmetry in the gravitational theory. We show that Jordan frame and Einstein frame formulation of the theory are physically nonequivalent.

First-order phase transitions can leave relic pockets of false vacua and their particles, that manifest as macroscopic Dark Matter. We compute one predictive model: a gauge theory with a dark quark relic heavier than the confinement scale. During the first-order phase transition to confinement, dark quarks remain in the false vacuum and get compressed, forming Fermi balls that can undergo gravitational collapse to stable dark dwarfs (bound states analogous to white dwarfs) near the Chandrasekhar limit, or to primordial black holes.