In this paper a characteristics-based open boundary condition (CBC) is proposed for the magnetohydrodynamic (MHD) system of equations. The algorithm is carefully designed and implemented in the context of a high-order flux reconstruction (FR) scheme under the Generalized Lagrange Multiplier (GLM)-MHD system of equations. It is implemented by adding the contribution of the characteristic equation directly to the corrected flux term in the FR scheme dispensing with solving time-dependent characteristic equations along boundary faces. The CBC method is shown to be more accurate and robust than commonly used zero normal derivative (ZND) and approximate Riemann solver boundary conditions (ARBC) in solving 1D, 2D, and 3D test problems. The CBC method is successfully applied to simulate challenging problems of magnetic reconnection for which other options failed to get stable results over long-period time integration.

We explore the multiphase structure of the circumgalactic medium (CGM) probed by synthetic spectra through a cosmological zoom-in galaxy formation simulation. We employ a Bayesian method for modelling a combination of absorption lines to derive physical properties of absorbers with a formal treatment of detections, including saturated systems, and non-detections in a uniform manner. We find that in the lines of sight passing through localized density structures, absorption lines of low, intermediate and high ions are present in the spectrum and overlap in velocity space. Low, intermediate and high ions can be combined to derive the mass-weighted properties of a density-varying peak, although the ions are not co-spatial within the structure. By contrast, lines of sight that go through the hot halo only exhibit detectable HI and high ions. In such lines of sight, the absorption lines are typically broad due to the complex velocity fields across the entire halo. We show that the derived gas density, temperature, and metallicity match closely the corresponding HI mass-weighted averages along the LOS. We also show that when the data quality allows, our Bayesian technique allows one to recover the underlying physical properties of LOS by incorporating both detections and non-detections. It is especially useful to include non-detections, of species such as NV or NeVIII, when the number of detections of strong absorbers, such as HI and OVI, is smaller than the number of model parameters (density, temperature, and metallicity).

Ubiquitous detections of cold/warm gas around galaxies indicate that the circumgalactic medium (CGM) is multiphase and dynamic. Recent state-of-the-art cosmological galaxy simulations have generally underproduced the column density of cold halo gas. We argue that this may be due to a mismatch of spatial resolution in the circumgalactic space and the relevant physical scales at which the cold gas operates. Using semi-analytic calculations and a set of magnetohydrodynamic (MHD) simulations, we present a multiphase model of the gaseous halos around galaxies, the circumgalactic mist (CGm). The CGm model is based on the idea that the observed cold halo gas may be a composite of cold, dense and small cloudlets embedded in a hot diffuse halo, resembling terrestrial clouds and mist. We show that the resulting cold gas from thermal instabilities conforms to a characteristic column density of $N_{\rm H}\approx 10^{17}\rm{cm^{-2}}$ as predicted by the $c_s t_{\rm cool}$ ansatz. The model implies a large number of cold clumps in the inner galactic halo with a small volume filling factor but large covering fraction. The model also naturally gives rise to spatial extents and differential covering fractions of cold, warm and hot gas. To self-consistently model the co-evolution of the CGM and star formation within galaxies, future simulations must address the mismatch of the spatial resolution and characteristic scale of cold gas.

We investigate the $\gamma$-ray and X-ray properties of the flat spectrum radio quasar (FSRQ) \4c50 at redshift $z= 1.517$. The {\it Fermi} Large Area Telescope (LAT) data indicate that this source was in an active state since 2013 July. During the active period, the $\gamma$-ray flux increased by more than a factor of three and two distinct flares were detected with the variability timescale as short as several hours. The $\gamma$-ray spectra can be well fitted by a log-parabola. From the fitting, we find a correlation between the peak energy and spectral curvature for the $\gamma$-ray spectra, which is the first time seen in $\gamma$-ray emission from a blazar. The {\it Swift} X-ray Telescope (XRT) data show that the source was variable at X-ray energies, but no evidence shows flux or spectral changes related to the $\gamma$-ray activity. The broad-band X-ray spectrum obtained with {\it Swift} XRT and {\it NuSTAR} is well described by a broken PL model, with an extremely hard spectrum ($\Gamma_{1} \sim 0.1$) below the break energy, $E_{\rm break} \sim 2.1 {\rm keV}$, and $\Gamma_{2} \sim 1.5$ above the break energy. The spectral steepening below $\sim 3$ keV is likely due to the low energy cut-off in the energy distribution of the photon-emitting electron population. Both the $\gamma$-ray and X-ray emission appear harder when brighter. The broad-band spectral energy distribution (SED) is constructed for the source, and we provide a model fit to the SED. Our modeling suggests that the emission region should be outside the broad line region, and the properties of the region indicate a jet with 42\% of the Eddington power during the active state.

