Periodic quasars are candidates for binary supermassive black holes (BSBHs) efficiently emitting low frequency gravitational waves. Recently, $\sim$150 candidates were identified from optical synoptic surveys. However, they may be false positives caused by stochastic quasar variability given the few cycles covered (typically 1.5). To independently test the binary hypothesis, we search for evidence of truncated or gapped circumbinary accretion discs (CBDs) in their spectral energy distributions (SEDs). Our work is motivated by CBD simulations that predict flux deficits as cutoffs from central cavities opened by secondaries or notches from minidiscs around both BHs. We find that candidate periodic quasars show SEDs similar to those of control quasars matched in redshift and luminosity. While seven of 138 candidates show a blue cutoff in the IR-optical-UV SED, six of which may represent CBDs with central cavities, the red SED fraction is similar to that in control quasars, suggesting no correlation between periodicity and SED anomaly. Alternatively, dust reddening may cause red SEDs. The fraction of extremely radio-loud quasars, e.g., blazars (with $R>100$), is tentatively higher than that in control quasars (at 2.5$\sigma$). Our results suggest that, assuming most periodic candidates are robust, IR-optical-UV SEDs of CBDs are similar to those of accretion discs of single BHs, if the periodicity is driven by BSBHs; the higher blazar fraction may signal precessing radio jets. Alternatively, most current candidate periodic quasars identified from few-cycle light curves may be false positives. Their tentatively higher blazar fraction and lower Eddington ratios may both be caused by selection biases.

Without an existing large scale coherent magnetic field in the early Universe, Population III (PopIII) stars would likely rotate at or near break-up speed. In this work, focusing on the accretion phase of PopIII stars, we investigate the possibility of generating a coherent magnetic field through large scale dynamo processes, as well as the corresponding field saturation level. Using results from hydrodynamic simulations, we demonstrate that primordial accretion disks are turbulent with a Shakura-Sunyaev disk parameter $\alpha_{ss} \gtrsim 10^{-3}$, and evidence helical turbulence with a dynamo number $\vert D_{\alpha \Omega} \vert \gg 10$. The presence of helical turbulence at these levels allows large scale dynamo modes to grow, and the saturation level is determined by the amount of net helicity remaining in the dynamo-active regions (a.k.a. the quenching problem). %We demonstrate that the magnetic field can reach approximate equipartition, with $B/B_{\rm eq} \sim 3$, indicating that the dynamo quenching problem could be alleviated through an accretion flow. We demonstrate that, if the accretion could successfully alleviate the quenching problem, the magnetic field can reach approximate equipartition with $B/B_{\rm eq} \sim 3$.

Under rapid cooling from molecular hydrogen, the accretion disks around Population III (PopIII) stars are believed to fragment, resulting in multiple accreting cores. In this paper, we build a theoretical framework for calculating the optical depth of H$_2$ ro-vibrational line cooling based on the vertical structure in these accretion disks. Applying this physically motivated prescription for the optical depth, we find that cooling in the inner disk with $r \lesssim 10 {\rm\ AU}$ is attenuated significantly as a result of high surface density; $PdV$ heating becomes more efficient than cooling, which prevents fragmentation in the inner disk. Despite this, cooling becomes dynamically important in the outer disk, favoring fragmentation. We argue that most of the resultant fragments are initially at the outer disk, and that any surviving fragment has to migrate slower than the disk-scale photo-evaporation process. Since type I migration is fast, migration slows down as a result of gap-opening in the disk structure. Two possible processes for gap-opening are studied: (1) through a massive, densely-cored ($\rho \gtrsim 10^{-8} {\rm\ g\ cm^{-3}}$) clump able to radiate away the excess gravitational potential energy, and (2) through a fast-growing central star, with $\dot{M} \gtrsim 2 \times 10^{-3} \, M_\odot {\rm\ yr^{-1}}$, whose gravity dominates the star-disk system and favors gap opening.

