GREX-PLUS (Galaxy Reionization EXplorer and PLanetary Universe Spectrometer) is a mission candidate for a JAXA's strategic L-class mission to be launched in the 2030s. Its primary sciences are two-fold: galaxy formation and evolution and planetary system formation and evolution. The GREX-PLUS spacecraft will carry a 1.2 m primary mirror aperture telescope cooled down to 50 K. The two science instruments will be onboard: a wide-field camera in the 2-8 $\mu$m wavelength band and a high resolution spectrometer with a wavelength resolution of 30,000 in the 10-18 $\mu$m band. The GREX-PLUS wide-field camera aims to detect the first generation of galaxies at redshift $z>15$. The GREX-PLUS high resolution spectrometer aims to identify the location of the water ``snow line'' in proto-planetary disks. Both instruments will provide unique data sets for a broad range of scientific topics including galaxy mass assembly, origin of supermassive blackholes, infrared background radiation, molecular spectroscopy in the interstellar medium, transit spectroscopy for exoplanet atmosphere, planetary atmosphere in the Solar system, and so on.

In this paper, we investigate the impact of higher-order distortions on the precise measurement of weak gravitational lensing shear and flexion. We begin by defining generalized higher-order distortions and outlining methods for measuring them. Then, using several lens models, we examine how these distortions affect shear and flexion measurements. Our results show that neglecting higher-order distortions can introduce systematic errors of a few percent in both shear and flexion measurements, indicating that these effects cannot be ignored. Although the strength of these errors depends on factors such as lensing strength and the size of background sources, we demonstrate that simultaneous measurement of higher-order distortions can reduce the systematic errors to below 1% in most cases.

Euclid will provide deep NIR imaging to $\sim$26.5 AB magnitude over $\sim$59 deg$^2$ in its deep and auxiliary fields. The Cosmic DAWN survey complements the deep Euclid data with matched depth multiwavelength imaging and spectroscopy in the UV--IR to provide consistently processed Euclid selected photometric catalogs, accurate photometric redshifts, and measurements of galaxy properties to a redshift of $z\sim 10$. In this paper, we present an overview of the survey, including the footprints of the survey fields, the existing and planned observations, and the primary science goals for the combined data set.

The 2017 observing campaign of the Event Horizon Telescope (EHT) delivered the first very long baseline interferometry (VLBI) images at the observing frequency of 230 GHz, leading to a number of unique studies on black holes and relativistic jets from active galactic nuclei (AGN). In total, eighteen sources were observed: the main science targets, Sgr A* and M87 along with various calibrators. We investigated the morphology of the sixteen AGN in the EHT 2017 data set, focusing on the properties of the VLBI cores: size, flux density, and brightness temperature. We studied their dependence on the observing frequency in order to compare it with the Blandford-Königl (BK) jet model. We modeled the source structure of seven AGN in the EHT 2017 data set using linearly polarized circular Gaussian components and collected results for the other nine AGN from dedicated EHT publications, complemented by lower frequency data in the 2-86 GHz range. Then, we studied the dependences of the VLBI core flux density, size, and brightness temperature on the frequency measured in the AGN host frame. We compared the observations with the BK jet model and estimated the magnetic field strength dependence on the distance from the central black hole. Our results indicate a deviation from the standard BK model, particularly in the decrease of the brightness temperature with the observing frequency. Either bulk acceleration of the jet material, energy transfer from the magnetic field to the particles, or both are required to explain the observations.