We introduce a Bayesian approach for modeling Voigt profiles in absorption spectroscopy and its implementation in the python package, BayesVP, publicly available at https://github.com/cameronliang/BayesVP. The code fits the absorption line profiles within specified wavelength ranges and generates posterior distributions for the column density, Doppler parameter, and redshifts of the corresponding absorbers. The code uses publicly available efficient parallel sampling packages to sample posterior and thus can be run on parallel platforms. BayesVP supports simultaneous fitting for multiple absorption components in high-dimensional parameter space. We provide other useful utilities in the package, such as explicit specification of priors of model parameters, continuum model, Bayesian model comparison criteria, and posterior sampling convergence check.

We present the first global 3D simulations of thermal convection in the oblate envelopes of rapidly-rotating solar-type stars. This has been achieved by exploiting the capabilities of the new Compressible High-ORder Unstructured Spectral difference (CHORUS) code. We consider rotation rates up to 85\% of the critical (breakup) rotation rate, which yields an equatorial radius that is up to 17\% larger than the polar radius. This substantial oblateness enhances the disparity between polar and equatorial modes of convection. We find that the convection redistributes the heat flux emitted from the outer surface, leading to an enhancement of the heat flux in the polar and equatorial regions. This finding implies that lower-mass stars with convective envelopes may not have darker equators as predicted by classical gravity darkening arguments. The vigorous high-latitude convection also establishes elongated axisymmetric circulation cells and zonal jets in the polar regions. Though the overall amplitude of the surface differential rotation, $\Delta \Omega$, is insensitive to the oblateness, the oblateness does limit the fractional kinetic energy contained in the differential rotation to no more than 61\%. Furthermore, we argue that this level of differential rotation is not enough to have a significant impact on the oblateness of the star.

We analyze circumgalactic medium (CGM) in a suite of high-resolution cosmological re-simulations of a Milky-Way size galaxy and show that CGM properties are quite sensitive to details of star formation--feedback loop modelling. The simulation that produces a realistic late-type galaxy, fails to reproduce existing observations of the CGM. In contrast, simulation that does not produce a realistic galaxy has the predicted CGM in better agreement with observations. This illustrates that properties of galaxies and properties of their CGM provide strong ${\it complementary}$ constraints on the processes governing galaxy formation. Our simulations predict that column density profiles of ions are well described by an exponential function of projected distance $d$: $N \propto e^{-d/h_s}$. Simulations thus indicate that the sharp drop in absorber detections at larger distances in observations does not correspond to a "boundary" of an ion, but reflects the underlying steep exponential column density profile. Furthermore, we find that ionization energy of ions is tightly correlated with the scale height $h_s$: $h_s \propto E_{\rm ion}^{0.74}$. At $z \approx 0$, warm gas traced by low-ionization species (e.g., Mg II and C IV) has $ h_s \approx 0.03-0.07 R_{\rm vir}$, while higher ionization species (O VI and Ne VIII) have $h_s \approx 0.32-0.45R_{\rm vir}$. Finally, the scale heights of ions in our simulations evolve slower than the virial radius for $z\leq 2$, but similarly to the halo scale radius, $r_s$. Thus, we suggest that the column density profiles of galaxies at different redshifts should be scaled by $r_s$ rather than the halo virial radius.

We present a measurement of the cosmic microwave background (CMB) gravitational lensing potential using data from the first two seasons of observations with SPTpol, the polarization-sensitive receiver currently installed on the South Pole Telescope (SPT). The observations used in this work cover 100 deg$^2$ of sky with arcminute resolution at 150 GHz. Using a quadratic estimator, we make maps of the CMB lensing potential from combinations of CMB temperature and polarization maps. We combine these lensing potential maps to form a minimum-variance (MV) map. The lensing potential is measured with a signal-to-noise ratio of greater than one for angular multipoles between $100< L <250$. This is the highest signal-to-noise mass map made from the CMB to date and will be powerful in cross-correlation with other tracers of large-scale structure. We calculate the power spectrum of the lensing potential for each estimator, and we report the value of the MV power spectrum between $100< L <2000$ as our primary result. We constrain the ratio of the spectrum to a fiducial $\Lambda$CDM model to be $A_{\rm MV}=0.92 \pm 0.14 {\rm\, (Stat.)} \pm 0.08 {\rm\, (Sys.)}$. Restricting ourselves to polarized data only, we find $A_{\rm POL}=0.92 \pm 0.24 {\rm\, (Stat.)} \pm 0.11 {\rm\, (Sys.)}$. This measurement rejects the hypothesis of no lensing at $5.9 \sigma$ using polarization data alone, and at $14 \sigma$ using both temperature and polarization data.

