Papers for Thursday, Oct 14 2021

We review the current status of general relativistic studies for coalescences of black hole--neutron star binaries. First, high-precision computations of black hole--neutron star binaries in quasiequilibrium circular orbits are summarized, focusing on the quasiequilibrium sequences and the mass-shedding limit. Next, the current status of numerical-relativity simulations for the merger of black hole--neutron star binaries is described. We summarize our understanding for the merger process, tidal disruption and its criterion, properties of the merger remnant and ejected material, gravitational waveforms, and gravitational-wave spectra. We also discuss expected electromagnetic counterparts to black hole--neutron star coalescences.

Multi-dimensional (magneto-)hydrodynamical simulations of physical processes in stellar interiors depend on a multitude of uncalibrated free parameters, which set the spatial and time scales of their computations. We aim to provide an asteroseismic calibration of the wave and convective Rossby numbers, and of the stiffness at the interface between the convective core and radiative envelope of intermediate-mass stars. We deduce these quantities for rotating dwarfs from the observed properties of their identified gravity and gravito-inertial modes. We rely on near-core rotation rates and asteroseismic models of 26 B- and 37 F-type dwarf pulsators derived from 4-year Kepler space photometry, high-resolution spectroscopy and Gaia astrometry in the literature to deduce their convective and wave Rossby numbers. We compute the stiffness at the convection/radiation interface from the inferred maximum buoyancy frequency at the interface and the convective turnover frequency in the core. We use those asteroseismically inferred quantities to make predictions of convective penetration levels, local flux levels of gravito-inertial waves triggered by the convective core, and of the cores' potential rotational and magnetic states. Our sample of 63 gravito-inertial mode pulsators covers near-core rotation rates from almost zero up to the critical rate. The frequencies of their identified modes lead to models with stiffness values between $10^{2.69}$ and $10^{3.60}$ for the B-type pulsators, while those of F-type stars cover the range from $10^{3.47}$ to $10^{4.52}$. The convective Rossby numbers derived from the maximum convective diffusion coefficient in the convective core, based on mixing length theory and a value of the mixing length coefficient relevant for these pulsators, vary between $10^{-2.3}$ and $10^{-0.8}$ for B-type stars and $10^{-3}$ and $10^{-1.5}$ for F-type stars. (abridged)

Optically luminous quasars at $z > 5$ are important probes of super-massive black hole (SMBH) formation. With new and future radio facilities, the discovery of the brightest low-frequency radio sources in this epoch would be an important new probe of cosmic reionization through 21-cm absorption experiments. In this work, we systematically study the low-frequency radio properties of a sample of 115 known spectroscopically confirmed $z>5$ quasars using the second data release of the Low Frequency Array (LOFAR) Two Metre Sky survey (LoTSS-DR2), reaching noise levels of $\sim$80 $\mu$Jy beam$^{-1}$ (at 144 MHz) over an area of $\sim5720$ deg$^2$. We find that 41 sources (36%) are detected in LoTSS-DR2 at $>2 \sigma$ significance and we explore the evolution of their radio properties (power, spectral index, and radio loudness) as a function of redshift and rest-frame ultra-violet properties. We obtain a median spectral index of $-0.29^{+0.10}_{-0.09}$ by stacking 93 quasars using LoTSS-DR2 and Faint Images of the Radio Sky at Twenty Centimetres (FIRST) data at 1.4 GHz, in line with observations of quasars at $z<3$. We compare the radio loudness of the high-$z$ quasar sample to a lower-$z$ quasar sample at $z\sim2$ and find that the two radio loudness distributions are consistent with no evolution, although the low number of high-z quasars means that we cannot rule out weak evolution. Furthermore, we make a first order empirical estimate of the $z=6$ quasar radio luminosity function, which is used to derive the expected number of high-$z$ sources that will be detected in the completed LoTSS survey. This work highlights the fact that new deep radio observations can be a valuable tool in selecting high-$z$ quasar candidates for follow-up spectroscopic observations by decreasing contamination of stellar dwarfs and reducing possible selection biases introduced by strict colour cuts.

We examine Suzaku, XMM-Newton, and Chandra observations of the Abell 2029/2033 system to investigate the nature of a bridge of X-ray emission joining the two galaxy clusters. By modelling the contributions from the outskirts of the two clusters, and excluding the emission from the southern infalling group and the background group LOS9, we find a significant excess of X-ray emission between the two clusters at the level of 6.5-7.0$\sigma$, depending on the choice of model, that cannot be explained by the overlap of the clusters. This excess component to the surface brightness is consistent with being emission from a filament with roughly 1.0 Mpc wide. The derived emission measure for the gas associated with the filament yields an average gas density of $3.7^{+1.0}_{-0.7} \times 10^{-5}$ cm$^{-3}$, corresponding roughly to 160 times the mean baryon density of the Universe. The Suzaku X-ray spectrum of the excess emission indicates that it is significantly colder ($1.4_{-0.5}^{+0.7}$ keV) than the cluster outskirts emission from the two clusters ($\sim$ 5 keV), statistically consistent with the temperature expected from the hottest and densest parts of the warm-hot intergalactic medium (WHIM). The geometry, density, and temperature are similar to those found from X-ray studies of the Abell 222/223 filament.

We present the deepest Spitzer/IRAC $3.6$, $4.5$, $5.8$ and $8.0\mu$m wide-area mosaics yet over the GOODS-N and GOODS-S fields as part of the GOODS Re-ionization Era wide-Area Treasury from Spitzer (GREATS) project. We reduced and mosaicked in a self-consistent way observations taken by the 11 different Spitzer/IRAC programs over the two GOODS fields from 12 years of Spitzer cryogenic and warm mission data. The cumulative depth in the $3.6\mu$m and $4.5\mu$m bands amounts to $\sim 4260$ hr, $\sim 1220$ hr of which are new very deep observations from the GREATS program itself. In the deepest area, the full-depth mosaics reach $\gtrsim200$ hr over an area of $\sim100$ arcmin$^2$, corresponding to a sensitivity of $\sim29$ AB magnitude at $3.6\mu$m ($1\sigma$ for point sources). Archival cryogenic $5.8\mu$m and $8.0\mu$m band data (a cumulative 976 hr) are also included in the release. The mosaics are projected onto the tangential plane of CANDELS/GOODS at a $0.3''$ pixel$^{-1}$ scale. This paper describes the methodology enabling, and the characteristics of, the public release of the mosaic science images, the corresponding coverage maps in the four IRAC bands, and the empirical Point-Spread Functions (PSFs). These PSFs enable mitigation of the source blending effects by taking into account the complex position-dependent variation in the IRAC images. The GREATS data products are in the Infrared Science Archive (IRSA). We also release the deblended $3.6$-to-$8.0\mu$m photometry for $9192$ Lyman-Break galaxies at $z\sim3.5-10$. GREATS will be the deepest mid-infrared imaging until JWST and, as such, constitutes a major resource for characterizing early galaxy assembly.

We perform 3D SPH simulations of warped, non-coplanar gravitationally unstable discs to show that as the warp propagates through the self-gravitating disc, it heats up the disc rendering it gravitationally stable. Thus losing their spiral structure and appearing completely axisymmetric. In their youth, protoplanetary discs are expected to be massive and self-gravitating, which results in non-axisymmetric spiral structures. However recent observations of young protoplanetary discs with ALMA have revealed that discs with large-scale spiral structure are rarely observed in the midplane. Instead, axisymmetric discs with some also having ring & gap structures are more commonly observed. Our work invloving warps, non-coplanar disc structures that are expected to commonly occur in young discs, potentially resolves this discrepancy between observations and theoretical predictions. We demonstrate that they are able to suppress the large-scale spiral structure of self-gravitating protoplanetary discs.

We investigate the origin of the measured over-abundance of alkali metals in the atmospheres of hot gas giants, relative to both their host stars and their atmospheric water abundances. We show that formation exterior to the water snow line followed by inward disc-driven migration results in excess accretion of oxygen-poor, refractory-rich material from within the snow-line. This naturally leads to enrichment of alkali metals in the planetary atmosphere relative to the bulk composition of its host star but relative abundances of water that are similar to the stellar host. These relative abundances cannot be explained by in situ formation which places the refractory elements in the planetary deep interior rather than the atmosphere. We therefore suggest that the measured compositions of the atmospheres of hot Jupiters are consistent with significant migration for at least a subset of hot gas giants. Our model makes robust predictions about atmospheric composition that can be confirmed with future data from JWST and Ariel.

