We investigate a cogenesis scenario within the twin Higgs setup which can naturally explain the nature of dark matter, the cosmic coincidence puzzle, little hierarchy problem, leptogenesis and the tiny neutrino masses. Three heavy Majorana neutrinos are introduced to the standard model sector and the twin sector respectively, which explain the tiny neutrino masses and generate the lepton asymmetry and the twin lepton asymmetry at the same time. The twin cogenesis scenario is general and applies to any viable twin Higgs model without hard $\mathbb{Z}_2$ breaking and evading the $\Delta N_{\rm eff}$ constraint. We demonstrate twin cogenesis in two models: fraternal twin Higgs model, and neutrino-philic twin two Higgs doublet model, a newly proposed model to lift the twin neutrino masses with spontaneous $\mathbb{Z}_2$ breaking. The MeV scale dark photon ensures the energy in the twin sector as well as the symmetric component of twin sector particles can be depleted. The lightest twin baryons are the dark matter candidates with masses approximately 5.5~GeV, which explain naturally the amount of dark matter and visible matter in the Universe are of the same order.

Over 54 years of hourly mean value of solar wind velocity from 27 Nov. 1963 to 31 Dec. 2017 are used to investigate characteristics of the rotation period of solar wind through auto-correlation analysis. Solar wind of high velocity is found to rotate faster than low-velocity wind, while its rotation rate increases with velocity increasing, but in contrast for solar wind of low velocity, its rotation rate decreases with velocity increasing. Our analysis shows that solar wind of a higher velocity statistically possesses a faster rotation rate for the entire solar wind. The yearly rotation rate of solar wind velocity does not follow the Schwable cycle, but it is significantly negatively correlated to yearly sunspot number when it leads by 3 years. Physical explanations are proposed to these findings.

We investigate the prospects of detecting gravitational waves from coalescing massive black hole binaries in the Universe with the TianQin observatory, a space-based gravitational wave interferometer proposed to be launched in the 2030s. To frame the scientific scope of the mission, in this paper, we carry out a preliminary estimation of the signal-to-noise ratio, detection rate and parameter estimation precision of massive black hole binaries detectable by TianQin. In order to make our results as robust as possible, we consider several models of the growth history of massive black holes, exploring the effect of some key astrophysical prescriptions as well the impact of the employed computational methods. In the most optimistic model, TianQin can detect as many as approximately 60 mergers per year. If TianQin detects a merger at redshift of 15, it will be capable of estimating its luminosity distance to within an accuracy of 10%; for a nearby event at redshift approximately 2, TianQin can issue early warnings 24 hours before coalescence, with a timing accuracy of around three hours and a sky localization ability of approximately 80 deg$^2$, thus enabling multi-messenger observations.

In response to the need for the Astro2020 Decadal Survey to explicitly engage early career astronomers, the National Academies of Sciences, Engineering, and Medicine hosted the Early Career Astronomer and Astrophysicist Focus Session (ECFS) on October 8-9, 2018 under the auspices of Committee of Astronomy and Astrophysics. The meeting was attended by fifty six pre-tenure faculty, research scientists, postdoctoral scholars, and senior graduate students, as well as eight former decadal survey committee members, who acted as facilitators. The event was designed to educate early career astronomers about the decadal survey process, to solicit their feedback on the role that early career astronomers should play in Astro2020, and to provide a forum for the discussion of a wide range of topics regarding the astrophysics career path. This white paper presents highlights and themes that emerged during two days of discussion. In Section 1, we discuss concerns that emerged regarding the coming decade and the astrophysics career path, as well as specific recommendations from participants regarding how to address them. We have organized these concerns and suggestions into five broad themes. These include (sequentially): (1) adequately training astronomers in the statistical and computational techniques necessary in an era of "big data", (2) responses to the growth of collaborations and telescopes, (3) concerns about the adequacy of graduate and postdoctoral training, (4) the need for improvements in equity and inclusion in astronomy, and (5) smoothing and facilitating transitions between early career stages. Section 2 is focused on ideas regarding the decadal survey itself, including: incorporating early career voices, ensuring diverse input from a variety of stakeholders, and successfully and broadly disseminating the results of the survey.