Little Red Dots (LRDs) are compact, point-like sources characterized by their red color and broad Balmer lines, which have been debated to be either dominated by active galactic nuclei (AGN) or dusty star-forming galaxies (DSFGs). Here we report two LRDs (ID9094 and ID2756) at z$_{\rm spec}$>7, recently discovered in the JWST FRESCO GOODS-North field. Both satisfy the "v-shape" colors and compactness criteria for LRDs and are identified as Type-I AGN candidates based on their broad H$\beta$ emission lines (full width at half maximum: 2280$\pm$490 km/s for ID9094 and 1070$\pm$240 km/s for ID2756) and narrow [OI] lines ($\sim$ 300-400 km/s). To investigate their nature, we conduct deep NOEMA follow-up observations targeting the [CII] 158${\rm \mu m}$ emission line and the 1.3 mm dust continuum. We do not detect [CII] or 1.3 mm continuum emission for either source. Notably, in the scenario that the two LRDs were DSFGs, we would expect significant detections: $>16\sigma$ for [CII] and $>3\sigma$ for the 1.3 mm continuum of ID9094, and $>5\sigma$ for [CII] of ID2756. Using the 3$\sigma$ upper limits of [CII] and 1.3 mm, we perform two analyses: (1) UV-to-FIR spectral energy distribution (SED) fitting with and without AGN components, and (2) comparison of their properties with the L$_{[CII]}$-SFR$_{tot}$ empirical relation. Both analyses are consistent with a scenario where AGN activity may contribute to the observed properties, though a dusty star-forming origin cannot be fully ruled out. Our results highlight the importance of far-infrared observations for studying LRDs, a regime that remains largely unexplored.

We present new high resolution imaging of a light-scattering dust ring and halo around the young star HD 35841. Using spectroscopic and polarimetric data from the Gemini Planet Imager in H-band (1.6 microns), we detect the highly inclined (i=85 deg) ring of debris down to a projected separation of ~12 au (~0.12") for the first time. Optical imaging from HST/STIS shows a smooth dust halo extending outward from the ring to >140 au (>1.4"). We measure the ring's scattering phase function and polarization fraction over scattering angles of 22-125 deg, showing a preference for forward scattering and a polarization fraction that peaks at ~30% near the ansae. Modeling of the scattered-light disk indicates that the ring spans radii of ~60-220 au, has a vertical thickness similar to that of other resolved dust rings, and contains grains as small as 1.5 microns in diameter. These models also suggest the grains have a low porosity, are more likely to consist of carbon than astrosilicates, and contain significant water ice. The halo has a surface brightness profile consistent with that expected from grains pushed by radiation pressure from the main ring onto highly eccentric but still bound orbits. We also briefly investigate arrangements of a possible inner disk component implied by our spectral energy distribution models, and speculate about the limitations of Mie theory for doing detailed analyses of debris disk dust populations.

Red giants undergo dramatic and complex structural transformations as they evolve. Angular momentum is transported between the core and envelope during this epoch, a poorly understood process. Here, we infer envelope and core rotation rates from Kepler observations of $\sim$1517 red giants. While many measurements are consistent with the existing studies, our investigation reveals systematic changes in the envelope-to-core rotation ratio and we report the discovery of anomalies such as clump stars with rapidly rotating cores, and red giants with envelopes rotating faster than their cores. We propose binary interactions as a possible mechanism by which some of these cores and envelopes are spun up. These results pose challenges to current theoretical expectations and can have major implications for compact remnants born from stellar cores.

Solar flares are often accompanied by filament/prominence eruptions ($\sim10^{4}$ K and $\sim 10^{10-11}$ cm$^{-3}$), sometimes leading to coronal mass ejections (CMEs) that directly affect the Earth's environment. `Superflares' are found on some active solar-type (G-type main-sequence) stars, but the association of filament eruptions/CMEs has not been established. Here we show that our optical spectroscopic observation of the young solar-type star EK Draconis reveals the evidence for a stellar filament eruption associated with a superflare. This superflare emitted a radiated energy of $2.0\times10^{33}$ erg, and blue-shifted hydrogen absorption component with a large velocity of $-510$ km s$^{-1}$ was observed shortly after. The temporal changes in the spectra greatly resemble those of solar filament eruptions. Comparing this eruption with solar filament eruptions in terms of the length scale and velocity strongly suggests that a stellar CME occurred. The erupted filament mass of $1.1\times10^{18}$ g is 10 times larger than those of the largest solar CMEs. The massive filament eruption and an associated CME provide the opportunity to evaluate how they affect the environment of young exoplanets/young Earth and stellar mass/angular-momentum evolution.