Cosmic rays (CRs) are an important component in the interstellar medium (ISM), but their effect on the dynamics of the disk-halo interface (< 10 kpc from the disk) is still unclear. We study the influence of CRs on the gas above the disk with high-resolution FIRE-2 cosmological simulations of late-type Lstar galaxies at redshift around zero. We compare runs with and without CR feedback (with constant anisotropic diffusion around 3e29 cm^2/s and streaming). Our simulations capture the relevant disk halo interactions, including outflows, inflows, and galactic fountains. Extra-planar gas in all of the runs satisfies dynamical balance, where total pressure balances the weight of the overlying gas. While the kinetic pressure from non-uniform motion (>1-kpc scale) dominates in the midplane, thermal and bulk pressures (or CR pressure if included) take over at large heights. We find that with CR feedback, (1) the warm (1e4K) gas is slowly accelerated by CRs; (2) the hot (> 5e5K) gas scale height is suppressed; (3) the warm-hot (2e4-5e5K) medium becomes the most volume-filling phase in the disk-halo interface. We develop a novel conceptual model of the near-disk gas dynamics in low-redshift Lstar galaxies: With CRs, the disk-halo interface is filled with CR-driven warm winds and hot super-bubbles that are propagating into the CGM with a small fraction falling back to the disk. Without CRs, most outflows from hot superbubbles are trapped by the existing hot halo and gravity, so typically they form galactic fountains.

Cluster strong lensing cosmography is a promising probe of the background geometry of the Universe and several studies have emerged, thanks to the increased quality of observations using space and ground-based telescopes. For the first time, we use a sample of five cluster strong lenses to measure the values of cosmological parameters and combine them with those from classical probes. In order to assess the degeneracies and the effectiveness of strong-lensing cosmography in constraining the background geometry of the Universe, we adopt four cosmological scenarios. We find good constraining power on the total matter density of the Universe ($\Omega_{\rm m}$) and the equation of state of the dark energy parameter $w$. For a flat $w$CDM cosmology, we find $\Omega_{\rm m} = 0.30_{-0.11}^{+0.09}$ and $w=-1.12_{-0.32}^{+0.17}$ from strong lensing only. Interestingly, we show that the constraints from the Cosmic Microwave Background (CMB) are improved by factors of 2.5 and 4.0 on $\Omega_{\rm m}$ and $w$, respectively, when combined with our posterior distributions in this cosmological model. In a scenario where the equation of state of dark energy evolves with redshift, the strong lensing constraints are compatible with a cosmological constant (i.e. $w=-1$). In a curved cosmology, our strong lensing analyses can accommodate a large range of values for the curvature of the Universe of $\Omega_{\rm k}=0.28_{-0.21}^{+0.16}$. In all cosmological scenarios, we show that our strong lensing constraints are complementary and in good agreement with measurements from the CMB, baryon acoustic oscillations and Type Ia supernovae. Our results show that cluster strong lensing cosmography is a potentially powerful probe to be included in the cosmological analyses of future surveys.

Most isolated neutron stars have been discovered thanks to the detection of their pulsed non-thermal emission, at wavelengths spanning from radio to gamma-rays. However, if the beamed non-thermal radiation does not intercept our line of sight or it is too faint or absent, isolated neutron stars can also be detected through their thermal emission, which peaks in the soft X-ray band and is emitted nearly isotropically. In the past thirty years, several thermally-emitting isolated neutron stars have been discovered thanks to X-ray all-sky surveys, observations targeted at the center of supernova remnants, or as serendipitous X-ray sources. Distinctive properties of these relatively rare X-ray sources are very soft spectra and high ratios of X-ray to optical flux. The recently released 4XMM-DR10 catalog contains more than half a million X-ray sources detected with the XMM-Newton telescope in the 0.2-10 keV range in observations carried out from 2000 to 2019. Based on a study of the spectral properties of these sources and on cross-correlations with catalogs of possible counterparts, we have carried out a search of isolated neutron stars, finding four potential candidates. The spectral and long-term variability analysis of these candidates, using also Chandra and Swift-XRT data, allowed us to point out the most interesting sources deserving further multiwavelength investigations.

We present deep ALMA dust continuum observations for a sample of luminous ($M_{\rm UV} < -22$) star-forming galaxies at $z \simeq 7$. We detect five of the six sources in the far-infrared (FIR), providing key constraints on the obscured star-formation rate (SFR) and the infrared-excess-$\beta$ (IRX-$\beta$) relation without the need for stacking. Despite the galaxies showing blue rest-frame UV slopes ($\beta \simeq -2$) we find that 35-75 percent of the total SFR is obscured. We find the IRX-$\beta$ relation derived for these $z \simeq 7$ sources is consistent with that found for local star-burst galaxies. Using our relatively high-resolution (FWHM $\simeq 0.7\,{\rm arcsec}$) observations we identify a diversity of dust morphologies in the sample. We find both compact emission that appears offset relative to the unobscured components and extended dust emission that is co-spatial with the rest-frame UV light. In the majority of the sources we detect strong rest-frame UV colour gradients (with up to $\Delta \beta \simeq 0.7$-$1.4$) as probed by the multi-band UltraVISTA ground-based data. The observed redder colours are spatially correlated with the location of the FIR detection. Our results show that, even in bright Lyman-break galaxies at $z \simeq 7$, the majority of the star-formation is actually occurring within the faintest components in the rest-frame UV, which have an obscured fraction of $f_{\rm obs} \ge 0.8$. As well as demonstrating the importance of dust obscured star-formation within the Epoch of Reionization, these observations provide an exciting taster of the rich spatially resolved datasets that will be obtained from JWST and high-resolution ALMA follow-up at these redshifts.

We investigate the [X/Mg] abundances of 16 elements for 82,910 Galactic disk stars from GALAH+ DR3. We fit the median trends of low-Ia and high-Ia populations with a two-process model, which describes stellar abundances in terms of a prompt core-collapse and delayed Type-Ia supernova component. For each sample star, we fit the amplitudes of these two components and compute the residual $\Delta$[X/H] abundances from this two-parameter fit. We find RMS residuals $\lesssim 0.07$ dex for well-measured elements and correlated residuals among some elements (such as Ba, Y, and Zn) that indicate common enrichment sources. From a detailed investigation of stars with large residuals, we infer that roughly $40\%$ of the large deviations are physical and $60\%$ are caused by problematic data such as unflagged binarity, poor wavelength solutions, and poor telluric subtraction. As one example of a population with distinctive abundance patterns, we identify 15 stars that have 0.3-0.6 dex enhancements of Na but normal abundances of other elements from O to Ni and positive average residuals of Cu, Zn, Y, and Ba. We measure the median elemental residuals of 14 open clusters, finding systematic $\sim0.1-0.4$ dex enhancements of O, Ca, K, Y, and Ba and $\sim0.2$ dex depletion of Cu in young clusters. Finally, we present a restricted three-process model where we add an asymptotic giant branch star (AGB) component to better fit Ba and Y. With the addition of the third process, we identify a population of stars, preferentially young, that have much higher AGB enrichment than expected from their SNIa enrichment.

The streaming instability is a fundamental process that can drive dust-gas dynamics and ultimately planetesimal formation in protoplanetary discs. As a linear instability, it has been shown that its growth with a distribution of dust sizes can be classified into two distinct regimes, fast- and slow-growth, depending on the dust-size distribution and the total dust-to-gas density ratio $\epsilon$. Using numerical simulations of an unstratified disc, we bring three cases in different regimes into nonlinear saturation. We find that the saturation states of the two fast-growth cases are similar to its single-species counterparts. The one with maximum dimensionless stopping time $\tau_\mathrm{s,max}=0.1$ and $\epsilon=2$ drives turbulent vertical dust-gas vortices, while the other with $\tau_\mathrm{s,max}=2$ and $\epsilon=0.2$ leads to radial traffic jams and filamentary structures of dust particles. The dust density distribution for the former is flat in low densities, while the one for the latter has a low-end cutoff. By contrast, the one slow-growth case results in a virtually quiescent state. Moreover, we find that in the fast-growth regime, significant dust segregation by size occurs, with large particles moving towards dense regions while small particles remain in the diffuse regions, and the mean radial drift of each dust species is appreciably altered from the (initial) drag-force equilibrium. The former effect may skew the spectral index derived from multi-wavelength observations and change the initial size distribution of a pebble cloud for planetesimal formation. The latter along with turbulent diffusion may influence the radial transport and mixing of solid materials in young protoplanetary discs.

