We investigated the properties of AGN environments, particularly environments where the association of luminous galaxies (LGs) is found within 4~Mpc from AGNs with redshifts of 0.8 -- 1.1. For comparison, three additional AGN environments, (namely, AGNs of all types, type~1 AGNs with X-ray and/or radio detection, and type~2 AGNs) and an environment of blue $M_{*}$, characteristic luminosity of the Schechter function, galaxies were investigated. The cross-correlation function with the surrounding galaxies was measured and compared between the AGN and blue galaxy samples. We also compared the distributions of color, absolute magnitude, and stellar mass of the galaxies around such target objects. The properties of clusters detected using surrounding galaxies selected based on a photometric redshift were examined and compared for different samples. The target AGNs were drawn from the Million Quasars (MILLIQUAS) catalog, and the blue galaxies were drawn from six redshift survey catalogs (SDSS, WiggleZ, DEEP2, VVDS, VIPERS, and PRIMUS). The galaxies used as a measure of the environment around the targets are drawn from S18a internal data released by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We found that, among the five AGN and blue galaxy samples considered, the environment of AGN-LG pairs is the most enriched with luminous galaxies. We also found an enhancement in the number of mass-selected clusters in the AGN-LG pair sample against those in the other samples. The results obtained in this study indicate that existence of multiple clusters is the major driver in the association of AGNs and LGs, rather than a single large-mass dark matter halo hosting the AGN.

Solar eruptions are the most powerful drivers of space weather. To understand their cause and nature, it is crucial to know how the coronal magnetic field evolves before eruption. Here we study the formation process of a relatively large-scale magnetic flux rope (MFR) in active region NOAA~12371 that erupts with a major flare and coronal mass ejection on 2015 June 21. A data-driven numerical magnetohydrodynamic model is employed to simulate three-dimensional coronal magnetic field evolution of one-day duration before the eruption. Comparison between the observed features and our modeled magnetic field discloses how the pre-eruption MFR forms. Initially, the magnetic field lines were weakly twisted as being simple sheared arcades. Then a long MFR was formed along the polarity inversion line due to the complex photospheric motion, which is mainly shearing rather than twisting. The presence of the MFR is evidenced by a coherent set of magnetic field lines with twist number above unity. Below the MFR a current sheet is shown in the model, suggesting that tether-cutting reconnection plays a key role in the MFR formation. The MFR's flux grows as more and more field lines are twisted due to continuous injection of magnetic helicity by the photospheric motions. Meanwhile, the height of the MFR's axis increases monotonely from its formation. By an analysis of the decay index of its overlying field, we suggest that it is because the MFR runs into the torus instability regime and becomes unstable that finally triggers the eruption.

Magnetic flux ropes (MFRs) are thought to be the central structure of solar eruptions, and their ideal MHD instabilities can trigger the eruption. Here we performed a study of all the MFR configurations that lead to major solar flares, either eruptive or confined, from 2011 to 2017 near the solar disk center. The coronal magnetic field is reconstructed from observed magnetograms, and based on magnetic twist distribution, we identified the MFR, which is defined as a coherent group of magnetic field lines winding an axis with more than one turn. It is found that 90% of the events possess pre-flare MFRs, and their three-dimensional structures are much more complex in details than theoretical MFR models. We further constructed a diagram based on two parameters, the magnetic twist number which controls the kink instability (KI), and the decay index which controls the torus instability (TI). It clearly shows lower limits for TI and KI thresholds, which are $n_{\rm crit} = 1.3$ and $|T_w|_{\rm crit} = 2$, respectively, as all the events above $n_{\rm crit}$ and nearly 90% of the events above $|T_w|_{\rm crit}$ erupted. Furthermore, by such criterion, over 70% of the events can be discriminated between eruptive and confined flares, and KI seems to play a nearly equally important role as TI in discriminating between the two types of flare. There are more than half of events with both parameters below the lower limits, and 29% are eruptive. These events might be triggered by magnetic reconnection rather than MHD instabilities.