In 2017, the next generation Very Large Array (ngVLA) Science Advisory Council, together with the international astronomy community, developed a set of five Key Science Goals (KSGs) to inform, prioritize and refine the technical capabilities of a future radio telescope array for high angular resolution operation from 1.2 - 116 GHz with 10 times the sensitivity of the Jansky VLA and ALMA. The resulting KSGs, which require observations at centimeter and millimeter wavelengths that cannot be achieved by any other facility, represent a small subset of the broad range of astrophysical problems that the ngVLA will be able address. This document presents an update to the original ngVLA KSGs, taking account of new results and progress in the 7+ years since their initial presentation, again drawing on the expertise of the ngVLA Science Advisory Council and the broader community in the ngVLA Science Working Groups. As the design of the ngVLA has also matured substantially in this period, this document also briefly addresses initial expectations for ngVLA data products and processing that will be needed to achieve the KSGs. The original ngVLA KSGs endure as outstanding problems of high priority. In brief, they are: (1) Unveiling the Formation of Solar System Analogues; (2) Probing the Initial Conditions for Planetary Systems and Life with Astrochemistry; (3) Charting the Assembly, Structure, and Evolution of Galaxies from the First Billion Years to the Present; (4) Science at the Extremes: Pulsars as Laboratories for Fundamental Physics; (5) Understanding the Formation and Evolution of Stellar and Supermassive Black Holes in the Era of Multi-Messenger Astronomy.

The current progress in the detection of terrestrial type exoplanets has opened a new avenue in the characterization of exoplanetary atmospheres and in the search for biosignatures of life with the upcoming ground-based and space missions. To specify the conditions favorable for the origin, development and sustainment of life as we know it in other worlds, we need to understand the nature of astrospheric, atmospheric and surface environments of exoplanets in habitable zones around G-K-M dwarfs including our young Sun. Global environment is formed by propagated disturbances from the planet-hosting stars in the form of stellar flares, coronal mass ejections, energetic particles, and winds collectively known as astrospheric space weather. Its characterization will help in understanding how an exoplanetary ecosystem interacts with its host star, as well as in the specification of the physical, chemical and biochemical conditions that can create favorable and/or detrimental conditions for planetary climate and habitability along with evolution of planetary internal dynamics over geological timescales. A key linkage of (astro) physical, chemical, and geological processes can only be understood in the framework of interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary and Earth sciences. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets will significantly expand the current definition of the habitable zone to the biogenic zone and provide new observational strategies for searching for signatures of life. The major goal of this paper is to describe and discuss the current status and recent progress in this interdisciplinary field and to provide a new roadmap for the future development of the emerging field of exoplanetary science and astrobiology.

CASA, the Common Astronomy Software Applications package, is the primary data processing software for the Atacama Large Millimeter/submillimeter Array (ALMA) and NSF's Karl G. Jansky Very Large Array (VLA), and is frequently used also for other radio telescopes. The CASA software can process data from both single-dish and aperture-synthesis telescopes, and one of its core functionalities is to support the data reduction and imaging pipelines for ALMA, VLA and the VLA Sky Survey (VLASS). CASA has recently undergone several exciting new developments, including an increased flexibility in Python (CASA 6), support of Very Long Baseline Interferometry (VLBI), performance gains through parallel imaging, data visualization with the new Cube Analysis Rendering Tool for Astronomy (CARTA), enhanced reliability and testing, and modernized documentation. These proceedings of the 2019 Astronomical Data Analysis Software & Systems (ADASS) conference give an update of the CASA project, and detail how these new developments will enhance user experience of CASA.