The early 1970s saw a new and surprising feature in the composition of solar energetic particles (SEPs), resonant enhancements up to 10,000-fold in the ratio 3He/4He that could even make 3He dominant over H in rare events. It was soon learned that these events also had enhancements in the abundances of heavier elements, such as a factor of ~10 enhancements in Fe/O, which was later seen to be part of a smooth increase in enhancements vs. mass-to-charge ratio A/Q from H to Pb, rising by a factor of ~1000. These events were also associated with streaming 10 - 100 keV electrons that produce type III radio bursts. In recent years we have found these "impulsive" SEP events to be accelerated in islands of magnetic reconnection from plasma temperatures of 2 - 3 MK on open field lines in solar jets. Similar reconnection on closed loops traps the energy of the particles to produce hot (>10 MK), bright flares. Sometimes impulsive SEP intensities are boosted by shock waves when the jets launch fast coronal mass ejections. No single theory yet explains both the sharp resonance in 3He and the smooth increase up to heavier elements; two processes seem to occur. Sometimes the efficient acceleration even exhausts the rare 3He in the source region, limiting its fluence.

Using the Atacama Large Millimetre/submillimeter Array (ALMA) and the Karl G. Jansky Very Large Array (VLA), we observed the Extended Green Object (EGO) G19.01$-$0.03 with sub-arcsecond resolution from 1.05 mm to 5.01 cm wavelengths. Our $\sim0.4''\sim1600$ AU angular resolution ALMA observations reveal a velocity gradient across the millimetre core MM1, oriented perpendicular to the previously known bipolar molecular outflow, that is consistently traced by 20 lines of 8 molecular species with a range of excitation temperatures, including complex organic molecules (COMs). Kinematic modelling shows the data are well described by models that include a disc in Keplerian rotation and infall, with an enclosed mass of $40-70 \mathrm{M}_{\odot}$ (within a 2000 AU outer radius) for a disc inclination angle of $i=40^{\circ}$, of which $5.4-7.2 \mathrm{M}_{\odot}$ is attributed to the disc. Our new VLA observations show that the 6.7 GHz Class II methanol masers associated with MM1 form a partial ellipse, consistent with an inclined ring, with a velocity gradient consistent with that of the thermal gas. The disc-to-star mass ratio suggests the disc is likely to be unstable and may be fragmenting into as-yet-undetected low mass stellar companions. Modelling the centimetre--millimetre spectral energy distribution of MM1 shows the ALMA 1.05 mm continuum emission is dominated by dust, whilst a free-free component, interpreted as a hypercompact HII region, is required to explain the VLA $\sim$5 cm emission. The high enclosed mass derived for a source with a moderate bolometric luminosity ($\sim$10$^{4} \mathrm{L}_{\odot}$) suggests that the MM1 disc may feed an unresolved high-mass binary system.

The cosmic large-scale structures of the Universe are mainly the result of the gravitational instability of initially small density fluctuations in the dark-matter distribution. Dark matter appears to be initially cold and behaves as a continuous and collisionless medium on cosmological scales, with evolution governed by the gravitational Vlasov--Poisson equations. Cold dark matter can accumulate very efficiently at focused locations, leading to a highly non-linear filamentary network with extreme matter densities. Traditionally, investigating the non-linear Vlasov--Poisson equations was typically reserved for massively parallelised numerical simulations. Recently, theoretical progress has allowed us to analyse the mathematical structure of the first infinite densities in the dark-matter distribution by elementary means. We review related advances, as well as provide intriguing connections to classical plasma problems, such as the beam-plasma instability.

The study of stellar surfaces can reveal information about the chemical composition, interior structure, and magnetic properties of stars. It is also critical to the detection and characterization of extrasolar planets, in particular those targeted in extreme precision radial velocity (EPRV) searches, which must contend with stellar variability that is often orders of magnitude stronger than the planetary signal. One of the most successful methods to map the surfaces of stars is Doppler imaging, in which the presence of inhomogeneities is inferred from subtle line shape changes in high resolution stellar spectra. In this paper, we present a novel, efficient, and closed-form solution to the problem of Doppler imaging of stellar surfaces. Our model explicitly allows for incomplete knowledge of the local (rest frame) stellar spectrum, allowing one to learn differences from spectral templates while simultaneously mapping the stellar surface. It therefore works on blended lines, regions of the spectrum where line formation mechanisms are not well understood, or stars whose spots have intrinsically different spectra from the rest of the photosphere. We implement the model within the open source starry framework, making it fast, differentiable, and easy to use in both optimization and posterior inference settings. As a proof-of-concept, we use our model to infer the surface map of the brown dwarf WISE 1049-5319B, finding close agreement with the solution of Crossfield et al. (2014). We also discuss Doppler imaging in the context of EPRV studies and describe an interpretable spectral-temporal Gaussian process for stellar spectral variability that we expect will be important for EPRV exoplanet searches.

Our view of the variety of stellar structures pervading the local Milky Way has been transformed by the application of clustering algorithms to the Gaia catalog. In particular, several stellar streams have been recently discovered that are comprised of hundreds to thousands of stars and span several hundred parsecs. We analyze one such structure, Theia 456, a low-density stellar stream extending nearly 200 pc and 20$^{\circ}$ across the sky. By supplementing Gaia astrometric data with spectroscopic metallicities from LAMOST and photometric rotation periods from the Zwicky Transient Facility (ZTF) and the Transiting Exoplanet Survey Satellite (TESS), we establish Theia 456's radial velocity coherence, and we find strong evidence that members of Theia 456 have a common age ($\simeq$175 Myr), common dynamical origin, and formed from chemically homogeneous pre-stellar material ([Fe/H] = $-$0.07 dex). Unlike well-known stellar streams in the Milky Way, which are in its halo, Theia 456 is firmly part of the thin disk. If our conclusions about Theia 456 can be applied to even a small fraction of the remaining $\simeq$8300 independent structures in the Theia catalog, such low-density stellar streams may be ubiquitous. We comment on the implications this has for the nature of star-formation throughout the Galaxy.

Recent studies with IceCube have shown signs of a time-integrated flux of astrophysical neutrinos from point-like sources such as TXS 0506+056 and NGC 1068. Time-variability of this neutrino emission from TXS 0506+056 has been studied extensively by assuming a temporal profile of the possible flare(s) or searching for temporal neutrino correlation with other electromagnetic counterparts. However, experimental evidence of the temporal profile of an astrophysical neutrino signal, besides the TXS 0506+056 source, remains lacking. In this study, we present a new KS-test based method for investigating time-variability. This new method complements the existing time-dependent search methods with a test for arbitrary time-variability, independent of an assumed temporal profile or electromagnetic counterpart. Additionally, this method provides a diagnostic tool for characterizing point-like source candidates in IceCube by distinguishing variable from steady neutrino emission and we show results of applying this method to a small catalog of candidate blazars.

We examine SDO/EVE data to better understand solar flare irradiance, and how that irradiance may vary for large events. We measure scaling laws relating GOES flare classes to irradiance in 21 lines measured with SDO/EVE, formed across a wide range of temperatures, and find that this scaling depends on the line formation temperature. We extrapolate these irradiance values to large events, exceeding X10. In order to create full spectra, however, we need a physical model of the irradiance. We present the first results of a new physical model of solar flare irradiance, NRLFLARE, that sums together a series of flare loops to calculate the spectral irradiance ranging from the X-rays through the far ultraviolet (~ 0 to 1250 Angstroms), constrained by GOES/XRS observations. We test this model against SDO/EVE data. The model spectra and time evolution compares well in high temperature emission, but cooler lines show large discrepancies. We speculate that the discrepancies are likely due to both a non-uniform cross section of the flaring loops as well as opacity effects. We then show that allowing the cross-sectional area to vary with height significantly improves agreement with observations, and is therefore a crucial parameter needed to accurately model the intensity of spectral lines, particularly in the transition region from 4.7 < log T < 6.0.