We investigate a global extreme-ultraviolet (EUV) wave associated with a coronal mass ejection (CME)-driven shock on 2017 September 10. The EUV wave is transmitted by north- and south-polar coronal holes (CHs), which is observed by the Solar Dynamics Observatory (SDO) and Solar Terrestrial Relations Observatory A (STEREO-A) from opposite sides of the Sun. We obtain key findings on how the EUV wave interacts with multiple coronal structures, and on its connection with the CME-driven shock: (1) the transmitted EUV wave is still connected with the shock that is incurvated to the Sun, after the shock has reached the opposite side of the eruption; (2) the south CH transmitted EUV wave is accelerated inside an on-disk, low-density region with closed magnetic fields, which implies that an EUV wave can be accelerated in both open and closed magnetic field regions; (3) part of the primary EUV wavefront turns around a bright point (BP) with a bipolar magnetic structure when it approaches a dim, low-density filament channel near the BP; (4) the primary EUV wave is diffused and apparently halted near the boundaries of remote active regions (ARs) that are far from the eruption, and no obvious AR related secondary waves are detected; (5) the EUV wave extends to an unprecedented scale of ~360{\deg} in latitudes, which is attributed to the polar CH transmission. These results provide insights into the effects of coronal density and magnetic field distributions on the evolution of an EUV wave, and into the connection between the EUV wave and the associated CME-driven shock.

We report the result of optical identifications of FIRST radio sources with the Hyper Suprime-Cam Subaru Strategic Program survey (HSC-SSP). The positional cross-match within 1" between the FIRST and HSC-SSP catalogs (i ~< 26) produced more than 3600 optical counterparts in the 156 deg^2 of the HSC-SSP field. The matched counterparts account for more than 50% of the FIRST sources in the search field, which substantially exceed previously reported fractions of SDSS counterparts (i ~< 22) of ~30%. Among the matched sample, 9% are optically unresolved sources such as radio-loud quasars. The optically faint (i > 21) radio galaxies (RGs) show that the fitting linear function of the 1.4 GHz source counts has a slope that is flatter than that of the bright RGs, while optically faint radio quasars show a slope steeper than that of bright radio quasars. The optically faint RGs show a flat slope in the i-band number counts down to 24 mag, implying either less massive or distant radio-active galactic nuclei (AGNs) beyond 24 mag. The photometric redshift and the comparison of colors with the galaxy models show that most of the matched RGs are distributed at redshifts from 0 to 1.5. The optically faint sample includes the high radio-loudness sources that are not seen in the optically bright sample. Such sources are located at redshift z > 1. This study gives ~1500 radio AGNs lying at the optically faint end and high-redshift regime not probed by previous searches.

We investigate how a weak coronal mass ejection (CME) launched on 2016 October 8 without obvious signatures in the low corona produced a relatively intense geomagnetic storm. Remote sensing observations from SDO, STEREO and SOHO and in situ measurements from WIND are employed to track the CME from the Sun to the Earth. Using a graduated cylindrical shell (GCS) model, we estimate the propagation direction and the morphology of the CME near the Sun. CME kinematics are determined from the wide-angle imaging observations of STEREO A and are used to predict the CME arrival time and speed at the Earth. We compare ENLIL MHD simulation results with in situ measurements to illustrate the background solar wind where the CME was propagating. We also apply a Grad--Shafranov technique to reconstruct the flux rope structure from in situ measurements in order to understand the geo-effectiveness associated with the CME magnetic field structure. Key results are obtained concerning how a weak CME can generate a relatively intense geomagnetic storm: (1) there were coronal holes at low latitudes, which could produce high speed streams (HSSs) to interact with the CME in interplanetary space; (2) the CME was bracketed between a slow wind ahead and a HSS behind, which enhanced the southward magnetic field inside the CME and gave rise to the unexpected geomagnetic storm.