New CCD photometric light curves of short period (P=0.285d) eclipsing binary RW Dor are presented. The observations performed with the PROMPT-8 robotic telescope at CTIO in Chile from March 2015 to March 2017. The other eclipse timings were obtained from the 2.15-m JS telescope at CASLEO, San Juan, Argentina in December 2011. By light-curve analysis, it is found that RW Dor is a W-type shallow contact binary with a fill-out factor $f \sim 11\%$ and high mass ratio $q \sim 1.587$ (1/q = 0.63), where the hotter component is the less massive one ($M_1 \sim 0.52M_{\odot}$ and $M_2 \sim 0.82M_{\odot}$). For orbital period investigation, the new fifteen eclipse times and those in previous published were compiled. Based on $O-C$ analysis with very weak evidence suggests that a long-term period decrease with a rate of $\mathrm{d}P/\mathrm{d}t = -9.61\times10^{-9}$ d $\textrm{yr}^{-1}$ is superimposed on a cyclic variation ($A_3$ = 0.0054 days and $P_3$ = 49.9 yrs). The long-term period decrease can be interpreted as mass transfer from the more massive component to the less massive one or combine with the angular momentum loss (AML) via magnetic braking. In addition, with the marginal contact phase, high mass ratio (1/q $>$ 0.4) and the long-term period decrease, all suggest that RW Dor is a newly formed contact binary via a Case A mass transfer and it will evolve into a deeper normal contact binary. If the cyclic change is correct, the light-travel time effect via the presence of a cool third body will be more plausible to explain for this.

The "gravitational million-body problem," to model the dynamical evolution of a self-gravitating, collisional N-body system with ~10^6 particles over many relaxation times, remains a major challenge in computational astrophysics. Unfortunately, current techniques to model such systems suffer from severe limitations. A direct N-body simulation with more than 10^5 particles can require months or even years to complete, while an orbit-sampling Monte Carlo approach cannot adequately model the dynamics in a dense cluster core, particularly in the presence of many black holes. We have developed a new technique combining the precision of a direct N-body integration with the speed of a Monte Carlo approach. Our Rapid And Precisely Integrated Dynamics code, the RAPID code, statistically models interactions between neighboring stars and stellar binaries while integrating directly the orbits of stars or black holes in the cluster core. This allows us to accurately simulate the dynamics of the black holes in a realistic globular cluster environment without the burdensome N^2 scaling of a full N-body integration. We compare RAPID models of idealized globular clusters to identical models from the direct N-body and Monte Carlo methods. Our tests show that RAPID can reproduce the half-mass radii, core radii, black hole ejection rates, and binary properties of the direct N-body models far more accurately than a standard Monte Carlo integration while remaining significantly faster than a full N-body integration. With this technique, it will be possible to create more realistic models of Milky Way globular clusters with sufficient rapidity to explore the full parameter space of dense stellar clusters.

New photometric observations of NSVS 01286630 were performed and two sets of four-color (B, V , Rc, Ic) light curves (LCs) were obtained. Using the 2013 version of the Wilson-Devinney (W-D) code, we analyzed these data. The photometric solutions reveal that NSVS 01286630 is an active detached eclipsing binary (EB) with a high orbit inclination. Remarkably, the temperature of the primary component (the hotter star) is higher than the secondary one, but, the value of mass ratio q(M2/M1) for NSVS 0128663 is more than 1, which can be explained that the surface of the secondary component of NSVS 01286630 is covered with big cool-star-spots. Based on our new CCD mid-eclipse times and the data published until now, the variations of the mid-eclipse times were reanalyzed in detail using a weighted least-squares method. It is discovered that the (O-C) diagram of the system shows a cyclic oscillation with a period of 3.61 yr and an amplitude of 0.001 days. The cyclic variation may be caused by the light travel time effect (LTTE) due to the presence of a third companion, whose mass we calculated is about 0.11 solar mass. The third body may affect the orbital evolution of the central binary system by transferring the angular momentum.

1SWASP\,J162117.36$+$441254.2 was originally classified as an EW-type binary with a period of 0.20785\,days. However, it was detected to have undergone a stellar outburst on June 3, 2016. Although the system was latter classified as a cataclysmic variable (CV) and the event was attributed as a dwarf-nova outburst, the physical reason is still unknown. This binary has been monitored photometrically since April 19, 2016 and many light curves were obtained before, during and after the outburst. Those light and color curves observed before the outburst indicate that the system is a special CV. The white dwarf is not accreting material from the secondary and there are no accretion disks surrounding the white dwarf. By comparing the light curves obtained from April 19 to September 14, 2016, it was found that magnetic activity of the secondary is associated with the outburst. We show strong evidence that the $L_1$ region on the secondary was heavily spotted before and after the outburst and thus quench the mass transfer, while the outburst is produced by a sudden mass accretion of the white dwarf. These results suggest that J162117 is a good astrophysical laboratory to study stellar magnetic activity and its influences on CV mass transfer and mass accretion.