The r-process nucleosynthesis in core-collapse supernovae (CC-SNe) is studied, with a focus on the explosion scenario induced by rotation and strong magnetic fields. Nucleosynthesis calculations are conducted based on magneto-hydrodynamical explosion models with a wide range of parameters for initial rotation and magnetic fields. The explosion models are classified in two different types: i.e., prompt-magnetic-jet and delayed-magnetic-jet, for which the magnetic fields of proto-neutron stars (PNSs) during collapse and the core-bounce are strong and comparatively moderate, respectively. Following the hydrodynamical trajectories of each explosion model, we confirmed that r-processes successfully occur in the prompt-magnetic-jets, which produce heavy nuclei including actinides. On the other hand, the r-process in the delayed-magnetic-jet is suppressed, which synthesizes only nuclei up to the second peak ($A \sim 130$). Thus, the r-process in the delayed-magnetic-jets could explain only "weak r-process" patterns observed in metal-poor stars rather than the "main r-process", represented by the solar abundances. Our results imply that core-collapse supernovae are possible astronomical sources of heavy r-process elements if their magnetic fields are strong enough, while weaker magnetic explosions may produce "weak r-process" patterns ($A \lesssim 130$). We show the potential importance and necessity of magneto-rotational supernovae for explaining the galactic chemical evolution, as well as abundances of r-process enhanced metal-poor stars. We also examine the effects of the remaining uncertainties in the nature of PNSs due to weak interactions that determine the final neutron-richness of ejecta. Additionally, we briefly discuss radioactive isotope yields in primary jets (e.g., $^{56}$Ni), with relation to several optical observation of SNe and relevant high-energy astronomical phenomena.

The HR 8799 system uniquely harbors four young super-Jupiters whose orbits can provide insights into the system's dynamical history and constrain the masses of the planets themselves. Using the Gemini Planet Imager (GPI), we obtained down to one milliarcsecond precision on the astrometry of these planets. We assessed four-planet orbit models with different levels of constraints and found that assuming the planets are near 1:2:4:8 period commensurabilities, or are coplanar, does not worsen the fit. We added the prior that the planets must have been stable for the age of the system (40 Myr) by running orbit configurations from our posteriors through $N$-body simulations and varying the masses of the planets. We found that only assuming the planets are both coplanar and near 1:2:4:8 period commensurabilities produces dynamically stable orbits in large quantities. Our posterior of stable coplanar orbits tightly constrains the planets' orbits, and we discuss implications for the outermost planet b shaping the debris disk. A four-planet resonance lock is not necessary for stability up to now. However, planet pairs d and e, and c and d, are each likely locked in two-body resonances for stability if their component masses are above $6~M_{\rm{Jup}}$ and $7~M_{\rm{Jup}}$, respectively. Combining the dynamical and luminosity constraints on the masses using hot-start evolutionary models and a system age of $42 \pm 5$~Myr, we found the mass of planet b to be $5.8 \pm 0.5~M_{\rm{Jup}}$, and the masses of planets c, d, and e to be $7.2_{-0.7}^{+0.6}~M_{\rm{Jup}}$ each.

We investigated the ability of the Euclid telescope to detect galaxy-scale gravitational lenses. To do so, we performed a systematic visual inspection of the $0.7\,\rm{deg}^2$ Euclid Early Release Observations data towards the Perseus cluster using both the high-resolution $I_{\scriptscriptstyle\rm E}$ band and the lower-resolution $Y_{\scriptscriptstyle\rm E}$, $J_{\scriptscriptstyle\rm E}$, $H_{\scriptscriptstyle\rm E}$ bands. Each extended source brighter than magnitude 23 in $I_{\scriptscriptstyle\rm E}$ was inspected by 41 expert human classifiers. This amounts to $12\,086$ stamps of $10^{\prime\prime}\,\times\,10^{\prime\prime}$. We found $3$ grade A and $13$ grade B candidates. We assessed the validity of these $16$ candidates by modelling them and checking that they are consistent with a single source lensed by a plausible mass distribution. Five of the candidates pass this check, five others are rejected by the modelling, and six are inconclusive. Extrapolating from the five successfully modelled candidates, we infer that the full $14\,000\,{\rm deg}^2$ of the Euclid Wide Survey should contain $100\,000^{+70\,000}_{-30\,000}$ galaxy-galaxy lenses that are both discoverable through visual inspection and have valid lens models. This is consistent with theoretical forecasts of $170\,000$ discoverable galaxy-galaxy lenses in Euclid. Our five modelled lenses have Einstein radii in the range 0.\!\!^{\prime\prime}68\,<\,\theta_\mathrm{E}\,<1.\!\!^{\prime\prime}24, but their Einstein radius distribution is on the higher side when compared to theoretical forecasts. This suggests that our methodology is likely missing small-Einstein-radius systems. Whilst it is implausible to visually inspect the full Euclid dataset, our results corroborate the promise that Euclid will ultimately deliver a sample of around $10^5$ galaxy-scale lenses.