We have measured the clustering of galaxies around active galactic nuclei (AGN) for which single-epoch virial masses of the super-massive black hole (SMBH) are available to investigate the relation between the large scale environment of AGNs and the evolution of SMBHs. The AGN samples used in this work were derived from the Sloan Digital Sky Survey (SDSS) observations and the galaxy samples were from 240~deg$^{2}$ S15b data of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). The investigated redshift range is 0.6--3.0, and the masses of the SMBHs lie in the range $10^{7.5}$--$10^{10}$~M$_{\odot}$. The absolute magnitude of the galaxy samples reaches to $M_{\lambda 310}$ $\sim$ $-$18 at rest frame wavelength 310~nm for the low-redshift end of the samples. More than 70\% of the galaxies in the analysis are blue. We found a significant dependence of the cross-correlation length on redshift, which primarily reflects the brightness dependence of the galaxy clustering. At the lowest redshifts the cross-correlation length increases from 7~$h^{-1}$~Mpc around $M_{\lambda 310}$$ = $$-19$~mag to $>$10~$h^{-1}$~Mpc beyond $M_{\lambda 310}$$ = $$-20$~mag. No significant dependence of the cross-correlation length on BH mass was found for whole galaxy samples dominated by blue galaxies, while there was an indication of BH mass dependence in the cross-correlation with red galaxies. These results provides a picture of the environment of AGNs studied in this paper being enriched with blue starforming galaxies, and a fraction of the galaxies are being evolved to red galaxies along with the evolution of SMBHs in that system.