Second-order mean-motion resonances lead to an interesting phenomenon in the sculpting of the period ratio distribution due to their shape and width in period-ratio/eccentricity space. As the osculating periods librate in resonance, the time-averaged period ratio approaches the exact commensurability. The width of second-order resonances increases with increasing eccentricity, and thus more eccentric systems have a stronger peak at commensurability when averaged over sufficient time. The libration period is short enough that this time-averaging behavior is expected to appear on the timescale of the Kepler mission. Using N-body integrations of simulated planet pairs near the 5:3 and 3:1 mean-motion resonances, we investigate the eccentricity distribution consistent with the planet pairs observed by Kepler. This analysis, an approach independent from previous studies, shows no statistically significant peak at the 3:1 resonance and a small peak at the 5:3 resonance, placing an upper limit on the Rayleigh scale parameter, $\sigma$, of the eccentricity of the observed Kepler planets at $\sigma=0.245$ (3:1) and $\sigma=0.095$ (5:3) at 95% confidence, consistent with previous results from other methods.

Continuous observations play an important role in the studies of solar variability. While such observations can be achieved from space with almost 100% duty cycle, it is difficult to accomplish very high duty cycle from the ground. In this context, we assess the duty cycle that has been achieved from the ground by analyzing the observations of a six station network of identical instruments, Global Oscillation Network Group (GONG). We provide a detailed analysis of the duty cycle using GONG observations spanning over 18 years. We also discuss duty cycle of individual sites and point out various factors that may impact individual site or network duty cycle. The mean duty cycle of the network is 93%, however it reduces by about 5% after all images pass through the stringent quality-control checks. The standard deviations in monthly and yearly duty cycle values are found to be 1.9% and 2.2%, respectively. These results provide a baseline that can be used in the planning of future ground-based networks.

The most intriguing question of modern astronomy is the question of our Universe formation. The Hubble diagram analysis with Type Ia Supernovae (SNe Ia) is widely used to estimate the cosmological parameters with high accuracy. The cosmological measurements allow us to better understand the early stages of the Universe evolution. Nevertheless, the redshift on the SN Ia Hubble diagram includes the contribution from the unknown peculiar velocities that decreases the accuracy of such measurements. We consider the contamination of peculiar velocities of the host galaxies for those SNe that exploded in galaxy clusters. For this purpose, we use the Pantheon cosmological sample and offer the procedure to minimise the influence of such uncertainties.

We aim at building a method that efficiently identifies young low-mass stars and brown dwarfs from low-resolution near-infrared spectra, by studying gravity-sensitive features and their evolution with age. We built a dataset composed of all publicly available ($\sim$2800) near-infrared spectra of dwarfs with spectral types between M0 and L3. First, we investigate methods for the derivation of the spectral type and extinction using comparison to spectral templates, and various spectral indices. Then, we examine gravity-sensitive spectral indices and apply machine learning methods, in order to efficiently separate young ($\lesssim$10 Myr) objects from the field. Using a set of six spectral indices for spectral typing, including two newly defined ones (TLI-J and TLI-K), we are able to achieve a precision below 1 spectral subtype across the entire spectral type range. We define a new gravity-sensitive spectral index (TLI-g) that consistently separates young from field objects, showing a performance superior to other indices from the literature. Even better separation between the two classes can be achieved through machine learning methods which use the entire NIR spectra as an input. Moreover, we show that the H- and K-bands alone are enough for this purpose. Finally, we evaluate the relative importance of different spectral regions for gravity classification as returned by the machine learning models. We find that the H-band broad-band shape is the most relevant feature, followed by the FeH absorption bands at 1.2 $\mu m$ and 1.24 $\mu m$ and the KI doublet at 1.24 $\mu m$.

Quadruply lensed quasars are visible only when the source quasar lies within the diamond caustic of the lensing galaxy. This condition creates a Malmquist-like selection effect in the population of observed quadruply lensed quasars, increasing the true caustic area. The inferred areas are therefore biased low by an amount proportional to the square of the uncertainty in the natural logarithm of the area, $\sigma_{\ln{A}}$. The inferred time delays, $\tau$, vary with the caustic area of the lensing galaxy as $\sqrt{A}$. As a consequence, model time delays are biased low, leading to underestimates of the Hubble constant. Classical Malmquist bias also plays a role, with inferred magnifications biased low. As magnification is strongly anti-correlated with the area of the caustic, this classical Malmquist bias has an effect opposite to, but smaller than, the effect of the area bias. The overall effect is estimated to be small (median bias of less than 1%), but has the potential to be larger. Current posterior plots do not include caustic area, which limits the capability to accurately estimate the bias. With anticipated increases in sample sizes (in particular with the completion of LSST), correction for these biases becomes increasingly important.

Magnetic flux ropes are the centerpiece of solar eruptions. Direct measurements for the magnetic field of flux ropes are crucial for understanding the triggering and energy release processes, yet they remain heretofore elusive. Here we report microwave imaging spectroscopy observations of an M1.4-class solar flare that occurred on 2017 September 6, using data obtained by the Expanded Owens Valley Solar Array. This flare event is associated with a partial eruption of a twisted filament observed in H{\alpha} by the Goode Solar Telescope at the Big Bear Solar Observatory. The extreme ultraviolet (EUV) and X-ray signatures of the event are generally consistent with the standard scenario of eruptive flares, with the presence of double flare ribbons connected by a bright flare arcade. Intriguingly, this partial eruption event features a microwave counterpart, whose spatial and temporal evolution closely follow the filament seen in H{\alpha} and EUV. The spectral properties of the microwave source are consistent with nonthermal gyrosynchrotron radiation. Using spatially resolved microwave spectral analysis, we derive the magnetic field strength along the filament spine, which ranges from 600-1400 Gauss from its apex to the legs. The results agree well with the non-linear force-free magnetic model extrapolated from the pre-flare photospheric magnetogram. We conclude that the microwave counterpart of the erupting filament is likely due to flare-accelerated electrons injected into the filament-hosting magnetic flux rope cavity following the newly reconnected magnetic field lines.

The unusual multiwavelength lightcurves of GRB 101225A are revisited by assuming that it is from an off-axis GRB powered by a newborn magnetar. We show that its optical afterglow lightcurve is fitted with the forward shock model by parameterizing its jet structure as a Gaussian function with a half opening angle of the jet core as $1.67^{\rm o}$. The derived initial Lorentz factor ($\Gamma_0$) is 120, and the viewing angle to the jet axis is $\theta_v=3.7^{\rm o}$. Tentative QPO signatures of $P=488$ seconds and $P=250\sim 300$ seconds are found with a confidence level of 90\% by analysing its X-ray flares observed in the time interval of $[4900,\ 7500]$ seconds. Its global gamma-ray/X-ray lightcurve and the QPO signatures are represented with the magnetar dipole radiation (DR) model by considering the magnetar precession motion, assuming that the magnetar spindown is dominated by the GW emission. The bulk Lorentz factor of the DR ejecta is limited to 8, being much lower than $\Gamma_0$. Comparing GRB 101225A with the extremely off-axis GRB 170817A, we suspect that the nature of the two-component jet in GRB 170817A is a combination of a co-axial GRB jet and a DR ejecta. GRB 101225A would be among the brightest ones of the CDF-S XT2 like X-ray transient population driven by newborn magnetars. Discussion on detectability of its gravitational wave emission is also presented.

We present a large-scale study of stellar rotation for T Tauri stars in the Orion Star-Forming Complex. We use the projected rotational velocity ($v\sin(i)$) estimations reported by the APOGEE-2 collaboration as well as individual masses and ages derived from the position of the stars in the HR diagram, considering Gaia-EDR3 parallaxes and photometry plus diverse evolutionary models. We find an empirical trend for $v\sin(i)$ decreasing with age for low-mass stars ($0.4 M_{\odot} < M_{\ast} < 1.2 M_{\odot}$). Our results support the existence of a mechanism linking $v\sin(i)$ to the presence of accreting protoplanetary disks, responsible for regulating stellar rotation in timescales of about 6 Myr, which is the timescale in which most of the T Tauri stars lose their inner disk. Our results provide important constraints to models of rotation in the early phases of evolution of young stars and their disks.