We present high-resolution and spatially-matched observations with JWST and ALMA of a starburst galaxy (PACS-830) at $z=1.46$. The NIRCam observations mainly trace the stellar light while the CO ($J$=5--4) observations map the dense molecular gas at kpc scales. Both datasets reveal the morphology to be that of a gas/dust rich bulge with two extending arms, together resembling a grand-design spiral galaxy. The more pronounced arm contributes 21 $\pm$ 6\% of the total CO emission. These results demonstrate that starburst activity at high redshift can be triggered, without undergoing a highly disruptive major merger. We assess the strength and distribution of star formation using two tracers: (1) Polycyclic Aromatic Hydrocarbons (PAHs) emission detected at $8~\mu$m ($L_8$) with a MIRI/F1800W image, and (2) $L_\mathrm{IR}$, inferred from the CO ($J$=5--4) map. The spatial profiles of the $L_\mathrm{IR}$ and $L_8$ are dissimilar, thus leading to a significant deficit of mid-IR ($L_8$) emission in the nucleus. We hypothesize that this is due to the destruction of PAH molecules by the intense ionizing radiation field or decreased emission in the photodissociation region, as seen in nearby star-forming regions and consistent with the galaxy-wide properties of distant starbursts. This study reveals spatial variations in the $L_8$ to $L_\mathrm{IR}$ ratio for the first time at $z>1$, in agreement with expectations from theory. Our analysis underscores the pivotal role of joint high-resolution observations with JWST and ALMA in discerning the different phases of the interstellar medium (ISM) and revealing internal physics in galaxy substructures.

Red-shifted components of chromospheric emission lines in the hard X-ray impulsive phase of solar flares have recently been studied through their 30 s evolution with the high resolution of IRIS. Radiative-hydrodynamic flare models show that these redshifts are generally reproduced by electron-beam generated chromospheric condensations. The models produce large ambient electron densities, and the pressure broadening of hydrogen Balmer series should be readily detected in observations. To accurately interpret upcoming spectral data of flares with the DKIST, we incorporate non-ideal, non-adiabatic line broadening profiles of hydrogen into the RADYN code. These improvements allow time-dependent predictions for the extreme Balmer line wing enhancements in solar flares. We study two chromospheric condensation models, which cover a range of electron beam fluxes ($1-5 \times 10^{11}$ erg s$^{-1}$ cm$^{-2}$) and ambient electron densities ($1 - 60 \times 10^{13}$ cm$^{-3}$) in the flare chromosphere. Both models produce broadening and redshift variations within 10 s of the onset of beam heating. In the chromospheric condensations, there is enhanced spectral broadening due to large optical depths at H$\alpha$, H$\beta$, and H$\gamma$, while the much lower optical depth of the Balmer series H12$-$H16 provides a translucent window into the smaller electron densities in the beam-heated layers below the condensation. The wavelength ranges of typical DKIST/ViSP spectra of solar flares will be sufficient to test the predictions of extreme hydrogen wing broadening and accurately constrain large densities in chromospheric condensations.