We examine the clustering of quasars over a wide luminosity range, by utilizing 901 quasars at $\overline{z}_{\rm phot}\sim3.8$ with $-24.73<M_{\rm 1450}<-22.23$ photometrically selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) S16A Wide2 date release and 342 more luminous quasars at $3.4<z_{\rm spec}<4.6$ having $-28.0<M_{\rm 1450}<-23.95$ from the Sloan Digital Sky Survey (SDSS) that fall in the HSC survey fields. We measure the bias factors of two quasar samples by evaluating the cross-correlation functions (CCFs) between the quasar samples and 25790 bright $z\sim4$ Lyman Break Galaxies (LBGs) in $M_{\rm 1450}<-21.25$ photometrically selected from the HSC dataset. Over an angular scale of \timeform{10.0"} to \timeform{1000.0"}, the bias factors are $5.93^{+1.34}_{-1.43}$ and $2.73^{+2.44}_{-2.55}$ for the low and high luminosity quasars, respectively, indicating no luminosity dependence of quasar clustering at $z\sim4$. It is noted that the bias factor of the luminous quasars estimated by the CCF is smaller than that estimated by the auto-correlation function (ACF) over a similar redshift range, especially on scales below \timeform{40.0"}. Moreover, the bias factor of the less-luminous quasars implies the minimal mass of their host dark matter halos (DMHs) is $0.3$-$2\times10^{12}h^{-1}M_{\odot}$, corresponding to a quasar duty cycle of $0.001$-$0.06$.

We present the luminosity function of z=4 quasars based on the Hyper Suprime-Cam Subaru Strategic Program Wide layer imaging data in the g, r, i, z, and y bands covering 339.8 deg^2. From stellar objects, 1666 z~4 quasar candidates are selected by the g-dropout selection down to i=24.0 mag. Their photometric redshifts cover the redshift range between 3.6 and 4.3 with an average of 3.9. In combination with the quasar sample from the Sloan Digital Sky Survey in the same redshift range, the quasar luminosity function covering the wide luminosity range of M1450=-22 to -29 mag is constructed. It is well described by a double power-law model with a knee at M1450=-25.36+-0.13 mag and a flat faint-end slope with a power-law index of -1.30+-0.05. The knee and faint-end slope show no clear evidence of redshift evolution from those at z~2. The flat slope implies that the UV luminosity density of the quasar population is dominated by the quasars around the knee, and does not support the steeper faint-end slope at higher redshifts reported at z>5. If we convert the M1450 luminosity function to the hard X-ray 2-10keV luminosity function using the relation between UV and X-ray luminosity of quasars and its scatter, the number density of UV-selected quasars matches well with that of the X-ray-selected AGNs above the knee of the luminosity function. Below the knee, the UV-selected quasars show a deficiency compared to the hard X-ray luminosity function. The deficiency can be explained by the lack of obscured AGNs among the UV-selected quasars.

We report the detection of a double planetary system around the evolved intermediate-mass star HD 47366 from precise radial-velocity measurements at Okayama Astrophysical Observatory, Xinglong Station, and Australian Astronomical Observatory. The star is a K1 giant with a mass of 1.81+-0.13M_sun, a radius of 7.30+-0.33R_sun, and solar metallicity. The planetary system is composed of two giant planets with minimum mass of 1.75^{+0.20}_{-0.17}Mjup and 1.86^{+0.16}_{-0.15}Mjup, orbital period of 363.3^{+2.5}_{-2.4} d and 684.7^{+5.0}_{-4.9} d, and eccentricity of 0.089^{+0.079}_{-0.060} and 0.278^{+0.067}_{-0.094}, respectively, which are derived by a double Keplerian orbital fit to the radial-velocity data. The system adds to the population of multi-giant-planet systems with relatively small orbital separations, which are preferentially found around evolved intermediate-mass stars. Dynamical stability analysis for the system revealed, however, that the best-fit orbits are unstable in the case of a prograde configuration. The system could be stable if the planets were in 2:1 mean-motion resonance, but this is less likely considering the observed period ratio and eccentricity. A present possible scenario for the system is that both of the planets have nearly circular orbits, namely the eccentricity of the outer planet is less than ~0.15, which is just within 1.4sigma of the best-fit value, or the planets are in a mutually retrograde configuration with a mutual orbital inclination larger than 160 degree.

We have obtained deep J & Ks-band images centered on a bright radio quiet QSO UM402 (z_{em}=2.856) using IRCS camera and AO systems on Subaru Telescope, as well as retrieved WFC3/F140W archive images. A faint galaxy (m_{k}=23.32 +/- 0.05 in the Vega system) that lies ~ 2.4" north of the QSO sightline has been clearly resolved in all three deep high resolution datasets, and appears as an irregular galaxy with two close components in the Ks-band images (separation ~ 0.3"). Given the small impact parameter (b=19.6 kpc, at z_{lls}=2.531), as well as the red color of (J-Ks)_{vega} ~ 1.6, it might be a candidate galaxy giving rise to the Lyman Limit system absorption at z_{abs}=2.531 seen in the QSO spectrum. After carefully subtracting the PSF from the QSO images, the host galaxy of this bright radio quiet QSO at z ~ 3 was marginally revealled. We placed a low limit of the host component of m_{k} ~ 23.3 according to our analyses.

We report the detection of a substellar companion orbiting an evolved intermediate-mass ($M_\star=2.4\,M_\odot$) star HD 14067 (G9 III) using precise Doppler technique. Radial velocities of this star can be well fitted either by a periodic Keplerian variation with a decreasing linear velocity trend (P=1455 days, $K_1=92.2$ m s$^{-1}$, $e=0.533$, and $\dot{\gamma}=-22.4$ m s$^{-1}$ yr$^{-1}$) or a single Keplerian orbit without linear trend (P=2850 days, $K_1=100.1$ m s$^{-1}$, and $e=0.697$). The minimum mass ($m_2\sin{i}=7.8\,M_{\rm J}$ for the model with a linear trend, or $m_2\sin{i}=9.0\,M_{\rm J}$ for the model without a linear trend) suggests a long-period giant planet around an evolved intermediate-mass star. The eccentricity of the orbit is among the highest known for planets ever detected around evolved stars.