We find that a fuzzy dark matter background and the mG scale magnetic field in the galactic center can give rise to a radiation with a very large energy release. The frequency of the radiation field is the same as the frequency of the oscillating axion-like background field. We show that there is an energy transfer between the fuzzy dark matter sector and the electromagnetic sector because of the presence of the generated radiation field and the galactic magnetic field. The energy release rate of radiation is found to be very slow in comparison with the energy of fuzzy dark matter but could be significant comparing with the energy of galactic magnetic field in the source region. Using this example, we show that the fuzzy dark matter together with a large scale magnetic field is possible to give rise to fruitful physics.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

We study Earth matter effect in oscillation of supernovae neutrinos. We show that detecting Earth matter effect gives an independent measurement of spectra of supernovae neutrinos, i.e. the flavor difference of the spectra of supernovae neutrinos. We study the effect of energy resolution and angular resolution of final electron or positron on detecting the signal of Earth matter effect. We show that varying the widths of energy bins in analysis can change the signal strength of Earth matter effect and the statistical fluctuation. A reasonable choice of energy bins can both suppress the statistical fluctuation and make out a good signal strength relative to the statistical fluctuation. Neutrino detectors with good energy resolution and good angular resolution are therefore preferred so that there are more freedom to vary energy bins and to optimize the signal of Earth matter effect in analyzing events of supernovae neutrinos.

The first V Rc Ic bands light curves of MQ UMa are presented and analyzed using the Wilson-Devinney (W-D) program. It is discovered that MQ UMa is an A-subtype contact binary with a high fill-out (f=82 %) and a low mass ratio (q = 0:195), which indicates that it is in the late evolutionary stage of latetype tidal-locked binary stars. The mass of the primary and secondary stars are estimated and the evolutionary status of the two components are placed in the H-R diagram. The W-D solutions also indicates that MQ UMa may be a triple system with an F5V type additional component. A sample of sixteen high fillout, low mass ratio overcontact binaries are collected and their possible evolution scenarios are discussed. Using the five times of light minimum recently observed together with those collected from literatures, the authors find that the observedcalculate (O-C) curve exhibits a cyclic period variation. The cyclic period change also reveals the presence of a tertiary component, which may play an important role in the formation and evolution of this binary system by drawing angular momentum from the central system.

The leading theory for the formation of the Earth's moon invokes a collision between a Mars-sized body and the proto-Earth to produce a disk of orbiting material that later condenses to form the Moon. Here we study the early evolution of the protolunar disk. First, we show that the disk opacity is large and cooling is therefore inefficient (t_{cool} \Omega >> 1). In this regime angular momentum transport in the disk leads to steady heating unless \alpha < (t_{cool} \Omega)^{-1} << 1. Following earlier work by Charnoz and Michaut, and Carballido et al., we show that once the disk is completely vaporized it is well coupled to the magnetic field. We consider a scenario in which turbulence driven by magnetic fields leads to a brief, hot phase where the disk is geometrically thick, wit h strong turbulent mixing. The disk cools by spreading until it decouples from the field. We point out that approximately half the accretion energy is dissipated in the boundary layer where the disk meets the Earth's surface. This creates high entropy material close to the Earth, driving convection and mixing. Finally, a hot, magnetized disk could drive bipolar outflows that remove mass and angular momentum from the Earth-Moon system.

Current and future surveys of large-scale cosmic structure are associated with a massive and complex datastream to study, characterize, and ultimately understand the physics behind the two major components of the 'Dark Universe', dark energy and dark matter. In addition, the surveys also probe primordial perturbations and carry out fundamental measurements, such as determining the sum of neutrino masses. Large-scale simulations of structure formation in the Universe play a critical role in the interpretation of the data and extraction of the physics of interest. Just as survey instruments continue to grow in size and complexity, so do the supercomputers that enable these simulations. Here we report on HACC (Hardware/Hybrid Accelerated Cosmology Code), a recently developed and evolving cosmology N-body code framework, designed to run efficiently on diverse computing architectures and to scale to millions of cores and beyond. HACC can run on all current supercomputer architectures and supports a variety of programming models and algorithms. It has been demonstrated at scale on Cell- and GPU-accelerated systems, standard multi-core node clusters, and Blue Gene systems. HACC's design allows for ease of portability, and at the same time, high levels of sustained performance on the fastest supercomputers available. We present a description of the design philosophy of HACC, the underlying algorithms and code structure, and outline implementation details for several specific architectures. We show selected accuracy and performance results from some of the largest high resolution cosmological simulations so far performed, including benchmarks evolving more than 3.6 trillion particles.