We present a novel and powerful Compressible High-ORder Unstructured Spectral-difference (CHORUS) code for simulating thermal convection and related fluid dynamics in the interiors of stars and planets. The computational geometries are treated as rotating spherical shells filled with stratified gas. The hydrodynamic equations are discretized by a robust and efficient high-order Spectral Difference Method (SDM) on unstructured meshes. The computational stencil of the spectral difference method is compact and advantageous for parallel processing. CHORUS demonstrates excellent parallel performance for all test cases reported in this paper, scaling up to 12,000 cores on the Yellowstone High-Performance Computing cluster at NCAR. The code is verified by defining two benchmark cases for global convection in Jupiter and the Sun. CHORUS results are compared with results from the ASH code and good agreement is found. The CHORUS code creates new opportunities for simulating such varied phenomena as multi-scale solar convection, core convection, and convection in rapidly-rotating, oblate stars.

We adopted an unstructured hydrodynamical solver CharLES to the problem of global convection in the Sun. With the aim to investigate the properties of solar turbulent convection and reproduce differential rotation pattern. We performed simulations in two spherical shells, with 1.3 and 10 million cells. In the first, coarse mesh, the solution does not reproduce realistic convection, and is dominated by numerical effects. In the second mesh, thermal conduction leads to cooling of bottom layers, that could not be compensated by solar irradiance. More simulations in the 10M cells mesh should be performed to investigate the influence of transport coefficients and numerical effects. Our estimate of the code performance suggests, that realistic simulations in even finer grids could be performed for reasonable computational cost.

We present measurements of $E$-mode polarization and temperature-$E$-mode correlation in the cosmic microwave background (CMB) using data from the first season of observations with SPTpol, the polarization-sensitive receiver currently installed on the South Pole Telescope (SPT). The observations used in this work cover 100~\sqdeg\ of sky with arcminute resolution at $150\,$GHz. We report the $E$-mode angular auto-power spectrum ($EE$) and the temperature-$E$-mode angular cross-power spectrum ($TE$) over the multipole range $500 < \ell \leq5000$. These power spectra improve on previous measurements in the high-$\ell$ (small-scale) regime. We fit the combination of the SPTpol power spectra, data from \planck\, and previous SPT measurements with a six-parameter \LCDM cosmological model. We find that the best-fit parameters are consistent with previous results. The improvement in high-$\ell$ sensitivity over previous measurements leads to a significant improvement in the limit on polarized point-source power: after masking sources brighter than 50\,mJy in unpolarized flux at 150\,GHz, we find a 95\% confidence upper limit on unclustered point-source power in the $EE$ spectrum of $D_\ell = \ell (\ell+1) C_\ell / 2 \pi < 0.40 \ \mu{\mbox{K}}^2$ at $\ell=3000$, indicating that future $EE$ measurements will not be limited by power from unclustered point sources in the multipole range $\ell < 3600$, and possibly much higher in $\ell.$

(Abridged) This paper presents an absorption-line study of the multiphase circumgalactic medium (CGM) based on observations of Lya, CII, CIV, SiII, SiIII, and SiIV absorption transitions in the vicinities of 195 galaxies at redshift z<0.176. The galaxy sample is established based on a cross-comparison between public galaxy and QSO survey data and is characterized by a median redshift of <z>=0.041, a median projected distance of <d>=362 kpc to the sightline of the background QSO, and a median stellar mass of log(M_star/M_sun) = 9.7 \pm 1.1. Comparing the absorber features identified in the QSO apectra with known galaxy properties has led to strong constraints for the CGM absorption properties at z<~0.176. First, abundant hydrogen gas is observed out to d~500 kpc, well beyond the dark matter halo radius Rh of individual galaxies, with a mean covering fraction of ~60%. In contrast, no heavy elements are detected at d>~0.7 Rh from either low-mass dwarfs or high-mass galaxies. The lack of detected heavy elements in low- and high-ionization states suggests that either there exists a chemical enrichment edge at d~0.7 Rh or gaseous clumps giving rise to the observed absorption lines cannot survive at these large distances. Considering all galaxies at d>Rh leads to a strict upper limit for the covering fraction of heavy elements of ~3% (at a 95% confidence level) over d=(1-9) Rh. At d<Rh, differential covering fraction between low- and high-ionization gas is observed, suggesting that the CGM becomes progressively more ionized from d<0.3 Rh to larger distances. Comparing CGM absorption observations at low and high redshifts shows that at a fixed-fraction of Rh the CGM exhibits stronger mean absorption at z=2.2 than at z~0. We discuss possible pseudo-evolution of the CGM as a result of misrepresentation of halo radius.