Binary neutron star mergers (NSMs) have been confirmed as one source of the heaviest observable elements made by the rapid neutron-capture (r-) process. However, modeling NSM outflows -- from the total ejecta masses to their elemental yields -- depends on the unknown nuclear equation of state (EOS) that governs neutron-star structure. In this work, we derive a phenomenological EOS by assuming that NSMs are the dominant sources of the heavy-element material in metal-poor stars with r-process abundance patterns. We start with a population synthesis model to obtain a population of merging neutron star binaries and calculate their EOS-dependent elemental yields. Under the assumption that these mergers were responsible for the majority of r-process elements in the metal-poor stars, we find parameters representing the EOS for which the theoretical NSM yields reproduce the derived abundances from observations of metal-poor stars. For our proof-of-concept assumptions, we find an EOS that is slightly softer than, but still in agreement with, current constraints, e.g., by the Neutron Star Interior Composition Explorer (NICER), with $R_{1.4}=12.25\pm 0.03$~km and $M_{\textrm TOV}$ of $2.17\pm 0.03$~M$_\odot$(statistical uncertainties, neglecting modeling systematics).

Baryonic matter can be accreted onto primordial back holes (PBHs) formed in the early Universe. The radiation from accreting PBHs is capable of altering the evolution of the intergalactic medium (IGM), leaving marks on the global 21 cm signal in the dark ages. For accreting PBHs with mass $M_{\rm PBH}=10^{3}(10^{4})~M_{\odot}$ and mass fraction $f_{\rm PBH}=10^{-1}(10^{-3})$, the brightness temperature deviation $\Delta \delta T_{b}$ reaches $\sim 18~(26)~\rm mK$ at redshift $z\sim 90$ ($\nu \sim 16~\rm MHz$), and the gradient of the brightness temperature $d\delta T_{b}/d\nu$ reaches $ \sim 0.8~(0.5)~\rm mK~MHz^{-1}$ at frequency $\nu\sim 28~\rm MHz$ ($z\sim 50$). For larger PBHs with higher mass fraction, the brightness temperature deviation is larger in the redshift range $z\sim 30-300$ ($\nu\sim 5-46~\rm MHz$), and the gradient is lower at the frequency range $\nu \sim 20-60~\rm MHz$ ($z\sim 23-70$). It is impossible to detect these low frequency radio signals from the Earth due to the influences of the Earth's ionosphere. However, after taking care of the essential factors properly e.g. the foreground and interference, future radio telescope in lunar orbit or on the farside surface of the Moon has a chance to detect the global 21 cm signals impacted by accreting PBHs and distinguish them from the standard model.

In recent years star formation has been discovered in the Milky Way's warp. These stars formed in the warp (warp stars) must eventually settle into the plane of the disc. We use an $N$-body$+$smooth particle hydrodynamics model of a warped galaxy to study how warp stars settle into the disc. By following warp stars in angular momentum space, we show that they first tilt to partially align with the main disc in a time scale of $\sim1$ Gyr. Then, once differential precession halts this process, they phase mix into an axisymmetric distribution on a time scale of $\sim 6$ Gyr. The warp stars end up contaminating the geometric thick disc. Because the warp in our fiducial simulation is growing, the {\it warp stars} settle to a distribution with a negative vertical age gradient as younger stars settle further from the mid-plane. While vertically extended, warp star orbits are still nearly circular and they are therefore subject to radial migration, with a net movement inwards. As a result warp stars can be found throughout the disc. The density distribution of a given population of warp stars evolves from a torus to an increasingly centrally filled-in density distribution. Therefore we argue that, in the Milky Way, warp stars should be found in the Solar Neighbourhood. Moreover, settled warp stars may constitute part of the young flaring population seen in the Milky Way's outskirts.

The tensions between the values of Hubble constant obtained from the early and the late Universe data pose a significant challenge to modern cosmology. Possible modifications of the flat homogeneous isotropic cosmological {\Lambda}CDM model are considered, in which the Universe contains dark energy, cold baryonic matter and dark matter. They are based on general relativity and satisfy two requirements: (1) the value of the Hubble constant, calculated from the value of the Hubble parameter at the recombination by the formulas of flat {\Lambda}CDM model, should be equal to 92% of the one based on low-reshift observations; (2) deviations from the {\Lambda}CDM model should not lead to effects that contradict astronomical observations and estimations obtained thereof. The analysis showed that there are few opportunities for choice. Either we should consider DM with negative pressure -\r{ho}_{dm}c^2<<p_{dm}<0, which weakly affects the evolution of the Universe and the observed manifestations of DM, or we should admit the mechanism of generation of new matter, for example, by a decay of DE.

In this article, we compare the properties of two coronal mass ejections (CMEs) that show similar source region characteristics but different evolutionary behavior in the later phases. We discuss the two events in terms of their near-Sun characteristics, interplanetary evolution, and geo-effectiveness. We carefully analyzed the initiation and propagation parameters of these events to establish the precise CME-ICME connection and their near-Earth consequences. The First event was associated with poor geo-magentic storm disturbance index (Dst $\approx$-20 nT) while the second event is associated with intense geomagnetic storm of DST $\approx$-119 nT. The configuration of the sunspots in the active regions and their evolution are observed by Helioseismic and Magnetic Imager (HMI). For source region imaging, we rely on data obtained from Atmospheric Imaging Assembly (AIA) on board Solar Dynamics Observatory (SDO) and H$\alpha$ filtergrams from Solar Tower Telescope at Aryabhatta Research Institute of Observational Sciences (ARIES). For both the CMEs, flux rope eruptions from the source region triggered flares of similar intensities ($\approx$M1). At the solar source region of the eruptions, we observed circular ribbon flare (CRF) for both the cases, suggesting fan-spine magnetic configuration in the active region corona.

The rotation state of small asteroids is affected by the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect, which is a net torque caused by solar radiation directly reflected and thermally reemitted from the surface. Due to this effect, the rotation period slowly changes, which can be most easily measured in light curves because the shift in the rotation phase accumulates over time quadratically. We collected archived light curves and carried out new photometric observations for asteroids (10115) 1992 SK, (1620) Geographos, and (1685) Toro. We applied the method of light curve inversion to fit observations with a convex shape model. The YORP effect was modeled as a linear change of the rotation frequency $\upsilon \equiv \mathrm{d}\omega / \mathrm{d}t$ and optimized together with other spin and shape parameters. We detected the acceleration $\upsilon = (8.3 \pm 0.6) \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$ of the rotation for asteroid (10115) 1992 SK. This observed value agrees well with the theoretical value of YORP-induced spin-up computed for our shape and spin model. For (1685) Toro, we obtained $\upsilon = (3.3 \pm 0.3) \times 10^{-9}\,\mathrm{rad}\,\mathrm{d}^{-2}$, which confirms an earlier tentative YORP detection. For (1620) Geographos, we confirmed the previously detected YORP acceleration and derived an updated value of $\upsilon$ with a smaller uncertainty. We also included the effect of solar precession into our inversion algorithm, and we show that there are hints of this effect in Geographos' data. The detected change of the spin rate of (10115) 1992 SK has increased the total number of asteroids with YORP detection to ten. In all ten cases, the $\mathrm{d}\omega / \mathrm{d}t$ value is positive, so the rotation of these asteroids is accelerated. It is unlikely to be just a statistical fluke, but it is probably a real feature that needs to be explained.

The present work proposes a new formalism for the inner regions of a neutrino-dominated accretion flows (NDAFs) by considering the self-gravity, where the neutrino opacity is high enough to make neutrinos trapped becoming a dominant factor in the transportation of energy and angular momentum over the magneto rotational instability. We investigate the possibility of gravitational instability and fragmentation to model the highly variable structure of the prompt emission in gamma-ray bursts (GRBs). The results lead us to introduce the gravitational instability, in these inner regions, as a source of a new viscosity which is of the same functional form as that of the $\beta$-prescription of viscosity. Such a consideration brings about fragmentation in the unstable inner disk. In addition, we find the consequent clumpy structure of this area capable to account for the temporal variability of GRB's light curve, especially for the lower choices of the parameter $\beta$, $\sim 10^{-5}$. Finally, we predict the formation of gravitational waves through the migration of fragments before being tidally disrupted. These waves appear to be detectable via a range of current and future detectors from LIGO to Cosmic Explorer.

Protoplanetary discs are complex dynamical systems where several processes may lead to the formation of ring-like structures and planets. These discs are flared following a profile where the vertical scale height increases with radius. In this work, we investigate the role of this disc flaring geometry on the formation of rings and holes. We combine a flattening law change with X-ray and FUV photoevaporative winds. We have used a semi-analytical 1D viscous {\alpha} approach, presenting the evolution of the disc mass and mass rate in a grid of representative systems. Our results show that changing the profile of the flared disc may favour the formation of ring-like features resembling those observed in real systems at the proper evolutionary times, with proper disc masses and accretion rate values. However, these features seem to be short-lived and further enhancements are still needed for better matching all the features seen in real systems.