New multi-color light curves of the very short period K-type eclipsing binary V1799 Ori were obtained and analyzed with the W-D code. The photometric solutions reveal that the system is a W-type shallow-contact binary with a mass ratio of $q=1.335(\pm0.005)$ and a degree of contact about $f = 3.5(\pm1.1)\%$. In general, the results are in good agreement with which is reported by Samec. The remarkable O'Connell effects in the light curves are well explained by employing star spots on the binary surface, which confirms that the system is active at present. Several new times of light minimum were obtained. All the available times of light minimum were collected, along with the recalculated and new obtained. Applying a least-squares method to the constructed O-C diagram, a new ephemeris was derived for V1799 Ori. The orbital period is found to show a continuous weak increase at a rate of $1.8(\pm0.6)\times10^{-8}$ days$\cdot$yr$^{-1}$. The extreme-shallow contact, together with the period increase, suggests that the binary may be at a critical stage predicted by the TRO theory. \keywords{binaries : close -- binaries : eclipsing -- stars: individual (V1799 Ori)

We present a new parallel code for computing the dynamical evolution of collisional N-body systems with up to N~10^7 particles. Our code is based on the the Henon Monte Carlo method for solving the Fokker-Planck equation, and makes assumptions of spherical symmetry and dynamical equilibrium. The principal algorithmic developments involve optimizing data structures, and the introduction of a parallel random number generation scheme, as well as a parallel sorting algorithm, required to find nearest neighbors for interactions and to compute the gravitational potential. The new algorithms we introduce along with our choice of decomposition scheme minimize communication costs and ensure optimal distribution of data and workload among the processing units. The implementation uses the Message Passing Interface (MPI) library for communication, which makes it portable to many different supercomputing architectures. We validate the code by calculating the evolution of clusters with initial Plummer distribution functions up to core collapse with the number of stars, N, spanning three orders of magnitude, from 10^5 to 10^7. We find that our results are in good agreement with self-similar core-collapse solutions, and the core collapse times generally agree with expectations from the literature. Also, we observe good total energy conservation, within less than 0.04% throughout all simulations. We analyze the performance of the code, and demonstrate near-linear scaling of the runtime with the number of processors up to 64 processors for N=10^5, 128 for N=10^6 and 256 for N=10^7. The runtime reaches a saturation with the addition of more processors beyond these limits which is a characteristic of the parallel sorting algorithm. The resulting maximum speedups we achieve are approximately 60x, 100x, and 220x, respectively.

We report here the discovery of an optical flare observed in R band from the red-dwarf eclipsing binary CU Cnc whose component stars are at the upper boundary of full convection (M1=0.43 and M2=0.4M0, M0 is the solar mass). The amplitude of the flare is the largest among those detected in R band (~0.52mag) and the duration time is about 73 minutes. As those observed on the Sun, quasi-periodic oscillations were seen during and after the flare. Three more R-band flares were found by follow up monitoring. In total, this binary was monitored photometrically by using R filter for 79.9 hours, which reveals a R-band flare rate about 0.05 flares per hour. These detections together with other strong chromospheric and coronal activities, i.e., very strong H_alpha and H_beta emission features and an EUV and X-ray source, indicate that it has very strong magnetic activity. Therefore, the apparent faintness (~1.4 magnitude in V) of CU Cnc compared with other single red dwarfs of the same mass can be plausibly explained by the high coverage of the dark spots.

Using the precise times of mid-egress of the eclipsing polar HU Aqr, we discovered that this polar is orbited by two or more giant planets. The two planets detected so far have masses of at least 5.9 and 4.5\,M_{Jup}. Their respective distances from the polar are 3.6 AU and 5.4 AU with periods of 6.54 and 11.96 years, respectively. The observed rate of period decrease derived from the downward parabolic change in O-C curve is a factor 15 larger than the value expected for gravitational radiation. This indicates that it may be only a part of a long-period cyclic variation, revealing the presence of one more planet. It is interesting to note that the two detected circumbinary planets follow the Titus-Bode law of solar planets with n=5 and 6. We estimate that another 10 years of observations will reveal the presence of the predicted third planet.