Absorption-line spectroscopy of multiply-lensed QSOs near a known foreground galaxy provides a unique opportunity to go beyond the traditional one-dimensional application of QSO probes and establish a crude three-dimensional (3D) map of halo gas around the galaxy that records the line-of-sight velocity field at different locations in the gaseous halo. Two intermediate-redshift galaxies are targeted in the field around the quadruply-lensed QSO HE0435-1223 at redshift z=1.689, and absorption spectroscopy along each of the lensed QSOs is carried out in the vicinities of these galaxies. One galaxy is a typical, star-forming L* galaxy at z=0.4188 and projected distance of rho=50 kpc from the lensing galaxy. The other is a super-L* barred spiral at z=0.7818 and rho=33 kpc. Combining known orientations of the quadruply-lensed QSO to the two foreground galaxies with the observed MgII absorption profiles along individual QSO sightlines has for the first time led to spatially resolved kinematics of tenuous halo gas on scales of 5-10 kpc at z>0.2. A MgII absorber is detected in every sightline observed through the halos of the two galaxies, and the recorded absorber strength is typical of what is seen in previous close QSO--galaxy pair studies. While the multi-sightline study confirms the unity covering fraction of MgII absorbing gas at rho < 50 kpc from star-forming disks, the galaxies also present two contrasting examples of complex halo gas kinematics. Different models, including a rotating disk, collimated outflows, and gaseous streams from either accretion or tidal/ram-pressure stripping, are considered for comparisons with the absorption-line observations, and infalling streams/stripped gas of width >~ 10 kpc are found to best describe the observed gas kinematics across multiple sightlines.

Gravitational lensing of the cosmic microwave background generates a curl pattern in the observed polarization. This "B-mode" signal provides a measure of the projected mass distribution over the entire observable Universe and also acts as a contaminant for the measurement of primordial gravity-wave signals. In this Letter we present the first detection of gravitational lensing B modes, using first-season data from the polarization-sensitive receiver on the South Pole Telescope (SPTpol). We construct a template for the lensing B-mode signal by combining E-mode polarization measured by SPTpol with estimates of the lensing potential from a Herschel-SPIRE map of the cosmic infrared background. We compare this template to the B modes measured directly by SPTpol, finding a non-zero correlation at 7.7 sigma significance. The correlation has an amplitude and scale-dependence consistent with theoretical expectations, is robust with respect to analysis choices, and constitutes the first measurement of a powerful cosmological observable.

We collect the 2LAC and MOJAVE quasi-simultaneous data to investigate the radio-gamma connection of blazars. The cross sample contains 166 sources. The statistic analysis based on this sample confirms positive correlations between these two bands, but the correlations become weaker as the gamma-ray energy increases. The statistic results between various parameters show negative correlations of gamma-ray photon spectral index with gamma-ray loudness for both FSRQs and BL Lacertae objects, positive correlations of gamma-ray variability index with the gamma-ray loudness for FSRQs, a negative correlation of the gamma-ray variability index with the gamma-ray photon spectral index for FSRQs, and negative correlations of gamma-ray photon spectral index with gamma-ray luminosity for FSRQs. These results suggest that the gamma-ray variability may be due to changes inside the gamma-ray emission region like the injected power, rather than changes in the photon density of the external radiation fields, and the variability amplitude tends to be larger as the gamma-rays are closer to the high energy peak of spectral energy distribution. No correlation of variability index found for BL Lacertae objects implies that variability behavior may differ below and above the peak energy.

It is well known that there are only two low-frequency-peaked BL Lac objects (LBLs: BL Lacertae and S5 0716+714) and one flat spectrum radio quasar (FSRQ: 3C 279) among more than 30 active galactic nuclei (AGNs) with detected TeV emissions. We study the spectral energy distribution (SED) of a famous LBL OJ 287, whose light curve has a 12-y period. Using a homogeneous one-zone synchrotron + synchrotron-self Compton model, we model the quasi-simultaneous broad-band SED of OJ 287. With some reasonable assumptions, we extrapolate the model to the high state of OJ 287 and predict its {\gamma}-ray emissions. Taking into account the absorption of {\gamma}-ray by the extragalactic background light (EBL), we find that the TeV emission of OJ 287 in high state is slightly higher than the sensitivity of H.E.S.S. The study on SEDs of OJ 287 has implications to unveil the origin of jet activity during its 12-y period and the properties of EBL.

The Taiwanese-American Occultation Survey (TAOS) will detect objects in the Kuiper Belt, by measuring the rate of occultations of stars by these objects, using an array of three to four 50cm wide-field robotic telescopes. Thousands of stars will be monitored, resulting in hundreds of millions of photometric measurements per night. To optimize the success of TAOS, we have investigated various methods of gathering and processing the data and developed statistical methods for detecting occultations. In this paper we discuss these methods. The resulting estimated detection efficiencies will be used to guide the choice of various operational parameters determining the mode of actual observation when the telescopes come on line and begin routine observations. In particular we show how real-time detection algorithms may be constructed, taking advantage of having multiple telescopes. We also discuss a retrospective method for estimating the rate at which occultations occur.