Machine Learning (ML) is the study of computer algorithms that improve automatically through experience. It is divided into supervised learning, where the computer is presented with examples of entries, and the goal is to learn a general rule that maps inputs to outputs, and unsupervised learning, where no label is provided to the learning algorithm, leaving it alone to find structures. Deep learning is a branch of machine learning based on artificial neural networks, which are an interconnected group of nodes, inspired by a simplification of neurons in a brain. In asteroid dynamics, machine learning methods have been recently used to identify members of asteroid families, small bodies images in astronomical fields, and to identify resonant arguments images of asteroids in three-body resonances, among other applications. Here we will conduct a full review of available literature in the field, and classify it in terms of metrics recently used by other authors to assess the state of the art of applications of machine learning in other astronomical sub-fields. While applications of machine learning to Solar System bodies, an area that includes imaging and spectrophotometry of small bodies, have already reached a state classified as progressing, with more established research communities and methodologies, and articles where the use of ML lead to discovery of new celestial objects or features, or to new insights in the area, ML applied to asteroid dynamics is still in the emerging phase, with smaller groups, methodologies still not well-established, and fewer papers producing new discoveries or insights. The still pioneering nature of research in the area means that much room of improvement is still available, and that several new enticing scientific discoveries may surge by applications of these new techniques in the next years.

Extensive photometric surveys are and will continue producing massive amounts of data on small bodies. Usually, these data will be sparsely obtained at arbitrary (and unknown)rotational phases. Therefore, new methods to process such data need to be developed to make the most of those large catalogs. We aim to produce a method to create phase curves of small bodies considering the uncertainties introduced by the nominal errors in the magnitudes and the effect introduced by rotational variations. We use the SLOAN Moving Objects Catalog data as a benchmark to construct phase curves of all small bodies in there, in u', g', r', i', and z' filters. We will obtain from the phase curves the absolute magnitudes and set up with them the absolute colors, which are the colors of the asteroids not affected by changes in phase angle. We select objects with $\geq3$ observations taken in at least one filter and spanned over a minimum of 5 degrees in phase angle. We developed a method that combines Monte Carlo simulations and Bayesian inference to estimate the absolute magnitudes using the HG$_{12}^*$ photometric system. We obtained almost 15\,000 phase curves, about 12,000 including all five filters. The absolute magnitudes and absolute colors are compatible with previously published data, supporting our method. The method we developed is fully automatic and well suited to be run on large amounts of data. Moreover, it includes the nominal uncertainties in the magnitudes and the whole distribution of possible rotational states of the objects producing, possibly, less precise values, i.e., larger uncertainties, but more accurate, i.e., closer to the actual value. To our knowledge, this work is the first to include the effect of rotational variations in such a manner.

Ultra-fast outflows (UFOs) have been detected in the high-quality X-ray spectra of a number of active galactic nuclei (AGN) with fairly high accretion rates and are thought to significantly contribute to the AGN feedback. After a decade of dedicated study, their launching mechanisms and structure are still not well understood, but variability techniques may provide useful constraints. In this work, therefore, we perform a flux-resolved X-ray spectroscopy on a highly accreting and variable NLS1 AGN, 1H 0707-495, using all archival XMM-Newton observations to study the structure of the UFO. We find that the wind spectral lines weaken at higher luminosities, most likely due to an increasing ionization parameter as previously found in a few similar sources. Instead, the velocity is anticorrelated with the luminosity, which is opposite to the trend observed in the NLS1 IRAS 13224-3809. Furthermore, the detection of the emission lines, which are not observed in IRAS 13224-3809, indicates a wind with a larger opening angle in 1H 0707-495, presumably due to a higher accretion rate. The emitting gas is found to remain broadly constant with the luminosity. We describe the variability of the wind with a scenario where the strong radiation extends the launch radius outwards and shields the outer emitting gas, similarly to super-Eddington compact objects, although other possible explanations are discussed. Our work provides several hints for a multi-phase outflow in 1H 0707-495.

We present results of detailed investigation of the poorly studied X-ray pulsar 2S 1845$-$024 based on the data obtained with $NuSTAR$ observatory during the type I outburst in 2017. Neither pulse phase-averaged, nor phase-resolved spectra of the source show evidence for a cyclotron absorption feature. We also used the data obtained from other X-ray observatories ($Swift$, $XMM-Newton$ and $Chandra$) to study the spectral properties as a function of orbital phase. The analysis revealed a high hydrogen column density for the source reaching $\sim$10$^{24}$ cm$^{-2}$ around the periastron. Using high-quality $Chandra$ data we were able to obtain an accurate localization of 2S 1845$-$024 at R.A. = 18$^{h}$48$^{m}$16$^{s}$.8 and Dec. = $-$2$^{\circ}$25'25".1 (J2000) that allowed us to use infrared (IR) data to roughly classify the optical counterpart of the source as an OB supergiant at the distance of $\gtrsim$15 kpc.

The propagation of cosmic rays can be described as a diffusive motion in most galactic environments. High-energy gamma-rays measured by Fermi have allowed inference of a gradient in the cosmic-ray density and spectral energy behavior in the Milky Way, which is not predicted by models. Here, a turbulence-dependent diffusion model is used to probe different types of cosmic-ray diffusion tensors. Crucially, it is demonstrated that the observed gradients can be explained through turbulence-dependent energy-scaling of the diffusion tensor.

The differences between Narrow Line Seyfert 1 galaxies (NLS1s) and Broad Line AGNs (BLAGNs) are not completely understood; it is thought that they may have different inclinations and/or physical characteristics. The FWHM(Hb)-luminosities correlations are found for NLS1s and their origin is the matter of debate. Here we investigated the spectroscopic parameters and their correlations considering a dusty, cone model of AGN. We apply a simple conical dust distribution (spreading out of broad line region, BLR), assuming that the observed surface of the model is in a good correlation with MIR emission. The dusty cone model in combination with a BLR provides the possibility to estimate luminosity dependence on the cone inclination. The FWHM(Hb)-luminosities correlations obtained from model in comparison with observational data show similarities which may indicate the influence of AGN inclination and structure to this correlation. An alternative explanation for FWHM(Hb)-luminosities correlations is the selection effect by the black hole mass. These FWHM(Hb)-luminosities correlations may be related to the starburst in AGNs, as well. The distinction between spectral properties of the NLS1s and BLAGNs could be caused by multiple effects: beside physical differencies between NLS1s and BLAGNs (NLS1s have lighter black hole mass than BLAGNs), inclination of the conical AGN geometry may have important role as well, where NLS1s may be seen in lower inclination angles.

The magnetic activity of the Sun changes with the solar cycle. Similar cycles are found in other stars as well, but their details are not known to a similar degree. Characterising stellar magnetic cycles is important for the understanding of the stellar and solar dynamos that are driving the magnetic activity. We present spectropolarimetric observations of five young, solar-type stars, and compare them to previous observations, with the aim to identify and characterise stellar equivalents of the solar cycle. We use Zeeman-Doppler imaging (ZDI) to map the surface magnetic field and brightness of our targets. The magnetic field is decomposed into spherical harmonic expansions, from which we report the strengths of the axisymmetric vs. non-axisymmetric, and poloidal vs. toroidal components, and compare them to the Rossby numbers of the stars. We present five new ZDI-maps of young, solar-type stars from Dec 2017. Of special interest is the case of V1358 Ori, that has gone through a polarity reversal between our observations and earlier ones. A less evident polarity reversal might also have occurred in HD 35296. There is a preference for more axisymmetric field, and possibly more toroidal field, for the more active stars with lower Rossby number, but a larger sample should be studied to draw any strong conclusions from this. For most of the individual stars, the amounts of toroidal and poloidal field have stayed on similar levels as in earlier observations. We find evidence for a magnetic polarity reversal having occurred in V1358 Ori. \c{hi}1 Ori could be an interesting target for future observations, with a possible short magnetic cycle of a few years. The correlation between the brightness maps and the magnetic field is mostly poor, which could indicate the presence of small-scale magnetic features of different polarities, that cancel each other out, and are not resolved in our maps.