The spectra of SN 2003dh, identified in the afterglow of GRB030329, are modeled using radiation transport codes. It is shown that SN 2003dh had a high explosion kinetic energy ($\sim 4 \times 10^{52}$ erg in spherical symmetry), making it one of the most powerful hypernovae observed so far, and supporting the case for association between hypernovae and Gamma Ray Bursts. However, the light curve derived from fitting the spectra suggests that SN 2003dh was not as bright as SN 1998bw, ejecting only $\sim 0.35\Msun$ of \Nifs. The spectra of SN 2003dh resemble those of SN 1998bw around maximum, but later they look more like those of the less energetic hypernova SN 1997ef. The spectra and the inferred light curve can be modeled adopting a density distribution similar to that used for SN 1998bw at $v > 25,000$\kms but more like that of SN 1997ef at lower velocities. The mass of the ejecta is $\sim 8\Msun$, somewhat less than in the other two hypernovae. The progenitor must have been a massive star ($M \sim 35-40\Msun$), as for other hypernovae. The need to combine different one-dimensional explosion models strongly indicates that SN 2003dh was an asymmetric explosion.

We investigate environmental effects on evolution of bright cluster galaxies in a Lambda-dominated cold dark matter universe using a combination of dissipationless N-body simulations and a semi-analytic galaxy formation model. We incorporate effects of ram-pressure stripping (RPS) and minor merger-induced small starburst (minor burst) into our model. By considering minor burst, observed morphology-radius relation is successfully reproduced. When we do not consider minor burst, the RPS hardly increases the intermediate B/T population. In addition, the RPS and minor burst are not important for colours or star formation rates of galaxies in the cluster core if star formation time-scale is properly chosen, because the star formation is sufficiently suppressed by consumption of the cold gas. We also find that SF in bulge-dominated galaxies is mainly terminated by starburst induced by major mergers in all environments.

We report on the multi-epoch observations of H2O maser emission in star forming region OH 43.8-0.1 carried out with VLBI Exploration of Radio Astrometry (VERA). The large-scale maser distributions obtained by single-beam VLBI mapping reveal new maser spots scattered in area of 0.7 x 1.0 arcsec, in addition to a `shell-like' structure with a scale of 0.3 x 0.5 arcsec which was previously mapped by Downes et al.(1979). Proper motions are also obtained for 43 spots based on 5-epoch monitoring with a time span of 281 days. The distributions of proper motions show a systematic outflow in the north-south direction with an expansion velocity of ~8 km/s, and overall distributions of maser spots as well as proper motions are better represented by a bipolar flow plus a central maser cluster with a complex structure, rather than a shell with a uniform expansion such as those found in Cep A R5 and W75N VLA2. The distance to OH 43.8-0.1 is also estimated based on the statistical parallax, yielding D = 2.8 +/- 0.5 kpc. This distance is consistent with a near kinematic distance and rules out a far kinematics distance (~9 kpc), and the LSR velocity of OH 43.8-0.1 combined with the distance provides a constraint on the flatness of the galactic rotation curve, that there is no systematic difference in rotation speeds at the Sun and at the position of OH 43.8-0.1, which is located at the galacto-centric radius of ~6.3 kpc.

We try to identify the sources and systematics of the uncertainty of the masses of the central black holes in active galactic nuclei (AGNs) which have been obtained by reverberation-mapping methods. We characterize the broad H$\beta$ emission-line profiles by the ratio of their full-width at half maximum (FWHM) to their line dispersion, i.e., the second moment of the line profile. We use this parameter to separate the reverberation-mapped AGNs into two populations, the first with narrower H$\beta$ lines with relatively extended wings, and the second with broader lines that are relatively flat-topped. The first population is characterized by higher Eddington ratios than the second. Within each population, we calibrate the black-hole mass scale by comparison of the reverberation-based mass with that predicted by the bulge velocity dispersion. We also use the distribution of ratios of the reverberation-based mass to the velocity-dispersion mass prediction in a comparison with a ``generalized thick disk'' model in order to see if inclination can plausibly account for the observed distribution. We find that the line dispersion is a less biased parameter in general than FWHM for black hole mass estimation, although we show that it is possible to empirically correct for the bias introduced by using FWHM to characterize the emission-line width. We also argue that inclination effects are apparent in a small subset of objects with the narrowest emission lines. Our principal conclusion is that the H$\beta$ profile is sensitive primarily to Eddington ratio, but that inclination effects play a role in some cases. (abridged)