We searched the orbital period changes in 182 EA-type (including the 101 Algol systems used by \cite{hal89}), 43 EB-type and 53 EW-type binaries with known both the mass ratio and the spectral type of their secondary components. We reproduced and improved the same diagram as Hall's (1989) according to the new collected data. Our plots do not support the conclusion derived by \cite{hal89} that all cases of cyclical period changes are restricted to binaries having the secondary component with spectral types later than F5. The presence of period changes also among stars with secondary component of early type indicates that the magnetic activity is one cause, but not the only one, for the period variation. It is discovered that cyclic period changes, likely due to the presence of a third body are more frequent in EW-type binaries among close binaries. Therefore, the most plausible explanation of the cyclical period changes is the LTTE via the presence of a third body. By using the century-long historical record of the times of light minimum, we analyzed the cyclical period change in the Algol binary WW Dra. It is found that the orbital period of the binary shows a $\sim112.2 \textbf{\textrm{yr}}$ cyclic variation with an amplitude of $\sim0.1977\textbf{\textrm{days}}$. The cyclic oscillation can be attributed to the LTTE via a third body with a mass no less than $6.43 M_{\odot}$. However, no spectral lines of the third body were discovered indicating that it may be a candidate black hole. The third body is orbiting the binary at a distance shorter than 14.4 AU and it may play an important role in the evolution of this system.

We derive a series of moment equations describing the motion and flavor transformation of neutrinos in supernova. We find a particular series of moments of neutrino density matrix in supernova. The emission angle distribution of neutrinos is described by this series of moments. We expand the equation of neutrinos using these moments and obtain moment equations. We find that these moments have very good property of convergence and the infinite series of equations can be truncated to equations with a small set of moments. Using a small set of moment equations the required computational power is reduced by about two orders of magnitude compared to that in multi-angle simulation. The study on non-linear flavor transformation of neutrinos is substantially simplified using these equations. Two flavor system of neutrinos is also considered and new equations describing the flavor polarization vectors of neutrinos are found.

HS 0705+6700 is a short-period (P=2.3 hours), close binary containing a hot sdB-type primary and a fully convective secondary. We have monitored this eclipsing binary for more than 2 years and as a result, 32 times of light minimum were obtained. Based on our new eclipse times together with these compiled from the literature, it is discovered that the O-C curve of HS 0705+6700 shows a cyclic variation with a period of 7.15 years and a semiamplitude of 92.4 s. The periodic change was analyzed for the light-travel time effect that may be due to the presence of a tertiary companion. The mass of the third body is determined to be M3 sin i = 0.0377 (+/-0.0043) Msun when a total mass of 0.617 Msun for HS 0705+6700 is adopted. For orbital inclinations i >= 32.8, the mass of the tertiary component would be below the stable hydrogen-burning limit of M3~0.072 Msun, and thus it would be a brown dwarf. The third body is orbiting the sdB-type binary at a distance shorter than 3.6 astronomical units (AU). HS 0705+6700 was formed through the evolution of a common envelope after the primary becomes a red giant. The detection of a sub-stellar companion in HS 0705+6700 system at this distance from the binary could give some constraints on stellar evolution in such systems and the interactions between red giants and their companions.

Up to now, most stellar-mass black holes were discovered in X-ray emitting binaries, in which the black holes are formed through a common-envelope evolu tion. Here we give evidence for the presence of a massive black hole candidate as a tertiary companion in the massive eclipsing binary V Puppis. We found that the orbital period of this short-period binary (P=1.45 days) shows a periodic variation while it undergoes a long-term increase. The cyclic period oscillation can be interpreted by the light-travel time effect via the presence of a third body with a mass no less than 10.4 solar mass. However, no spectral lines of the third body were discovered indicating that it is a massive black hole candidate. The black hole candidate may correspond to the weak X-ray source close to V Puppis discovered by Uhuru, Copernicus, and ROSAT satellites produced by accreting materials from the massive binary via stellar wind. The circumstellar matter with many heavy elements around this binary may be formed by the supernova explosion of the progenitor of the massive black hole. All of the observations suggest that a massive black hole is orbiting the massive close binary V Puppis with a period of 5.47 years. Meanwhile, we found the central close binary is undergoing slow mass transfer from the secondary to the primary star on a nuclear time scale of the secondary component, revealing that the system has passed through a rapid mass-transfer stage.

Solar neutrino physics enters a stage of precision measurements. In this connection we present a precise analytic description of the neutrino conversion in the context of LMA MSW solution of the solar neutrino problem. Using the adiabatic perturbation theory we derive an analytic formula for the $\nu_e$ survival probability which takes into account the non-adiabatic corrections and the regeneration effect inside the Earth. The probability is averaged over the neutrino production region. We find that the non-adiabatic corrections are of the order $10^{-9}-10^{-7}$. Using the formula for the Earth regeneration effect we discuss features of the zenith angle dependence of the $\nu_e$ flux. In particular, we show that effects of small structures at the surface of the Earth can be important.