Some Quantum Gravity (QG) theories, aiming at unifying general relativity and quantum mechanics, predict an energy-dependent modified dispersion relation for photons in vacuum leading to a Violation of Lorentz Invariance (LIV). One way to test these theories is to monitor TeV photons time-of-flight emitted by distant, highly energetic and highly variable astrophysical sources such as flaring active galactic nuclei. Only one time-lag detection was reported so far. We have recently shown however that significant intrinsic time-lags should arise from in situ blazar emission processes at TeV energies and should consequently interfere with LIV searches. In this contribution we will review how intrinsic time delays and LIV-induced propagation effects can simultaneously impact blazars' observed spectral energy distributions and lightcurves. Using a time-dependent approach, we provide predictions on both contributions for various cases in the frame of a standard one zone synchrotron-self-Compton (SSC) model. We will also introduce hints and methods on how to disentangle intrinsic time delays from extrinsic ones in order to highlight LIV effects.

The inner solar system possesses a unique orbital structure in which there are no planets inside the Mercury orbit and the mass is concentrated around the Venus and Earth orbits. The origins of these features still remain unclear. We propose a novel concept that the building blocks of the inner solar system formed at the dead-zone inner edge in the early phase of the protosolar disk evolution, where the disk is effectively heated by the disk accretion. First, we compute the dust evolution in a gas disk with a dead zone and obtain the spatial distribution of rocky planetesimals. The disk is allowed to evolve both by a viscous diffusion and magnetically-driven winds. We find that the rocky planetesimals are formed in concentrations around $\sim$ 1 au with a total mass comparable to the mass of the current inner solar system in the early phase of the disk evolution within $\lesssim0.1$ Myr. Based on the planetesimal distribution and the gas disk structure, we subsequently perform \textit{N}-body simulations of protoplanets to investigate the dynamical configuration of the planetary system. We find that the protoplanets can grow into planets without significant orbital migration because of the rapid clearing of the inner disk by the magnetically-driven disk winds. Our model can explain the origins of the orbital structure of the inner solar system. Several other features such as the rocky composition can also be explained by the early formation of rocky planetesimals.

Direction dependent calibration of widefield radio interferometers estimates the systematic errors along multiple directions in the sky. This is necessary because with most systematic errors that are caused by effects such as the ionosphere or the receiver beam shape, there is significant spatial variation. Fortunately, there is some deterministic behavior of these variations in most situations. We enforce this underlying smooth spatial behavior of systematic errors as an additional constraint onto spectrally constrained direction dependent calibration. Using both analysis and simulations, we show that this additional spatial constraint improves the performance of multi-frequency direction dependent calibration.

Following a dynamical encounter with Sgr A*, binaries in the Galactic Centre (GC) can be tidally separated and one member star ejected as a hyper-velocity star (HVS) with a velocity beyond the escape speed of the Milky Way. As GC-born objects located in more observationally accessible regions of the sky, HVSs offer insight into the stellar population in the inner parsecs of the Milky Way. We perform a suite of simulations ejecting HVSs from the GC, exploring how detectable HVS populations depend on assumptions concerning the GC stellar population, focusing on those that are both gravitationally unbound from the Galaxy and which would appear in current and/or future data releases from the \textit{Gaia} space mission with precise astrometry and measured radial velocities. We show that predictions are sensitive to two parameters in particular: the shape of the stellar initial mass function (IMF) in the GC and the ejection rate of HVSs. The absence of confident HVS candidates in \textit{Gaia} Data Release 2 excludes scenarios in which the HVS ejection rate is $\gtrsim3\times10^{-2} \, \mathrm{yr^{-1}}$. Stricter constraints will be placed on these parameters when more HVS candidates are unearthed in future \textit{Gaia} data releases -- assuming recent determinations of the GC IMF shape, one confident HVS \textit{at minimum} is expected in \textit{Gaia} DR3 and DR4 as long as the HVS ejection rate is greater than $\sim 10^{-3} \, \mathrm{yr^{-1}}$ and $\sim10^{-5} \, \mathrm{yr^{-1}}$, respectively.

Primordial Black Holes (PBH) can form in the early universe and might comprise a significant fraction of the dark matter. Interestingly, they are accompanied by the generation of Gravitational Wave (GW) signals and they could contribute to the merger events currently observed by the LIGO/Virgo Collaboration (LVC). In this thesis, we study the PBH scenario, addressing various properties at the formation epoch and the computation of abundance beyond the Gaussian paradigm, while also developing the theoretical description of PBH evolution through accretion and mergers, with particular focus on modelling their GW signatures. In a second part, we compare the primordial scenario with current GW data, seizing the possible contribution of PBH binaries to LVC signals and forecasting the potential of future GW detectors, such as Einstein Telescope and LISA, to detect mergers of primordial binaries and the stochastic GW background induced at second order by the PBH formation mechanism.

We investigate the effect of damped oscillations on a nearly flat inflationary potential and the features they produce in the power-spectrum and bi-spectrum. We compare the model with the Planck data using Plik unbinned and CamSpec clean likelihood and we are able to obtain noticeable improvement in fit compared to the power-law $\Lambda$CDM model. We are able to identify three plausible candidates each for the two likelihoods used. We find that the best-fit to Plik and CamSpec likelihoods match closely to each other. The improvement comes from various possible outliers at the intermediate to small scales. We also compute the bi-spectrum for the best-fits. At all limits, the amplitude of bi-spectrum, $f_{NL}$ is oscillatory in nature and its peak value is determined by the amplitude and frequency of the oscillations in the potential, as expected. We find that the bi-spectrum consistency relation strictly holds at all scales in all the best-fit candidates.

It is well established that the transverse MHD waves are ubiquitous in the solar corona. One of the possible mechanisms for heating both open (e.g. coronal holes) and closed (e.g. coronal loops) magnetic field regions of the solar corona is due to the MHD wave-driven turbulence. In this work, we have studied the variation in the filling factor of overdense structures in the solar corona due to the generation of the transverse MHD wave-driven turbulence. Using 3D MHD simulations, we estimate the density filling factor of an open magnetic structure by calculating the fraction of the volume occupied by the overdense plasma structures to the entire volume of the simulation domain. Next, we perform forward modeling and generate synthetic spectra of Fe XIII 10749 \AA\ and 10800 \AA\ density sensitive line pairs using FoMo. Using the synthetic images, we again estimate the filling factors. The estimated filling factors obtained from both methods are in reasonable agreement. Also, our results match fairly well with the observations of filling factors in coronal holes and loops. Our results show that the generation of turbulence increases the filling factor of the solar corona.

Solar campfires are fine-scale heating events, recently observed by Extreme Ultraviolet Imager (EUI), onboard Solar Orbiter. Here we use EUI 174\AA\ images, together with EUV images from SDO/AIA, and line-of-sight magnetograms from SDO/HMI to investigate the magnetic origin of 52 randomly selected campfires in the quiet solar corona. We find that (i) the campfires are rooted at the edges of photospheric magnetic network lanes; (ii) most of the campfires reside above the neutral line between majority-polarity magnetic flux patch and a merging minority-polarity flux patch, with a flux cancelation rate of $\sim$10$^{18}$Mx hr$^{-1}$; (iii) some of the campfires occur repeatedly from the same neutral line; (iv) in the large majority of instances, campfires are preceded by a cool-plasma structure, analogous to minifilaments in coronal jets; and (v) although many campfires have `complex' structure, most campfires resemble small-scale jets, dots, or loops. Thus, `campfire' is a general term that includes different types of small-scale solar dynamic features. They contain sufficient magnetic energy ($\sim$10$^{26}$-10$^{27}$ erg) to heat the solar atmosphere locally to 0.5--2.5MK. Their lifetimes range from about a minute to over an hour, with most of the campfires having a lifetime of $<$10 minutes. The average lengths and widths of the campfires are 5400$\pm$2500km and 1600$\pm$640km, respectively. Our observations suggest that (a) the presence of magnetic flux ropes may be ubiquitous in the solar atmosphere and not limited to coronal jets and larger-scale eruptions that make CMEs, and (b) magnetic flux cancelation is the fundamental process for the formation and triggering of most campfires.