We use the Fisher information matrix method to calculate the parameter estimation accuracy of inspiraling supermassive black holes binaries for TianQin, a space-borne laser interferometric detector aimed at detecting gravitational waves in the millihertz frequency band. The `restricted' post-Newtonian waveform in which third order post-Newtonian (3PN) phase including spin effects (spin-orbit $\beta$ and spin-spin $\sigma$) and first-order eccentricity contribution is employed. Monte Carlo simulations using $10^3$ binaries for mass pairs with component masses in the range of $({10^5},{10^7}){M_ \odot }$ and cosmological redshift $z=0.5$ show that the medians of the root-mean-square error distributions for the chirp mass $M_c$ and symmetric mass ratio $\eta$ are in the range of $\sim 0.02\% - 0.7\% $ and $\sim 4\% - 8\% $, respectively. The luminosity distance $D_L$ can be determined to be $\sim 1\% - 3\% $, and the angular resolution of source $\Delta \Omega $ is better than 12 deg$^2$. The corresponding results for $z=1.0$ and $2.0$, which are deteriorated with the decreasing of the signal-to-noise ratio, have also been given. We show that adding spin parameters degrades measurement accuracy of the mass parameters (${M_c}$, $\eta$), and the time and the orbital phase of coalescence ($t_c$, $\phi _c$); the inclusion of the first-order eccentricity correction to the phase worsens the estimation accuracy comparing with the circular cases. We also show the effects of post-Newtonian order on parameter estimation accuracy by comparing the results based on second order and third order post-Newtonian phases. Moreover, we calculate the horizon distance of supermassive black hole binaries for TianQin.

Coronal heating is a big question for modern astronomy. Daily measurement of 985 solar spectral irradiances (SSIs) at the spectral intervals 1-39 nm and 116-2416 nm during March 1 2003 to October 28 2017 is utilized to investigate characteristics of solar rotation in the solar atmosphere by means of the Lomb \,-\, Scargle periodogram method to calculate their power spectra. The rotation period of coronal plasma is obtained to be 26.3 days, and that of the solar atmosphere at the bottom of the photosphere modulated by magnetic structures is 27.5 days. Here we report for the first time that unexpectedly the coronal atmosphere is found to rotate faster than the underlying photosphere. When time series of SSIs are divided into different cycles, and the ascending and descending periods of a solar cycle, rotation rate in the corona is also found to be larger than that in the photosphere, and this actually gives hidden evidence: it is small-scale magnetic activity that heats the corona.

Solar chromosphere and coronal heating is a big question for astrophysics. Daily measurement of 985 solar spectral irradiances (SSIs) at the spectral intervals 1-39 nm and 116-2416 nm during March 1 2003 to October 28 2017 is utilized to investigate phase relation respectively with daily sunspot number, the Mount Wilson Sunspot Index, and the Magnetic Plage Strength Index. All SSIs which form in the whole heated region: the upper photosphere, chromosphere, transition region, and corona are found to be significantly more correlated to weak magnetic activity than to strong magnetic activity, and to dance in step with weak magnetic activity. All SSIs which form in the low photosphere (the unheated region), which indicate the "energy" leaked from the solar subsurface are found to be more related to strong magnetic activity instead and in anti-phase with weak magnetic activity. In the upper photosphere and chromosphere, strong magnetic activity should lead SSI by about a solar rotation, also displaying that weak magnetic activity should take effect on heating there. It is thus small-scale weak magnetic activity that effectively heats the upper solar atmosphere.

TianQin is a space-based laser interferometric gravitational wave detector aimed at detecting gravitational waves at low frequencies (0.1 mHz -- 1 Hz). It is formed by three identical drag-free spacecrafts in an equilateral triangular constellation orbiting around the Earth. The distance between each pair of spacecrafts is approximately $1.7 \times 10^{5} ~\rm{km}$. The spacecrafts are interconnected by infrared laser beams forming up to three Michelson-type interferometers. The detailed mission design and the study of science objectives for the TianQin project depend crucially on the orbit and the response of the detector. In this paper, we provide the analytic expressions for the coordinates of the orbit for each spacecraft in the heliocentric-ecliptic coordinate system to the leading orders. This enables a sufficiently accurate study of science objectives and data analysis, and serves as a first step to further orbit design and optimization. We calculate the response of a single Michelson detector to plane gravitational waves in arbitrary waveform which is valid in the full range of the sensitive frequencies. It is then used to generate the more realistic sensitivity curve of TianQin. We apply this model on a reference white-dwarf binary as a proof of principle.

Using Keck/HIRES spectra {\Delta}v ~ 7 km/s, we analyze forbidden lines of [O I] 6300 {\AA}, [O I] 5577 {\AA} and [S II] 6731 {\AA} from 33 T Tauri stars covering a range of disk evolutionary stages. After removing a high velocity component (HVC) associated with microjets, we study the properties of the low velocity component (LVC). The LVC can be attributed to slow disk winds that could be magnetically (MHD) or thermally (photoevaporative) driven. Both of these winds play an important role in the evolution and dispersal of protoplanetary material. LVC emission is seen in all 30 stars with detected [O I] but only in 2 out of eight with detected [S II] , so our analysis is largely based on the properties of the [O I] LVC. The LVC itself is resolved into broad (BC) and narrow (NC) kinematic components. Both components are found over a wide range of accretion rates and their luminosity is correlated with the accretion luminosity, but the NC is proportionately stronger than the BC in transition disks. The FWHM of both the BC and NC correlates with disk inclination, consistent with Keplerian broadening from radii of 0.05 to 0.5 AU and 0.5 to 5 AU, respectively. The velocity centroids of the BC suggest formation in an MHD disk wind, with the largest blueshifts found in sources with closer to face-on orientations. The velocity centroids of the NC however, show no dependence on disk inclination. The origin of this component is less clear and the evidence for photoevaporation is not conclusive.

G-band bright points (GBPs) are thought to be the foot-points of magnetic flux tubes. The aim of this paper is to investigate the relation between the diffusion regimes of GBPs and the associated longitudinal magnetic field strengths. Two high resolution observations of different magnetized environments were acquired with the Hinode/Solar Optical Telescope. Each observation was recorded simultaneously with G-band filtergrams and Narrow-band Filter Imager (NFI) Stokes I and V images. GBPs are identified and tracked automatically, and then categorized into several groups by their longitudinal magnetic field strengths, which are extracted from the calibrated NFI magnetograms using a point-by-point method. The Lagrangian approach and the distribution of diffusion indices approach are adopted separately to explore the diffusion regime of GBPs for each group. It is found that the values of diffusion index and diffusion coefficient both decrease exponentially with the increasing longitudinal magnetic field strengths whichever approach is used. The empirical formulas deduced from the fitting equations are proposed to describe these relations. Stronger elements tend to diffuse more slowly than weak elements, independently of the magnetic flux of the surrounding medium. This may be because the magnetic energy of stronger elements is not negligible compared with the kinetic energy of the gas, and therefore the flows cannot perturb them so easily.Yang

We propose a framework in which visible matter interacts with matter from a hidden sector through mass mixings of St\"uckelberg U(1) gauge fields. In contrast to other Z' mediation scenarios, our setup has the added appealing features that (i) the choice of Z's can be significantly broadened without necessarily introducing unwanted exotic matter and (ii) there can be sizable tree-level interactions between the visible and hidden sectors. String theory embeddings of this scenario and their phenomenological features are briefly discussed.

We consider a gauged $U(1)_x$ extension of the standard model and of the minimal supersymmetric standard model where the dark matter fields are charged under $U(1)_x$ and carry lepton number while the standard model fields and fields of the minimal supersymmetric standard model are neutral under $U(1)_x$. We consider leptogenesis in this class of models with all fundamental interactions having no violation of lepton number, and the total $B-L$ in the universe vanishes. Such leptogenesis leads to equal and opposite lepton numbers in the visible sector and in the dark sector, and thus also produces asymmetric dark matter. Part of the lepton numbers generated in the leptonic sector subsequently transfer to the baryonic sector via sphaleron interactions. The stability of the dark particles is protected by the $U(1)_x$ gauge symmetry. A kinetic mixing between the $U(1)_x$ and the $U(1)_Y$ gauge bosons allows for dissipation of the symmetric component of dark matter. The case when $U(1)_x$ is $U(1)_{B-L}$ is also discussed for the supersymmetric case. This case is particularly interesting in that we have a gauged $U(1)_{B-L}$ which ensures the conservation of $B-L$ with an initial condition of a vanishing $B-L$ in the universe. Phenomenological implications of the proposed extensions are discussed, which include implications for electroweak physics, neutrino masses and mixings, and lepton flavor changing processes such as $\ell_i \to \ell_j \gamma$. We also briefly discuss the direct detection of the dark matter in the model.

We consider the possibility that some primordial fields decay purely into the dark sector creating asymmetric dark matter. This asymmetry is subsequently transmuted into leptons and baryons. Within this paradigm we compute the amount of asymmetric dark matter created from the out of equilibrium decays of the primordial fields with CP violating Yukawa couplings. The dark matter asymmetry is then transferred to the visible sector by the asymmetry transfer equation and generates an excess of $B-L$. Baryogenesis occurs via sphaleron processes which conserve $B-L$ but violate $B+L$. A mechanism for the annihilation of the symmetric component of dark matter is also discussed. The model leads to multi-component dark matter consisting of both bosonic and fermionic components.

In accreting young stars one of the prominent spectral features in the near infrared is the Paschen and Brackett series in emission. We examine hydrogen line ratios for 16 classical T Tauri stars from SpeX spectra and and assess the trends with veiling and accretion. The observed line ratios are compared to two theoretical models for line formation: (1) Baker and Menzel's (1938) Case B for radiative ionization and recombination and (2) a set of local line excitation calculations designed to replicate the conditions in T Tauri winds and magnetic accretion columns (Kwan & Fischer 2011). While the comparison between Case B and observed line ratios implies a wide range in electron density and temperature among the hydrogen line formation regions in T Tauri stars, the predictions of the local line excitation models give consistent results across multiple diagnostics. Under the assumptions of the local line excitation calculations, we find that n_H in the hydrogen line formation region is constrained to 2x10^10 -- 2x10^11 cm^-3, where stars with higher accretion rates have densities at the higher end of this range. Because of uncertainties in extinction, temperature is not well delineated but falls within the range expected for collisional excitation to produce the line photons. We introduce new diagnostics for assessing extinction based on near infrared hydrogen line ratios from the local line excitation calculations.

Solar-cycle related variation of differential rotation is investigated through analyzing the rotation rates of magnetic fields, distributed along latitudes and varying with time at the time interval of August 1976 to April 2008. More pronounced differentiation of rotation rates is found to appear at the ascending part of a Schwabe cycle than at the descending part on an average. The coefficient $B$ in the standard form of differential rotation, which represents the latitudinal gradient of rotation, may be divided into three parts within a Schwabe cycle. Part one spans from the start to the $4^{th}$ year of a Schwabe cycle, within which the absolute $B$ is approximately a constant or slightly fluctuates. Part two spans from the $4^{th}$ to the $7^{th}$ year, within which the absolute $B$ decreases. Part three spans from the $7^{th}$ year to the end, within which the absolute $B$ increases. Strong magnetic fields repress differentiation of rotation rates, so that rotation rates show less pronounced differentiation, but weak magnetic fields seem to just reflect differentiation of rotation rates. The solar-cycle related variation of solar differential rotation is inferred to the result of both the latitudinal migration of the surface torsional pattern and the repression of strong magnetic activity to differentiation of rotation rates.

We discuss the possibility of cogenesis generating the ratio of baryon asymmetry to dark matter in a Stueckelberg U(1) extension of the standard model and of the minimal supersymmetric standard model. For the U(1) we choose $L_{\mu}-L_{\tau}$ which is anomaly free and can be gauged. The dark matter candidate arising from this extension is a singlet of the standard model gauge group but is charged under $L_{\mu}-L_{\tau}$. Solutions to the Boltzmann equations for relics in the presence of asymmetric dark matter are discussed. It is shown that the ratio of the baryon asymmetry to dark matter consistent with the current WMAP data, i.e., the cosmic coincidence, can be successfully explained in this model with the depletion of the symmetric component of dark matter from resonant annihilation via the Stueckelberg gauge boson. For the extended MSSM model it is shown that one has a two component dark matter picture with asymmetric dark matter being the dominant component and the neutralino being the subdominant component (i.e., with relic density a small fraction of the WMAP cold dark matter value). Remarkably, the subdominant component can be detected in direct detection experiments such as SuperCDMS and XENON-100. Further, it is shown that the class of Stueckelberg models with a gauged $L_{\mu}-L_{\tau}$ will produce a dramatic signature at a muon collider with the $\sigma(\mu^+\mu^-\to \mu^+\mu^-,\tau^+\tau^-)$ showing a detectable $Z'$ resonance while $\sigma(\mu^+\mu^-\to e^+e^-)$ is devoid of this resonance. Asymmetric dark matter arising from a $U(1)_{B-L}$ Stueckelberg extension is also briefly discussed. Finally, in the models we propose the asymmetric dark matter does not oscillate and there is no danger of it being washed out from oscillations.

Long-term variations of solar differential rotation and sunspot activity are investigated through re-analyzing the data on parameters of the differential rotation law obtained by Makarov, Tlatov, and Callebaut (1997), Javaraiah, Bertello, and Ulrich (2005a, b), and Javaraiah et al. (2009). Our results indicate that the solar surface rotation rate at the Equator (indicated by the A parameter of the standard solar rotation law) shows a secular decrease since cycle 12 onwards, given by about $1\,-\,1.5\times10^{-3}$($deg\ day^{-1} year^{-1}$). The B parameter of the standard differential rotation law seems to also show a secular decrease since cycle 12 onwards, but of weak statistical significance. The rotation rate averaged on latitudes ($0^{o}\,--\,40^{o}$) does not show a secular trend of statistical significance. Moreover, the average sunspot area shows a secular increase of statistical significance since cycle 12 onwards, while a negative correlation is found between the level of sunspot activity (indicated by the average sunspot area) and the solar equatorial rotation in the long run.

We have developed a method to equip homodyne interferometers with the capability to operate with constant high sensitivity over many fringes for continuous real-time tracking. The method can be considered as an extension of the "J_1...J_4" methods, and its enhancement to deliver very sensitive angular measurements through Differential Wavefront Sensing is straightforward. Beam generation requires a sinusoidal phase modulation of several radians in one interferometer arm. On a stable optical bench, we have demonstrated a long-term sensitivity over thousands of seconds of 0.1 mrad/sqrt[Hz] that correspond to 20 pm/sqrt[Hz] in length, and 10 nrad/sqrt[Hz] in angle at millihertz frequencies.

In order to probe the mechanism of variations of the Solar Constant on the inter-solar-cycle scale, total solar irradiance (TSI, the so-called Solar Constant) in the time interval of 7 November 1978 to 20 September 2010 is decomposed into three components through the empirical mode decomposition and time-frequency analyses. The first component is the rotation signal, counting up to 42.31% of the total variation of TSI, which is understood to be mainly caused by large magnetic structures, including sunspot groups. The second is an annual-variation signal, counting up to 15.17% of the total variation, the origin of which is not known at this point in time. Finally, the third is the inter-solar-cycle signal, counting up to 42.52%, which are inferred to be caused by the network magnetic elements in quiet regions, whose magnetic flux ranges from $(4.27-38.01)\times10^{19}$ Mx.

Sunspot groups observed by Royal Greenwich Observatory/US Air Force/NOAA from May 1874 to November 2008 and the Carte Synoptique solar filaments from March 1919 to December 1989 are used to investigate the relative phase shift of the paired wings of butterfly diagrams of sunspot and filament activities. Latitudinal migration of sunspot groups (or filaments) does asynchronously occur in the northern and southern hemispheres, and there is a relative phase shift between the paired wings of their butterfly diagrams in a cycle, making the paired wings spatially asymmetrical on the solar equator. It is inferred that hemispherical solar activity strength should evolve in a similar way within the paired wings of a butterfly diagram in a cycle, making the paired wings just and only keep the phase relationship between the northern and southern hemispherical solar activity strengths, but a relative phase shift between the paired wings of a butterfly diagram should bring about an almost same relative phase shift of hemispheric solar activity strength.

We present results on the compact steep-spectrum quasar 3C 48 from observations with the VLBA, MERLIN and EVN at multiple radio frequencies. In the 1.5-GHz VLBI images, the radio jet is characterized by a series of bright knots. The active nucleus is embedded in the southernmost VLBI component A, which is further resolved into two sub-components A1 and A2 at 4.8 and 8.3 GHz. A1 shows a flat spectrum and A2 shows a steep spectrum. The most strongly polarized VLBI components are located at component C $\sim$0.25 arcsec north of the core. The polarization angles at C show gradual changes across the jet width at all observed frequencies, indicative of a gradient in the emission-weighted intrinsic polarization angle across the jet and possibly a systematic gradient in the rotation measure; moreover, the percentage of polarization increases near the curvature at C, likely consistent with the presence of a local jet-ISM interaction and/or changing magnetic-field directions. The hot spot B shows a higher rotation measure, and has no detected proper motion. These facts provide some evidence for a stationary shock in the vicinity of B. Comparison of the present VLBI observations with those made 8.43 years ago suggests a proper motion of $\beta_{app}=3.7\pm0.4 c$ for A2 to the north. The apparent superluminal motion suggests that the relativistic jet plasma moves at a velocity of $\gtrsim0.96 c$ if the jet is viewed at an inclination angle less than $20\degr$. A simple precessing jet model and a hydrodynamical isothermal jet model with helical-mode Kelvin-Helmholtz instabilities are used to fit the oscillatory jet trajectory of 3C 48 defined by the bright knots.

We present multi--frequency radio observational results of the quasar 3C~48. The observations were carried out with the Very Large Array (VLA) at five frequencies of 0.33, 1.5, 4.8, 8.4, and 22.5 GHz, and with the Multi--Element Radio Linked Interferometer Network (MERLIN) at the two frequencies of 1.6 and 5 GHz. The source shows a one--sided jet to the north within 1\arcsec, which then extends to the northeast and becomes diffuse. Two bright components (N2 and N3), containing most of the flux density are present in the northern jet. The spectral index of the two components is $\alpha_{N2}\sim-0.99\pm0.12$ and $\alpha_{N3}\sim-0.84\pm0.23$ ($S\propto\nu^{\alpha}$). Our images show the presence of an extended structure surrounding component N2, suggestive of strong interaction between the jet and the interstellar medium (ISM) of the host galaxy. A steep--spectrum component, labelled as S, located 0.25 ardsec southwest to the flat--spectrum component which could be the core of 3C 48, is detected at a significance of $>15\sigma$. Both the location and the steepness of the spectrum of component S suggest the presence of a counter--jet in 3C 48.

We present images of NRAO530 observed with the EVN (VLBI) at 5 GHz, the MERLIN at 1 .6 and 5 GHz, and the VLA at 5 and 8 GHz showing the complex morphology on scales from pc to kpc. The VLBI image shows a core-jet structure indicating a somehow oscillation trajectory on a scale of 30 mas, north to the strongest compact component (core). A core-jet structure extended to several hundreds mas at about P.A. -50 deg and a distant component located 11 arcsec west to the core are detected in both the MERLIN and the VLA observations. An arched structure of significant emission between the core and the distant component is also revealed in both the MERLIN image at 1.6 cm and the VLA images at 8.4 and 5 GHz. The core component shows a flat spectrum with alpha = -0.02 (S proportional to the frequency power -alpha) while alpha = 0.8 for the distant component. The steep spectrum of the distant component and the detection of the arched emission suggests that the western distant component is a lobe or a hot-spot powered by the nucleus of NRAO530. A patch of diffuse emission, 12 arcsec nearly east (P.A. = 70 deg) to the core component, is also observed with the VLA at 5 GHz, suggesting a presence of a counter lobe in the source.