We present the discovery of two nearly identically-sized sub-Neptune transiting planets orbiting HD 63935, a bright ($V=8.6$ mag), sun-like ($T_{eff}=5560K$) star at 49 pc. TESS identified the first planet, HD 63935 b (TOI-509.01), in Sectors 7 and 34. We identified the second signal (HD 63935 c) in Keck HIRES and Lick APF radial velocity data as part of our followup campaign. It was subsequently confirmed with TESS photometry in Sector 34 as TOI-509.02. Our analysis of the photometric and radial velocity data yields a robust detection of both planets with periods of $9.0600 \pm 0.007$ and $21.40 \pm 0.0019$ days, radii of $2.99 \pm 0.14$ and $2.90 \pm 0.13$ $R_\oplus$, and masses of $10.8 \pm 1.8$ and $11.1 \pm 2.4$ $M_\oplus$. We calculate densities for planets b and c consistent with a few percent of the planet mass in hydrogen/helium envelopes. We also describe our survey's efforts to choose the best targets for JWST atmospheric followup. These efforts suggest that HD 63935 b will have the most clearly visible atmosphere of its class. It is the best target for transmission spectroscopy (ranked by Transmission Spectroscopy Metric, a proxy for atmospheric observability) in the so-far uncharacterized parameter space comprising sub-Neptune-sized (2.6 $R_\oplus$ $<$ $R_p$ $<$ 4 $R_\oplus$), moderately-irradiated (100 $F_\oplus$ $<$ $F_p$ $<$ 1000 $F_\oplus$) planets around G-stars. Planet c is also a viable target for transmission spectroscopy, and given the indistinguishable masses and radii of the two planets, the system serves as a natural laboratory for examining the processes that shape the evolution of sub-Neptune planets.

The HI 21-cm optical depth ($\tau_b$) can be considerably large as the kinetic and spin temperature of the inter-galactic medium (IGM) is expected to be very low during cosmic dawn. It will be particularly higher at regions with HI over-density. We revisit the validity of the widely used linearized equation for estimating the HI 21-cm differential brightness temperature ($T_b$) which assumes $\tau_b << 1$ and approximates $[1-\exp({-\tau_b})]$ as $\tau_b$. We consider two scenarios, one without any additional cooling mechanism or radio background (referred as the standard scenario) and the other (referred as the excess-cooling} scenario) assumes the EDGES-like absorption profile and an excess cooling mechanism. We find that given a measured global absorption signal, consistent with the standard (excess-cooling) scenario, the linearized equation overestimates the spin temperature by $\sim 5\%(10\%)$. Further, using numerical simulations, we study the impact that the large optical depth has on various signal statistics. We observe that the variance, skewness and kurtosis, calculated at simulation resolution ($\sim 0.5 h^{-1} \, {\rm Mpc}$), are over-predicted up to $\sim 30\%$, $30\%$ and $15\%$ respectively for the standard and up to $\sim 90\%$, $50\%$ and $50\%$ respectively for the excess-cooling scenario. Moreover, we find that the probability distribution function of $T_b$ is squeezed and becomes more Gaussian in shape if no approximation is made. The spherically averaged HI power spectrum is overpredicted by up to $\sim 25 \%$ and $80\%$ at all scales for the standard and excess-cooling scenarios respectively.

Dark matter particles may bind with nuclei if there exists an attractive force of sufficient strength. We show that a dark photon mediator of mass $\sim (10 - 100)$ MeV that kinetically mixes with Standard Model electromagnetism at the level of $\sim 10^{-3}$ generates keV-scale binding energies between dark matter and heavy elements, while forbidding the ability to bind with light elements. In underground direct detection experiments, the formation of such bound states liberates keV-scale energy in the form of electrons and photons, giving rise to mono-energetic electronic signals with a time-structure that may contain daily and seasonal modulations. We show that data from liquid-xenon detectors provides exquisite sensitivity to this scenario, constraining the galactic abundance of such dark particles to be at most $\sim 10^{-18} - 10^{-12}$ of the galactic dark matter density for masses spanning $\sim (1 - 10^5)$ GeV. However, an exponentially small fractional abundance of these dark particles is enough to explain the observed electron recoil excess at XENON1T.

The next decade is expected to see the launch of one or more space based gravitational wave detectors: the European lead Laser Interferometer Space Antenna (LISA); and one or more Chinese mission concepts, Taiji and TianQin. One of the primary scientific targets for these missions are the mergers of black holes with masses between $10^3 M_\odot$ and $10^8 M_\odot$. These systems may produce detectable electromagnetic signatures in additional to gravitational waves due to the presence of gas in mini-disks around each black hole, and a circumbinary disk surrounding the system. The electromagnetic emission may occur before, during and after the merger. In order to have the best chance of capturing all phases of the emission it is imperative that the gravitational wave signals can be detected in low latency, and used to produce reliable estimates for the sky location and distance to help guide the search for counterparts. The low latency detection also provides a starting point for the ``global fit'' of the myriad signals that are simultaneously present in the data. Here a low latency analysis pipeline is presented that is capable of analyzing months of data in just a few hours using a laptop from the last decade. The problem of performing a global fit is avoided by whitening out the bright foreground produced by nearby galactic binaries. The performance of the pipeline is illustrated using simulated data from the LISA Data Challenge.

We consider the Multiverse as an ensemble of universes. Using standard statistical physics analysis we get that the Cosmological Constant (CC) is exponentially small. The small and finite CC is achieved without any anthropic reasoning. We then quantize the CC. The quantization allows a precise summation of the possible contributions and using the measured value of the CC yields a prediction on the temperature of the Multiverse that we define. Furthermore, quantization allows the interpretation of a single Universe as a superposition of different eigenstates with different energy levels rather than the existence of an actual Multiverse.

The standard model is a remarkably consistent and complete quantum field theory but its coupling to gravity and the Higgs field remain problematic, as reflected in the cosmological constant problem, the Weyl anomaly, and the hierarchy puzzle. We point out that 36 conformally-coupled dimension-zero scalar fields can simultaneously cancel the vacuum energy and both terms in the Weyl anomaly, if the Higgs and graviton fields are emergent. The cancellation is highly non-trivial: given the standard model gauge group $SU(3)\times SU(2)\times U(1)$, it requires precisely 48 Weyl fermions, i.e., three generations of standard model fermions, including right-handed neutrinos. The dimension-zero scalars have a four-derivative Lagrangian, usually taken to imply vacuum instability. However, using the Euclidean inner product natural in the context of our recent proposal arXiv:2109.06204, we find no negative norm or negative energy states. Hence the vacuum is stable. Moreover, the scalars possess a scale invariant power spectrum extending to long wavelengths, suggesting a new explanation for the primordial scalar perturbations in cosmology, without the need for inflation. These intriguing results, spanning a vast range of scales, suggest dimension-zero scalars may play a key role in fundamental physics. We discuss how the Higgs and graviton fields might emerge in this context.

Similarly to light, gravitational waves can be gravitationally lensed as they propagate near massive astrophysical objects such as galaxies, stars, or black holes. In recent years, forecasts have suggested a reasonable chance of strong gravitational-wave lensing detections with the LIGO-Virgo-Kagra detector network at design sensitivity. As a consequence, methods to analyze lensed detections have seen rapid development. However, the impact of higher-order modes on the lensing analyses is still under investigation. In this work, we show that the presence of higher-order modes enables the identification of the individual images types for the observed gravitational-wave events when two lensed images are detected, which would lead to unambiguous identification of lensing. In addition, we show that we can analyze higher-order-mode content with greater accuracy with strongly lensed gravitational wave events.

We consider the possibility that dark matter (DM) only interacts with the Standard Model leptons, but not quarks at tree-level, and analyze the future lepton collider prospects of such leptophilic DM in the mono-photon and mono-$Z$ (both leptonic and hadronic) channels. Adopting a model-independent effective field theory framework, we consider all possible dimension-6 operators of scalar-pseudoscalar (S-P), vector-axialvector (V-A) and tensor-axialtensor (T-AT) types for a fermionic DM and derive the collider sensitivities on the effective cut-off scale $\Lambda$ as a function of the DM mass. As a concrete example, we take the beam configurations of the International Linear Collider with $\sqrt s=1$ TeV and 1000 fb$^{-1}$ integrated luminosity, including the effect of beam polarization, and show that it can probe leptophilic DM at $3\sigma$ level up to $\Lambda$ values of 4.8 TeV, 6.5 TeV and 5.3 TeV for the S-P, V-A and T-AT type operators respectively. This is largely complementary to the direct and indirect searches for leptophilic DM, and can potentially provide the best-ever sensitivity in the low-mass DM regime.