ASASSN-18am/SN 2018gk is a newly discovered member of the rare group of luminous, hydrogen-rich supernovae (SNe) with a peak absolute magnitude of $M_V \approx -20$ mag that is in between normal core-collapse SNe and superluminous SNe. These SNe show no prominent spectroscopic signatures of ejecta interacting with circumstellar material (CSM), and their powering mechanism is debated. ASASSN-18am declines extremely rapidly for a Type II SN, with a photospheric-phase decline rate of $\sim6.0$ mag (100d)$^{-1}$. Owing to the weakening of HI and the appearance of HeI in its later phases, ASASSN-18am is spectroscopically a Type IIb SN with a partially stripped envelope. However, its photometric and spectroscopic evolution show significant differences from typical SNe IIb. Using a radiative diffusion model, we find that the light curve requires a high synthesised $^{56}$Ni mass of $M_{\rm Ni} \sim0.4\rm M_\odot$ and ejecta with high kinetic energy $E_{\rm kin}= (5-9)\times10^{51} $ erg. Introducing a magnetar central engine still requires $M_{\rm Ni} \sim0.3\, M_\odot$ and $E_{\rm kin}= 3\times10^{51} $ erg. The high $^{56}$Ni mass is consistent with strong iron-group nebular lines in its spectra, which are also similar to several SNe Ic-BL with high $^{56}$Ni yields. The earliest spectrum shows "flash ionization" features, from which we estimate a mass-loss rate of $ \dot{M}\approx 2\times10^{-4} \, \rm M_\odot~yr^{-1} $. This wind density is too low to power the luminous light curve by ejecta-CSM interaction. We measure expansion velocities as high as $ 17,000 \rm~km~s^{-1}$ for H$\alpha$, which is remarkably high compared to other SNe II. We estimate an oxygen core mass of 1.7-3.1 $\rm M_\odot$ using the [OI] luminosity measured from a nebular-phase spectrum, implying a progenitor with a zero-age main sequence mass of 19-24 $\rm M_\odot$.

We present observations of ASASSN-19dj, a nearby tidal disruption event (TDE) discovered in the post-starburst galaxy KUG 0810+227 by the All-Sky Automated Survey for Supernovae (ASAS-SN) at a distance of d $\simeq$ 98 Mpc. We observed ASASSN-19dj from $-$21 to 392 days relative to peak UV/optical emission using high-cadence, multi-wavelength spectroscopy and photometry. From the ASAS-SN $g$-band data, we determine that the TDE began to brighten on 2019 February 6.8 and for the first 25 days the rise was consistent with a flux $\propto$ $t^2$ power-law. ASASSN-19dj peaked in the UV/optical on 2019 March 6.5 (MJD = 58548.5) at a bolometric luminosity of $L = (6.2 \pm 0.2) \times 10^{44} \text{ erg s}^{-1}$. Initially remaining roughly constant in X-rays and slowly fading in the UV/optical, the X-ray flux increased by over an order of magnitude $\sim$225 days after peak, resulting from the expansion of the X-ray emitting surface. The late-time X-ray emission is well-fit by a blackbody with an effective radius of $\sim 1 \times 10^{12} \text{ cm}$ and a temperature of $\sim 6 \times 10^{5} \text{ K}$. Analysis of Catalina Real-Time Transient Survey images reveals a nuclear outburst roughly 14.5 years earlier with a smooth decline and a luminosity of $L_V$ $\geq$ $1.4 \times 10^{43}$ erg s$^{-1}$, although the nature of the flare is unknown. ASASSN-19dj occurred in the most extreme post-starburst galaxy yet to host a TDE, with Lick H$\delta_{A}$ = $7.67 \pm 0.17$ \AA.

We present observations of the unusually luminous Type II supernova (SN) 2016gsd. With a peak absolute magnitude of V = $-$19.95 $\pm$ 0.08, this object is one of the brightest Type II SNe, and lies in the gap of magnitudes between the majority of Type II SNe and the superluminous SNe. Its light curve shows little evidence of the expected drop from the optically thick phase to the radioactively powered tail. The velocities derived from the absorption in H$\alpha$ are also unusually high with the blue edge tracing the fastest moving gas initially at 20000 km s$^{-1}$, and then declining approximately linearly to 15000 km s$^{-1}$ over $\sim$100 d. The dwarf host galaxy of the SN indicates a low-metallicity progenitor which may also contribute to the weakness of the metal lines in its spectra. We examine SN 2016gsd with reference to similarly luminous, linear Type II SNe such as SNe 1979C and 1998S, and discuss the interpretation of its observational characteristics. We compare the observations with a model produced by the JEKYLL code and find that a massive star with a depleted and inflated hydrogen envelope struggles to reproduce the high luminosity and extreme linearity of SN 2016gsd. Instead, we suggest that the influence of interaction between the SN ejecta and circumstellar material can explain the majority of the observed properties of the SN. The high velocities and strong H$\alpha$ absorption present throughout the evolution of the SN may imply a circumstellar medium configured in an asymmetric geometry.

The structural property of the magnetic field in flare-bearing solar active regions (ARs) is one of the key aspects for understanding and forecasting solar flares. In this paper, we make a comparative analysis on the chirality and magnetic configurations associated with two X-class two-ribbon flares happening in AR 10930 and AR 11158. The photospheric magnetic fields of the two ARs were observed by space-based instruments, and the corresponding coronal magnetic fields were calculated based on the nonlinear force-free field model. The analysis shows that the electric current in the two ARs was distributed mostly around the main polarity inversion lines (PILs) where the flares happened, and the magnetic chirality (indicated by the signs of force-free factor $\alpha$) along the main PILs is opposite for the two ARs, i.e., left-handed ($\alpha<0$) for AR 10930 and right-handed ($\alpha>0$) for AR 11158. It is found that, for both the flare events, a prominent magnetic connectivity (featured by co-localized strong $\alpha$ and strong current density distributions) was formed along the main PIL before flare and was totally broken after flare eruption. The two branches of the broken magnetic connectivity, combined with the prominent magnetic connectivity before flare, compose the opposite magnetic configurations in the two ARs owing to their opposite chirality, i.e., Z-shaped configuration in AR 10930 with left-handed chirality and inverse Z-shaped configuration in AR 11158 with right-handed chirality. It is speculated that two-ribbon flares can be generally classified to these two magnetic configurations by chirality in the flare source regions of ARs.

We investigate the interior structure, perturbations, and the associated quasinormal modes of a quantum black hole model recently proposed by Bodendorfer, Mele, and M\"unch (BMM). Within the framework of loop quantum gravity, the quantum parameters in the BMM model are introduced through polymerization, consequently replacing the Schwarzschild singularity with a spacelike transition surface. By treating the quantum geometry corrections as an `effective' matter contribution, we first prove the violation of energy conditions (in particular the null energy condition) near the transition surface and then investigate the required junction conditions on it. In addition, we study the quasinormal modes of massless scalar field perturbations, electromagnetic perturbations, and axial gravitational perturbations in this effective model. As expected, the quasinormal spectra deviate from their classical counterparts in the presence of quantum corrections. Interestingly, we find that the quasinormal frequencies of perturbations with different spins share the same qualitative tendency with respect to the change of the quantum parameters in this model.

The proposed ARIANNA-200 neutrino detector, located at sea-level on the Ross Ice Shelf, Antarctica, consists of 200 autonomous and independent detector stations separated by 1 kilometer in a uniform triangular mesh, and serves as a pathfinder mission for the future IceCube-Gen2 project. The primary science mission of ARIANNA-200 is to search for sources of neutrinos with energies greater than 10^17 eV, complementing the reach of IceCube. An ARIANNA observation of a neutrino source would provide strong insight into the enigmatic sources of cosmic rays. ARIANNA observes the radio emission from high energy neutrino interactions in the Antarctic ice. Among radio based concepts under current investigation, ARIANNA-200 would uniquely survey the vast majority of the southern sky at any instant in time, and an important region of the northern sky, by virtue of its location on the surface of the Ross Ice Shelf in Antarctica. The broad sky coverage is specific to the Moore's Bay site, and makes ARIANNA-200 ideally suited to contribute to the multi-messenger thrust by the US National Science Foundation, Windows on the Universe - Multi-Messenger Astrophysics, providing capabilities to observe explosive sources from unknown directions. The ARIANNA architecture is designed to measure the angular direction to within 3 degrees for every neutrino candidate, which too plays an important role in the pursuit of multi-messenger observations of astrophysical sources.

We present nearly 500 days of observations of the tidal disruption event ASASSN-18pg, spanning from 54 days before peak light to 441 days after peak light. Our dataset includes X-ray, UV, and optical photometry, optical spectroscopy, radio observations, and the first published spectropolarimetric observations of a TDE. ASASSN-18pg was discovered on 2018 July 11 by the All-Sky Automated Survey for Supernovae (ASAS-SN) at a distance of $d=78.6$ Mpc, and with a peak UV magnitude of $m\simeq14$ it is both one of the nearest and brightest TDEs discovered to-date. The photometric data allow us to track both the rise to peak and the long-term evolution of the TDE. ASASSN-18pg peaked at a luminosity of $L\simeq2.2\times10^{44}$ erg s$^{-1}$, and its late-time evolution is shallower than a flux $\propto t^{-5/3}$ power-law model, similar to what has been seen in other TDEs. ASASSN-18pg exhibited Balmer lines and spectroscopic features consistent with Bowen fluorescence prior to peak which remained detectable for roughly 225 days after peak. Analysis of the two-component H$\alpha$ profile indicates that, if they are the result of reprocessing of emission from the accretion disk, the different spectroscopic lines may be coming from regions between $\sim10$ and $\sim60$ light-days from the black hole. No X-ray emission is detected from the TDE and there is no evidence of a jet or strong outflow detected in the radio. Our spectropolarimetric observations give no strong evidence for significant asphericity in the emission region, with the emission region having an axis ratio of at least $\sim0.65$.

Solar filaments are dark structures on the solar disk, with an elongated spine and several barbs extending out from the spine. When appearing above the solar limb, a filament is called a prominence, with several feet extending down to the solar surface. It was generally thought that filament barbs are simply the prominence feet veering away from the spine and down to the solar surface. However, it was recently noticed that there might be another dynamic type of barbs, which were proposed to be due to filament thread longitudinal oscillation. If this is the case, the dynamic barbs would not extend down to the solar surface. With the quadrature observations of a filament barb on 2011 June 5 from the {\it Solar Dynamics Observatory} and the {\it STEREO} satellites, we confirm that the filament barb is due to filament thread longitudinal oscillations. Viewed from the side, the filament barb looks like an appendix along the spine of the prominence, and does not extend down to the solar surface as a foot.

In this paper, we investigate complex-valued Euclidean wormholes in the Starobinsky inflation. Due to the properties of the concave inflaton potential, the classicality condition at both ends of the wormhole can be satisfied, as long as the initial condition of the inflaton field is such that it is located sufficiently close to the hilltop. We compare the probabilities of classicalized wormholes with the Hartle-Hawking compact instantons and conclude that the Euclidean wormholes are probabilistically preferred than compact instantons, if the inflation lasts more than 50 e-foldings. Our result assumes that the Euclidean path integral is the correct effective description of quantum gravity. This opens a new window for various future investigations that can be either confirmed or refuted by future experiments.

We report observations of the hydrogen-deficient supernova (SN) 2019bkc/ATLAS19dqr. With B- and r-band decline between peak and 10 days post peak of Delta m_10(B)=5.24+/-0.07 mag and Delta m_10(r)=3.85+/-0.10$ mag, respectively, SN 2019bkc is the most rapidly declining SN I discovered so far. While its closest matches are the rapidly declining SN 2005ek and SN 2010X, the light curves and spectra of SN 2019bkc show some unprecedented characteristics. SN 2019bkc appears "hostless," with no identifiable host galaxy near its location, although it may be associated with the galaxy cluster MKW1 at z = 0.02. We evaluate a number of existing models of fast-evolving SNe, and we find that none of them can satisfactorily explain all aspects of SN 2019bkc observations.

We present photometric and spectroscopic follow-up observations of the highly luminous Type Ibn supernova ASASSN-14ms, which was discovered on UT 2014-12-26.61 at $m_V \sim 16.5$. With a peak absolute $V$-band magnitude brighter than $-20.5$, a peak bolometric luminosity of $1.7 \times 10^{44}$ ergs s$^{-1}$, and a total radiated energy of $2.1 \times 10^{50}$ ergs, ASASSN-14ms is one of the most luminous Type Ibn supernovae yet discovered. In simple models, the most likely power source for this event is a combination of the radioactive decay of $^{56}$Ni and $^{56}$Co at late times and the interaction of supernova ejecta with the progenitor's circumstellar medium at early times, although we cannot rule out the possibility of a magnetar-powered light curve. The presence of a dense circumstellar medium is indicated by the intermediate-width He I features in the spectra. The faint ($m_g \sim 21.6$) host galaxy SDSS J130408.52+521846.4 has an oxygen abundance below $12+\log(O/H) \lesssim 8.3$, a stellar mass of $M_* \sim 2.6 \times 10^8 M_{\odot}$, and a star formation rate of $\textrm{SFR} \sim 0.02$ $M_{\odot}$ yr$^{-1}$.

We present photometric and spectroscopic observations of the unusual Type Ia supernova ASASSN-18tb, including a series of SALT spectra obtained over the course of nearly six months and the first observations of a supernova by the Transiting Exoplanet Survey Satellite (TESS). We confirm a previous observation by Kollmeier et al. (2019) showing that ASASSN-18tb is the first relatively normal Type Ia supernova to exhibit clear broad ($\sim1000$ km s$^{-1}$) H$\alpha$ emission in its nebular phase spectra. We find that this event is best explained as a sub-Chandrasekhar mass explosion with $M_{Ni} \approx 0.3\; \rm{M}_\odot$. Despite the strong H$\alpha$ signature at late times, we find that the early rise of the supernova shows no evidence for deviations from a single-component power-law and is best fit with a moderately shallow power-law of index $1.69\pm0.04$. We find that the H$\alpha$ luminosity remains approximately constant after its initial detection at phase +37 d, and that the H$\alpha$ velocity evolution does not trace that of the Fe~III$~\lambda4660$ emission. These suggest that the H$\alpha$ emission arises from circumstellar medium (CSM) rather than swept up material from a non-degenerate companion. However, ASASSN-18tb is strikingly different from other known CSM-interacting Type Ia supernovae in a number of significant ways. Those objects typically show an H$\alpha$ luminosity two orders of magnitude higher than what is seen in ASASSN-18tb, pushing them away from the empirical light-curve relations that define "normal" Type Ia supernovae. Conversely, ASASSN-18tb exhibits a fairly typical light curve and luminosity for an underluminous or transitional SN Ia, with $M_R \approx -18.1$ mag. Moreover, ASASSN-18tb is the only SN Ia showing H$\alpha$ from CSM interaction to be discovered in an early-type galaxy.

Pure CO ice has been irradiated with electrons of energy in the range $150-1000$~eV with the Interstellar Energetic-Process System (IEPS). The main products of irradiation are carbon chains C$_n$ ($n=3$, 5, 6, 8, 9, 10, 11, 12), suboxides, C$_n$O ($n=2$, 3, 4, 5, 6, 7), and C$_n$O$_2$ ($n=1$, 3, 4, 5, 7) species. \ce{CO2} is by far the most abundant reaction product in all the experiments. The destruction cross-section of CO peaks at about 250 eV, decreases with the energy of the electrons and is more than one order of magnitude higher than for gas-phase CO ionization. The production cross-section of carbon dioxide has been also derived and is characterized by the competition between chemistry and desorption. Desorption of CO and of new species during the radiolysis follows the electron distribution in the ice. Low energy electrons having short penetration depths induce significant desorption. Finally, as the ice thickness approaches the electron penetration depth the abundance of the products starts to saturate. Implications on the atmospheric photochemistry of cold planets hosting surface CO ices are also discussed.

We present optical photometry and spectroscopy of the Type Ia supernova SN2018cqj/ATLAS18qtd. The supernova exploded in an isolated region at $\sim 65$~kpc from the S0 galaxy IC~550 at $z=0.0165$ ($D\approx 74$~Mpc) and has a redshift consistent with a physical association to this galaxy. Multicolor photometry show that SN2018cqj/ATLAS18qtd is a low-luminosity ($M_{B_{max}}\approx -17.9$ mag), fast-declining Type Ia with color stretch $s_{BV} \approx 0.6$ and $B$-band decline rate $\Delta m_{15}(B) \approx 1.77$ mag. Two nebular-phase spectra obtained as part of the 100IAS survey at +193 and +307 days after peak show the clear detection of a narrow H$\alpha$ line in emission that is resolved in the first spectrum with $\rm FWHM \approx 1200$ km s$^{-1}$ and $L_{H\alpha} \approx 3.8\times 10^{37}$ erg s$^{-1}$. The detection of a resolved H$\alpha$ line with a declining luminosity is broadly consistent with recent models where hydrogen is stripped from the non-degenerate companion in a single-degenerate progenitor system. However, the amount of hydrogen consistent with the luminosities of the H$\alpha$ line would be $\sim 10^{-3}$ M$_{\odot}$, significantly less than theoretical model predictions in the classical single-degenerate progenitor systems. SN2018cqj/ATLAS18qtd is the second low-luminosity, fast-declining Type Ia SN after SN2018fhw/ASASSN-18tb that shows narrow H$\alpha$ in emission in its nebular-phase spectra.

The Askaryan Radio Array is an ultra-high energy (UHE, $>10^{17}$\,eV) neutrino detector designed to observe neutrinos by searching for the radio waves emitted by the relativistic byproducts of neutrino-nucleon interactions in Antarctic ice. In this paper, we present constraints on the diffuse flux of ultra-high energy neutrinos between $10^{16}-10^{21}$\,eV resulting from a search for neutrinos in two complementary analyses, both analyzing four years of data (2013-2016) from the two deep stations (A2, A3) operating at that time. We place a 90\% CL upper limit on the diffuse neutrino flux at $10^{18}$ eV of $EF(E)=3.9\times10^{-16}$~$\textrm{cm}^2$$\textrm{s}^{-1}$$\textrm{sr}^{-1}$. The analysis leverages more than quadruple the livetime of the previous ARA limit, and at this energy is the most sensitive reported by an in-ice radio neutrino detector by more than a factor of four. Looking forward, ARA has yet to analyze more than half of its archived data, and expects that dataset to double in size in the next three years as livetime accumulates in the full five-station array. With this additional livetime, we project ARA will have world-leading sensitivity to neutrinos above $10^{17}$ eV by 2022.

We report the lens mass and distance measurements of the nearby microlensing event TCP J05074264+2447555. We measure the microlens parallax vector ${\pi}_{\rm E}$ using Spitzer and ground-based light curves with constraints on the direction of lens-source relative proper motion derived from Very Large Telescope Interferometer (VLTI) GRAVITY observations. Combining this ${\pi}_{\rm E}$ determination with the angular Einstein radius $\theta_{\rm E}$ measured by VLTI GRAVITY observations, we find that the lens is a star with mass $M_{\rm L} = 0.495 \pm 0.063~M_{\odot}$ at a distance $D_{\rm L} = 429 \pm 21~{\rm pc}$. We find that the blended light basically all comes from the lens. The lens-source proper motion is $\mu_{\rm rel,hel} = 26.55 \pm 0.36~{\rm mas\,yr^{-1}}$, so with currently available adaptive-optics (AO) instruments, the lens and source can be resolved in 2021. This is the first microlensing event whose lens mass is unambiguously measured by interferometry + satellite parallax observations, which opens a new window for mass measurements of isolated objects such as stellar-mass black holes.

We show that standard candles can provide some valuable information about the density contrast, which could be particularly important at redshifts where other observations are not available. We use an inversion method to reconstruct the local radial density profile from luminosity distance observations assuming background cosmological parameters obtained from large scale observations. Using type Ia Supernovae% (SNe) , Cepheids and the cosmological parameters from the Planck mission we reconstruct the radial density profiles along two different directions of the sky. We compare these profiles to other density maps obtained from luminosity density, in particular Keenan et al. 2013 and the \TMPP galaxy catalogue. The method independently confirms the existence of inhomogeneities, could be particularly useful to correctly normalize density maps from galaxy surveys with respect to the average density of the Universe, and could clarify the apparent discrepancy between local and large scale estimations of the Hubble constant. When better observational supernovae data will be available, the accuracy of the reconstructed density profiles will improve and will allow to further investigate the existence of structures whose size is beyond the reach of galaxy surveys.

The past decade has welcomed the emergence of cosmic neutrinos as a new messenger to explore the most extreme environments of the universe. The discovery measurement of cosmic neutrinos, announced by IceCube in 2013, has opened a new window of observation that has already resulted in new fundamental information that holds the potential to answer key questions associated with the high-energy universe, including: what are the sources in the PeV sky and how do they drive particle acceleration; where are cosmic rays of extreme energies produced, and on which paths do they propagate through the universe; and are there signatures of new physics at TeV-PeV energies and above? The planned advancements in neutrino telescope arrays in the next decade, in conjunction with continued progress in broad multimessenger astrophysics, promise to elevate the cosmic neutrino field from the discovery to the precision era and to a survey of the sources in the neutrino sky. The planned detector upgrades to the IceCube Neutrino Observatory, culminating in IceCube-Gen2 (an envisaged $400M facility with anticipated operation in the next decade, described in this white paper) are the cornerstone that will drive the evolution of neutrino astrophysics measurements.

In this catalog we compile information for all supernovae discovered by the All-Sky Automated Survey for SuperNovae (ASAS-SN) as well as all other bright ($m_{peak}\leq17$), spectroscopically confirmed supernovae found in 2017, totaling 308 supernovae. We also present UV through near-IR magnitudes gathered from public databases of all host galaxies for the supernovae in the sample. We perform statistical analyses of our full bright supernova sample, which now contains 949 supernovae discovered since 2014 May 1, including supernovae from our previous catalogs. This is the fourth of a series of yearly papers on bright supernovae and their hosts from the ASAS-SN team, and this work presents updated data and measurements, including light curves, redshifts, classifications, and host galaxy identifications, that supersede information contained in any previous publications.

We present the discovery of ASASSN-18jd (AT 2018bcb), a luminous optical/UV/X-ray transient located in the nucleus of the galaxy 2MASX J22434289$-$1659083 at $z=0.1192$. Swift UVOT photometry shows the UV SED of the transient to be well modeled by a slowly shrinking blackbody with temperature $T \sim 2.5 \times 10^{4} \rm ~K$, a maximum observed luminosity of $L_\text{max} = 4.5^{+0.6}_{-0.3} \times 10^{44} \rm ~erg ~s^{-1}$, and a total radiated energy of $E = 9.6^{+1.1}_{-0.6} \times 10^{51} \rm ~erg$. X-ray data from Swift XRT and XMM-Newton show a transient, variable X-ray flux with blackbody and power-law components. Optical spectra show strong, roughly constant broad Balmer emission as well as transient features attributable to He II, N III-V, O III, and coronal Fe. While ASASSN-18jd shares similarities with Tidal Disruption Events (TDEs), it is also similar to the "rapid turn-on" events seen in quiescent galaxies and in faint Active Galactic Nuclei (AGNs).

RNO is the mid-scale discovery instrument designed to make the first observation of neutrinos from the cosmos at extreme energies, with sensitivity well beyond current instrument capabilities. This new observatory will be the largest ground-based neutrino telescope to date, enabling the measurement of neutrinos above $10^{16}$ eV, determining the nature of the astrophysical neutrino flux that has been measured by IceCube at higher energies, similarly extending the reach of multi-messenger astrophysics to the highest energies, and enabling investigations of fundamental physics at energies unreachable by particle accelerators on Earth.

We present the results of the monitoring campaign of the Type IIn supernova (SN) 2018cnf (aka ASASSN-18mr). It was discovered about 10 days before the maximum light (on MJD = 58293.4+-5.7 in the V band, with MV = -18.13+-0.15 mag). The multiband light curves show an immediate post-peak decline with some minor luminosity fluctuations, followed by a flattening starting about 40 days after maximum. The early spectra are relatively blue and show narrow Balmer lines with P Cygni profiles. Additionally, Fe II, O I, He I and Ca II are detected. The spectra show little evolution with time, with intermediate-width features becoming progressively more prominent, indicating stronger interaction of the SN ejecta with the circumstellar medium. The inspection of archival images from the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) survey has revealed a variable source at the SN position, with a brightest detection in December 2015 at Mr = -14.66+-0.17 mag. This was likely an eruptive phase from the massive progenitor star started from at least mid-2011, and that produced the circumstellar environment within which the star exploded as a Type IIn SN. The overall properties of SN 2018cnf closely resemble those of transients such as SN 2009ip. This similarity favours a massive hypergiant, perhaps a luminous blue variable, as progenitor for SN 2018cnf.

We report on analysis of englacial radio-frequency (RF) pulser data received over horizontal baselines of 1--5 km, based on broadcasts from two sets of transmitters deployed to depths of up to 1500 meters at the South Pole. First, we analyze data collected usingtwo RF bicone transmitters 1400 meters below the ice surface, and frozen into boreholes drilled for the IceCube experiment in 2011. Additionally, in Dec., 2018, a fat-dipole antenna, fed by one of three high-voltage (~1 kV), fast (~(1-5 ns)) signal generators was lowered into the 1700-m deep icehole drilled for the South Pole Ice Core Experiment (SPICE), approximately 3 km from the geographic South Pole. Signals from transmitters were recorded on the five englacial multi-receiver ARA stations, with receiver depths between 60--200 m. We confirm the long, >1 km RF electric field attenuation length, test our observed signal arrival timing distributions against models, and measure birefringent asymmetries at the 0.15% level.

Longitudinal oscillations of solar filament have been investigated via numerical simulations continuously, but mainly in one dimension (1D), where the magnetic field line is treated as a rigid flux tube. Whereas those one-dimensional simulations can roughly reproduce the observed oscillation periods, implying that gravity is the main restoring force for filament longitudinal oscillations, the decay time in one-dimensional simulations is generally longer than in observations. In this paper, we perform a two-dimensional (2D) non-adiabatic magnetohydrodynamic simulation of filament longitudinal oscillations, and compare it with the 2D adiabatic case and 1D adiabatic and non-adiabatic cases. It is found that, whereas both non-adiabatic processes (radiation and heat conduction) can significantly reduce the decay time, wave leakage is another important mechanism to dissipate the kinetic energy of the oscillating filament when the magnetic field is weak so that gravity is comparable to Lorentz force. In this case, our simulations indicate that the pendulum model might lead to an error of ~100% in determining the curvature radius of the dipped magnetic field using the longitudinal oscillation period when the gravity to Lorentz force ratio is close to unity.

A Monte Carlo simulation program for the radio detection of Ultra High Energy (UHE) neutrino interactions in the Antarctic ice as viewed by the Antarctic Impulsive Transient Antenna (ANITA) is described in this article. The program, icemc, provides an input spectrum of UHE neutrinos, the parametrization of the Askaryan radiation generated by their interaction in the ice, and the propagation of the radiation through ice and air to a simulated model of the third and fourth ANITA flights. This paper provides an overview of the icemc simulation, descriptions of the physics models used and of the ANITA electronics processing chain, data/simulation comparisons to validate the predicted performance, and a summary of the impact of published results.

We present a comprehensive dataset of optical and near-infrared photometry and spectroscopy of type~Ia supernova (SN) 2016hnk, combined with integral field spectroscopy (IFS) of its host galaxy, MCG -01-06-070, and nearby environment. Properties of the SN local environment are characterized by means of single stellar population synthesis applied to IFS observations taken two years after the SN exploded. SN 2016hnk spectra are compared to other 1991bg-like SNe Ia, 2002es-like SNe Ia, and Ca-rich transients. In addition, abundance stratification modelling is used to identify the various spectral features in the early phase spectral sequence and the dataset is also compared to a modified non-LTE model previously produced for the sublumnious SN 1999by. SN 2016hnk is consistent with being a sub-luminous (M$_{\rm B}=-16.7$ mag, s$_{\rm BV}$=0.43$\pm$0.03), highly reddened object. IFS of its host galaxy reveals both a significant amount of dust at the SN location, as well as residual star formation and a high proportion of old stellar populations in the local environment compared to other locations in the galaxy, which favours an old progenitor for SN 2016hnk. Inspection of a nebular spectrum obtained one year after maximum contains two narrow emission lines attributed to the forbidden [Ca II] $\lambda\lambda$7291,7324 doublet with a Doppler shift of 700 km s$^{-1}$. Based on various observational diagnostics, we argue that the progenitor of SN 2016hnk was likely a near Chandrasekhar-mass ($M_{\rm Ch}$) carbon-oxygen white dwarf that produced 0.108 $M_\odot$ of $^{56}$Ni. Our modeling suggests that the narrow [Ca II] features observed in the nebular spectrum are associated with $^{48}$Ca from electron capture during the explosion, which is expected to occur only in white dwarfs that explode near or at the $M_{\rm Ch}$ limit.

The Askaryan Radio Array (ARA) is an ultra-high energy (UHE) neutrino telescope at the South Pole consisting of an array of radio antennas aimed at detecting the Askaryan radiation produced by neutrino interactions in the ice. Currently, the experiment has five stations in operation that have been deployed in stages since 2012. This contribution focuses on the development of a search for a diffuse flux of neutrinos in two ARA stations (A2 and A3) from 2013-2016. A background of $\sim 0.01-0.02$ events is expected in one station in each of two search channels in horizontal- and vertical-polarizations. The expected new constraints on the flux of ultra-high energy neutrinos based on four years of analysis with two stations improve on the previous limits set by ARA by a factor of about two. The projected sensitivity of ARA's five-station dataset is beginning to be competitive with other neutrino telescopes at high energies near $10^{10.5}\,$GeV.

We report comprehensive multi-wavelength observations of a peculiar Type Ia-like supernova ("SN Ia-pec") ASASSN-15pz. ASASSN-15pz is a spectroscopic "twin" of SN 2009dc, a so-called "Super-Chandrasekhar-mass" SN, throughout its evolution, but it has a peak luminosity M_B,peak = -19.69 +/- 0.12 mag that is \approx 0.6 mag dimmer and comparable to the SN 1991T sub-class of SNe Ia at the luminous end of the normal width-luminosity relation. The synthesized Ni56 mass of M_Ni56 = 1.13 +/- 0.14 M_sun is also substantially less than that found for several 2009dc-like SNe. Previous well-studied 2009dc-like SNe have generally suffered from large and uncertain amounts of host-galaxy extinction, which is negligible for ASASSN-15pz. Based on the color of ASASSN-15pz, we estimate a host extinction for SN 2009dc of E(B-V)_host=0.12 mag and confirm its high luminosity (M_B, peak[2009dc] \approx -20.3 mag). The 2009dc-like SN population, which represents ~1% of SNe Ia, exhibits a range of peak luminosities, and do not fit onto the tight width-luminosity relation. Their optical light curves also show significant diversity of late-time (>~ 50 days) decline rates. The nebular-phase spectra provide powerful diagnostics to identify the 2009dc-like events as a distinct class of SNe Ia. We suggest referring to these sources using the phenomenology-based "2009dc-like SN Ia-pec" instead of "Super-Chandrasekhar SN Ia," which is based on an uncertain theoretical interpretation.

We present the results of our photometric and spectroscopic follow-up of the intermediate-luminosity optical transient AT 2017jfs. At peak, the object reaches an absolute magnitude of Mg=-15.46+-0.15 mag and a bolometric luminosity of 5.5x10^41 erg/s. Its light curve has the double-peak shape typical of Luminous Red Novae (LRNe), with a narrow first peak bright in the blue bands, while the second peak is longer lasting and more luminous in the red and near-infrared (NIR) bands. During the first peak, the spectrum shows a blue continuum with narrow emission lines of H and Fe II. During the second peak, the spectrum becomes cooler, resembling that of a K-type star, and the emission lines are replaced by a forest of narrow lines in absorption. About 5 months later, while the optical light curves are characterized by a fast linear decline, the NIR ones show a moderate rebrightening, observed until the transient disappeared in solar conjunction. At these late epochs, the spectrum becomes reminiscent of that of M-type stars, with prominent molecular absorption bands. The late-time properties suggest the formation of some dust in the expanding common envelope or an IR echo from foreground pre-existing dust. We propose that the object is a common-envelope transient, possibly the outcome of a merging event in a massive binary, similar to NGC4490-2011OT1.

In order to understand the nature of the accelerating expansion of the late-time universe, it is important to experimentally determine whether dark energy is a cosmological constant or dynamical in nature. If dark energy already exists prior to inflation, which is a reasonable assumption, then one expects that a dynamical dark energy would leave some footprint in the anisotropy of the late-time accelerated expansion. To demonstrate the viability of this notion, we invoke the quintessence field with the exponential potential as one of the simplest dynamical dark energy models allowed by observations. We investigate the effects of its quantum fluctuations (the physical origin of the perturbation being isocurvature) generated during inflation and having fully positive correlation with the primordial curvature perturbations, and estimate the anisotropy of the cosmic expansion so induced. We show that the primordial amplitude of quantum fluctuations of quintessence field {\delta\phi}_P can be related to the tensor-to-scalar ratio r, and we calculate the perturbed luminosity distance to first order and the associated luminosity distance power spectrum, which is an estimator of anisotropicity of late-time accelerated expansion. We find that the gravitational potential at large scales and late times is less decayed in QCDM compared to that in {\Lambda}CDM so that the smaller the redshift and multipole, the more relative deficit of power in QCDM. Our results of luminosity distance power spectrum also show the similar conclusions of suppression as that of the previous investigation regarding the effect of quantum fluctuations of quintessence field on the CMB temperature anisotropies.

In this study, we present the observations of extreme-ultraviolet (EUV) waves associated with an M6.5 flare on 2013 April 11. The event was observed by Solar Dynamics Observatory (SDO) in different EUV channels. The flare was also associated with a halo CME and type II radio bursts. We observed both fast and slow components of the EUV wave. The speed of the fast component, which is identified as a fast-mode MHD wave, varies in the range from 600 to 640 km s^-1 , whereas the speed of the slow-component is ~140 km s^-1 . We observed an unusual phenomenon that, as the fast-component EUV wave passes through two successive magnetic quasi-separatrix layers (QSLs), two stationary wave fronts are formed locally. We propose that part of the outward-propagating fast-mode EUV wave is converted into slow-mode magnetohydrodynamic waves, which are trapped in local magnetic field structures, forming successive stationary fronts. Along the other direction, the fast-component EUV wave also creates oscillations in a coronal loop lying ~225 Mm away from the flare site. We have computed the energy of the EUV wave to be of the order of 10^20 J.

We present a detailed study of the 2017 eruption of the classical nova ASASSN-17pf (LMCN 2017-11a), which is located in the Large Magellanic Cloud, including data from AAVSO, ASAS-SN, SALT, SMARTS, SOAR, and the Neil Gehrels \textit{Swift} Observatory. The optical light-curve is characterized by multiple maxima (flares) on top of a slowly evolving light-curve (with a decline time, $t_2>$ 100 d). The maxima correlate with the appearance of new absorption line systems in the optical spectra characterized by increasing radial velocities. We suggest that this is evidence of multiple episodes of mass-ejection with increasing expansion velocities. The line profiles in the optical spectra indicate very low expansion velocities (FWHM $\sim$ 190 km s$^{-1}$), making this nova one of the slowest expanding ever observed, consistent with the slowly evolving light-curve. The evolution of the colors and spectral energy distribution show evidence of decreasing temperatures and increasing effective radii for the pseudo-photosphere during each maximum. The optical and infrared light-curves are consistent with dust formation 125 days post-discovery. We speculate that novae showing several optical maxima have multiple mass-ejection episodes leading to shocks that may drive $\gamma$-ray emission and dust formation.

The recent discovery of a staggering diversity of planets beyond the Solar System has brought with it a greatly expanded search space for habitable worlds. The Kepler exoplanet survey has revealed that most planets in our interstellar neighborhood are larger than Earth and smaller than Neptune. Collectively termed super-Earths and mini-Neptunes, some of these planets may have the conditions to support liquid water oceans, and thus Earth-like biology, despite differing in many ways from our own planet. In addition to their quantitative abundance, super-Earths are relatively large and are thus more easily detected than true Earth twins. As a result, super-Earths represent a uniquely powerful opportunity to discover and explore a panoply of fascinating and potentially habitable planets in 2020 - 2030 and beyond.

We report distinctly double-peakedH-alpha and H-beta emission lines in the late-time, nebular-phase spectra (>~200 d) of the otherwise normal at early phases (<~ 100 d) Type IIP supernova ASASSN-16at (SN 2016X). Such distinctly double-peaked nebular Balmer lines have never been observed for a Type II SN. The nebular-phase Balmer emission is driven by the radioactive Co56 decay, so the observed line-profile bifurcation suggests a strong bipolarity in the Ni56 distribution or in the line-forming region of the inner ejecta. The strongly bifurcated blue- and red-shifted peaks are separated by ~3x10^3 km/s and are roughly symmetrically positioned with respect to the host-galaxy rest frame, implying that the inner ejecta are composed of two almost detached blobs. The red peak progressively weakens relative to the blue peak, and disappears in the 740 d spectrum. One possible reason for the line-ratio evolution is increasing differential extinction from continuous formation of dust within the envelope, which is also supported by the near-infrared flux excess that develops after ~100 d.

We analyze the 476 SN Ia host galaxies from the All-Sky Automated Survey for Supernova (ASAS-SN) Bright Supernova Catalogs to determine the observed relative Type Ia supernova (SN) rates as a function of luminosity and host galaxy properties. We find that the luminosity distribution of the SNe Ia in our sample is reasonably well described by a Schechter function with a faint-end slope $\alpha \approx 1.5$ and a knee $M_{\star} \approx -18.0$. Our specific SN Ia rates are consistent with previous results but extend to far lower host galaxy masses. We find an overall rate that scales as $(M_{\star}/10^{10} M_{\odot})^{\alpha}$ with $\alpha \approx -0.5$. This shows that the specific SN Ia rate continues rising towards lower masses even in galaxies as small as $\log(M_{\star} / M_{\odot}) \lesssim 7.0$, where it is enhanced by a factor of $\sim10-20$ relative to host galaxies with stellar masses $\sim10^{10}M_{\odot}$. We find no strong dependence of the specific SN Ia rate on the star formation activity of the host galaxies, but additional observations are required to improve the constraints on the star formation rates.

The ANtarctic Impulsive Transient Antenna (ANITA) NASA long-duration balloon payload completed its fourth flight in December 2016, after 28 days of flight time. ANITA is sensitive to impulsive broadband radio emission from interactions of ultra-high-energy neutrinos in polar ice (Askaryan emission). We present the results of two separate blind analyses searching for signals from Askaryan emission in the data from the fourth flight of ANITA. The more sensitive analysis, with a better expected limit, has a background estimate of $0.64^{+0.69}_{-0.45}$ and an analysis efficiency of $82\pm2\%$. The second analysis has a background estimate of $0.34^{+0.66}_{-0.16}$ and an analysis efficiency of $71\pm6\%$. Each analysis found one event in the signal region, consistent with the background estimate for each analysis. The resulting limit further tightens the constraints on the diffuse flux of ultra-high-energy neutrinos at energies above $10^{19.5}$ eV.

As part of the 100IAS survey, a program aimed to obtain nebular-phase spectra for a volume-limited and homogeneous sample of Type Ia supernovae (SNe Ia), we observed ASASSN-18tb (SN 2018fhw) at 139 days past maximum light. ASASSN-18tb was a fast-declining, sub-luminous event that fits well within the observed photometric and spectroscopic distributions of the SN Ia population. We detect a prominent H$\alpha$ emission line of $L_{{\rm H}\alpha}=2.2\pm0.2\times10^{38}$ ergs s$^{-1}$ with FWHM $\approx1100$ km s$^{-1}$ in the nebular-phase spectrum of this SN Ia. High luminosity H$\alpha$ emission ($L_{{\rm H}\alpha}\gtrsim 10^{40}$ ergs~s$^{-1}$) has previously been discovered in a rare class of SNe Ia-like objects showing CSM interactions (SNe Ia-CSM). They predominantly belong to over-luminous ($M_{\rm max}<-19$ mag in optical) 1991T-like SNe Ia and are exclusively found in star-forming galaxies. By contrast, ASASSN-18tb is a sub-luminous SN Ia ($M_{B, {\rm max}}\sim -17.7$ mag) found in an early-type galaxy dominated by old stellar populations. We discuss possible origins for the observed hydrogen. Out of 75 SNe Ia for which we have so far obtained nebular spectra in 100IAS, no other SN shows a $\sim 1000 \,{\rm km s^{-1}}$ H$\alpha$ emission line with comparable line luminosity as ASASSN-18tb, emphasizing the rarity of such emission in the nebular phase. Based on preliminary results from our survey, the rate for ASASSN-18tb-like nebular H$\alpha$ emission could be as high as $\sim 10\%$ level among sub-luminous SNe Ia.

Recently, the ANITA collaboration reported on two upward-going extensive air shower events consistent with a primary particle that emerges from the surface of the ice. These events may be of $\nu_\tau$ origin, in which the neutrino interacts within the Earth to produce a $\tau$ lepton that emerges from the Earth, decays in the atmosphere, and initiates an extensive air shower. In this paper we estimate an upper bound on the ANITA acceptance to a diffuse $\nu_\tau$ flux detected via $\tau$-lepton-induced air showers within the bounds of Standard Model (SM) uncertainties. By comparing this estimate with the acceptance of Pierre Auger Observatory and IceCube and assuming SM interactions, we conclude that a $\nu_\tau$ origin of these events would imply a neutrino flux at least two orders of magnitude above current bounds.

We report the first unambiguous quasi-periodic large-scale extreme-ultraviolet (EUV) wave or shock that was detected by the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory. During the whip-like unwinding eruption of a small filament on 2012 April 24, multiple consecutive large-scale wavefronts emanating from AR11467 were observed simultaneously along the solar surface and a closed transequatorial loop system. In the meantime, an upward propagating dome-shaped wavefront was also observed, whose initial speed and deceleration are about 1392 km/s and 1.78 km/s^2, respectively. Along the solar surface, the quasi-peridoic wavefronts had a period of about 163 +/- 21 seconds and propagated at a nearly constant speed of 747 +/- 26 km/s; they interacted with active region AR11469 and launched a sympathetic upward propagating secondary EUV wave. The wavefronts along the loop system propagated at a speed of 897 km/s, and they were reflected back at the southern end of the loop system at a similar speed. In addition to the propagating waves, a standing kink wave was also present in the loop system simultaneously. Periodicity analysis reveals that the period of the wavefronts was consistent with that of the unwinding helical structures of the erupting filament. Based on these observational facts, we propose that the observed quasi-periodic EUV wavefronts were most likely excited by the periodic unwinding motion of the filament helical structures. In addition, two different seismological methods are applied to derive the magnetic field strength of the loop system, and for the first time the reliability of these inversion techniques are tested with the same magnetic structure.

Using data collected by the Askaryan Radio Array (ARA) experiment at the South Pole, we have used long-baseline propagation of radio-frequency signals to extract information on the radio-frequency index-of-refraction in South Polar ice. Owing to the increasing ice density over the upper 150--200 meters, rays are observed along two, nearly parallel paths, one of which is direct and a second which refracts through an inflection point, with differences in both arrival time and arrival angle that can be used to constrain the neutrino properties. We also observe indications, for the first time, of radio-frequency ice birefringence for signals propagating along predominantly horizontal trajectories, corresponding to an asymmetry of order 0.1% between the ordinary and extra-ordinary paths, numerically compatible with previous measurements of birefringent asymmetries for vertically-propagating radio-frequency signals at South Pole. Taken together, these effects offer the possibility of redundantly measuring the range from receiver to a neutrino interaction in Antarctic ice, if receiver antennas are deployed at shallow (25 m<z<100 m) depths. Such range information is essential in determining both the neutrino energy, as well as the incident neutrino direction.

The combination of magnetohydrodynamic (MHD) simulation and multi-wavelength observations is an effective way to study mechanisms of magnetic flux rope eruption. We develop a data-driven MHD model using the zero-$\beta$ approximation. The initial condition is provided by nonlinear force-free field derived by the magneto-frictional method based on vector magnetic field observed by the Helioseismic and Magnetic Imager (HMI) aboard the Solar Dynamics Observatory (SDO). The bottom boundary uses observed time series of the vector magnetic field and the vector velocity derived by the Differential Affine Velocity Estimator for Vector Magnetograms (DAVE4VM). We apply the data-driven model to active region 11123 observed from 06:00 UT on 2011 November 11 to about 2 hours later. The evolution of the magnetic field topology coincides with the flare ribbons observed in the 304 and 1600 \AA $ $ wavebands by the Atmospheric Imaging Assembly. The morphology, propagation path, and propagation range of the flux rope are comparable with the observations in 304 \AA . We also find that a data-constrained boundary condition, where the bottom boundary is fixed to the initial values, reproduces a similar simulation result. This model can reproduce the evolution of a magnetic flux rope in its dynamic eruptive phase.

Supernova (SN) 2018oh (ASASSN-18bt) is the first spectroscopically-confirmed type Ia supernova (SN Ia) observed in the $Kepler$ field. The $Kepler$ data revealed an excess emission in its early light curve, allowing to place interesting constraints on its progenitor system (Dimitriadis et al. 2018, Shappee et al. 2018b). Here, we present extensive optical, ultraviolet, and near-infrared photometry, as well as dense sampling of optical spectra, for this object. SN 2018oh is relatively normal in its photometric evolution, with a rise time of 18.3$\pm$0.3 days and $\Delta$m$_{15}(B)=0.96\pm$0.03 mag, but it seems to have bluer $B - V$ colors. We construct the "uvoir" bolometric light curve having peak luminosity as 1.49$\times$10$^{43}$erg s$^{-1}$, from which we derive a nickel mass as 0.55$\pm$0.04M$_{\odot}$ by fitting radiation diffusion models powered by centrally located $^{56}$Ni. Note that the moment when nickel-powered luminosity starts to emerge is +3.85 days after the first light in the Kepler data, suggesting other origins of the early-time emission, e.g., mixing of $^{56}$Ni to outer layers of the ejecta or interaction between the ejecta and nearby circumstellar material or a non-degenerate companion star. The spectral evolution of SN 2018oh is similar to that of a normal SN Ia, but is characterized by prominent and persistent carbon absorption features. The C II features can be detected from the early phases to about 3 weeks after the maximum light, representing the latest detection of carbon ever recorded in a SN Ia. This indicates that a considerable amount of unburned carbon exists in the ejecta of SN 2018oh and may mix into deeper layers.

We present an exquisite, 30-min cadence Kepler (K2) light curve of the Type Ia supernova (SN Ia) 2018oh (ASASSN-18bt), starting weeks before explosion, covering the moment of explosion and the subsequent rise, and continuing past peak brightness. These data are supplemented by multi-color Pan-STARRS1 and CTIO 4-m DECam observations obtained within hours of explosion. The K2 light curve has an unusual two-component shape, where the flux rises with a steep linear gradient for the first few days, followed by a quadratic rise as seen for typical SNe Ia. This "flux excess" relative to canonical SN Ia behavior is confirmed in our $i$-band light curve, and furthermore, SN 2018oh is especially blue during the early epochs. The flux excess peaks 2.14$\pm0.04$ days after explosion, has a FWHM of 3.12$\pm0.04$ days, a blackbody temperature of $T=17,500^{+11,500}_{-9,000}$ K, a peak luminosity of $4.3\pm0.2\times10^{37}\,{\rm erg\,s^{-1}}$, and a total integrated energy of $1.27\pm0.01\times10^{43}\,{\rm erg}$. We compare SN 2018oh to several models that may provide additional heating at early times, including collision with a companion and a shallow concentration of radioactive nickel. While all of these models generally reproduce the early K2 light curve shape, we slightly favor a companion interaction, at a distance of $\sim$$2\times10^{12}\,{\rm cm}$ based on our early color measurements, although the exact distance depends on the uncertain viewing angle. Additional confirmation of a companion interaction in future modeling and observations of SN 2018oh would provide strong support for a single-degenerate progenitor system.

On 2018 Feb. 4.41, the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered ASASSN-18bt in the K2 Campaign 16 field. With a redshift of z=0.01098 and a peak apparent magnitude of B_{max}=14.31, ASASSN-18bt is the nearest and brightest SNe Ia yet observed by the Kepler spacecraft. Here we present the discovery of ASASSN-18bt, the K2 light curve, and pre-discovery data from ASAS-SN and the Asteroid Terrestrial-impact Last Alert System (ATLAS). The K2 early-time light curve has an unprecedented 30-minute cadence and photometric precision for an SN~Ia light curve, and it unambiguously shows a ~4 day nearly linear phase followed by a steeper rise. Thus, ASASSN-18bt joins a growing list of SNe Ia whose early light curves are not well described by a single power law. We show that a double-power-law model fits the data reasonably well, hinting that two physical processes must be responsible for the observed rise. However, we find that current models of the interaction with a non-degenerate companion predict an abrupt rise and cannot adequately explain the initial, slower linear phase. Instead, we find that existing, published models with shallow 56Ni are able to span the observed behavior and, with tuning, may be able to reproduce the ASASSN-18bt light curve. Regardless, more theoretical work is needed to satisfactorily model this and other early-time SNe~Ia light curves. Finally, we use Swift X-ray non-detections to constrain the presence of circumstellar material (CSM) at much larger distances and lower densities than possible with the optical light curve. For a constant density CSM these non-detections constrain rho<4.5 * 10^5 cm^-3 at a radius of 4 *10^15 cm from the progenitor star. Assuming a wind-like environment, we place mass-loss limits of Mdot< 8 * 10^-6 M_sun yr^-1 for v_w=100 km s^-1, ruling out some symbiotic progenitor systems.

We employ VLTI GRAVITY to resolve, for the first time, the two images generated by a gravitational microlens. The measurements of the image separation \theta_{-,+}=3.78 +/- 0.05 mas, and hence the Einstein radius \theta_E =1.87 +/- 0.03 mas, are precise. This demonstrates the robustness of the method, provided that the source is bright enough for GRAVITY (K <~ 10.5) and the image separation is of order or larger than the fringe spacing. When \theta_E is combined with a measurement of the "microlens parallax" \pi_E, the two will together yield the lens mass and lens-source relative parallax and proper motion. Because the source parallax and proper motion are well measured by Gaia, this means that the lens characteristics will be fully determined, whether or not it proves to be luminous. This method can be a powerful probe of dark, isolated objects, which are otherwise quite difficult to identify, much less characterize. Our measurement contradicts Einstein's (1936) prediction that "the luminous circle [i.e., microlensed image] cannot be distinguished" from a star.

The Planck data on cosmic microwave background indicates that the Starobinsky-type model with concave inflation potential is favored over the convex-type chaotic inflation. Is there any reason for that? Here we argue that if our universe began with a Euclidean wormhole, then the Starobinsky-type inflation is probabilistically favored. It is known that for a more generic choice of parameters than that originally assumed by Hartle and Hawking, the Hartle-Hawking wave function is dominated by Euclidean wormholes, which can be interpreted as the creation of two classical universes from nothing. We show that only one end of the wormhole can be classicalized for a convex potential, while both ends can be classicalized for a concave potential. The latter is therefore more probable.

Ultra-high energy neutrinos are detectable through impulsive radio signals generated through interactions in dense media, such as ice. Subsurface in-ice radio arrays are a promising way to advance the observation and measurement of astrophysical high-energy neutrinos with energies above those discovered by the IceCube detector ($\geq$1 PeV) as well as cosmogenic neutrinos created in the GZK process ($\geq$100 PeV). Here we describe the $\textit{NuPhase}$ detector, which is a compact receiving array of low-gain antennas deployed 185 m deep in glacial ice near the South Pole. Signals from the antennas are digitized and coherently summed into multiple beams to form a low-threshold interferometric phased array trigger for radio impulses. The NuPhase detector was installed at an Askaryan Radio Array (ARA) station during the 2017/18 Austral summer season. $\textit{In situ}$ measurements with an impulsive, point-source calibration instrument show a 50% trigger efficiency on impulses with voltage signal-to-noise ratios (SNR) of $\le$2.0, a factor of $\sim$1.8 improvement in SNR over the standard ARA combinatoric trigger. Hardware-level simulations, validated with $\textit{in situ}$ measurements, predict a trigger threshold of an SNR as low as 1.6 for neutrino interactions that are in the far field of the array. With the already-achieved NuPhase trigger performance included in ARASim, a detector simulation for the ARA experiment, we find the trigger-level effective detector volume is increased by a factor of 1.8 at neutrino energies between 10 and 100 PeV compared to the currently used ARA combinatoric trigger. We also discuss an achievable near term path toward lowering the trigger threshold further to an SNR of 1.0, which would increase the effective single-station volume by more than a factor of 3 in the same range of neutrino energies.

These proceedings address a recent publication by the ANITA collaboration of four upward- pointing cosmic-ray-like events observed in the first flight of ANITA. Three of these events were consistent with stratospheric cosmic-ray air showers where the axis of propagation does not inter- sect the surface of the Earth. The fourth event was consistent with a primary particle that emerges from the surface of the ice suggesting a possible {\tau}-lepton decay as the origin of this event. These proceedings follow-up on the modeling and testing of the hypothesis that this event was of {\tau} neutrino origin.

Previously we identified a new class of early galaxy that we estimate contributes up to 30\% of the ionizing photons responsible for reionization. These are low mass halos in the range $M_h =10^{6.5}-10^{8} M_{\odot}$ that have been chemically enriched by supernova ejecta from prior Pop III star formation. Despite their low star formation rates, these Metal Cooling halos (MCs) are significant sources of ionizing radiation, especially at the onset of reionization, due to their high number density and ionizing escape fractions. Here we present a fully-coupled radiation hydrodynamic simulation of reionization that includes these MCs as well the more massive hydrogen atomic line cooling halos. Our method is novel: we perform halo finding inline with the radiation hydrodynamical simulation, and assign escaping ionizing fluxes to halos using a probability distribution function (PDF) measured from the galaxy-resolving Renaissance Simulations. The PDF captures the mass dependence of the ionizing escape fraction as well as the probability that a halo is actively forming stars. With MCs, reionization starts earlier than if only halos of $10^8 M_{\odot}$ and above are included, however the redshift when reionization completes is only marginally affected as this is driven by more massive galaxies. Because star formation is intermittent in MCs, the earliest phase of reionization exhibits a stochastic nature, with small H II regions forming and recombining. Only later, once halos of mass $\sim 10^9 M_{\odot}$ and above begin to dominate the ionizing emissivity, does reionization proceed smoothly in the usual manner deduced from previous studies. This occurs at $z\approx 10$ in our simulation.

Magnetic reconnection is a fundamental physical process in various astrophysical, space, and laboratory environments. Many pieces of evidence for magnetic reconnection have been uncovered. However, its specific processes that could be fragmented and turbulent have been short of direct observational evidence. Here, we present observations of a super-hot current sheet during SOL2017-09-10T X8.2-class solar flare that display the fragmented and turbulent nature of magnetic reconnection. As bilateral plasmas converge toward the current sheet, significant plasma heating and non-thermal motions are detected therein. Two oppositely directed outflow jets are intermittently expelled out of the fragmenting current sheet, whose intensity shows a power-law distribution in spatial frequency domain. The intensity and velocity of the sunward outflow jets also display a power-law distribution in temporal frequency domain. The length-to-width ratio of the current sheet is estimated to be larger than theoretical threshold of and thus ensures occurrence of tearing mode instability. The observations therefore suggest fragmented and turbulent magnetic reconnection occurring in the long stretching current sheet.

The Askaryan Radio Array (ARA) reports an observation of radio emission coincident with the "Valentine's Day" solar flare on Feb. 15$^{\rm{th}}$, 2011 in the prototype "Testbed" station. We find $\sim2000$ events that passed our neutrino search criteria during the 70 minute period of the flare, all of which reconstruct to the location of the sun. A signal analysis of the events reveals them to be consistent with that of bright thermal noise correlated across antennas. This is the first natural source of radio emission reported by ARA that is tightly reconstructable on an event-by-event basis. The observation is also the first for ARA to point radio from individual events to an extraterrestrial source on the sky. We comment on how the solar flares, coupled with improved systematic uncertainties in reconstruction algorithms, could aid in a mapping of any above-ice radio emission, such as that from cosmic-ray air showers, to astronomical locations on the sky.

In this paper, we present the observations of two successive fast-mode extreme ultraviolet (EUV) wave events observed on 2016 July 23. Both fast-mode waves were observed by the Atmospheric Imaging Assembly (AIA) instrument on board the Solar Dynamics Observatory (SDO) satellite, with a traveling speed of ~ 675 and 640 km/s, respectively. These two wave events were associated with two filament eruptions and two GOES M-class solar flares from the NOAA active region 12565, which was located near the western limb. The EUV waves mainly move toward the south direction. We observed the interaction of the EUV waves with a helmet streamer further away in the south. When either or one of the EUV waves penetrates into the helmet streamer, a slowly propagating wave with a traveling speed of ~ 150 km/s is observed along the streamer. We suggest that the slowly-moving waves are slow-mode waves, and interpret this phenomenon as the magnetohydrodynamic (MHD) wave mode conversion from the fast mode to the slow mode. Besides, we observed several stationary fronts in the north and south of the source region.

We report a luminous Type II supernova, ASASSN-15nx, with a peak luminosity of M_V=-20 mag, that is between typical core-collapse supernovae and super-luminous supernovae. The post-peak optical light curves show a long, linear decline with a steep slope of 2.5 mag/100 days (i.e., an exponential decline in flux), through the end of observations at phase ~260 days. In contrast, the light curves of hydrogen rich supernovae (SNe II-P/L) always show breaks in their light curves at phase ~100 days, before settling onto Co56 radioactive decay tails with a decline rate of about 1 mag/100 days. The spectra of ASASSN-15nx do not exhibit the narrow emission-line features characteristic of Type IIn SNe, which can have a wide variety of light-curve shapes usually attributed to strong interactions with a dense circumstellar medium (CSM). ASASSN-15nx has a number of spectroscopic peculiarities, including a relatively weak and triangularly-shaped H-alpha emission profile with no absorption component. The physical origin of these peculiarities is unclear, but the long and linear post-peak light curve without a break suggests a single dominant powering mechanism. Decay of a large amount of Ni56 (M_Ni56 = 1.6 +/- 0.2 M_sun) can power the light curve of ASASSN-15nx, and the steep light-curve slope requires substantial gamma-ray escape from the ejecta, which is possible given a low-mass hydrogen envelope for the progenitor. Another possibility is strong CSM interactions powering the light curve, but the CSM needs to be sculpted to produce the unique light-curve shape and to avoid producing SN IIn-like narrow emission lines.

The Antarctic Impulsive Transient Antenna (ANITA), a NASA long-duration balloon payload, searches for radio emission from interactions of ultra-high-energy neutrinos in polar ice. The third flight of ANITA (ANITA-III) was launched in December 2014 and completed a 22-day flight. We present the results of three analyses searching for Askaryan radio emission of neutrino origin. In the most sensitive of the analyses, we find one event in the signal region on an expected a priori background of $0.7^{+0.5}_{-0.3}$. Though consistent with the background estimate, the candidate event remains compatible with a neutrino hypothesis even after additional post-unblinding scrutiny.

We present nebular-phase spectra of the Type Ia supernova (SN Ia) 2016brx, a member of the 1991bg-like subclass that lies at the faint end of the SN Ia luminosity function. Nebular spectra are available for only three other 1991bg-like SNe, and their Co line centers are all within <~ 500 km/s of each other. In contrast, the nebular Co line center of SN 2016brx is blue-shifted by >1500 km/s compared to them and by ~1200 km/s compared to the rest frame. This is a significant shift relative to the narrow nebular line velocity dispersion of <~ 2000 km/s of these SNe. The large range of nebular line shifts implies that the Ni56 in the ejecta of SN 1991bg-like events is off-center by ~1000 km/s rather than universally centrally confined as previously suggested. With the addition of SN 2016brx, the Co nebular line shapes of 1991bg-like objects appear to connect with the brighter SNe Ia that show double-peak profiles, hinting at a continuous distribution of line profiles among SNe Ia. One class of models to produce both off-center and bi-modal Ni56 distributions is collisions of white dwarfs with unequal and equal masses.

We investigate the abiotic production of oxygen and its photochemical byproduct ozone through water vapor photolysis in moist atmospheres of temperate terrestrial exoplanets. The amount of water vapor available for photolysis in the middle atmosphere of a planet can be limited by an atmospheric cold-trap, the formation of which largely depends on the amount of non-condensable gases. We study this effect using a photochemical model coupled to a 1D radiative-convective equilibrium model in atmospheres with N2, CO2 and H2O as the main constituents. We find that in atmospheres with a low N2 inventory, water vapor mixing ratios in the middle atmosphere can be over two orders of magnitude higher compared to atmospheres with an Earth-like N2 inventory. Without a strong surface sink, the non-condensable oxygen can build up rapidly, drying out the upper atmosphere. With a moderate surface sink, the planet can approach a steady state with significant oxygen mixing ratios in which oxygen production is balanced by surface uptake. We use a radiative transfer model to study the spectroscopic fingerprint of these atmospheres in transit observations. Spectral signatures of abiotic oxygen and ozone can be of comparable magnitude as in spectra of Earth seen as an exoplanet. Middle atmospheric water vapor is unlikely to be a usable indicator of the abiotic origin of oxygen because of the influence of oxygen on the water vapor distribution. This suggests that atmospheric oxygen and ozone cannot be used as binary bioindicators and their interpretation will likely require atmospheric and planetary models.

With the rapid development of telescopes, both temporal cadence and the spatial resolution of observations are increasing. This in turn generates vast amount of data, which can be efficiently searched only with automated detections in order to derive the features of interest in the observations. A number of automated detection methods and algorithms have been developed for solar activities, based on the image processing and machine learning techniques. In this paper, after briefly reviewing some automated detection methods, we describe our efficient and versatile automated detection method for solar filaments. It is able not only to recognize filaments, determine the features such as the position, area, spine, and other relevant parameters, but also to trace the daily evolution of the filaments. It is applied to process the full disk H-alpha data observed in nearly three solar cycles, and some statistic results are presented.

We report on an upward traveling, radio-detected cosmic-ray-like impulsive event with characteristics closely matching an extensive air shower. This event, observed in the third flight of the Antarctic Impulsive Transient Antenna (ANITA), a NASA-sponsored long-duration balloon payload, is consistent with a similar event reported in a previous flight. These events may be produced by the atmospheric decay of an upward-propagating $\tau$-lepton produced by a $\nu_{\tau}$ interaction, although their relatively steep arrival angles create tension with the standard model (SM) neutrino cross section. Each of the two events have $a~posteriori$ background estimates of $\lesssim 10^{-2}$ events. If these are generated by $\tau$-lepton decay, then either the charged-current $\nu_{\tau}$ cross section is suppressed at EeV energies, or the events arise at moments when the peak flux of a transient neutrino source was much larger than the typical expected cosmogenic background neutrinos.

Solar prominences are subject to all kinds of perturbations during their lifetime, and frequently demonstrate oscillations. The study of prominence oscillations provides an alternative way to investigate their internal magnetic and thermal structures as the oscillation characteristics depend on their interplay with the solar corona. Prominence oscillations can be classified into longitudinal and transverse types. We perform three-dimensional ideal magnetohydrodynamic simulations of prominence oscillations along a magnetic flux rope, with the aim to compare the oscillation periods with those predicted by various simplified models and to examine the restoring force. We find that the longitudinal oscillation has a period of about 49 minutes, which is in accordance with the pendulum model where the field-ligned component of gravity serves as the restoring force. In contrast, the horizontal transverse oscillation has a period of about 10 minutes and the vertical transverse oscillation has a period of about 14 minutes, and both of them can be nicely fitted with a two-dimensional slab model. We also find that the magnetic tension force dominates most of the time in transverse oscillations, except for the first minute when magnetic pressure overwhelms.

The balloon-borne HiCal radio-frequency (RF) transmitter, in concert with the ANITA radio-frequency receiver array, is designed to measure the Antarctic surface reflectivity in the RF wavelength regime. The amplitude of surface-reflected transmissions from HiCal, registered as triggered events by ANITA, can be compared with the direct transmissions preceding them by O(10) microseconds, to infer the surface power reflection coefficient $\cal{R}$. The first HiCal mission (HiCal-1, Jan. 2015) yielded a sample of 100 such pairs, resulting in estimates of $\cal{R}$ at highly-glancing angles (i.e., zenith angles approaching $90^\circ$), with measured reflectivity for those events which exceeded extant calculations. The HiCal-2 experiment, flying from Dec., 2016-Jan., 2017, provided an improvement by nearly two orders of magnitude in our event statistics, allowing a considerably more precise mapping of the reflectivity over a wider range of incidence angles. We find general agreement between the HiCal-2 reflectivity results and those obtained with the earlier HiCal-1 mission, as well as estimates from Solar reflections in the radio-frequency regime. In parallel, our calculations of expected reflectivity have matured; herein, we use a plane-wave expansion to estimate the reflectivity R from both a flat, smooth surface (and, in so doing, recover the Fresnel reflectivity equations) and also a curved surface. Multiplying our flat-smooth reflectivity by improved Earth curvature and surface roughness corrections now provides significantly better agreement between theory and the HiCal 2a/2b measurements.

Hydrogen-poor superluminous supernovae (SLSNe-I) have been predominantly found in low-metallicity, star-forming dwarf galaxies. Here we identify Gaia17biu/SN 2017egm as an SLSN-I occurring in a "normal" spiral galaxy (NGC 3191) in terms of stellar mass (several times 10^10 M_sun) and metallicity (roughly Solar). At redshift z=0.031, Gaia17biu is also the lowest redshift SLSN-I to date, and the absence of a larger population of SLSNe-I in dwarf galaxies of similar redshift suggests that metallicity is likely less important to the production of SLSNe-I than previously believed. With the smallest distance and highest apparent brightness for an SLSN-I, we are able to study Gaia17biu in unprecedented detail. Its pre-peak near-ultraviolet to optical color is similar to that of Gaia16apd and among the bluest observed for an SLSN-I while its peak luminosity (M_g = -21 mag) is substantially lower than Gaia16apd. Thanks to the high signal-to-noise ratios of our spectra, we identify several new spectroscopic features that may help to probe the properties of these enigmatic explosions. We detect polarization at the ~0.5% level that is not strongly dependent on wavelength, suggesting a modest, global departure from spherical symmetry. In addition, we put the tightest upper limit yet on the radio luminosity of an SLSN-I with <5.4x10^26 erg/s/Hz (at 10 GHz), which is almost a factor of 40 better than previous upper limits and one of the few measured at an early stage in the evolution of an SLSN-I. This limit largely rules out an association of this SLSNe-I with known populations of gamma-ray burst (GRB) like central engines.

White-light flares are the flares with emissions visible in the optical continuum. They are thought to be rare and pose the most stringent requirements in energy transport and heating in the lower atmosphere. Here we present a nearly circular white-light flare on 2015 March 10 that was well observed by the Optical and Near-infrared Solar Eruption Tracer and Solar Dynamics Observatory. In this flare, there appear simultaneously both impulsive and gradual white-light kernels. The generally accepted thick-target model would be responsible for the impulsive kernels but not sufficient to interpret the gradual kernels. Some other mechanisms including soft X-ray backwarming or downward-propagating Alfven waves, acting jointly with electron beam bombardment, provide a possible interpretation. However, the origin of this kind of white-light kernels is still an open question that induces more observations and researches in the future to decipher it.

We present observations of the rise and peak of the Type IIn supernova SN 2017hcc/ATLAS17lsn obtained by the All-Sky Automated Survey for Supernovae (ASAS-SN) and Swift UVOT. The light curve of SN 2017hcc/ATLAS17lsn peaks at $V\simeq 13.7$ mag, which from the estimated redshift of the host galaxy ($z=0.0168$, $D\simeq 73$ Mpc) implies an absolute peak magnitude $M_{V,peak} \simeq -20.7$ mag. The near-UV to optical spectral energy distribution of SN 2017hcc/ATLAS17lsn from Swift UVOT is consistent with a hot, but cooling blackbody with $\rm T_{bb}\simeq 16500$ K on Oct. 28.4 and $\rm T_{bb} \simeq 11700$ K on Nov. 19.6. The estimated peak bolometric luminosity $L_{bol, peak}\simeq 1.3\times 10^{44}$ erg s$^{-1}$ makes SN2017hcc/ATLAS17lsn one of the most luminous Type IIn supernovae studied to date. From the bolometric light curve we constrain the risetime to be $\sim 27$ days and the total radiated energy of the event to date is $4\times 10^{50}$ erg.

The Eddington-inspired-Born-Infeld (EiBI) model is reformulated within the mimetic approach. In the presence of a mimetic field, the model contains non-trivial vacuum solutions which could be free of spacetime singularity because of the Born-Infeld nature of the theory. We study a realistic primordial vacuum universe and prove the existence of regular solutions, such as primordial inflationary solutions of de Sitter type or bouncing solutions. Besides, the linear instabilities present in the EiBI model are found to be avoidable for some interesting bouncing solutions in which the physical metric as well as the auxiliary metric are regular at the background level.

The vacuum, static, and spherically symmetric solutions in the mimetic Born-Infeld gravity are studied. The mimetic Born-Infeld gravity is a reformulation of the Eddington-inspired-Born-Infeld (EiBI) model under the mimetic approach. Due to the mimetic field, the theory contains non-trivial vacuum solutions different from those in Einstein gravity. We find that with the existence of the mimetic field, the spacelike singularity inside a Schwarzschild black hole could be altered to a lightlike singularity, even though the curvature invariants still diverge at the singularity. Furthermore, in this case, the maximal proper time for a timelike radially-infalling observer to reach the singularity is found to be infinite.

This catalog summarizes information for all supernovae discovered by the All-Sky Automated Survey for SuperNovae (ASAS-SN) and all other bright ($m_{peak}\leq17$), spectroscopically confirmed supernovae discovered in 2016. We then gather the near-IR through UV magnitudes of all host galaxies and the offsets of the supernovae from the centers of their hosts from public databases. We illustrate the results using a sample that now totals 668 supernovae discovered since 2014 May 1, including the supernovae from our previous catalogs, with type distributions closely matching those of the ideal magnitude limited sample from Li et al. (2011). This is the third of a series of yearly papers on bright supernovae and their hosts from the ASAS-SN team.

We present and discuss here the observations of a small long duration GOES B- class flare associated with a quiescent filament eruption, a global EUV wave and a CME on 2011 May 11. The event was well observed by the Solar Dynamics Observatory (SDO), GONG H{\alpha}, STEREO and HiRAS spectrograph. As the fil- ament erupted, ahead of the filament we observed the propagation of EIT wave fronts, as well as two flare ribbons on both sides of the polarity inversion line (PIL) on the solar surface. The observations show the co-existence of two types of EUV waves, i.e., a fast and a slow one. A type II radio burst with up to the third harmonic component was also associated with this event. The evolution of pho- tospheric magnetic field showed flux emergence and cancellation at the filament site before its eruption.

The Antarctic Impulsive Transient Antenna (ANITA) is a NASA long-duration balloon experiment with the primary goal of detecting ultra-high-energy ($>10^{18}\,\mbox{eV}$) neutrinos via the Askaryan Effect. The fourth ANITA mission, ANITA-IV, recently flew from Dec 2 to Dec 29, 2016. For the first time, the Tunable Universal Filter Frontend (TUFF) boards were deployed for mitigation of narrow-band, anthropogenic noise with tunable, switchable notch filters. The TUFF boards also performed second-stage amplification by approximately 45 dB to boost the $\sim\,\mu\mbox{V-level}$ radio frequency (RF) signals to $\sim$ mV-level for digitization, and supplied power via bias tees to the first-stage, antenna-mounted amplifiers. The other major change in signal processing in ANITA-IV is the resurrection of the $90^{\circ}$ hybrids deployed previously in ANITA-I, in the trigger system, although in this paper we focus on the TUFF boards. During the ANITA-IV mission, the TUFF boards were successfully operated throughout the flight. They contributed to a factor of 2.8 higher total instrument livetime on average in ANITA-IV compared to ANITA-III due to reduction of narrow-band, anthropogenic noise before a trigger decision is made.

By far cosmology is one of the most exciting subject to study, even more so with the current bulk of observations we have at hand. These observations might indicate different kinds of doomsdays, if dark energy follows certain patterns. Two of these doomsdays are the Little Rip (LR) and Little Sibling of the Big Rip (LSBR). In this work, aside from proving the unavoidability of the LR and LSBR in the Eddington-inspired-Born-Infeld (EiBI) scenario, we carry out a quantum analysis of the EiBI theory with a matter field, which, from a classical point of view would inevitably lead to a universe that ends with either LR or LSBR. Based on a modified Wheeler-DeWitt equation, we demonstrate that such fatal endings seems to be avoidable.

The primary science goal of the NASA-sponsored ANITA project is measurement of ultra-high energy neutrinos and cosmic rays, observed via radio-frequency signals resulting from a neutrino- or cosmic ray- interaction with terrestrial matter (atmospheric or ice molecules, e.g.). Accurate inference of the energies of these cosmic rays requires understanding the transmission/reflection of radio wave signals across the ice-air boundary. Satellite-based measurements of Antarctic surface reflectivity, using a co-located transmitter and receiver, have been performed more-or-less continuously for the last few decades. Satellite-based reflectivity surveys, at frequencies ranging from 2--45 GHz and at near-normal incidence, yield generally consistent reflectivity maps across Antarctica. Using the Sun as an RF source, and the ANITA-3 balloon borne radio-frequency antenna array as the RF receiver, we have also measured the surface reflectivity over the interval 200-1000 MHz, at elevation angles of 12-30 degrees, finding agreement with the Fresnel equations within systematic errors. To probe low incidence angles, inaccessible to the Antarctic Solar technique and not probed by previous satellite surveys, a novel experimental approach ("HiCal-1") was devised. Unlike previous measurements, HiCal-ANITA constitute a bi-static transmitter-receiver pair separated by hundreds of kilometers. Data taken with HiCal, between 200--600 MHz shows a significant departure from the Fresnel equations, constant with frequency over that band, with the deficit increasing with obliquity of incidence, which we attribute to the combined effects of possible surface roughness, surface grain effects, radar clutter and/or shadowing of the reflection zone due to Earth curvature effects.

A class of dynamical dark energy models is constructed through an extended version of fermion fields corresponding to phantom dark ghost spinors, which are spin one half with mass dimension one. We find that if these spinors interact with torsion fields in a homogeneous and isotropic universe, then it does not imply any future dark energy singularity or any abrupt event, though the fermion has a negative kinetic energy. In fact, the equation of state of this dark energy model will asymptotically approach the value $w=-1$ from above without crossing the phantom divide and inducing therefore a de Sitter state. Consequently, we expect the model to be stable because no real phantom fields will be created. At late time, the torsion fields will vanish as the corresponding phantom dark ghost spinors dilute. As would be expected, intuitively, this result is unaffected by the presence of cold dark matter although the proof is not as straightforward as in general relativity.

We present the ultraviolet (UV) spectroscopic evolution of a tidal disruption event (TDE) for the first time. After the discovery of the nearby TDE iPTF16fnl, we obtained a series of observations with the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST). The dominant emission features closely resemble those seen in the UV spectra of the TDE ASASSN-14li and are also similar to those of N-rich quasars. However, there is significant evolution in the shape and central wavelength of the line profiles over the course of our observations, such that at early times the lines are broad and redshifted, while at later times the lines are significantly narrower and peak near the wavelengths of their corresponding atomic transitions. Like ASASSN-14li, but unlike N-rich quasars, iPTF16fnl shows neither MgII$\lambda 2798$\AA\ nor CIII]$\lambda 1909$\AA\ emission features. We also present optical photometry and spectroscopy, which suggest that the complex HeII profiles observed in the optical spectra of many TDEs are in part due to the presence of NIII and CIII Wolf-Rayet features, which can potentially serve as probes of the far-UV when space-based observations are not possible. Finally, we use Swift XRT and UVOT observations to place strong limits on the X-ray emission and determine the characteristic temperature, radius, and luminosity of the emitting material. We find that iPTF16fnl is subluminous and evolves more rapidly than other optically discovered TDEs.

Precise (sub-meter level) real-time navigation using a space-capable single-frequency global positioning system (GPS) receiver and ultra-rapid (real-time) ephemerides from the international global navigation satellite systems service is proposed for low-Earth-orbiting (LEO) satellites. The C/A code and L1 carrier phase measurements are combined and single-differenced to eliminate first-order ionospheric effects and receiver clock offsets. A random-walk process is employed to model the phase ambiguities in order to absorb the time-varying and satellite-specific higher-order measurement errors as well as the GPS clock correction errors. A sequential Kalman filter which incorporates the known orbital dynamic model is developed to estimate orbital states and phase ambiguities without matrix inversion. Real flight data from the single-frequency GPS receiver onboard China's SJ-9A small satellite are processed to evaluate the orbit determination accuracy. Statistics from internal orbit assessments indicate that three-dimensional accuracies of better than 0.50 m and 0.55 mm/s are achieved for position and velocity, respectively.

The large-amplitude longitudinal oscillations of solar filaments have been observed and explored for more than ten years. Previous studies are mainly based on the one-dimensional rigid flux tube model with a single magnetic dip. However, it is noticed that there might be two magnetic dips, and hence two threads, along one magnetic field line. Following the previous work, we intend to investigate the kinematics of the filament longitudinal oscillations when two threads are magnetically connected, which is done by solving one-dimensional radiative hydrodynamic equations with the numerical code MPI-AMRVAC. Two different types of perturbations are considered, and the difference from previous works resulting from the filament thread-thread interaction is investigated. We find that even with the inclusion of the thread-thread interaction, the oscillation period is modified weakly, by at most 20% compared to the traditional pendulum model with one thread. However, the damping timescale is significantly affected by the thread-thread interaction. Hence, we should take it into account when applying the consistent seismology to the filaments where two threads are magnetically connected.

Given its peak luminosity and early-time spectra, ASASSN-15lh was classified as the most luminous supernova (SN) ever discovered (Dong et al. 2016). Here we report a UV rebrightening of ASASSN-15lh observed with Swift during our follow-up campaign. The rebrightening began at $t \simeq$ 90 days (observer frame) after the primary peak and was followed by a $\sim 120$-day long plateau in the bolometric luminosity, before starting to fade again at $t\simeq 210$ days. ASASSN-15lh rebrightened in the Swift UV bands by $\Delta m_{UVW2} \simeq -1.75$ mag, $\Delta m_{UVM2} \simeq -1.25$ mag, and $\Delta m_{UVW1} \simeq -0.8$ mag, but did not rebrighten in the optical bands. Throughout its initial decline, subsequent rebrightening, and renewed decline, the spectra did not show evidence of interactions between the ejecta and circumstellar medium (CSM) such as narrow emission lines. There are hints of weak H$\alpha$ emission at late-times, but Margutti et al. (2016) have shown that it is narrow line emission consistent with star formation in the host nucleus. By fitting a blackbody we find that during the rebrightening the effective photospheric temperature increased from $T_{BB} \simeq 11000$ K to $T_{BB} \simeq 18000$ K. Over the $\sim$ 550 days since its detection, ASASSN-15lh has radiated $\sim 1.7 - 1.9 \times 10^{52}$ ergs. Although its physical nature remains uncertain, the evolution of ASASSN-15lh's photospheric radius, its radiated energy, and the implied event rate, are all more similar to those of H-poor superluminous supernovae (SLSNe-I) than to tidal disruption events (TDEs).

This manuscript presents information for all supernovae discovered by the All-Sky Automated Survey for SuperNovae (ASAS-SN) during 2015, its second full year of operations. The same information is presented for bright ($m_V\leq17$), spectroscopically confirmed supernovae discovered by other sources in 2015. As with the first ASAS-SN bright supernova catalog, we also present redshifts and near-UV through IR magnitudes for all supernova host galaxies in both samples. Combined with our previous catalog, this work comprises a complete catalog of 455 supernovae from multiple professional and amateur sources, allowing for population studies that were previously impossible. This is the second of a series of yearly papers on bright supernovae and their hosts from the ASAS-SN team.

Autonomous orbit determination via integration of epoch-differenced gravity gradients and starlight refraction is proposed in this paper for low-Earth-orbiting satellites operating in GPS-denied environments. The starlight refrac-tion can compensate for the significant along-track position error using solely gravity gradients and benefit from the integration in view of accuracy improvement in radial and cross-track position estimates. The between-epoch dif-ferencing of gravity gradients is employed to eliminate slowly varying measurement biases and noises near the orbit revolution frequency. The refraction angle is directly used as measurement and its Jacobian matrix is derived from an implicit observation equation. An information fusion filter based on sequential extended Kalman filter is devel-oped for the orbit determination. Truth-model simulations are used to test the performance of the algorithm and the effects of differencing intervals and orbital heights are analyzed. A semi-simulation study using actual gravity gra-dient data from the Gravity field and steady-state Ocean Circulation Explorer (GOCE) combined with simulated starlight refraction measurements is further conducted and a three-dimensional position accuracy of better than 100 m is achieved.

We report on a search for ultra-high-energy (UHE) neutrinos from gamma-ray bursts (GRBs) in the data set collected by the Testbed station of the Askaryan Radio Array (ARA) in 2011 and 2012. From 57 selected GRBs, we observed no events that survive our cuts, which is consistent with 0.12 expected background events. Using NeuCosmA as a numerical GRB reference emission model, we estimate upper limits on the prompt UHE GRB neutrino fluence and quasi-diffuse flux from $10^{7}$ to $10^{10}$ GeV. This is the first limit on the prompt UHE GRB neutrino quasi-diffuse flux above $10^{7}$ GeV.

It is important to study the fine structures of solar filaments with high-resolution observations since it can help us understand the magnetic and thermal structures of the filaments and their dynamics. In this paper, we study a newly-formed filament located inside the active region NOAA 11762, which was observed by the 1.6 m New Solar Telescope (NST) at Big Bear Solar Observatory (BBSO) from 16:40:19 UT to 17:07:58 UT on 2013 June 5. As revealed by the H$\alpha$ filtergrams, cool material is seen to be injected into the filament spine with a speed of 5--10 km s$^{-1}$. At the source of the injection, brightenings are identified in the chromosphere, which is accompanied by magnetic cancellation in the photosphere, implying the importance of magnetic reconnection in replenishing the filament with plasmas from the lower atmosphere. Counter-streamings are detected near one endpoint of the filament, with the plane-of-the-sky speed being 7--9 km s$^{-1}$ in the H$\alpha$ red-wing filtergrams and 9--25 km s$^{-1}$ in the blue-wing filtergrams. The observations are indicative of that this active region filament is supported by a sheared arcade without magnetic dips, and the counter-streamings are due to unidirectional flows with alternative directions, rather than due to the longitudinal oscillations of filament threads as in many other filaments.

A novel instantaneous GNSS ambiguity resolution approach which makes use of only single-frequency carrier phase measurements for ultra-short baseline attitude determination is proposed. The Monte Carlo sampling method is employed to obtain the probability density function of ambiguities from a quaternion-based GNSS-attitude model and the LAMBDA method strengthened with a screening mechanism is then utilized to fix the integer values. Experimental results show that 100% success rate could be achieved for ultra-short baselines.

Solar filaments can be formed via chromospheric evaporation followed by condensation in the corona or by the direct injection of cool plasma from the chromosphere to the corona. In this paper, with high-resolution H$\alpha$ data observed by the 1.6 m New Solar Telescope of the Big Bear Solar Observatory on 2015 August 21, we confirmed that an active-region filament is maintained by the continuous injection of cold chromospheric plasma. We find that the filament is rooted along a bright ridge in H$\alpha$, which corresponds to the intersection of a magnetic quasi-separatrix layer with the solar surface. This bright ridge consists of many small patches and the sizes of these patches are comparable to the width of the filament threads. It is found that upflows originate from the brighter patches of the ridge, whereas the downflows move toward the weaker patches of the ridge. The whole filament is composed of two opposite directional streams, implying that longitudinal oscillations are not the only cause of the counterstreamings, and unidirectional siphon flows with alternative directions are another possibility.

We present ground-based and Swift photometric and spectroscopic observations of the tidal disruption event (TDE) ASASSN-15oi, discovered at the center of 2MASX J20390918-3045201 ($d\simeq216$ Mpc) by the All-Sky Automated Survey for SuperNovae (ASAS-SN). The source peaked at a bolometric luminosity of $L\simeq1.3\times10^{44}$ ergs s$^{-1}$ and radiated a total energy of $E\simeq6.6\times10^{50}$ ergs over the first $\sim3.5$ months of observations. The early optical/UV emission of the source can be fit by a blackbody with temperature increasing from $T\sim2\times10^4$ K to $T\sim4\times10^4$ K while the luminosity declines from $L\simeq1.3\times10^{44}$ ergs s$^{-1}$ to $L\simeq2.3\times10^{43}$ ergs s$^{-1}$, requiring the photosphere to be shrinking rapidly. The optical/UV luminosity decline during this period is most consistent with an exponential decline, $L\propto e^{-(t-t_0)/\tau}$, with $\tau \simeq46.5$ days for $t_0\simeq57241.6$ (MJD), while a power-law decline of $L\propto (t-t_0)^{-\alpha}$ with $t_0\simeq57212.3$ and $\alpha=1.62$ provides a moderately worse fit. ASASSN-15oi also exhibits roughly constant soft X-ray emission that is significantly weaker than the optical/UV emission. Spectra of the source show broad helium emission lines and strong blue continuum emission in early epochs, although these features fade rapidly and are not present $\sim3$ months after discovery. The early spectroscopic features and color evolution of ASASSN-15oi are consistent with a TDE, but the rapid spectral evolution is unique among optically-selected TDEs.

It has been revealed that the magnetic topology in the solar atmosphere displays hemispheric preference, i.e., negative/positive helicity in the northern/southern hemisphere, respectively. However, the strength of the hemispheric rule and its cyclic variation are controversial. In this paper, we apply a new method based on filament drainage to 571 erupting filaments from 2010 May to 2015 December in order to determine the filament chirality and its hemispheric preference. It is found that $91.6\%$ of our sample of erupting filaments follow the hemispheric rule of helicity sign. It is found that the strength of the hemispheric preference of the quiescent filaments decreases slightly from $\sim$$97\%$ in the rising phase to $\sim$$85\%$ in the declining phase of solar cycle 24, whereas the strength of the intermediate filaments keeps a high value around $96\pm4\%$ all the time. Only the active-region filaments show significant variations. Their strength of the hemispheric rule rises from $\sim$$63\%$ to $\sim$$95\%$ in the rising phase, and keeps a high value $82\pm5\%$ during the declining phase. Furthermore, during a half-year period around the solar maximum, their hemispheric preference totally vanishes. Besides, we also diagnose the magnetic configurations of the filaments based on our indirect method, and found that in our sample of erupting events, $89\%$ are inverse-polarity filaments with a flux rope magnetic configuration, whereas $11\%$ are normal-polarity filaments with a sheared arcade configuration.

For almost 20 years the physical nature of globally propagating waves in the solar corona (commonly called "EIT waves") has been controversial and subject to debate. Additional theories have been proposed over the years to explain observations that did not fit with the originally proposed fast-mode wave interpretation. However, the incompatibility of observations made using the Extreme-ultraviolet Imaging Telescope (EIT) onboard the Solar and Heliospheric Observatory with the fast-mode wave interpretation was challenged by differing viewpoints from the twin Solar Terrestrial Relations Observatory spacecraft and higher spatial/temporal resolution data from the Solar Dynamics Observatory. In this article, we reexamine the theories proposed to explain "EIT waves" to identify measurable properties and behaviours that can be compared to current and future observations. Most of us conclude that "EIT waves" are best described as fast-mode large-amplitude waves/shocks that are initially driven by the impulsive expansion of an erupting coronal mass ejection in the low corona.

Due to their low cost and low power consumption, single-frequency GPS receivers are considered suitable for low-cost space applications such as small satellite missions. Recently, requirements have emerged for real-time accurate orbit determination at sub-meter level in order to carry out onboard geocoding of high-resolution imagery, open-loop operation of altimeters and radio occultation. This study proposes an improved real-time kinematic positioning method for LEO satellites using single-frequency receivers. The C/A code and L1 phase are combined to eliminate ionospheric effects. The epoch-differenced carrier phase measurements are utilized to acquire receiver position changes which are further used to smooth the absolute positions. A kinematic Kalman filter is developed to implement kinematic orbit determination. Actual flight data from China small satellite SJ-9A are used to test the navigation performance. Results show that the proposed method outperforms traditional kinematic positioning method in terms of accuracy. A 3D position accuracy of 0.72 m and 0.79 m has been achieved using the predicted portion of IGS ultra-rapid products and broadcast ephemerides, respectively.

It is believed that type II radio bursts are generated by shock waves. In order to understand the generation conditions of type II radio bursts, in this paper, we analyze the physical parameters of a shock front. The type II radio burst we selected was observed by Siberian Solar Radio Telescope (SSRT) and Learmonth radio station and was associated with a limb CME occurring on 2014 January 8 observed by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO). The evolution of the CME in the inner corona presents a double-layered structure that propagates outward. We fit the outer layer of the structure with a partial circle and divide it into 7 directions from -45$^\circ$ to 45$^\circ$ with an angular separation of 15$^\circ$. We measure the outer layer speed along the 7 directions, and find that the speed in the direction of -15$^\circ$ with respect to the central direction is the fastest. We use the differential emission measure (DEM) method to calculate the physical parameters at the outer layer at the moment when the type II radio burst was initiated, including the temperature ($T$), emission measure ($EM$), temperature ratio ($T_{d}/T_{u}$), compression ratio ($X$), and Alfv\'en Mach number ($M_{A}$). We compare the quantities $X$ and $M_{A}$ to that obtained from band-splitting in the radio spectrum, and find that this type II radio burst is generated at a small region of the outer layer that is located at the sector in 45$^\circ$ direction. The results suggest that the generation of type II radio bursts (shock) requires larger values of $X$ and $M_{A}$ rather than simply a higher speed of the disturbance.

Outbursts on young stars are usually interpreted as accretion bursts caused by instabilities in the disk or the star-disk connection. However, some protostellar outbursts may not fit into this framework. In this paper, we analyze optical and near-infrared spectra and photometry to characterize the 2015 outburst of the probable young star ASASSN-15qi. The $\sim 3.5$ mag brightening in the $V$ band was sudden, with an unresolved rise time of less than one day. The outburst decayed exponentially by 1 mag for 6 days and then gradually back to the pre-outburst level after 200 days. The outburst is dominated by emission from $\sim10,000$ K gas. An explosive release of energy accelerated matter from the star in all directions, seen in a spectacular cool, spherical wind with a maximum velocity of 1000 km/s. The wind and hot gas both disappeared as the outburst faded and the source returned to its quiescent F-star spectrum. Nebulosity near the star brightened with a delay of 10-20 days. Fluorescent excitation of H$_2$ is detected in emission from vibrational levels as high as $v=11$, also with a possible time delay in flux increase. The mid-infrared spectral energy distribution does not indicate the presence of warm dust emission, although the optical photospheric absorption and CO overtone emission could be related to a gaseous disk. Archival photometry reveals a prior outburst in 1976. Although we speculate about possible causes for this outburst, none of the explanations are compelling.

An innovative orbit determination method which makes use of gravity gradients for Low-Earth-Orbiting satellites is proposed. The measurement principle of gravity gradiometry is briefly reviewed and the sources of measurement error are analyzed. An adaptive hybrid least squares batch filter based on linearization of the orbital equation and unscented transformation of the measurement equation is developed to estimate the orbital states and the measurement biases. The algorithm is tested with the actual flight data from the European Space Agency Gravity field and steady-state Ocean Circulation Explorer. The orbit determination results are compared with the GPS-derived orbits. The radial and cross-track position errors are on the order of tens of meters, whereas the along-track position error is over one order of magnitude larger. The gravity gradient based orbit determination method is promising for potential use in GPS-denied spacecraft navigation.

In this Note, a new approach to spacecraft positioning based on GGT inversion is presented. The gravity gradient tensor is initially measured in the gradiometer reference frame (GRF) and then transformed to the Earth-Centered Earth-Fixed (ECEF) frame via attitude information as well as Earth rotation parameters. Matrix Eigen-Decomposition is introduced to directly translate GGT into position based on the fact that the eigenvalues and eigenvectors of GGT are simplespecific functions of spherical coordinates of the observation position. without the need of an initial position. Unlike the strategy of inertial navigation aiding, no prediction or first guess of the spacecraft position is needed. The method makes use of the J2 gravity model, and is suitable for space navigation where higher frequency terrain contributions to the GGT signals can be neglected.

Spaceborne gravity gradients are proposed in this paper to provide autonomous orbit determination capabilities for near Earth satellites. The gravity gradients contain useful position information which can be extracted by matching the observations with a precise gravity model. The extended Kalman filter is investigated as the principal estimator. The stochastic model of orbital motion, the measurement equation and the model configuration are discussed for the filter design. An augmented state filter is also developed to deal with unknown significant measurement biases. Simulations are conducted to analyze the effects of initial errors, data-sampling periods, orbital heights, attitude and gradiometer noise levels, and measurement biases. Results show that the filter performs well with additive white noise observation errors. Degraded observability for the along-track position is found for the augmented state filter. Real flight data from the GOCE satellite are used to test the algorithm. Radial and cross-track position errors of less than 100 m have been achieved.

We report on four radio-detected cosmic-ray (CR) or CR-like events observed with the Antarctic Impulsive Transient Antenna (ANITA), a NASA-sponsored long-duration balloon payload. Two of the four were previously identified as stratospheric CR air showers during the ANITA-I flight. A third stratospheric CR was detected during the ANITA-II flight. Here we report on characteristics these three unusual CR events, which develop nearly horizontally, 20-30~km above the surface of the Earth. In addition, we report on a fourth steeply upward-pointing ANITA-I CR-like radio event which has characteristics consistent with a primary that emerged from the surface of the ice. This suggests a possible $\tau$-lepton decay as the origin of this event, but such an interpretation would require significant suppression of the Standard Model ${\tau}$-neutrino cross section.

Ultra-high energy neutrinos are interesting messenger particles since, if detected, they can transmit exclusive information about ultra-high energy processes in the Universe. These particles, with energies above $10^{16}\mathrm{eV}$, interact very rarely. Therefore, detectors that instrument several gigatons of matter are needed to discover them. The ARA detector is currently being constructed at South Pole. It is designed to use the Askaryan effect, the emission of radio waves from neutrino-induced cascades in the South Pole ice, to detect neutrino interactions at very high energies. With antennas distributed among 37 widely-separated stations in the ice, such interactions can be observed in a volume of several hundred cubic kilometers. Currently 3 deep ARA stations are deployed in the ice of which two have been taking data since the beginning of the year 2013. In this publication, the ARA detector "as-built" and calibrations are described. Furthermore, the data reduction methods used to distinguish the rare radio signals from overwhelming backgrounds of thermal and anthropogenic origin are presented. Using data from only two stations over a short exposure time of 10 months, a neutrino flux limit of $3 \cdot 10^{-6} \mathrm{GeV} / (\mathrm{cm^2 \ s \ sr})$ is calculated for a particle energy of 10^{18}eV, which offers promise for the full ARA detector.

The variations of the 3-D coronal magnetic fields associated with the X3.4-class flare of active region 10930 are studied in this paper. The coronal magnetic field data are reconstructed from the photospheric vector magnetograms obtained by the Hinode satellite and using the nonlinear force-free field extrapolation method developed in our previous work (He et al., 2011). The 3-D force-free factor $\alpha$, 3-D current density, and 3-D magnetic energy density are employed to analyze the coronal data. The distributions of $\alpha$ and current density reveal a prominent magnetic connectivity with strong negative $\alpha$ values and strong current density before the flare. This magnetic connectivity extends along the main polarity inversion line and is found to be totally broken after the flare. The distribution variation of magnetic energy density reveals the redistribution of magnetic energy before and after the flare. In the lower space of the modeling volume the increase of magnetic energy dominates, and in the higher space the decrease of energy dominates. The comparison with the flare onset imaging observation exhibits that the breaking site of the magnetic connectivity and site with the highest values of energy density increase coincide with the location of flare initial eruption. We conclude that a cramped positive $\alpha$ region appearing in the photosphere causes the breaking of the magnetic connectivity. A scenario for flare initial eruption is proposed in which the Lorentz force acting on the isolated electric current at the magnetic connectivity breaking site lifts the associated plasmas and causes the initial ejection.

After the {\em Solar and Heliospheric Observatory} ({\em SOHO}) was launched in 1996, the aboard Extreme Ultraviolet Imaging Telescope (EIT) observed a global coronal wave phenomenon, which was initially named "EIT wave" after the telescope. The bright fronts are immediately followed by expanding dimmings. It has been shown that the brightenings and dimmings are mainly due to plasma density increase and depletion, respectively. Such a spectacular phenomenon sparked long-lasting interest and debates. The debates were concentrated on two topics, one is about the driving source, and the other is about the nature of this wavelike phenomenon. The controversies are most probably because there may exist two types of large-scale coronal waves that were not well resolved before the {\em Solar Dynamics Observatory} ({\em SDO}) was launched: one is a piston-driven shock wave straddling over the erupting coronal mass ejection (CME), and the other is an apparently propagating front, which may correspond to the CME frontal loop. Such a two-wave paradigm was proposed more than 13 years ago, and now is being recognized by more and more colleagues. In this paper, we review how various controversies can be resolved in the two-wave framework and how important it is to have two different names for the two types of coronal waves.

The discovery of stationary "EIT waves" about 16 years ago posed a big challenge to the then favorite fast-mode wave model for coronal "EIT waves". It encouraged the proposing of various non-wave models, and played an important role in approaching the recent converging viewpoint, {\it i.e.} there are two types of EUV waves. However, it was recently discovered that a stationary wave front can also be generated when a fast-mode wave passes through a magnetic quasi-separatrix layer (QSL). In this paper, we perform a magnetohydrodynamic (MHD) numerical simulation of the interaction between a fast-mode wave and a magnetic QSL, and a stationary wave front is reproduced. The analysis of the numerical results indicates that near the plasma beta $\sim 1$ layer in front of the magnetic QSL, part of the fast-mode wave is converted to a slow-mode MHD wave, which is then trapped inside the magnetic loops, forming a stationary wave front. Our research implies that we have to be cautious in identifying the nature of a wave since there may be mode conversion during the propagation of the waves driven by solar eruptions.

For fifty years, cosmic-ray air showers have been detected by their radio emission. We present the first laboratory measurements that validate electrodynamics simulations used in air shower modeling. An experiment at SLAC provides a beam test of radio-frequency (RF) radiation from charged particle cascades in the presence of a magnetic field, a model system of a cosmic-ray air shower. This experiment provides a suite of controlled laboratory measurements to compare to particle-level simulations of RF emission, which are relied upon in ultra-high-energy cosmic-ray air shower detection. We compare simulations to data for intensity, linearity with magnetic field, angular distribution, polarization, and spectral content. In particular, we confirm modern predictions that the magnetically induced emission in a dielectric forms a cone that peaks at the Cherenkov angle and show that the simulations reproduce the data within systematic uncertainties.

We present and interpret the observations of extreme ultraviolet (EUV) waves associated with a filament eruption on 2011 May 11.The filament eruption also produces a small B-class two ribbon flare and a coronal mass ejection (CME). The event is observed by the Solar Dynamic Observatory (SDO) with high spatio-temporal resolution data recorded by Atmospheric Imaging Assembly (AIA). As the filament erupts, we observe two types of EUV waves (slow and fast) propagating outwards. The faster EUV wave has a propagation velocity of ~ 500 km/s and the slower EUV wave has an initial velocity of ~ 120 km/s. We report for the first time that not only the slower EUV wave stops at a magnetic separatrix to form bright stationary fronts, but also the faster EUV wave transits a magnetic separatrix, leaving another stationary EUV front behind.

If the Hartle-Hawking wave function is the correct boundary condition of our universe, the history of our universe will be well approximated by an instanton. Although this instanton should be classicalized at infinity, as long as we are observing a process of each history, we may detect a non-classicalized part of field combinations. When we apply it to a dark energy model, this non-classicalized part of fields can be well embedded to a quintessence and a phantom model, i.e., a quintom model. Because of the property of complexified instantons, the phantomness will be naturally free from a big rip singularity. This phantomness does not cause perturbative instabilities, as it is an effect emergent from the entire wave function. Our work may thus provide a theoretical basis for the quintom models, whose equation of state (EoS) can cross the cosmological constant boundary (CCB) phenomenologically.

The first flight of the Antarctic Impulsive Transient Antenna (ANITA) experiment recorded 16 radio signals that were emitted by cosmic-ray induced air showers. For 14 of these events, this radiation was reflected from the ice. The dominant contribution to the radiation from the deflection of positrons and electrons in the geomagnetic field, which is beamed in the direction of motion of the air shower. This radiation is reflected from the ice and subsequently detected by the ANITA experiment at a flight altitude of 36km. In this paper, we estimate the energy of the 14 individual events and find that the mean energy of the cosmic-ray sample is 2.9 EeV. By simulating the ANITA flight, we calculate its exposure for ultra-high energy cosmic rays. We estimate for the first time the cosmic-ray flux derived only from radio observations. In addition, we find that the Monte Carlo simulation of the ANITA data set is in agreement with the total number of observed events and with the properties of those events.

We present the first survey of dense gas towards Planck Galactic Cold Clumps (PGCCs). Observations in the J=1-0 transitions of HCO+ and HCN towards 621 molecular cores associated with PGCCs were performed using the Purple Mountain Observatory 13.7-m telescope. Among them, 250 sources have detection, including 230 cores detected in HCO+ and 158 in HCN. Spectra of the J=1-0 transitions from CO, 13CO, and C18O at the centers of the 250 cores were extracted from previous mapping observations to construct a multi-line data set. The significantly low detection rate of asymmetric double-peaked profiles, together with the well consistence among central velocities of CO, HCO+, and HCN spectra, suggests that the CO-selected Planck cores are more quiescent compared to classical star-forming regions. The small difference between line widths of C18O and HCN indicates that the inner regions of CO-selected Planck cores are not more turbulent than the exterior. The velocity-integrated intensities and abundances of HCO+ are positively correlated with those of HCN, suggesting these two species are well coupled and chemically connected. The detected abundances of both HCO+ and HCN are significantly lower than values in other low- to high-mass star-forming regions. The low abundances may be due to beam dilution. On the basis of the inspection of the parameters given in the PGCC catalog, we suggest that there may be about 1 000 PGCC objects having sufficient reservoir of dense gas to form stars.

We report the discovery of ASASSN-15lh (SN 2015L), which we interpret as the most luminous supernova yet found. At redshift z = 0.2326, ASASSN-15lh reached an absolute magnitude of M_{u,AB} = -23.5+/-0.1 and bolometric luminosity L_bol = (2.2+/-0.2)x 10^45 ergs s^-1, which is more than twice as luminous as any previously known supernova. It has several major features characteristic of the hydrogen-poor super-luminous supernovae (SLSNe-I), whose energy sources and progenitors are currently poorly understood. In contrast to most previously known SLSNe-I that reside in star-forming dwarf galaxies, ASASSN-15lh appears to be hosted by a luminous galaxy (M_K ~ -25.5) with little star formation. In the 4 months since first detection, ASASSN-15lh radiated (1.1+/- 0.2)x10^52 ergs, challenging the magnetar model for its engine.

There is an apparent power deficit relative to the $\Lambda$CDM prediction of the CMB spectrum at large scales, which, though not yet statistically significant, persists from WMAP to Planck data. Proposals that invoke some form of initial condition for the inflation have been made to address this apparent power suppression, albeit with conflicting conclusions. By studying the curvature perturbations of a scalar field in the FLRW universe parameterized by the equation of state parameter $w$, we find that the large-scale spectrum at the end of inflation reflects the super-horizon spectrum of the initial state. The large-scale spectrum is suppressed if the universe begins with the adiabatic vacuum in a super-inflation ($w < -1$) or positive-pressure ($w > 0$) era. In the latter case, there is however no causal mechanism to establish the initial adiabatic vacuum. On the other hand, as long as the universe begins with the adiabatic vacuum in an era with $-1 < w < 0$, even if there exists an intermediate positive-pressure era, the large-scale spectrum would be enhanced rather than suppressed. We further calculate the spectrum of a two-stage inflation model with a two-field potential and show that the result agrees with that obtained from the ad hoc single-field analysis.

In this paper, a modified Eddington-inspired-Born-Infeld (EiBI) theory with a pure trace term $g_{\mu\nu}R$ being added to the determinantal action is analysed from a cosmological point of view. It corresponds to the most general action constructed from a rank two tensor that contains up to first order terms in curvature. This term can equally be seen as a conformal factor multiplying the metric $g_{\mu\nu}$. This very interesting type of amendment has not been considered within the Palatini formalism despite the large amount of works on the Born-Infeld-inspired theory of gravity. This model can provide smooth bouncing solutions which were not allowed in the EiBI model for the same EiBI coupling. Most interestingly, for a radiation filled universe there are some regions of the parameter space that can naturally lead to a de Sitter inflationary stage without the need of any exotic matter field. Finally, in this model we discover a new type of cosmic "quasi-sudden" singularity, where the cosmic time derivative of the Hubble rate becomes very large but finite at a finite cosmic time.

We improve our filament automated detection method which was proposed in our previous works. It is then applied to process the full disk H$\alpha$ data mainly obtained by Big Bear Solar Observatory (BBSO) from 1988 to 2013, spanning nearly 3 solar cycles. The butterfly diagrams of the filaments, showing the information of the filament area, spine length, tilt angle, and the barb number, are obtained. The variations of these features with the calendar year and the latitude band are analyzed. The drift velocities of the filaments in different latitude bands are calculated and studied. We also investigate the north-south (N-S) asymmetries of the filament numbers in total and in each subclass classified according to the filament area, spine length, and tilt angle. The latitudinal distribution of the filament number is found to be bimodal. About 80% of all the filaments have tilt angles within [0{\deg}, 60{\deg}]. For the filaments within latitudes lower (higher) than 50{\deg} the northeast (northwest) direction is dominant in the northern hemisphere and the southeast (southwest) direction is dominant in the southern hemisphere. The latitudinal migrations of the filaments experience three stages with declining drift velocities in each of solar cycles 22 and 23, and it seems that the drift velocity is faster in shorter solar cycles. Most filaments in latitudes lower (higher) than 50{\deg} migrate toward the equator (polar region). The N-S asymmetry indices indicate that the southern hemisphere is the dominant hemisphere in solar cycle 22 and the northern hemisphere is the dominant one in solar cycle 23.

We present high-resolution observations of a quiescent solar prominence which was consisted of a vertical and a horizontal foot encircled by an overlying spine, and counter-streaming mass flows were ubiquitous in the prominence. While the horizontal foot and the spine were connecting to the solar surface, the vertical foot was suspended above the solar surface and supported by a semicircular bubble structure. The bubble first collapsed and then reformed at a similar height, finally, it started to oscillate for a long time. We find that the collapsing and oscillation of the bubble boundary were tightly associated with a flare-like feature located at the bottom of the bubble. Based on the observational results, we propose that the prominence should be composed of an overlying horizontal spine encircling a low-lying horizontal and a vertical foot, in which the horizontal foot was consisted of shorter field lines running partially along the spine and with the both ends connecting to the solar surface, while the vertical foot was consisted of piling-up dips due to the sagging of the spine fields and supported by a bipolar magnetic system formed by parasitic polarities (i.e., the bubble). The upflows in the vertical foot were possibly caused by the magnetic reconnection at the separator between the bubble and the overlying dips, which intruded into the persistent downflow field and formed the picture of counter-streaming mass flows. In addition, the counter-streaming flows in the horizontal foot were possibly caused by the imbalanced pressure at the both ends.

A magnetic flux rope structure is believed to exist in most coronal mass ejections (CMEs). However, it has been long debated whether the flux rope exists before eruption or is formed during eruption via magnetic reconnection. The controversy has been continuing because of our lack of routine measurements of the magnetic field in the pre-eruption structure, such as solar filaments. However, recently an indirect method was proposed to infer the magnetic field configuration based on the sign of helicity and the bearing direction of the filament barbs. In this paper, we apply this method to two erupting filament events, one on 2014 September 2 and the other on 2011 March 7, and find that the first filament is supported by a magnetic flux rope and the second filament is supported by a sheared arcade, i.e., the first one is an inverse-polarity filament and the second one is a normal-polarity filament. With the identification of the magnetic configurations in these two filaments, we stress that a flux rope is not a necessary condition for the pre-CME structure.

We investigate cusp singularities in f(R) gravity, especially for Starobinsky and Hu-Sawicki dark energy models. We illustrate that, by using double-null numerical simulations, a cusp singularity can be triggered by gravitational collapses. This singularity can be cured by adding a quadratic term, but this causes a Ricci scalar bump that can be observed by an observer outside the event horizon. Comparing with cosmological parameters, it seems that it would be difficult to see super-Planckian effects by astrophysical experiments. On the other hand, at once there exists a cusp singularity, it can be a mechanism to realize a horizon scale curvature singularity that can be interpreted by a firewall.

We propose a new scalar-tensor model which induces significant deviation from general relativity inside dense objects like neutron stars, while passing solar-system and terrestrial experiments, extending a model proposed by Damour and Esposito-Farese. Unlike their model, we employ a massive scalar field dubbed asymmetron so that it not only realizes proper cosmic evolution but also can account for the cold dark matter. In our model, asymmetron undergoes spontaneous scalarization inside dense objects, which results in reduction of the gravitational constant by a factor of order unity. This suggests that observational tests of constancy of the gravitational constant in high density phase are the effective ways to look into the asymmetron model.

X-ray flares have routinely been observed from the supermassive black hole, Sagittarius A$^\star$ (Sgr A$^\star$), at our Galactic center. The nature of these flares remains largely unclear, despite of many theoretical models. In this paper, we study the statistical properties of the Sgr A$^\star$ X-ray flares, by fitting the count rate (CR) distribution and the structure function (SF) of the light curve with a Markov Chain Monte Carlo (MCMC) method. With the 3 million second \textit{Chandra} observations accumulated in the Sgr A$^\star$ X-ray Visionary Project, we construct the theoretical light curves through Monte Carlo simulations. We find that the $2-8$ keV X-ray light curve can be decomposed into a quiescent component with a constant count rate of $\sim6\times10^{-3}~$count s$^{-1}$ and a flare component with a power-law fluence distribution $dN/dE\propto E^{-\alpha_{\rm E}}$ with $\alpha_{\rm E}=1.65\pm0.17$. The duration-fluence correlation can also be modelled as a power-law $T\propto E^{\alpha_{\rm ET}}$ with $\alpha_{\rm ET} < 0.55$ ($95\%$ confidence). These statistical properties are consistent with the theoretical prediction of the self-organized criticality (SOC) system with the spatial dimension $S = 3$. We suggest that the X-ray flares represent plasmoid ejections driven by magnetic reconnection (similar to solar flares) in the accretion flow onto the black hole.

The Ultra-Fast Flash Observatory pathfinder (UFFO-p) is a new space mission dedicated to detect Gamma-Ray Bursts (GRBs) and rapidly follow their afterglows in order to provide early optical/ultraviolet measurements. A GRB location is determined in a few seconds by the UFFO Burst Alert & Trigger telescope (UBAT) employing the coded mask imaging technique and the detector combination of Yttrium Oxyorthosilicate (YSO) scintillating crystals and multi-anode photomultiplier tubes. The results of the laboratory tests of UBAT's functionality and performance are described in this article. The detector setting, the pixel-to-pixel response to X-rays of different energies, the imaging capability for <50 keV X-rays, the localization accuracy measurements, and the combined test with the Block for X-ray and Gamma-Radiation Detection (BDRG) scintillator detector to check the efficiency of UBAT are all described. The UBAT instrument has been assembled and integrated with other equipment on UFFO-p and should be launched on board the Lomonosov satellite in late-2015.

Two-ribbon brightenings are one of the most remarkable characteristics of an eruptive solar flare and are often used for predicting the occurrence of coronal mass ejections (CMEs). Nevertheless, it was called in question recently whether all two-ribbon flares are eruptive. In this paper, we investigate a two ribbon-like white-light (WL) flare that is associated with a failed magnetic flux rope (MFR) eruption on 2015 January 13, which has no accompanying CME in the WL coronagraph. Observations by \textit{Optical and Near-infrared Solar Eruption Tracer} and \textit{Solar Dynamics Observatory} reveal that, with the increase of the flare emission and the acceleration of the unsuccessfully erupting MFR, two isolated kernels appear at the WL 3600 {\AA} passband and quickly develop into two elongated ribbon-like structures. The evolution of the WL continuum enhancement is completely coincident in time with the variation of \textit{Fermi} hard X-ray 26--50 keV flux. Increase of continuum emission is also clearly visible at the whole FUV and NUV passbands observed by \textit{Interface Region Imaging Spectrograph}. Moreover, in one WL kernel, the \ion{Si}{4}, \ion{C}{2}, and \ion{Mg}{2} h/k lines display significant enhancement and non-thermal broadening. However, their Doppler velocity pattern is location-dependent. At the strongly bright pixels, these lines exhibit a blueshift; while at moderately bright ones, the lines are generally redshifted. These results show that the failed MFR eruption is also able to produce a two-ribbon flare and high-energy electrons that heat the lower atmosphere, causing the enhancement of the WL and FUV/NUV continuum emissions and chromospheric evaporation.

We attempt to propose a method for automatically detecting the solar filament chirality and barb bearing. We first introduce the unweighted undirected graph concept and adopt the Dijkstra shortest-path algorithm to recognize the filament spine. Then, we use the polarity inversion line (PIL) shift method for measuring the polarities on both sides of the filament, and employ the connected components labeling method to identify the barbs and calculate the angle between each barb and the spine to determine the bearing of the barbs, i.e., left or right. We test the automatic detection method with H-alpha filtergrams from the Big Bear Solar Observatory (BBSO) H-alpha archive and magnetograms observed with the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO). Four filaments are automatically detected and illustrated to show the results. The barbs in different parts of a filament may have opposite bearings. The filaments in the southern hemisphere (northern hemisphere) mainly have left-bearing (right-bearing) barbs and positive (negative) magnetic helicity, respectively. The tested results demonstrate that our method is efficient and effective in detecting the bearing of filament barbs. It is demonstrated that the conventionally believed one-to-one correspondence between filament chirality and barb bearing is not valid. The correct detection of the filament axis chirality should be done by combining both imaging morphology and magnetic field observations.

Magnetic reconnection, a change of magnetic field connectivity, is a fundamental physical process in which magnetic energy is released explosively. It is responsible for various eruptive phenomena in the universe. However, this process is difficult to observe directly. Here, the magnetic topology associated with a solar reconnection event is studied in three dimensions (3D) using the combined perspectives of two spacecraft. The sequence of extreme ultraviolet (EUV) images clearly shows that two groups of oppositely directed and non-coplanar magnetic loops gradually approach each other, forming a separator or quasi-separator and then reconnecting. The plasma near the reconnection site is subsequently heated from $\sim$1 to $\ge$5 MK. Shortly afterwards, warm flare loops ($\sim$3 MK) appear underneath the hot plasma. Other observational signatures of reconnection, including plasma inflows and downflows, are unambiguously revealed and quantitatively measured. These observations provide direct evidence of magnetic reconnection in a 3D configuration and reveal its origin.

The Askaryan Radio Array (ARA) is an ultra-high energy ($>10^{17}$ eV) cosmic neutrino detector in phased construction near the South Pole. ARA searches for radio Cherenkov emission from particle cascades induced by neutrino interactions in the ice using radio frequency antennas ($\sim150-800$ MHz) deployed at a design depth of 200 m in the Antarctic ice. A prototype ARA Testbed station was deployed at $\sim30$ m depth in the 2010-2011 season and the first three full ARA stations were deployed in the 2011-2012 and 2012-2013 seasons. We present the first neutrino search with ARA using data taken in 2011 and 2012 with the ARA Testbed and the resulting constraints on the neutrino flux from $10^{17}-10^{21}$ eV.

Ellerman bombs (EBs) are tiny brightenings often observed near sunspots. The most impressive characteristic of the EB spectra is the two emission bumps in both wings of the H$\alpha$ and \ion{Ca}{II} 8542 {\AA} lines. High-resolution spectral data of three small EBs were obtained on 2013 June 6 with the largest solar telescope, the 1.6 meter New Solar Telescope (NST), at the Big Bear Solar Observatory. The characteristics of these EBs are analyzed. The sizes of the EBs are in the range of 0.3\arcsec\--0.8\arcsec\ and their durations are only 3--5 minutes. Our semi-empirical atmospheric models indicate that the heating occurs around the temperature minimum region with a temperature increase of 2700--3000 K, which is surprisingly higher than previously thought. The radiative and kinetic energies are estimated to be as high as 5$\times$10$^{25}$--3.0$\times$10$^{26}$ ergs despite the small size of these EBs. Observations of the magnetic field show that the EBs appeared just in a parasitic region with mixed polarities and accompanied by mass motions. Nonlinear force-free field extrapolation reveals that the three EBs are connected with a series of magnetic field lines associated with bald patches, which strongly implies that these EBs should be produced by magnetic reconnection in the solar lower atmosphere. According to the lightcurves and the estimated magnetic reconnection rate, we propose that there is a three phase process in EBs: pre-heating, flaring and cooling phases.

By combining the data of the Two Micron All Sky Survey, the \textit{Wide Field Infrared Survey Explorer} and the \textit{AKARI} satellite, we study the infrared colour properties of a sample of 2712 nearby radio-luminous galaxies (RLGs). These RLGs are divided into radio-loud (RL) active galactic nuclei (AGNs), mainly occurring at redshifts of $0.05<z<0.3$ and star-forming-dominated RLGs (SFGs), mainly occurring at redshifts of $0.01<z<0.15$. RL AGNs and SFGs are separately distributed in the ([3.4]$-$[4.6])$-$([4.6]$-$[12]) two-colour diagram, in which the RL AGNs display a double-core distribution, and the SFGs display a single-core distribution. SFGs have a redder [4.6]$-$[12] colour than RL AGNs due to the significant contribution from the dust component of SFGs. We find simple criteria of MIR colour separation between RL AGNs and SFGs such that: 95$\%$ of RL AGNs have [4.6]$-$[12] $<$ 3.0 and 94$\%$ of SFGs have [4.6]$-$[12] $>$ 3.0. We also analyse the MIR colours of RL AGNs divided into low- and high-excitation radio galaxies (LERGs and HERGs, respectively). The ([3.4]$-$[4.6])$-$([4.6]$-$[12]) diagram clearly shows separate distributions of LERGs and HERGs and a region of overlap, which suggests that LERGs and HERGs have different MIR properties. LERGs are responsible for the double-core distribution of RL AGNs on the ([3.4]$-$[4.6])$-$([4.6]$-$[12]) diagram. In addition, we also suggest 90$-$140$\mu$m band spectral index $\alpha(90,140)<-1.4$ as a criterion of selecting nearby active galaxies with non-thermal emissions at far-infrared wavelengths.

We report and analyze the observational evidence of global kink oscillations in a solar filament as observed in H alpha by National Solar Observatory (NSO)/Global Oscillation Network Group (GONG) instrument. An M1.1-class flare in active region 11692 on 2013 March 15 induced a global kink mode in the filament lying in the south-west of AR11692.We find periods of about 61 - 67 minutes and damping times of 92 - 117 minutes at three vertical slice positions chosen in and around the filament apex. We find that the waves are damped. From the observed global kink mode period and damping time scale using the theory of resonant absorption we perform prominence seismology. We estimate a lower cut-off value for the inhomogeneity length-scale to be around 0.34 - 0.44 times the radius of the filament cross-section.

Type II radio bursts are thought to be a signature of coronal shocks. In this paper, we analyze a short-lived type II burst that started at 07:40 UT on 2011 February 28. By carefully checking white-light images, we find that the type II radio burst is not accompanied by a coronal mass ejection, only with a C2.4 class flare and narrow jet. However, in the extreme-ultraviolet (EUV) images provided by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO), we find a wave-like structure that propagated at a speed of $\sim$ 600 km s$^{-1}$ during the burst. The relationship between the type II radio burst and the wave-like structure is in particular explored. For this purpose, we first derive the density distribution under the wave by the differential emission measure (DEM) method, which is used to restrict the empirical density model. We then use the restricted density model to invert the speed of the shock that produces the observed frequency drift rate in the dynamic spectrum. The inverted shock speed is similar to the speed of the wave-like structure. This implies that the wave-like structure is most likely a coronal shock that produces the type II radio burst. We also examine the evolution of the magnetic field in the flare-associated active region and find continuous flux emergence and cancellation taking place near the flare site. Based on these facts, we propose a new mechanism for the formation of the type II radio burst, i.e., the expansion of the strongly-inclined magnetic loops after reconnected with nearby emerging flux acts as a piston to generate the shock wave.

The Dvali-Gabadadze-Porrati brane-world model provides a possible approach to address the late-time cosmic acceleration. However, it has subsequently been pointed out that a ghost instability will arise on the self-accelerating branch. Here, we carefully investigate whether this ghost problem could be possibly cured by introducing the Gauss-Bonnet term in the five-dimensional bulk action, a natural generalization to the Dvali-Gabadadze-Porrati model. Our analysis is carried out for a background where a de Sitter brane is embedded in an anti-de Sitter bulk. Our result shows that the ghost excitations cannot be avoided even in this modified model.

We describe a successful effort to produce a laser comb around 1.55 $\mu$m in the astronomical H band using a method based on a line-referenced, electro-optical-modulation frequency comb. We discuss the experimental setup, laboratory results, and proof of concept demonstrations at the NASA Infrared Telescope Facility (IRTF) and the Keck-II telescope. The laser comb has a demonstrated stability of $<$ 200 kHz, corresponding to a Doppler precision of ~0.3 m/s. This technology, when coupled with a high spectral resolution spectrograph, offers the promise of $<$1 m/s radial velocity precision suitable for the detection of Earth-sized planets in the habitable zones of cool M-type stars.

The Eddington-inspired-Born-Infeld scenario (EiBI) can prevent the Big Bang singularity for a matter content whose equation of state is constant and positive. In a recent paper we showed that, on the contrary, it is impossible to smooth a big rip in the EiBI setup. In fact the situations are still different for other singularities. In this paper we show that a big freeze singularity in general relativity (GR) can in some cases be smoothed to a sudden or a type IV singularity under the EiBI scenario. Similarly, a sudden or a type IV singularity in GR can be replaced in some regions of the parameter space by a type IV singularity or a loitering behaviour, respectively, in the EiBI framework. Furthermore, we find that the auxiliary metric related to the physical connection usually has a smoother behaviour than that based on the physical metric. In addition, we show that bound structures close to a big rip or a little rip will be destroyed before the advent of the singularity and will remain bound close to a sudden, big freeze or type IV singularity. We then constrain the model following a cosmographic approach, which is well-known to be model-independent, for a given Friedmann-Lema\^itre-Robertson-Walker geometry. It turns out that among the various past or present singularities, the cosmographic analysis can pick up the physical region that determines the occurrence of a type IV singularity or a loitering effect in the past. Moreover, to determine which of the future singularities or doomsdays is more probable, observational constraints on higher order cosmographic parameters is required.

We study the impact of gravitational slingshot on the distribution of cold dark matter in early and modern era galaxies. Multiple gravitational encounters of a lower mass dark matter particle with massive baryonic astrophysical bodies would lead to an average energy gain for the dark matter, similar to second order Fermi acceleration. We calculate the average energy gain and model the integrated effect on the dark matter profile. We find that such slingshot effect was most effective in the early history of galaxies where first generation stars were massive, which smeared the dark matter distribution at the galactic center and flattened it from an initial cusp profile. On the other hand, slingshot is less effective after the high mass first generation stars and stellar remnants are no longer present. Our finding may help to resolve the cusp-core problem, and we discuss implications for the existing observation-simulation discrepancies and phenomena related to galaxy mergers.

We report the observation of the first two harmonics of the horizontally polarized kink waves excited in a coronal loop system lying at south-east of AR 11719 on 2013 April 11. The detected periods of the fundamental mode ($P_1$), its first overtone ($P_2$) in the northern half, and that in the southern one are $530.2 \pm 13.3$, $300.4 \pm 27.7$, and $334.7 \pm 22.1$ s, respectively. The periods of the first overtone in the two halves are the same considering uncertainties in the measurement. We estimate the average electron density, temperature, and length of the loop system as $(5.1 \pm 0.8) \times 10^8$ cm$^{-3}$, $0.65 \pm 0.06$ MK, and $203.8 \pm 13.8$ Mm, respectively. As a zeroth order estimation, the magnetic field strength, $B = 8.2 \pm 1.0$ G, derived by the coronal seismology using the fundamental kink mode matches with that derived by a potential field model. The extrapolation model also shows the asymmetric and nonuniform distribution of the magnetic field along the coronal loop. Using the amplitude profile distributions of both the fundamental mode and its first overtone, we observe that the antinode positions of both the fundamental mode and its first overtone shift towards the weak field region along the coronal loop. The results indicate that the density stratification and the temperature difference effects are larger than the magnetic field variation effect on the period ratio. On the other hand, the magnetic field variation has a greater effect on the eigenfunction of the first overtone than the density stratification does for this case.

The Born-Infeld determinantal gravity has been recently proposed as a way to smooth the Big Bang singularity. This theory is formulated on the Weitzenbock space-time and the teleparallel representation is used instead of the standard Riemannian representation. We find that although this theory is shown to be singularity-free for certain region of the parameter space in which the divergence of the Hubble rate at the high energy regime is substituted by a de-Sitter stage or a bounce in a Friedmann-Lemaitre-Robertson-Walker universe, cosmological singularities such as Big Rip, Big Bang, Big Freeze, and Sudden singularities can emerge in other regions of the configuration space of the theory. We also show that all these singular events exist even though the Universe is filled with a perfect fluid with a constant equation of state.

From an elevation of ~38 km, the balloon-borne ANtarctic Impulsive Transient Antenna (ANITA) is designed to detect the up-coming radio frequency (RF) signal resulting from a sub-surface neutrino-nucleon collision. Although no neutrinos have been discovered thus far, ANITA is nevertheless the only experiment to self-trigger on radio frequency emissions from cosmic-ray induced atmospheric air showers. In the majority of those cases, down-coming RF signals are observed via their reflection from the Antarctic ice sheet and back up to the ANITA interferometer. Estimating the energy scale of the incident cosmic rays therefore requires an estimate of the fractional power reflected at the air-ice interface. Similarly, inferring the energy of neutrinos interacting in-ice from observations of the upwards-directed signal refracting out to ANITA also requires consideration of signal coherence across the interface. By comparing the direct Solar RF signal intensity measured with ANITA to the surface-reflected Solar signal intensity, as a function of incident elevation angle relative to the surface {\Theta}, we estimate the power reflection coefficients R({\Theta}). We find general consistency between our average measurements and the values of R({\Theta}) expected from the Fresnel equations, separately for horizontal- vs. vertical-polarizations.

We present the first stereoscopic and Doppler observations of simultaneous transverse oscillations of a prominence and a filament and longitudinal oscillation of another filament launched by a single shock wave. Using H-alpha Doppler observations, we derive the three-dimensional oscillation velocities at different heights along the prominence axis. The results indicate that the prominence has a larger oscillation amplitude and damping time at higher altitude, but the periods at different heights are the same (i.e., 13.5 minutes). This suggests that the prominence oscillates like a linear vertical rigid body with one end anchored on the Sun. One of the filaments shows weak transverse oscillation after the passing of the shock, which is possibly due to the low altitude of the filament and the weakening (due to reflection) of the shock wave before the interaction. Large amplitude longitudinal oscillation is observed in the other filament after the passing of the shock wave. The velocity amplitude and period are about 26.8 km/s and 80.3 minutes, respectively. We propose that the orientation of a filament or prominence relative to the normal vector of the incoming shock should be an important factor for launching transverse or longitudinal filament oscillations. In addition, the restoring forces of the transverse prominence are most likely due to the coupling of gravity and magnetic tension of the supporting magnetic field, while that for the longitudinal filament oscillation is probably the resultant force of gravity and magnetic pressure.

The first galaxies in the Universe are the building blocks of all observed galaxies. We present scaling relations for galaxies forming at redshifts $z \ge 15$ when reionization is just beginning. We utilize the ``Rarepeak' cosmological radiation hydrodynamics simulation that captures the complete star formation history in over 3,300 galaxies, starting with massive Population III stars that form in dark matter halos as small as ~$10^6 M_\odot$. We make various correlations between the bulk halo quantities, such as virial, gas, and stellar masses and metallicities and their respective accretion rates, quantifying a variety of properties of the first galaxies up to halo masses of $10^9 M_\odot$. Galaxy formation is not solely relegated to atomic cooling halos with virial temperatures greater than $10^4$ K, where we find a dichotomy in galaxy properties between halos above and below this critical mass scale. Halos below the atomic cooling limit have a stellar mass -- halo mass relationship $\log M_\star \simeq 3.5 + 1.3\log(M_{\rm vir} / 10^7 M_\odot)$. We find a non-monotonic relationship between metallicity and halo mass for the smallest galaxies. Their initial star formation events enrich the interstellar medium and subsequent star formation to a median of $10^{-2} Z_\odot$ and $10^{-1.5} Z_\odot$, respectively, in halos of total mass $10^7 M_\odot$ that is then diluted by metal-poor inflows, well beyond Population III pre-enrichment levels of $10^{-3.5} Z_\odot$. The scaling relations presented here can be employed in models of reionization, galaxy formation and chemical evolution in order to consider these galaxies forming prior to reionization.

We show a new scheme of nonperturbative pair production by high energy photons ($\omega\gtrsim m$) in a strong external field is achievable at the next high intensity laser experiments. The pair momentum is boosted and for $\omega\gtrsim 1.2m$ the pair yield is increased when the external field is formed by two laser pulses converging at a small angle. These characteristics are nonperturbative in origin and related to the presence of magnetic field in addition to electric field. By enhancing the signal over perturbative backgrounds, these features allow the employment of above-threshold photons $\omega>2m$, which further increases the pair yield. We note the close relation of this photon-pair conversion mechanism to spontaneous pair creation, recommending it as an accessible stepping stone experiment using state-of-the-art or soon-to-be laser technology.

Prominences, or filaments, are a striking phenomenon in the solar atmosphere. Besides their own rich features and dynamics, they are related to many other activities, such as solar flares and coronal mass ejections (CMEs). In the past several years we have been investigating the prominence formation, oscillations, and eruptions through both data analysis and radiative hydrodynamic and magnetohydrodynamic (MHD) simulations. This paper reviews our progress on these topics, which includes: (1) With updated radiative cooling function, the coronal condensation becomes a little faster than previous work; (2) Once a seed condensation is formed, it can grow via siphon flow spontaneously even if the evaporation stops; (3) A scaling law was obtained to relate the length of the prominence thread to various parameters, indicating that higher prominences tend to have shorter threads, which is consistent with the fact that threads are long in active region prominences and short in quiescent prominences; (4) It was proposed that long-time prominence oscillations out of phase might serve as a precursor for prominence eruptions and CMEs; (5) An ensemble of oscillating prominence threads may explain the counter-streaming motion.

Coronal bright points (CBPs) are small-scale and long-duration brightenings in the lower solar corona. They are often explained in terms of magnetic reconnection. We aim to study the sub-structures of a CBP and clarify the relationship among the brightenings of different patches inside the CBP. The event was observed by the X-ray Telescope (XRT) aboard the Hinode spacecraft on 2009 August 22$-$23. The CBP showed repetitive brightenings (or CBP flashes). During each of the two successive CBP flashes, i.e., weak and strong flashes which are separated by $\sim$2 hr, the XRT images revealed that the CBP was composed of two chambers, i.e., patches A and B. During the weak flash, patch A brightened first, and patch B brightened $\sim$2 min later. During the transition, the right leg of a large-scale coronal loop drifted from the right side of the CBP to the left side. During the strong flash, patch B brightened first, and patch A brightened $\sim$2 min later. During the transition, the right leg of the large-scale coronal loop drifted from the left side of the CBP to the right side. In each flash, the rapid change of the connectivity of the large-scale coronal loop is strongly suggestive of the interchange reconnection. For the first time we found reciprocatory reconnection in the CBP, i.e., reconnected loops in the outflow region of the first reconnection process serve as the inflow of the second reconnection process.

Magnetic flux rope (MFR) is a coherent and helical magnetic field structure that is recently found probably to appear as an elongated hot-channel prior to a solar eruption. In this paper, we investigate the relationship between the hot-channel and associated prominence through analyzing a limb event on 2011 September 12. In the early rise phase, the hot-channel was cospatial with the prominence initially. It then quickly expanded, resulting in a separation of the top of the hot-channel from that of the prominence. Meanwhile, both of them experienced an instantaneous morphology transformation from a $\Lambda$ shape to a reversed-Y shape and the top of these two structures showed an exponential increase in height. These features are a good indication for the occurrence of the kink instability. Moreover, the onset of the kink instability is found to coincide in time with the impulsive enhancement of the flare emission underneath the hot-channel, suggesting that the ideal kink instability likely also plays an important role in triggering the fast flare reconnection besides initiating the impulsive acceleration of the hot-channel and distorting its morphology. We conclude that the hot-channel is most likely the MFR system and the prominence only corresponds to the cool materials that are collected in the bottom of the helical field lines of the MFR against the gravity.

We study cosmological perturbation theory within the framework of unimodular gravity. We show that the Lagrangian constraint on the determinant of the metric required by unimodular gravity leads to an extra constraint on the gauge freedom of the metric perturbations. Although the main equation of motion for the gravitational potential remains the same, the shift variable, which is gauge artifact in General Relativity, cannot be set to zero in unimodular gravity. This non-vanishing shift variable affects the propagation of photons throughout the cosmological evolution and therefore modifies the Sachs-Wolfe relation between the relativistic gravitational potential and the microwave temperature anisotropies. However, for adiabatic fluctuations the difference between the result in General Relativity and unimodular gravity is suppressed on large angular scales. Thus, no strong constraints on the theory can be derived.

We calculate the low red-shift Taylor expansion for the luminosity distance for an observer at the center of a spherically symmetric matter inhomogeneity with a non vanishing cosmological constant. We then test the accuracy of the formulas comparing them to the numerical calculation for different cases for both the luminosity distance and the radial coordinate. The formulas can be used as a starting point to understand the general non linear effects of a local inhomogeneity in presence of a cosmological constant, without making any special assumption about the inhomogeneity profile.

The brane-world scenario provides an intriguing possibility to explore the phenomenological cosmology implied by string/M theory. In this paper, we consider a modified Randall-Sundrum single brane model with two natural generalizations: a Gauss-Bonnet term in the five-dimensional bulk action as well as an induced gravity term in the four-dimensional brane action, which are the leading-order corrections to the Randall-Sundrum model. We study the influence of these combined effects on the evolution of the primordial gravitational waves generated during an extreme slow-roll inflation on the brane. The background, for the early inflationary era, is then modeled through a de Sitter brane embedded within an anti-de Sitter bulk. In this framework, we show that both effects tend to suppress the Randall-Sundrum enhancement of the amplitude of the tensor perturbations at relatively high energies. Moreover, the Gauss-Bonnet effect, relative to standard general relativity, will abruptly enhance the tensor-to-scalar ratio and break the standard consistency relation at high energies, which cannot be evaded by invoking the induced gravity effect, even though the induced gravity strength would mildly counterbalance these significant changes at high energies. We note that the brane-world model with or without the induced gravity effect fulfills the consistency relation. Finally, we discuss some implications of the observational constraints.

One longstanding mystery in bio-evolution since Darwin's time is the origin of the Cambrian explosion that happened around 540 million years ago (Mya), where an extremely rapid increase of species occurred. Here we suggest that a nearby GRB event ~500 parsecs away, which should occur about once per 5 Gy, might have triggered the Cambrian explosion. Due to a relatively lower cross section and the conservation of photon number in Compton scattering, a substantial fraction of the GRB photons can reach the sea level and would induce DNA mutations in organisms protected by a shallow layer of water or soil, thus expediting the bio-diversification. This possibility of inducing genetic mutations is unique among all candidate sources for major incidents in the history of bio-evolution. A possible evidence would be the anomalous abundance of certain nuclear isotopes with long half-lives transmuted by the GRB photons in geological records from the Cambrian period. Our notion also imposes constraints on the evolution of exoplanet organisms and the migration of panspermia.

The dynamics of a filament channel are observed with imaging and spectroscopic telescopes before and during the filament eruption on 2011 January 29. The extreme ultraviolet (EUV) spectral observations reveal that there is no EUV counterparts of the H$\alpha$ counter-streamings in the filament channel, implying that the ubiquitous H$\alpha$ counter-streamings found by previous research are mainly due to longitudinal oscillations of filament threads, which are not in phase between each other. However, there exist larger-scale patchy counter-streamings in EUV along the filament channel from one polarity to the other, implying that there is another component of uni-directional flow (in the range of $\pm$10 km s$^{-1}$) inside each filament thread in addition to the implied longitudinal oscillation. Our results suggest that the flow direction of the larger-scale patchy counter-streaming plasma in the EUV is related to the intensity of the plage or active network, with the upflows being located at brighter areas of the plage and downflows at the weaker areas. Besides, we propose a new method to determine the chirality of an erupting filament based on the skewness of the conjugate filament drainage sites. It suggests that the right-skewed drainage corresponds to sinistral chirality, whereas the left-skewed drainage corresponds to dextral chirality.

Recent years have seen tremendous progress in our understanding of the cosmos, which in turn points to even deeper questions to be further addressed. Concurrently the laser technology has undergone dramatic revolutions, providing exciting opportunity for science applications. History has shown that the symbiosis between direct observations and laboratory investigation is instrumental in the progress of astrophysics. We believe that this remains true in cosmology. Current frontier phenomena related to particle astrophysics and cosmology typically involve one or more of the following conditions: (1) extremely high energy events; (2) very high density, high temperature processes; (3) super strong field environments. Laboratory experiments using high intensity lasers can calibrate astrophysical observations, investigate underlying dynamics of astrophysical phenomena, and probe fundamental physics in extreme limits. In this article we give an overview of the exciting prospect of laser cosmology. In particular, we showcase its unique capability of investigating frontier cosmology issues such as cosmic accelerator and quantum gravity.

High-resolution H$\alpha$ observations indicate that filaments consist of an assembly of thin threads. In quiescent filaments, the threads are generally short, whereas in active region filaments, the threads are generally long. In order to explain these observational features, we performed one-dimensional radiative hydrodynamic simulations of filament formation along a dipped magnetic flux tube in the framework of the chromospheric evaporation-coronal condensation model. The geometry of a dipped magnetic flux tube is characterized by three parameters, i.e., the depth ($D$), the half-width ($w$), and the altitude ($h$) of the magnetic dip. The parameter survey in the numerical simulations shows that allowing the filament thread to grow in 5 days, the maximum length ($L_{th}$) of the filament thread increases linearly with $w$, and decreases linearly with $D$ and $h$. The dependence is fitted into a linear function $L_{th}=0.84w-0.88D-2.78h+17.31$ Mm. Such a relation can qualitatively explain why quiescent filaments have shorter threads and active region filaments have longer threads.

The luminosity distance can be used to determine the properties of large scale structure around the observer. To this purpose we develop a new inversion method to map luminosity distance to a LTB metric based on the use of the exact analytical solution for Einstein equations. The main advantages of this approach are an improved numerical accuracy and stability, an exact analytical setting of the initial conditions for the differential equations which need to be solved and the validity for any sign of the functions determining the LTB geometry. Given the fully analytical form of the differential equations, this method also simplifies the calculation of the red-shift expansion around the apparent horizon point where the numerical solution becomes unstable. We test the method by inverting the supernovae Ia luminosity distance function corresponding to the the best fit $\Lambda CDM$ model. We find that only a limited range of initial conditions is compatible with observations, or a transition from red to blue shift can occur at relatively low redshift. Despite LTB solutions without a cosmological constant have been shown not to be compatible with all different set of available observational data, those studies normally fit data assuming a special functional ansatz for the inhomogeneity profile, which often depend only on few parameters. Inversion methods on the contrary are able to fully explore the freedom in fixing the functions which determine a LTB solution. Another important possible application is not about LTB solutions as cosmological models, but rather as tools to study the effects on the observations made by a generic observer located in an inhomogeneous region of the Universe where a fully non perturbative treatment involving exact solutions of Einstein equations is required.

We perform a model independent study of the neutrino momentum distribution at freeze-out, treating the freeze-out temperature as a free parameter. Our results imply that measurement of neutrino reheating, as characterized by the measurement of the effective number of neutrinos $N_\nu$, amounts to the determination of the neutrino kinetic freeze-out temperature within the context of the standard model of particle physics where the number of neutrino flavors is fixed and no other massless (fractional) particles arise. At temperatures on the order of the neutrino mass, we show how cosmic background neutrino properties i.e. energy density, pressure, particle density, are modified in a physically consistent way as a function of neutrino mass and $N_\nu$.

The flavor composition of ultra-high energy cosmic neutrinos (UHECN) carries precious information about the physical properties of their sources, the nature of neutrino oscillations and possible exotic physics involved during the propagation. Since UHECN with different incoming directions would propagate through different amounts of matter in Earth and since different flavors of charged leptons produced in the neutrino-nucleon charged-current (CC) interaction would have different energy-loss behaviors in the medium, measurement of the angular distribution of incoming events by a neutrino observatory can in principle be employed to help determine the UHECN flavor ratio. In this paper we report on our investigation of the feasibility of such an attempt. Simulations were performed, where the detector configuration was based on the proposed Askaryan Radio Array (ARA) Observatory at the South Pole, to investigate the expected event-direction distribution for each flavor. Assuming $\nu_{\mu}$-$\nu_{\tau}$ symmetry and invoking the standard oscillation and the neutrino decay scenarios, the probability distribution functions (PDF) of the event directions are utilized to extract the flavor ratio of cosmogenic neutrinos on Earth. The simulation results are summarized in terms of the probability of flavor ratio extraction and resolution as functions of the number of observed events and the angular resolution of neutrino directions. We show that it is feasible to constrain the UHECN flavor ratio using the proposed ARA Observatory.

Mergers of binary neutron stars create conditions of supranuclear density $n\gtrsim n_{\rm nuc}\simeq 0.17 {\rm fm}^{-3}$ and moderate temperature $50\lesssim T \lesssim 90 {\rm MeV}$. These events thus probe a sensitive region of the density-temperature phase diagram of QCD matter. We study photon production by the QCD conformal anomaly for a signal of a possible transition to quark degrees of freedom during the merger. We discuss energy loss due to photon radiation as a cooling mechanism that is sensitive to the bulk viscosity and thermal conductivity of the quark matter.

Within the framework of RG improved inflationary cosmology motivated by asymptotically safe gravity, we study the dynamics of a scalar field which can be interpreted as the Higgs field. The background trajectories of this model can provide sufficient inflationary e-folds and a graceful exit to a radiation dominated phase. We study the possibility of generating primordial curvature perturbations through the Standard Model Higgs boson. This can be achieved under finely tuned parameter choices by making use of the modulated reheating mechanism. The primordial non-gaussianity is expected to be sizable in this model. Though tightly constrained by the newly released Planck CMB data, this model provides a potentially interesting connection between collider and early universe physics.

Recent years have seen dramatic progress in cosmology and particle astrophysics. So much so that anyone who dares to offer an overview would certainly risk him- or herself for being incomplete and biased at best, and even incorrect due to the author's limited expertise. It is with this understanding and excuse that I risk myself in offering this review. After a brief summary of Planck mission's first results, I highlight some selected theoretical and experimental advancement in dark energy, dark matter, and cosmic neutrinos research. It is hoped that with a glance through these exciting development, one would be convinced that we are now a step closer to the ultimate understanding of our universe, while major breakthroughs are still required.

There exists dark energy models that predict the occurrence of "little rip". At the point of little rip the Hubble rate and its cosmic time derivative approach infinity, which is quite similar to the big rip singularity except that the former happens at infinite future while the latter at a finite cosmic time; both events happen in the future and at high energies. In the case of the big rip, a combination of ultra-violet and infra-red effects can smooth its doomsday. We therefore wonder if the little rip can also be smoothed in a similar way. We address the ultra-violet and infra-red effects in general relativity through a brane-world model with a Gauss-Bonnet term in the bulk and an induced gravity term on the brane. We find that the little rip is transformed in this case into a sudden singularity, or a "big brake". Even though the big brake is smoother than the little rip in that the Hubble rate is finite at the event, the trade-off is that it takes place sooner, at a finite cosmic time. In our estimate, the big brake would happen at roughly 1300Gyr.

We revisit the method of using the photon-graviton conversion mechanism in the presence of the external magnetic field to probe small-scale primordial magnetic fields that may exist between the last scattering surface and present. Specifically, we investigate impacts on the conversion efficiency due to the presence of matter, including the plasma collective effect and the atomic polarizability. In general, these effects tend to reduce the conversion probability. Under this more realistic picture and based on the precision of COBE's measurement of CMB (cosmic microwave background) blackbody spectrum, we find an upper bound for the primordial magnetic field strength, B < 30G, at the time of recombination. Although at present the bound based on the photon-graviton conversion mechanism is not as tight as that obtained by the direct use of CMB temperature anisotropy, it nevertheless provides an important independent constraint on primordial magnetic fields and at epochs in addition to the recombination. The bound can be significantly improved if the CMB blackbody spectrum measurement becomes more precise in future experiments such as PIXIE.

Teleparallel theories of gravity have a long history. They include a special case referred to as the Teleparallel Equivalent of General Relativity (TEGR, aka GR$_{\|}$). Recently this theory has been generalized to f(T) gravity. Tight constraints from observations suggest that f(T) gravity is not as robust as initially hoped. This might hint at hitherto undiscovered problems at the theoretical level. In this work, we point out that a generic f(T) theory can be expected to have certain problems including superluminal propagating modes, the presence of which can be revealed by using the characteristic equations that govern the dynamics in f(T) gravity and/or the Hamiltonian structure of the theory via Dirac constraint analysis. We use several examples from simpler gauge field theories to explain how such superluminal modes could arise. We also point out problems with the Cauchy development of a constant time hypersurface in FLRW spacetime in f(T) gravity. The time evolution from a FLRW (and as a special case, Minkowski spacetime) initial condition is not unique.

A new multi-wavelength solar telescope, Optical and Near-infrared Solar Eruption Tracer (ONSET) of Nanjing University, was constructed, being fabricated by Nanjing Institute of Astronomical Optics & Technology and run in cooperation with Yunnan Astronomical Observatory. ONSET is able to observe the Sun in three wavelength windows: He {\small I} 10830 \AA, H$\alpha$, and white-light at 3600 \AA and 4250 \AA, which are selected in order to obtain the dynamics in the corona, chromosphere, and the photosphere simultaneously. Full-disk or partial-disk solar images with a field of 10 arcmin at three wavelengths can be obtained nearly simultaneously. It is designed to trace solar eruptions with high spatial and temporal resolutions. This telescope was installed at a new solar observing site near Fuxian Lake in Yunnan Province, southwest China. The site is located at E102N24, with an altitude of 1722 m. The seeing is stable and of high quality. We give a brief description of the scientific objectives and the basic structure of the ONSET. Some preliminary results are also presented.

We present a simple toy model corresponding to a network of frustrated topological defects of domain walls or cosmic strings that exist previous to the standard slow-roll inflationary era of the universe. Such a network (i) can produce a slower inflationary era than that of the standard scenario if it corresponds to a network of frustrated domain walls or (ii) can induce a vanishing universal acceleration; i.e., the universe would expand at a constant speed, if it corresponds to a network of frustrated cosmic strings red. Those features are phenomenologically modeled by a Chaplygin gas that can interpolate between a network of frustrated topological defects and a de Sitter-like or a power-law inflationary era. We show that this scenario can alleviate the quadruple anomaly of the cosmic microwave background spectrum. Using the method of the Bogoliubov coefficients, we obtain the spectrum of the gravitational waves as would be measured today for the whole range of frequencies. We comment on the possible detection of this spectrum by the planned detectors like BBO and DECIGO.

There exists the "entropy problem" of the early universe, that is, why did the universe begin with an extremely low entropy and how did it evolve into such high entropy at late times? It has been long believed that inflation cannot be the solution since it requires an extremely low entropy to ever occur. However, we point out that since the inflation is always accompanied with a horizon, the correct probability of inflation is associated with the quantum entanglement entropy, which should in principle be larger than what considered previously. This motivates us to reexamine the issue by computing the evolution of the cosmological entanglement entropy in the early universe. We invoke a toy model of nonlinear generalized Chaplygin gas (GCG), which has the advantage of providing a smooth and unitary transition between the inflation epoch and the radiation dominant era. We found that soon after the onset of the inflation, the total entanglement entropy rapidly decreases to a minimum, and it rises monotonically afterwards throughout the remainder of the inflation and the radiation epochs. This indicates that the universe does not need to begin with an extremely low entropy; its smallness can be naturally induced by the dynamics of inflation itself. We believe that our computation largely captures the essential feature of entropy evolution and can provide us insights beyond the toy model.

We show that the Stephenson-Kilmister-Yang (SKY) equation combined with Camenzind's matter current term naturally provides the cosmological constant and dark radiation as integration constants of the SKYC field equation. To characterize the property of the dark radiation, we develop a method to separate it from the ordinary radiation. We found a special property of Camenzind's matter current, namely that the solution space for radiation in fact belongs to that of the vacuum solution of SKY equation. We also found that his matter current does not obey the conservation condition suggested by Kilmister. Finally, we discuss the possible role of dark radiation emergent from the SKYC theory in recent cosmic-microwave-background observations and its implications to the inflation scenario.

We clarify in a quantitative way the impact that distinct chemical $T_c$ and kinetic $T_k$ freeze-out temperatures have on the reduction of the neutrino fugacity $\Upsilon_\nu$ below equilibrium, i.e. $\Upsilon_\nu<1$, and the increase of the neutrino temperature $T_\nu$ via partial reheating. We establish the connection between $\Upsilon_\nu$ and $T_k$ via the modified reheating relation $T_\nu(\Upsilon_\nu)/T_\gamma$, where $T_\gamma$ is the temperature of the background radiation. Our results demonstrate that one must introduce the chemical nonequilibrium parameter, i.e., the fugacity, $\Upsilon_\nu$, as an additional standard cosmological model parameter in the evaluation of CMB fluctuations as its value allows measurement of $T_k$.

It is found that both microflare-sized impulsive heating at one leg of the loop and a suddenly imposed velocity perturbation can propel the prominence to oscillate along the magnetic dip. An extensive parameter survey results in a scaling law, showing that the period of the oscillation, which weakly depends on the length and height of the prominence, and the amplitude of the perturbations, scales with $\sqrt{R/g_\odot}$, where $R$ represents the curvature radius of the dip, and $g_\odot$ is the gravitational acceleration of the Sun. This is consistent with the linear theory of a pendulum, which implies that the field-aligned component of gravity is the main restoring force for the prominence longitudinal oscillations, as confirmed by the force analysis. However, the gas pressure gradient becomes non-negligible for short prominences. The oscillation damps with time in the presence of non-adiabatic processes. Compared to heat conduction, the radiative cooling is the dominant factor leading to the damping. A scaling law for the damping timescale is derived, i.e., $\tau\sim l^{1.63} D^{0.66}w^{-1.21}v_{0}^{-0.30}$, showing strong dependence on the prominence length $l$, the geometry of the magnetic dip (characterized by the depth $D$ and the width $w$), and the velocity perturbation amplitude $v_0$. The larger the amplitude, the faster the oscillation damps. It is also found that mass drainage significantly reduces the damping timescale when the perturbation is too strong.

We developed a method to automatically detect and trace solar filaments in H\alpha\ full-disk images. The program is able not only to recognize filaments and determine their properties, such as the position, the area, the spine, and other relevant parameters, but also to trace the daily evolution of the filaments. The program consists of three steps: First, preprocessing is applied to correct the original images; Second, the Canny edge-detection method is used to detect filaments; Third, filament properties are recognized through the morphological operators. To test the algorithm, we applied it to the observations from the Mauna Loa Solar Observatory (MLSO), and the program is demonstrated to be robust and efficient. H\alpha\ images obtained by MLSO from 1998 to 2009 are analyzed, and a butterfly diagram of filaments is obtained. It shows that the latitudinal migration of solar filaments has three trends in the Solar Cycle 23: The drift velocity was fast from 1998 to the solar maximum; After the solar maximum, it became relatively slow. After 2006, the migration became divergent, signifying the solar minimum. About 60% filaments with the latitudes larger than $50^{\circ}$ migrate towards the polar regions with relatively high velocities, and the latitudinal migrating speeds in the northern and the southern hemispheres do not differ significantly in the Solar Cycle 23.

We introduce a new model of gamma ray burst (GRB) that explains its observed prompt signals, namely, its primary quasi-thermal spectrum and high energy tail. This mechanism can be applied to either assumption of GRB progenitor: coalescence of compact objects or hypernova explosion. The key ingredients of our model are: (1) The initial stage of a GRB is in the form of a relativistic quark-gluon plasma "lava"; (2) The expansion and cooling of this lava results in a QCD phase transition that induces a sudden gravitational stoppage of the condensed non-relativistic baryons and form a hadrosphere; (3) Acoustic shocks and Alfven waves (magnetoquakes) that erupt in episodes from the epicenter efficiently transport the thermal energy to the hadrospheric surface and induce a rapid detachment of leptons and photons from the hadrons; (4) The detached $e^+e^-$ and $\gamma$ form an opaque, relativistically hot leptosphere, which expands and cools to $T \sim mc^2$, or 0.5 MeV, where $e^+e^- \to 2\gamma$ and its reverse process becomes unbalanced, and the GRB photons are finally released; (5) The "mode-conversion" of Alfven waves into electromagnetic waves in the leptosphere provides a "snowplow" acceleration and deceleration that gives rise to both the high energy spectrum of GRB and the erosion of its thermal spectrum down to a quasi-thermal distribution. According to this model, the observed GRB photons should have a redshifted peak frequency at $E_p \sim \Gamma(1 + \beta/2)mc^2/(1 + z),$ where $\Gamma\sim {\cal{O}}(1)$ is the Lorentz factor of the bulk flow of the lava, which may be determined from the existing GRB data.

The origin of ultra-high energy cosmic rays (UHECR) has been widely regarded as one of the major questions in the frontiers of particle astrophysics. Gamma ray bursts (GRB), the most violent explosions in the universe second only to the Big Bang, have been a popular candidate site for UHECR productions. The recent IceCube report on the non-observation of GRB induced neutrinos therefore attracts wide attention. This dilemma requires a resolution: either the assumption of GRB as UHECR accelerator is to be abandoned or the expected GRB induced neutrino yield was wrong. It has been pointed out that IceCube has overestimated the neutrino flux at GRB site by a factor of $\sim 5$. In this paper we point out that, in addition to the issue of neutrino production at source, the neutrino oscillation and the possible neutrino decay during their flight from GRB to Earth should further reduce the detectability of IceCube, which is most sensitive to the muon-neutrino flavor as far as point-source identification is concerned. Specifically, neutrino oscillation will reduce the muon-neutrino flavor ratio from 2/3 per neutrino at GRB source to 1/3 on Earth, while neutrino decay, if exists and under the assumption of normal hierarchy of mass eigenstates, would result in a further reduction of muon-neutrino ratio to 1/8. With these in mind, we note that there have been efforts in recent years in pursuing other type of neutrino telescopes based on Askaryan effect, which can in principle observe and distinguish all three flavors with comparable sensitivities. Such new approach may therefore be complementary to IceCube in shedding more lights on this cosmic accelerator question.

The Swift Gamma-ray Burst (GRB) observatory responds to GRB triggers with optical observations in ~ 100 s, but cannot respond faster than ~ 60 s. While some ground-based telescopes respond quickly, the number of sub-60 s detections remains small. In mid- to late-2013, the Ultra-Fast Flash Observatory-Pathfinder is to be launched on the Lomonosov spacecraft to investigate early optical GRB emission. This pathfinder mission is necessarily limited in sensitivity and event rate; here we discuss a next generation rapid-response space observatory. We list science topics motivating our instruments, those that require rapid optical-IR GRB response, including: A survey of GRB rise shapes/times, measurements of optical bulk Lorentz factors, investigation of magnetic dominated (vs. non-magnetic) jet models, internal vs. external shock origin of prompt optical emission, the use of GRBs for cosmology, and dust evaporation in the GRB environment. We also address the impacts of the characteristics of GRB observing on our instrument and observatory design. We describe our instrument designs and choices for a next generation observatory as a second instrument on a low-earth orbit spacecraft, with a 120 kg instrument mass budget. Restricted to relatively modest mass and power, we find that a coded mask X-ray camera with 1024 cm2 of detector area could rapidly locate about 64 GRB triggers/year. Responding to the locations from the X-ray camera, a 30 cm aperture telescope with a beam-steering system for rapid (~ 1 s) response and a near-IR camera should detect ~ 29 GRB, given Swift GRB properties. Am additional optical camera would give a broadband optical-IR slope, allowing dynamic measurement of dust extinction at the source, for the first time.

We study non-linear massive gravity in the spherically symmetric context. Our main motivation is to investigate the effect of helicity-0 mode which remains elusive after analysis of cosmological perturbation around an open Friedmann-Lemaitre-Robertson-Walker (FLRW) universe. The non-linear form of the effective energy-momentum tensor stemming from the mass term is derived for the spherically symmetric case. Only in the special case where the area of the two sphere is not deviated away from the FLRW universe, the effective energy momentum tensor becomes completely the same as that of cosmological constant. This opens a window for discriminating the non-linear massive gravity from general relativity (GR). Indeed, by further solving these spherically symmetric gravitational equations of motion in vacuum to the linear order, we obtain a solution which has an arbitrary time-dependent parameter. In GR, this parameter is a constant and corresponds to the mass of a star. Our result means that Birkhoff's theorem no longer holds in the non-linear massive gravity and suggests that energy can probably be emitted superluminously (with infinite speed) on the self-accelerating background by the helicity-0 mode, which could be a potential plague of this theory.

Supernovae observations strongly support the presence of a cosmological constant, but its value, which we will call apparent, is normally determined assuming that the universe can be accurately described by a homogeneous model. Even in the presence of a cosmological constant we cannot exclude nevertheless the presence of a small local inhomogeneity which could affect the apparent value of the cosmological constant. Neglecting the presence of the inhomogeneity can in fact introduce a systematic misinterpretation of cosmological data, leading to the distinction between an apparent and the true value of the cosmological constant. But is such a difference distinguishable? Recently we set out to model the local inhomogeneity with a {\Lambda}LTB solution and computed the relation between the apparent and the true value of the cosmological constant. In this essay we reproduce the essence of our model with the emphasis on its physical implications.

We analyse multiwavelength observations of an M2.9/1N flare that occurred in AR NOAA 11112 on 16 October 2010. AIA 211 {\AA} EUV images reveal the presence of a faster coronal wave (decelerating from ~1390 to ~830 km/s) propagating ahead of a slower wave (decelerating from ~416 to ~166 km/s) towards the western limb. The dynamic radio spectrum from Sagamore Hill radio telescope shows the presence of metric type II radio burst, which reveals the presence of a coronal shock wave (speed~800 km/s). The speed of the faster coronal wave derived from AIA 211 {\AA} images is found to be comparable to the coronal shock speed. AIA 171 {\AA} high-cadence observations showed that a coronal loop, which was located at the distance of ~0.32 Rs to the west of the flaring region, started to oscillate by the end of the impulsive phase of the flare. The results indicate that the faster coronal wave may be the first driver of the transversal oscillations of coronal loop. As the slower wave passed through the coronal loop, the oscillations became even stronger. There was a plasmoid eruption observed in EUV and a white-light CME was recorded, having velocity of ~340-350 km/s. STEREO 195 {\AA} images show an EIT wave, propagating in the same direction of the lower-speed coronal wave observed in AIA, but decelerating from ~320 to ~254 km/s. These observations reveal the co-existence of both waves (i.e. coronal Moreton and EIT waves), and type II radio burst seems to be associated with the coronal Moreton wave.

We consider a generalization of the Randall-Sundrum single brane-world scenario (RS2). More precisely, the generalization is described through curvature corrections, corresponding to a Gauss-Bonnet term in the bulk and a Hilbert-Einstein term, as well as the strength of the induced gravity term, on the brane. We are mainly interested in analyzing the early inflationary era of the brane, which we model within the extreme slow-roll limit, i.e., under a de Sitter like brane inflation, where the inflaton field is confined on the brane. We compute the scalar perturbations in this model and compare our results with those previously obtained for the RS2 scenario with and without an induced gravity term on the brane or a Gauss-Bonnet term in the bulk. The amplitude of the scalar perturbations is decreased as compared with a pure RS2 model. In addition, the effect from the Gauss-Bonnet correction in an induced gravity brane-world model is to decrease the amplitude of the scalar perturbations, and a similar result is obtained for the induced gravity effect in a Gauss-Bonnet brane-world. In general, in the high energy limit the amplitude is highly suppressed by the Gauss-Bonnet effect. Finally, we constrain the model using the latest WMAP7 data.

We consider the implications of Lorentz-invariance violation (LIV) on cosmogenic neutrino observations, with particular focus on the constraints imposed on several well-developed models for ultra-high energy cosmogenic neutrino production by recent results from the Antarctic Impulsive Transient Antenna (ANITA) long-duration balloon payload, and Radio Ice Cherenkov Experiment (RICE) at the South Pole. Under a scenario proposed originally by Coleman and Glashow, each lepton family may attain maximum velocities that can exceed the speed of light, leading to energy-loss through several interaction channels during propagation. We show that future observations of cosmogenic neutrinos will provide by far the most stringent limit on LIV in the neutrino sector. We derive the implied level of LIV required to suppress observation of predicted fluxes from several mainstream cosmogenic neutrino models, and specifically those recently constrained by the ANITA and RICE experiments. We simulate via detailed Monte Carlo code the propagation of cosmogenic neutrino fluxes in the presence of LIV-induced energy losses. We show that this process produces several detectable effects in the resulting attenuated neutrino spectra, even at LIV-induced neutrino superluminality of (u_{\nu}-c)/c ~ 10^{-26}, about 13 orders of magnitude below current bounds.

We confront two types of phantom dark energy potential with observational data. The models we consider are the power-law potential, $V\propto {\phi}^{\mu}$, and the exponential potential, $V\propto \exp({\lambda}{\phi}/{M_P})$. We fit the models to the latest observations from SN-Ia, CMB and BAO, and obtain tight constraints on parameter spaces. Furthermore, we apply the goodness-of-fit and the information criteria to compare the fitting results from phantom models with that from the cosmological constant and the quintessence models presented in our previous work. The results show that the cosmological constant is statistically most preferred, while the phantom dark energy fits slightly better than the quintessence does.

We report quadrature observations of an extreme-ultraviolet (EUV) wave event on 2011 January 27 obtained by the Extreme Ultraviolet Imager (EUVI) onboard \emph{Solar Terrestrial Relations Observatory} (\emph{STEREO}), and the Atmospheric Imaging Assembly (AIA) onboard the \emph{Solar Dynamics Observatory} (\emph{SDO}). Two components are revealed in the EUV wave event. A primary front is launched with an initial speed of $\sim$440 km s$^{-1}$. It appears significant emission enhancement in the hotter channel but deep emission reduction in the cooler channel. When the primary front encounters a large coronal loop system and slows down, a secondary much fainter front emanates from the primary front with a relatively higher starting speed of $\sim$550 km s$^{-1}$. Afterwards the two fronts propagate independently with increasing separation. The primary front finally stops at a magnetic separatrix, while the secondary front travels farther before it fades out. In addition, upon the arrival of the secondary front, transverse oscillations of a prominence are triggered. We suggest that the two components are of different natures. The primary front belongs to a non-wave coronal mass ejection (CME) component, which can be reasonably explained with the field-line stretching model. The multi-temperature behavior may be caused by considerable heating due to the nonlinear adiabatic compression on the CME frontal loop. For the secondary front, most probably it is a linear fast-mode magnetohydrodynamic (MHD) wave that propagates through a medium of the typical coronal temperature. X-ray and radio data provide us with complementary evidence in support of the above scenario.

Microflares are small activities in solar low atmosphere, some are in the low corona, and others in the chromosphere. Observations show that some of the microflares are triggered by magnetic reconnection between emerging flux and a pre-existing background magnetic field. We perform 2.5D compressible resistive MHD simulations of magnetic reconnection with gravity considered. The background magnetic field is a canopy-type configuration which is rooted at the boundary of the solar supergranule. By changing the bottom boundary conditions in the simulation, new magnetic flux emerges up at the center of the supergranule and reconnects with the canopy-type magnetic field. We successfully simulate the coronal and chromospheric microflares, whose current sheets are located at the corona and the chromosphere, respectively. The microflare of coronal origin has a bigger size and a higher temperature enhancement than that of chromospheric origin. In the microflares of coronal origin, we also found a hot jet ($\sim$$1.8 \times 10^6$ K), which is probably related to the observational EUV/SXR jets, and a cold jet ($\sim$$10^4$ K), which is similar to the observational H$\alpha$/Ca surges, whereas there is only an H$\alpha$/Ca bright point in the microflares of chromospheric origin. The study of parameter dependence shows that the size and strength of the emerging magnetic flux are the key parameters which determine the height of the reconnection location, and further determine the different observational features of the microflares.

A new code, named MAP, is written in Fortran language for magnetohydrodynamics (MHD) calculation with the adaptive mesh refinement (AMR) and Message Passing Interface (MPI) parallelization. There are several optional numerical schemes for computing the MHD part, namely, modified Mac Cormack Scheme (MMC), Lax-Friedrichs scheme (LF) and weighted essentially non-oscillatory (WENO) scheme. All of them are second order, two-step, component-wise schemes for hyperbolic conservative equations. The total variation diminishing (TVD) limiters and approximate Riemann solvers are also equipped. A high resolution can be achieved by the hierarchical block-structured AMR mesh. We use the extended generalized Lagrange multiplier (EGLM) MHD equations to reduce the non-divergence free error produced by the scheme in the magnetic induction equation. The numerical algorithms for the non-ideal terms, e.g., the resistivity and the thermal conduction, are also equipped in the MAP code. The details of the AMR and MPI algorithms are described in the paper.

Filament longitudinal oscillations have been observed on the solar disk in H$\alpha$. We intend to find an example of the longitudinal oscillations of a prominence, where the magnetic dip can be seen directly, and examine what is the restoring force of such kind of oscillations. We carry out a multiwavelength data analysis of the active region prominence oscillations above the western limb on 2007 February 8. Besides, we perform a one-dimensional hydrodynamic simulation of the longitudinal oscillations. The high-resolution observations by Hinode/SOT indicate that the prominence, seen as a concave-inward shape in lower-resolution Extreme Ultraviolet (EUV) images, actually consists of many concave-outward threads, which is indicative of the existence of magnetic dips. After being injected into the dip region, a bulk of prominence material started to oscillate for more than 3.5 hours, with the period being 52 min. The oscillation decayed with time, with the decay timescale being 133 min. Our hydrodynamic simulation can well reproduce the oscillation period, but the damping timescale in the simulation is 1.5 times as long as the observations. The results clearly show the prominence longitudinal oscillations around the dip of the prominence and our study suggests that the restoring force of the longitudinal oscillations might be the gravity. Radiation and heat conduction are insufficient to explain the decay of the oscillations. Other mechanisms, such as wave leakage and mass accretion, have to be considered. The possible relation between the longitudinal oscillations and the later eruption of a prominence thread, as well as a coronal mass ejection (CME), is also discussed.

Gravitational waves are predicted by the general theory of relativity. It has been shown that gravitational waves have a nonlinear memory, displacing test masses permanently. This is called the Christodoulou memory. We proved that the electromagnetic field contributes at highest order to the nonlinear memory effect of gravitational waves, enlarging the permanent displacement of test masses. In experiments like LISA or LIGO which measure distances of test masses, the Christodoulou memory will manifest itself as a permanent displacement of these objects. It has been suggested to detect the Christodoulou memory effect using radio telescopes investigating small changes in pulsar's pulse arrival times. The latter experiments are based on present-day technology and measure changes in frequency. In the present paper, we study the electromagnetic Christodoulou memory effect and compute it for binary neutron star mergers. These are typical sources of gravitational radiation. During these processes, not only mass and momenta are radiated away in form of gravitational waves, but also very strong magnetic fields are produced and radiated away. Moreover, a large portion of the energy is carried away by neutrinos. We give constraints on the conditions, where the energy transported by electromagnetic radiation is of similar or slightly higher order than the energy radiated in gravitational waves or in form of neutrinos. We find that for coalescing neutron stars, large magnetic fields magnify the Christodoulou memory as long as the gaseous environment is sufficiently rarefied. Thus the observed effect on test masses of a laser interferometer gravitational wave detector will be enlarged by the contribution of the electromagnetic field. Therefore, the present results are important for the planned experiments.

We present a model that attempts to fuse the inflationary era and the subsequent radiation dominated era under a unified framework so as to provide a smooth transition between the two. The model is based on a modification of the generalized Chaplygin gas. We constrain the model observationally by mapping the primordial power spectrum of the scalar perturbations to the latest data of WMAP7. We compute as well the spectrum of the primordial gravitational waves as would be measured today.

Starting from a realistically sheared magnetic arcade connecting chromospheric, transition region to coronal plasma, we simulate the in-situ formation and sustained growth of a quiescent prominence in the solar corona. Contrary to previous works, our model captures all phases of the prominence formation, including the loss of thermal equilibrium, its successive growth in height and width to macroscopic dimensions, and the gradual bending of the arched loops into dipped loops, as a result of the mass accumulation. Our 2.5-dimensional, fully thermodynamically and magnetohydrodynamically consistent model mimics the magnetic topology of normal-polarity prominences above a photospheric neutral line, and results in a curtain-like prominence above the neutral line through which the ultimately dipped magnetic field lines protrude at a finite angle. The formation results from concentrated heating in the chromosphere, followed by plasma evaporation and later rapid condensation in the corona due to thermal instability, as verified by linear instability criteria. Concentrated heating in the lower atmosphere evaporates plasma from below to accumulate at the top of coronal loops and supply mass to the later prominence constantly. This is the first evaporation-condensation model study where we can demonstrate how the formed prominence stays in a force balanced state, which can be compared to the Kippenhahn-Schluter type magnetohydrostatic model, all in a finite low-beta corona.

In this paper we study in detail the effect of our recently proposed model of parity and charge-parity (PCP) violating varying alpha on the Cosmic Microwave Background (CMB) photon passing through the intra galaxy-cluster medium (ICM). The ICM is well known to be composed of magnetized plasma. According to our model, the polarization and intensity of the CMB would be affected when traversing through the ICM due to non-trivial scalar photon interactions. We have calculated the evolution of such polarization and intensity collectively, known as the stokes parameters of the CMB photon during its journey through the ICM and tested our results against the Sunyaev-Zel'dovich (SZ) measurement on Coma galaxy cluster. Our model contains a PCP violating parameter, $\beta$, and a scale of alpha variation $\omega$. Using the derived constrained on the photon-to-scalar conversion probability, ${\bar P}_{\gamma \rightarrow \phi}$, for Coma cluster in ref.[34] we found a contour plot in the ($\omega,\beta$) parameter plane. The $\beta =0$ line in this parameter space corresponds to well-studied Maxwell-dilaton type models which has lower bound on $\omega \gtrapprox 6.4 \times 10^{9}$ GeV. In general, as the absolute value of $\beta$ increases, lower bound on $\omega$ also increases. Our model in general predicts the modification of the CMB polarization with a non-trivial dependence on the parity violating coupling parameter $\beta$. However, it is unconstrained in this particular study. We show that this effect can in principle be detected in the future measurements on CMB polarization such that $\beta$ can also be constrained.

It has been shown that spectroscopy of transiting extrasolar planets can potentially provide a wealth of information about their atmospheres. Herein, we set up the inverse problem in spectroscopic retrieval. We use non-linear optimal estimation to retrieve the atmospheric state (pioneered for Earth sounding by Rodgers 1976, 2000). The formulation quantifies the the degrees-of-freedom and information content of the spectrum with respect to geophysical parameters; herein, we focus specifically on temperature and composition. First, we apply the technique to synthetic near infrared spectra, and explore the influence of spectral signal-to-noise ratio and resolution (the two important parameters when designing a future instrument) on the information content of the data. As expected, we find that the number of retrievable parameters increase with increasing signal-to-noise and resolution, although the gains quickly asymptote for large values. Second, we apply the methods to the previously studied atmosphere of HD198733b, and compare the results of our retrieval with those obtained by others.

It has been shown that spectroscopy of transiting extrasolar planets can potentially provide a wealth of information about their atmospheres. Herein, we set up the inverse problem in spectroscopic retrieval. We use non-linear optimal estimation to retrieve the atmospheric state (pioneered for Earth sounding by Rodgers 1976, 2000). The formulation quantifies the the degrees-of-freedom and information content of the spectrum with respect to geophysical parameters; herein, we focus specifically on temperature and composition. First, we apply the technique to synthetic near infrared spectra, and explore the influence of spectral signal-to-noise ratio and resolution (the two important parameters when designing a future instrument) on the information content of the data. As expected, we find that the number of retrievable parameters increase with increasing signal-to-noise and resolution, although the gains quickly asymptote for large values. Second, we apply the methods to the previously studied atmosphere of HD198733b, and compare the results of our retrieval with those obtained by others.

We propose the experiments on the collision of laser light and high intensity electromagnetic pulses generated by relativistic flying mirrors, with electron bunches produced by a conventional accelerator and with laser wake field accelerated electrons for studying extreme field limits in the nonlinear interaction of electromagnetic waves. The regimes of dominant radiation reaction, which completely changes the electromagnetic wave-matter interaction, will be revealed in the laser plasma experiments. This will result in a new powerful source of ultra short high brightness gamma-ray pulses. A possibility of the demonstration of the electron-positron pair creation in vacuum in a multi-photon processes can be realized. This will allow modeling under terrestrial laboratory conditions neutron star magnetospheres, cosmological gamma ray bursts and the Leptonic Era of the Universe.

Supernovae observations strongly support the presence of a cosmological constant, but its value, which we will call apparent, is normally determined assuming that the Universe can be accurately described by a homogeneous model. Even in the presence of a cosmological constant we cannot exclude nevertheless the presence of a small local inhomogeneity which could affect the apparent value of the cosmological constant. Neglecting the presence of the inhomogeneity can in fact introduce a systematic misinterpretation of cosmological data, leading to the distinction between an apparent and true value of the cosmological constant. We establish the theoretical framework to calculate the corrections to the apparent value of the cosmological constant by modeling the local inhomogeneity with a $\Lambda LTB$ solution. Our assumption to be at the center of a spherically symmetric inhomogeneous matter distribution correspond to effectively calculate the monopole contribution of the large scale inhomogeneities surrounding us, which we expect to be the dominant one, because of other observations supporting a high level of isotropy of the Universe around us. By performing a local Taylor expansion we analyze the number of independent degrees of freedom which determine the local shape of the inhomogeneity, and consider the issue of central smoothness, showing how the same correction can correspond to different inhomogeneity profiles. Contrary to previous attempts to fit data using large void models our approach is quite general. The correction to the apparent value of the cosmological constant is in fact present for local inhomogeneities of any size, and should always be taken appropriately into account both theoretically and observationally.

The standard flare model, which was proposed based on observations and magnetohydrodynamic theory, can successfully explain many observational features of solar flares. However, this model is just a framework, with many details awaiting to be filled in, including how reconnection is triggered. In this paper, we address an unanswered question: where do flare ribbons stop? With the data analysis of the 2003 May 29 flare event, we tentatively confirmed our conjecture that flare ribbons finally stop at the intersection of separatrices (or quasi-separatrix layer in a general case) with the solar surface. Once verified, such a conjecture can be used to predict the final size and even the lifetime of solar flares.

Coronal "EIT waves" appear as EUV bright fronts propagating across a significant part of the solar disk. The intriguing phenomenon provoked continuing debates on their nature and their relation with coronal mass ejections (CMEs). In this paper, we first summarize all the observational features of "EIT waves", which should be accounted for by any successful model. The theoretical models constructed during the past 10 years are then reviewed. Finally, the implication of the "EIT wave" research to the understanding of CMEs is discussed. The necessity is pointed out to revisit the nature of CME frontal loop.

We investigate the evolution of the linear cosmological perturbations in f(R) gravity, an alternative to dark energy for explaining the late-time cosmic acceleration. We numerically calculate the early-time evolution with an approximation we contrive to solve a problem that commonly appears when one solves the full evolution equations. With the approximate evolution equations we can fairly assess the effect of the gravity modification on the early-time evolution, thereby examining the validity of the general-relativity (GR) approximation that is widely used for the early universe. In particular, we compare the CMB photon density perturbation and the matter density perturbation obtained respectively by our approximation and the conventional GR approximation. We find that the effect of the gravity modification at early times in f(R) gravity may not be negligible. We conclude that to be self-consistent, in the f(R) theory one should employ the approximation presented in this paper instead of that of GR in the treatment of the early-time evolution.

We report on studies of the viability and sensitivity of the Askaryan Radio Array (ARA), a new initiative to develop a Teraton-scale ultra-high energy neutrino detector in deep, radio-transparent ice near Amundsen-Scott station at the South Pole. An initial prototype ARA detector system was installed in January 2011, and has been operating continuously since then. We report on studies of the background radio noise levels, the radio clarity of the ice, and the estimated sensitivity of the planned ARA array given these results, based on the first five months of operation. Anthropogenic radio interference in the vicinity of the South Pole currently leads to a few-percent loss of data, but no overall effect on the background noise levels, which are dominated by the thermal noise floor of the cold polar ice, and galactic noise at lower frequencies. We have also successfully detected signals originating from a 2.5 km deep impulse generator at a distance of over 3 km from our prototype detector, confirming prior estimates of kilometer-scale attenuation lengths for cold polar ice. These are also the first such measurements for propagation over such large slant distances in ice. Based on these data, ARA-37, the 200 km^2 array now under construction, will achieve the highest sensitivity of any planned or existing neutrino detector in the 10^{16}-10^{19} eV energy range.

We confront tracker field quintessence with observational data. The potentials considered in this paper include $V(\phi)\propto\phi^{-\alpha}$, $\exp(M_{p}/\phi)$, $\exp(M_{p}/\phi)-1$, $\exp(\beta M_{p}/\phi)$ and $\exp(\gamma M_{p}/\phi)-1$; while the data come from the latest SN Ia, CMB and BAO observations. Stringent parameter constraints are obtained. In comparison with the cosmological constant via information criteria, it is found that models with potentials $\exp(M_{p}/\phi)$, $\exp(M_{p}/\phi)-1$ and $\exp(\gamma M_{p}/\phi)-1$ are not supported by the current data.

We introduce a new thermochemical kinetics and photochemical model. We use high-temperature bidirectional reaction rates for important H, C, O and N reactions (most importantly for CH$_4$ to CO interconversion), allowing us to attain thermochemical equilibrium, deep in an atmosphere, purely kinetically. This allows ab initio chemical modeling of an entire atmosphere, from deep-atmosphere thermochemical equilibrium to the photochemically dominated regime. We use our model to explore the atmospheric chemistry of cooler ($T_{eff} < 10^3$ K) extrasolar giant planets. In particular, we choose to model the nearby hot Neptune GJ436b, the only planet in this temperature regime for which spectroscopic measurements and estimates of chemical abundances now exist. Recent {\it Spitzer} measurements with retrieval have shown that methane is driven strongly out of equilibrium and is deeply depleted on the dayside of GJ 436b, whereas quenched carbon monoxide is abundant. This is surprising because GJ 436b is cooler than many of the heavily irradiated hot Jovians and thermally favorable for CH$_4$, and thus requires an efficient mechanism for destroying it. We include realistic estimates of ultraviolet flux from the parent dM star GJ 436, to bound the direct photolysis and photosensitized depletion of CH$_4$. While our models indicate fairly rich disequilibrium conditions are likely in cooler exoplanets over a range of planetary metallicities, we are unable to generate the conditions for substantial CH$_4$ destruction. One possibility is an anomalous source of abundant H atoms between 0.01-1 bars (which attack CH$_4$), but we cannot as yet identify an efficient means to produce these hot atoms.

We report a spectroscopic analysis of an EIT wave event that occurred in active region 11081 on 2010 June 12 and was associated with an M2.0 class flare. The wave propagated near circularly. The south-eastern part of the wave front passed over an upflow region nearby a magnetic bipole. Using EIS raster observations for this region, we studied the properties of plasma dynamics in the wave front, as well as the interaction between the wave and the upflow region. We found a weak blueshift for the Fe XII {\lambda}195.12 and Fe XIII {\lambda}202.04 lines in the wave front. The local velocity along the solar surface, which is deduced from the line of sight velocity in the wave front and the projection effect, is much lower than the typical propagation speed of the wave. A more interesting finding is that the upflow and non-thermal velocities in the upflow region are suddenly diminished after the transit of the wave front. This implies a significant change of magnetic field orientation when the wave passed. As the lines in the upflow region are redirected, the velocity along the line of sight is diminished as a result. We suggest that this scenario is more in accordance with what was proposed in the field-line stretching model of EIT waves.

The Ultra-Fast Flash Observatory (UFFO) pathfinder is a space mission devoted to the measurement of Gamma-Ray Bursts (GRBs), especially their early light curves which will give crucial information on the progenitor stars and central engines of the GRBs. It consists of two instruments: the UFFO Burst Alert & Trigger telescope (UBAT) for the detection of GRB locations and the Slewing Mirror Telescope (SMT) for the UV/optical afterglow observations, upon triggering by UBAT. The UBAT employs a coded-mask {\gamma}/X-ray camera with a wide field of view (FOV), and is comprised of three parts: a coded mask, a hopper, and a detector module (DM). The UBAT DM consists of a LYSO scintillator crystal array, multi-anode photo multipliers, and analog and digital readout electronics. We present here the design and fabrication of the UBAT DM, as well as its preliminary test results.

Hundreds of gamma-ray burst (GRB) optical light curves have been measured since the discovery of optical afterglows. However, even after nearly 7 years of operation of the Swift Observatory, only a handful of measurements have been made soon (within a minute) after the gamma ray signal. This lack of early observations fails to address burst physics at short time scales associated with prompt emissions and progenitors. Because of this lack of sub-minute data, the characteristics of the rise phase of optical light curve of short-hard type GRB and rapid-rising GRB, which may account for ~30% of all GRB, remain practically unknown. We have developed methods for reaching sub-minute and sub-second timescales in a small spacecraft observatory. Rather than slewing the entire spacecraft to aim the optical instrument at the GRB position, we use rapidly moving mirror to redirect our optical beam. As a first step, we employ motorized slewing mirror telescope (SMT), which can point to the event within 1s, in the UFFO Pathfinder GRB Telescope onboard the Lomonosov satellite to be launched in Nov. 2011. UFFO's sub-minute measurements of the optical emission of dozens of GRB each year will result in a more rigorous test of current internal shock models, probe the extremes of bulk Lorentz factors, provide the first early and detailed measurements of fast-rise GRB optical light curves, and help verify the prospect of GRB as a new standard candle. We will describe the science and the mission of the current UFFO Pathfinder project, and our plan of a full-scale UFFO-100 as the next step.

It has been established that cold plasma condensations can form in a magnetic loop subject to localized heating of the footpoints. In this paper, we use grid-adaptive numerical simulations of the radiative hydrodynamic equations to parametrically investigate the filament formation process in a pre-shaped loop with both steady and finite-time chromospheric heating. Compared to previous works, we consider low-lying loops with shallow dips, and use a more realistic description for the radiative losses. We demonstrate for the first time that the onset of thermal instability satisfies the linear instability criterion. The onset time of the condensation is roughly \sim 2 hr or more after the localized heating at the footpoint is effective, and the growth rate of the thread length varies from 800 km hr-1 to 4000 km hr-1, depending on the amplitude and the decay length scale characterizing this localized chromospheric heating. We show how single or multiple condensation segments may form in the coronal portion. In the asymmetric heating case, when two segments form, they approach and coalesce, and the coalesced condensation later drains down into the chromosphere. With a steady heating, this process repeats with a periodicity of several hours. While our parametric survey confirms and augments earlier findings, we also point out that steady heating is not necessary to sustain the condensation. Once the condensation is formed, it can keep growing also when the localized heating ceases. Finally, we show that the condensation can survive continuous buffeting by perturbations resulting from the photospheric p-mode waves.

In this report we construct a phenomenological model in which the time variation of the fine structure constant, $\alpha$, is induced by a parity and charge-parity (PCP) violating interaction. Such a PCP violation in the photon sector has a distinct physical origin from that in the conventional models of this kind. We calculate the cosmological birefringence so induced in our model and show that it in turn produces a new non-vanishing multipole moment correlation between the temperature and the polarization anisotropies in the CMB spectrum. We have also calculated the amount of optical rotation due to a strong background magnetic field and the effect of our new PCP violating term on the variation of $\alpha$ during the cosmic evolution. We found that only in the radiation dominated era can the contribution of the new PCP violating term to the variation of $\alpha$ be non-vanishing.

We propose a phenomenological model for the early universe where there is a smooth transition between an early "quintessence" phase and a radiation dominated era. The matter content is modelled by an appropriately modified Chaplygin gas for the early universe. We constrain the model observationally by mapping the primordial power spectrum of the scalar perturbations to the latest data of WMAP7. We compute as well the spectrum of the primordial gravitational waves as would be measured today. We show that the high frequencies region of the spectrum depends on the free parameter of the model and most importantly this region of the spectrum can be within the reach of future gravitational waves detectors.

Matos, Guzman and Nunez proposed a model of galactic halo based on an exponential-potential scalar field that could induce a rotation curve that is constant for all radii. We demonstrate that with suitable boundary conditions, such scalar field dark matter (SDM) can not only produce the observed constant rotation curve at large radius but also give rise to the correct power-law scaling near the galactic core region. This solves the existing cusp-core problem faced by the conventional cold dark matter (CDM) model.

"EIT waves" are a globally propagating wavelike phenomenon. They were often interpreted as a fast-mode magnetoacoustic wave in the corona, despite various discrepancies between the fast-mode wave model and observations. To reconcile these discrepancies, we once proposed that "EIT waves" are apparent propagation of the plasma compression due to successive stretching of the magnetic field lines pushed by the erupting flux rope. According to this model, an "EIT wave" should be preceded by a fast-mode wave, which however was rarely observed. With the unprecedented high cadence and sensitivity of the {\em Solar Dynamics Observatory} ({\em SDO}) observations, we discern a fast-moving wave front with a speed of 560 km s$^{-1}$, ahead of an "EIT wave", which had a velocity of $\sim 190$ km s$^{-1}$, in the "EIT wave" event on 2010 July 27. The results, suggesting that "EIT waves" are not fast-mode waves, confirm the prediction of our fieldline stretching model for "EIT wave". In particular, it is found that the coronal Moreton wave was $\sim 3$ times faster than the "EIT wave" as predicted.

We set the first limits on the ultra-high energy (UHE) neutrino fluence at energies greater than 10^9 GeV from gamma-ray bursts (GRBs) based on data from the second flight of the ANtarctic Impulsive Transient Antenna (ANITA). During the 31 day flight of ANITA-II, 26 GRBs were recorded by Swift or Fermi. Of these, we analyzed the 12 GRBs which occurred during quiet periods when the payload was away from anthropogenic activity. In a blind analysis, we observe 0 events on a total background of 0.0044 events in the combined prompt window for all 12 low-background bursts. We also observe 0 events from the remaining 14 bursts. We place a 90% confidence level limit on the E^-4 prompt neutrino fluence of 2.5x10^17 GeV^3/cm^2 between 10^8 and 10^12 GeV from GRB090107A. This is the first reported limit on the UHE neutrino fluence from GRBs above 10^9 GeV, and the strongest limit above 10^8 GeV.

We have conducted a search for extended energy deposition trails left by ultra-relativistic magnetic monopoles interacting in Antarctic ice. The non-observation of any satisfactory candidates in the 31 days of accumulated ANITA-II flight data results in an upper limit on the diffuse flux of relativistic monopoles. We obtain a 90% C.L. limit of order 10^{-19}/(cm^2-s-sr) for values of Lorentz boost factor 10^{10}<gamma at the anticipated energy E=10^{16} GeV. This bound is stronger than all previously published experimental limits for this kinematic range.

We demonstrate that an extremely small but positive quantum correction, or the Casimir energy, to the cosmological constant can arise from a massive bulk fermion field in the Randall-Sundrum model. Specifically, a cosmological constant doubly descended from the Planck-electroweak hierarchy and as minute as the observed dark energy scale can be naturally achieved without fine-tuning of the bulk fermion mass. To ensure the stabilization of the system, we discuss two stabilization mechanisms under this setup. It is found that the Goldberger-Wise mechanism can be successfully introduced in the presence of a massive bulk fermion, without spoiling the smallness of the quantum correction.

Recently Disney et al. (2008) found a striking correlation among the five basic parameters that govern the galactic dynamics: R50, R90, Lr, Md, and MHI . They suggested that this is in conflict with the LCDM model, which predicts the hierarchical formation of cosmic structures from bottom up. In light of the importance of the issue, we performed a similar analysis on global parameters of galaxies with a significantly larger database and two additional parameters, LJ and RJ, of the near-infrared J band. We used databases from the Arecibo Legacy Fast Arecibo L-band Feed Array Survey for the atomic gas properties, the Sloan Digital Sky Survey for the optical properties, and the Two Micron All Sky Survey for the near-infrared properties, of the galaxies. We conducted principal component analysis (PCA) to find relations among these observational variables and confirmed that the five parameters in the work of Disney et al. are indeed correlated. The first principal component dominates the correlations among the five parameters and can explain 86% of the variation in the data. When color (g - i) is included, the first component still dominates and the color forms a second principal component that is almost independent of other parameters. The overall trend in our near-infrared PCA is very similar, except that color (i - J) seems even more decoupled from all other parameters. The dominance of the first principal component suggests that the structure of galaxies is governed by a single physical parameter. This confirms the observational results in Disney et al. However, based on the importance of the baryon physics in galaxy formation, we are not convinced that the hierarchical structure formation scenario and the notion of cold dark matter are necessarily flawed.

The radio approach for detecting the ultra-high energy cosmic neutrinos has become a mature field. The Cherenkov pulse in radio detection originates from the charge excess of particle showers due to Askaryan effect. The conventional way of calculating the Cherenkov pulse by making far- field approximation fails when the size of elongated showers become comparable with detection distance. We investigate the Cherenkov pulse in near-field by a numerical code based on the finite- difference time-domain (FDTD) method. Our study shows that the near-field radiation exhibits very different behaviors from the far-field one and therefore can be easily recognized. For ground array neutrino detectors, the near-field radiation would provide a unique signature for ultra high energy electromagnetic showers induced by the electron neutrino charge-current interaction. This can be useful in neutrino flavor identification.

With gravity, ionization, and radiation being considered, we perform 2.5D compressible resistive MHD simulations of chromospheric magnetic reconnection using the CIP-MOCCT scheme. The temperature distribution of the quiet-Sun atmospheric model VALC and the helium abundance (10%) are adopted. Our 2.5D MHD simulation reproduces qualitatively the temperature enhancement observed in chromospheric microflares. The temperature enhancement $\Delta T$ is demonstrated to be sensitive to the background magnetic field, whereas the total evolution time $\Delta t$ is sensitive to the magnitude of the anomalous resistivity. Moveover, we found a scaling law, which is described as $\Delta T/\Delta t \sim {n_H}^{-1.5} B^{2.1} {\eta_0}^{0.88}$. Our results also indicate that the velocity of the upward jet is much greater than that of the downward jet and the X-point may move up or down.

This is an erratum to our paper in Physical Review D82:022004,2010, corresponding to preprint: arXiv:1003.2961 .

We consider the production of primordial micro black holes (MBH) remnants in the early universe. These objects induce the universe to be in a matter-dominated era before the onset of inflation. Effects of such an epoch on the CMB power spectrum are discussed and computed both analytically and numerically. By comparison with the latest observational data from the WMAP collaboration, we find that our model is able to explain the quadrupole anomaly of the CMB power spectrum.

Accurate mass determination of clusters of galaxies is crucial if they are to be used as cosmological probes. However, there are some discrepancies between cluster masses determined based on gravitational lensing, and X-ray observations assuming strict hydrostatic equilibirium (i.e., the equilibrium gas pressure is provided entirely by thermal pressure). Cosmological simulations suggest that turbulent gas motions remaining from hierarchical structure formation may provide a significant contribution to the equilibrium pressure in clusters. We analyze a sample of massive clusters of galaxies drawn from high resolution cosmological simulations, and find a significant contribution (20%-45%) from non-thermal pressure near the center of relaxed clusters, and, in accord with previous studies, a minimum contribution at about 0.1 Rvir, growing to about 30%-45% at the virial radius, Rvir. Our results strongly suggest that relaxed clusters should have significant non-thermal support in their core region. As an example, we test the validity of strict hydrostatic equilibirium in the well-studied massive galaxy cluster Abell 1689 using the latest high resolution gravitational lensing and X-ray observations. We find a contribution of about 40% from non-thermal pressure within the core region of A1689, suggesting an alternate explanation for the mass discrepancy: the strict hydrostatic equilibirium is not valid in this region.

The kinematic properties of coronal mass ejections (CMEs) suffer from the projection effects, and it is expected that the real velocity should be larger and the real angular width should be smaller than the apparent values. Several attempts have been tried to correct the projection effects, which however led to a too large average velocity probably due to the biased choice of the CME events. In order to estimate the overall influence of the projection effects on the kinematic properties of the CMEs, we perform a forward modeling of the real distributions of the CME properties, such as the velocity, the angular width, and the latitude, by requiring their projected distributions to best match the observations. Such a matching is conducted by Monte Carlo simulations. According to the derived real distributions, it is found that (1) the average real velocity of all non-full-halo CMEs is about 514 km s$^{-1}$, and the average real angular width is about 33$^\circ$, in contrast to the corresponding apparent values of 418 km s$^{-1}$ and 42.7$^\circ$ in observations; (2) For the CMEs with the angular width in the range of $20^\circ- 120^\circ$, the average real velocity is 510 km s$^{-1}$ and the average real angular width is 43.4$^\circ$, in contrast to the corresponding apparent values of 392 km s$^{-1}$ and 52$^\circ$ in observations.

We investigate the cosmological and the local tests of the f(R) theory of modified gravity via the observations of (1) the cosmic expansion and (2) the cosmic structures and via (3) the solar-system experiments. To fit the possible cosmic expansion histories under consideration, for each of them we reconstruct f(R), known as "designer f(R)". We then test the designer f(R) via the cosmic-structure constraints on the metric perturbation ratio Psi/Phi and the effective gravitational coupling G_eff and via the solar-system constraints on the Brans-Dicke theory with the chameleon mechanism. We find that among the designer f(R) models specified by the CPL effective equation of state w_eff, only the model closely mimicking general relativity with a cosmological constant (LambdaCDM) can survive all the tests. Accordingly, these tests rule out the frequently studied "w_eff = -1" designer f(R) models which are distinct in cosmic structures from LambdaCDM. When considering only the cosmological tests, we find that the surviving designer f(R) models, although exist for a variety of w_eff, entail fine-tuning.

EIT waves are a wavelike phenomenon propagating outward from the coronal mass ejection (CME) source region, with expanding dimmings following behind. We present a spectroscopic study of an EIT wave/dimming event observed by Hinode/EIS. Although the identification of the wave front is somewhat affected by the pre-existing loop structures, the expanding dimming is well defined. We investigate the line intensity, width, and Doppler velocity for 4 EUV lines. In addition to the significant blue shift implying plasma outflows in the dimming region as revealed in previous studies, we find that the widths of all the 4 spectral lines increase at the outer edge of the dimmings. We illustrate that this feature can be well explained by the field line stretching model, which claims that EIT waves are apparently moving brightenings that are generated by the successive stretching of the closed field lines.

We comment on the general solution of the scalar field dark matter provided in the paper "Remarks on the spherical scalar field halo in galaxies" by Kamal K. Nandi, Ildar Valitov and Nail G. Migranov. The authors made a mistake in the general form of the tangential pressure profile p_t(r), which deviates from the correct profile, especially when r is small. Although this mistake does not alter significantly the value of w(r) when the integration constant D is small, we found that it does result in an overestimate of w(r) when D is large.

In 2009 Ho\v{r}ava proposed a power-counting renormalizable quantum gravity theory. Afterwards a term in the action that softly violates the detailed balance condition has been considered with the attempt of obtaining a more realistic theory in its IR-limit. This term is proportional to $\omega R^{(3)}$, where $\omega$ is a constant parameter and $R^{(3)}$ is the spatial Ricci scalar. In this paper we derive constraints on this IR-modified Ho\v{r}ava theory using the late-time cosmic accelerating expansion observations. We obtain a lower bound of $|\omega|$ that is nontrivial and depends on $\Lambda_W$, the cosmological constant of the three dimensional spatial action in the Ho\v{r}ava gravity. We find that to preserve the detailed balance condition, one needs to fine-tune $\Lambda_W$ such that $- 2.29\times 10^{-4}< (c^2 \Lambda_W)/(H^2_0 \currentDE) - 2 < 0 $, where $H_0$ and $\currentDE$ are the Hubble parameter and dark energy density fraction in the present epoch, respectively. On the other hand, if we do not insist on the detailed balance condition, then the valid region for $\Lambda_W$ is much relaxed to $-0.39< (c^2 \Lambda_W)/(H^2_0 \currentDE) - 2 < 0.12$. We find that although the detailed balance condition cannot be ruled out, it is strongly disfavored.

It has been suggested by Chen and Lai that the proper description of the large scale structure formation of the universe in the post-reionization era, which is conventionally characterized via gas hydrodynamics, should include the plasma collective effects in the formulation. Specifically, it is the combined pressure from the baryon thermal motions and the residual long-range electrostatic potentials resulted from the imperfect Debye shielding, that fights against the gravitational collapse. As a result, at small-scales the baryons would oscillate at the ion-acoustic, instead of the conventional neutral acoustic, frequency. In this paper we extend and improve the Chen-Lai formulation with the attention to the Landau damping of the ion-acoustic oscillations. Since T_e \sim T_i in the post-reionization era, the ion acoustic oscillations would inevitably suffer the Landau damping which severely suppresses the baryon density spectrum in the regimes of intermediate and high wavenumber k. To describe this Landau-damping phenomenon more appropriately, we find it necessary to modify the filtering wavenumber k_f in our analysis. It would be interesting if our predicted Landau damping of the ion-acoustic oscillations can be observed at high redshifts.

The Antarctic Impulsive Transient Antenna (ANITA) completed its second long-duration balloon flight in January 2009, with 31 days aloft (28.5 live days) over Antarctica. ANITA searches for impulsive coherent radio Cherenkov emission from 200 to 1200 MHz, arising from the Askaryan charge excess in ultra-high energy neutrino-induced cascades within Antarctic ice. This flight included significant improvements over the first flight in the payload sensitivity, efficiency, and a flight trajectory over deeper ice. Analysis of in-flight calibration pulses from surface and sub-surface locations verifies the expected sensitivity. In a blind analysis, we find 2 surviving events on a background, mostly anthropogenic, of 0.97+-0.42 events. We set the strongest limit to date for 1-1000 EeV cosmic neutrinos, excluding several current cosmogenic neutrino models.

We report the observation of sixteen cosmic ray events of mean energy of 1.5 x 10^{19} eV, via radio pulses originating from the interaction of the cosmic ray air shower with the Antarctic geomagnetic field, a process known as geosynchrotron emission. We present the first ultra-wideband, far-field measurements of the radio spectral density of geosynchrotron emission in the range from 300-1000 MHz. The emission is 100% linearly polarized in the plane perpendicular to the projected geomagnetic field. Fourteen of our observed events are seen to have a phase-inversion due to reflection of the radio beam off the ice surface, and two additional events are seen directly from above the horizon.

Motivated by the possible extension into a supersymmetric Randall-Sundrum (RS) model, we investigate the properties of the vacuum expectation value (VEV) of the stress-energy tensor for a quantized bulk Dirac spinor field in the RS geometry and compare it with that for a real scalar field. This is carried out via the Green function method based on first principles without invoking the degeneracy factor, whose validity in a warp geometry is a priori unassured. In addition, we investigate the local behavior of the Casimir energy near the two branes. One salient feature we found is that the surface divergences near the two branes have opposite signs. We argue that this is a generic feature of the fermionic Casimir energy density due to its parity transformation in the fifth dimension. Furthermore, we investigate the self-consistency of the RS metric under the quantum correction due to the stress-energy tensor. It is shown that the VEV of the stress-energy tensor and the classical one become comparable near the visible brane if k ~ M ~ M_Pl (the requirement of no hierarchy problem), where k is the curvature of the RS warped geometry and M the 5-dimensional Planck mass. In that case the self-consistency of RS model that includes bulk fields is in doubt. If, however, k <~ M, then an approximate self-consistency of the RS-type metric may still be satisfied.

Halo coronal mass ejections (CMEs) were found to be significantly faster than normal CMEs, which was a long-standing puzzle. In order to solve the puzzle, we first investigate the observed properties of 31 limb CMEs that display clearly loop-shaped frontal loops. The observational results show a strong tendency that slower CMEs are weaker in the white-light intensity. Then, we perform a Monte Carlo simulation of 20 000 artificial limb CMEs that have average velocity of $\sim$523 km s$^{-1}$. The Thomson scattering of these events is calculated when they are assumed to be observed as limb and halo events, respectively. It is found that the white-light intensity of many slow CMEs becomes remarkably reduced as they turn from being viewed as a limb event to as a halo event. When the intensity is below the background solar wind fluctuation, it is assumed that they would be missed by coronagraphs. The average velocity of "detectable" halo CMEs is $\sim$922 km s$^{-1}$, very close to the observed value. It also indicates that wider events are more likely to be recorded. The results soundly suggest that the higher average velocity of halo CMEs is due to that a majority of slow events and a part of narrow fast events carrying less material are so faint that they are blended with the solar wind fluctuations, and therefore are not observed.

"EIT waves" are a wavelike phenomenon propagating in the corona, which were initially observed in the extreme ultraviolet (EUV) wavelength by the EUV Imaging Telescope (EIT). Their nature is still elusive, with the debate between fast-mode wave model and non-wave model. In order to distinguish between these models, we investigate the relation between the EIT wave velocity and the local magnetic field in the corona. It is found that the two parameters show significant negative correlation in most of the EIT wave fronts, {\it i.e.}, EIT wave propagates more slowly in the regions of stronger magnetic field. Such a result poses a big challenge to the fast-mode wave model, which would predict a strong positive correlation between the two parameters. However, it is demonstrated that such a result can be explained by the fieldline stretching model, \emph{i.e.,} that "EIT waves" are apparently-propagating brightenings, which are generated by successive stretching of closed magnetic field lines pushed by the erupting flux rope during coronal mass ejections (CMEs).

Detection of molecules using infrared spectroscopy probes the conditions and compositions of exoplanet atmospheres. Water (H2O), methane (CH4), carbon dioxide (CO2), and carbon monoxide (CO) have been detected in two hot Jupiters. These previous results relied on space-based telescopes that do not provide spectroscopic capability in the 2.4 - 5.2 micron spectral region. Here we report ground-based observations of the dayside emission spectrum for HD 189733b between 2.0-2.4 micron and 3.1-4.1 micron, where we find a bright emission feature. Where overlap with space-based instruments exists, our results are in excellent agreement with previous measurements. A feature at ~3.25 micron is unexpected and difficult to explain with models that assume local thermodynamic equilibrium (LTE) conditions at the 1 bar to 1 x 10-6 bar pressures typically sampled by infrared measurements. The most likely explanation for this feature is that it arises from non-LTE emission from CH4, similar to what is seen in the atmospheres of planets in our own Solar System. These results suggest that non-LTE effects may need to be considered when interpreting measurements of strongly irradiated exoplanets.

Askaryan Radio Array (ARA), a large-scale radio Cherenkov observatory which scientists propose to develop in Antarctica, aims at discovering the origin and evolution of the cosmic accelerators that produce the highest energy cosmic rays by means of observing the ultra high energy (UHE) cosmogenic neutrinos. To optimize ARA's angular resolution of the incoming UHE neutrinos, which is essential for pointing pack to its source, the relation between the reconstruction capabilities of ARA and its design is studied. It is found that with the noise effect taken into account, in order to make this neutrino angular resolution as good as possible and detection efficiency as high as possible, the optimal choice for ARA geometry would be the station spacing of 1.6 km and the antenna spacing of 40 m.

The radio approach for detecting the ultra-high energy cosmic neutrinos has become a mature field. The Cherenkov signals in radio detection are originated from the charge excess of particle showers due to Askaryan effect. The conventional way of calculating the Cherenkov pulses by making Fraunhofer approximation fails when the sizes of the elongated showers become comparable with the detection distances. We present a calculation method of Cherenkov pulses based on the finite-difference time-domain (FDTD) method, and attain a satisfying effeciency via the GPU- acceleration. Our method provides a straightforward way of the near field calculation, which would be important for ultra high energy particle showers, especailly the electromagnetic showers induced by the high energy leptons produced in the neutrino charge current interactions.

Magnetic field in the solar lower atmosphere can be measured by the use of the Zeeman and Hanle effects. In contrast, the coronal magnetic field well above the solar surface, which directly controls various eruptive phenomena, can not be precisely measureed with the traditional techniques. Several attempts are being made to probe the coronal magnetic field, such as force-free extrapolation based on the photospheric magnetograms, gyroresonance radio emissions, and coronal seismology based on MHD waves in the corona. Compared to the waves trapped in the localized coronal loops, EIT waves are the only global-scale wave phenomenon, and thus are the ideal tool for the coronal global seismology. In this paper, we review the observations and modelings of EIT waves, and illustrate how they can be applied to probe the global magnetic field in the corona.

Hundreds of gamma-ray burst (GRB) UV-optical light curves have been measured since the discovery of optical afterglows, however, even after nearly 5 years of operation of the SWIFT observatory, only a handful of measurements have been made soon (within a minute) after the gamma ray signal. This lack of early observations fails to address burst physics at the short time scales associated with burst events and progenitors. Because of this lack of sub-minute data, the characteristics of the UV-optical light curve of short-hard type GRB and rapid-rising GRB, which may account for ~30% of all GRB, remain practically unknown. We have developed methods for reaching the sub-minute and the sub-second timescales in a small spacecraft observatory appropriate for launch on a microsatellite. Rather than slewing the entire spacecraft to aim the UV-optical instrument at the GRB position, we use rapidly moving mirrors to redirect our optical beam. Our collaboration has produced a unique MEMS (microelectromechanical systems) micromirror array which can point and settle on target in only 1 ms. This technology is proven, flying successfully as the MTEL (MEMS Telescope for Extreme Lightning) on the Tatiana-2 Spacecraft in September 2009 and as the KAMTEL on the International Space Station in April 2008. The sub-minute measurements of the UV-optical emission of dozens of GRB each year will result in a more rigorous test of current internal shock models, probe the extremes of bulk Lorentz factors, and provide the first early and detailed measurements of fast-rise and short type GRB UV-optical light curves.

Over 300 extrasolar planets (exoplanets) have been detected orbiting nearby stars. We now hope to conduct a census of all planets around nearby stars and to characterize their atmospheres and surfaces with spectroscopy. Rocky planets within their star's habitable zones have the highest priority, as these have the potential to harbor life. Our science goal is to find and characterize all nearby exoplanets; this requires that we measure the mass, orbit, and spectroscopic signature of each one at visible and infrared wavelengths. The techniques for doing this are at hand today. Within the decade we could answer long-standing questions about the evolution and nature of other planetary systems, and we could search for clues as to whether life exists elsewhere in our galactic neighborhood.

We comment on the calculation mistake in the paper "w and w' of scalar field models of dark energy" by Takeshi Chiba [arXiv:astro-ph/0510598], where w is the dark energy equation of state and w' is the time derivative of w in units of the Hubble time. The author made a mistake while rewriting the phantom equation of motion, which led to an incorrect generic bound for the phantom model and an incorrect bound for the tracker phantom model on the w-w' plane.

Classification of dark energy models in the plane of w and w', where w is the dark energy equation of state and w' its time-derivative in units of the Hubble time, has been studied in the literature. We take the current SN Ia, CMB and BAO data, invoke a widely used parametrization of the dark energy equation of state, and obtain the constraints on the w -w' plane. We find that dark energy models including the cosmological constant, phantom, non-phantom barotropic fluids, and monotonic up-rolling quintessence are ruled out at the 68.3% confidence level based on the current observational data. Down-rolling quintessence, including the thawing and the freezing models, is consistent with the current observations. All the above-mentioned models are still consistent with the data at the 95.4% confidence level.

Recently we proposed a new approach to the testing of dark energy models based on the observational data. In that work we focused particularly on quintessence models for demonstration and invoked a widely used parametrization of the dark energy equation of state. In this paper we take the more recent SN Ia, CMB and BAO data, invoke the same parametrization, and apply this method of consistency test to five categories of dark energy models, including the LCDM model, the generalized Chaplygin gas, and three quintessence models: exponential, power-law and inverse-exponential potentials. We find that the exponential potential of quintessence is ruled out at the 95.4% confidence level, while the other four models are consistent with data. This consistency test can be efficiently performed since for all models it requires the constraint of only a single parameter space that by choice can be easily accessed.

We propose a new approach to the consistency test of dark energy models with observations. To test a category of dark energy models, we suggest introducing a characteristic Q(z) that in general varies with the redshift z but in those models plays the role of a (constant) distinct parameter. Then, by reconstructing dQ(z)/dz from observational data and comparing it with zero we can assess the consistency between data and the models under consideration. For a category of models that passes the test, we can further constrain the distinct parameter of those models by reconstructing Q(z) from data. For demonstration, in this paper we concentrate on quintessence. In particular we examine the exponential potential and the power-law potential via a widely used parametrization of the dark energy equation of state, w(z) = w_0 + w_a z/(1+z), for data analysis. This method of the consistency test is particularly efficient because for all models we invoke the constraint of only a single parameter space that by choice can be easily accessed. The general principle of our approach is not limited to dark energy. It may also be applied to the testing of various cosmological models and even the models in other fields beyond the scope of cosmology.

Using the NICMOS instrument on the Hubble Space Telescope, we have measured the dayside spectrum of HD 209458b between 1.5--2.5 microns. The emergent spectrum is dominated by features due to the presence of methane (CH4) and water vapor (H2O), with smaller contributions from carbon dioxide (CO2). Combining this near-infrared spectrum with existing mid-infrared measurements shows the existence of a temperature inversion and confirms the interpretation of previous photometry measurements. We find a family of plausible solutions for the molecular abundance and detailed temperature profile. Observationally resolving the ambiguity between abundance and temperature requires either (1) improved wavelength coverage or spectral resolution of the dayside emission spectrum, or (2) a transmission spectrum where abundance determinations are less sensitive to the temperature structure.

More and more evidence indicates that "EIT waves" are strongly related to coronal mass ejections (CMEs). However, it is still not clear how the two phenomena are related to each other. We investigate a CME event on 1997 September 9, which was well observed by both EUV imaging telescope (EIT) and the high-cadence MK3 coronagraph at Mauna Loa Solar Observatory, and compare the spatial relation between the "EIT wave" fronts and the CME leading loops. It is found that "EIT wave" fronts are co-spatial with the CME leading loops, and the expanding EUV dimmings are co-spatial with the CME cavity. It is also found that the CME stopped near the boundary of a coronal hole, a feature common to observations of "EIT waves". It is suggested that "EIT waves"/dimmings are the EUV counterparts of the CME leading loop/cavity, based on which we propose that, as in the case of "EIT waves", CME leading loops are apparently-moving density enhancements that are generated by successive stretching (or opening-up) of magnetic loops.

We comment on the calculational mistake in the paper "Modeling galaxy halos using dark matter with pressure" by Somnath Bharadwaj and Sayan Kar. The authors made a mistake while calculating the metric, which led to an overestimate of the deflection angle of light passing through the halos for -1<w_r<-0.5 and an underestimate of the deflection angle for -0.5<w_r<0. In addition, solution for w_r>0 should not exist. Although the Bharadwaj-Kar solution should be corrected, it appears that the characteristics of the deflection angle under the supposed non-conventional non-ideal fluid equation of state for the dark matter halo remain sensitive to the impact parameter and may be verifiable through observations.

We discuss design considerations and simulation results for IceRay, a proposed large-scale ultra-high energy (UHE) neutrino detector at the South Pole. The array is designed to detect the coherent Askaryan radio emission from UHE neutrino interactions in the ice, with the goal of detecting the cosmogenic neutrino flux with reasonable event rates. Operating in coincidence with the IceCube neutrino detector would allow complete calorimetry of a subset of the events. We also report on the status of a testbed IceRay station which incorporates both ANITA and IceCube technology and will provide year-round monitoring of the radio environment at the South Pole.

One of the most tantalizing questions in astronomy and astrophysics, namely the origin and the evolution of the cosmic accelerators that produce the highest energy cosmic rays (UHECR), may be best addressed through the observation of ultra high energy (UHE) cosmogenic neutrinos. Neutrinos travel from their source undeflected by magnetic fields and unimpeded by interactions with the cosmic microwave background. At high energies, neutrinos could be detected in dense, radio frequency (RF) transparent media via the Askaryan effect. The abundant cold ice covering the geographic South Pole, with its exceptional RF clarity, has been host to several pioneering efforts to develop this approach, including RICE and ANITA. Building on the expertise gained in these efforts, and the infrastructure developed in the construction of the IceCube optical Cherenkov observatory, a low-cost array of radio frequency antenna stations could be deployed near the Pole to efficiently detect a significant number of UHE neutrinos with degree scale angular resolution within the next decade. Such an array, if installed in close proximity to IceCube, could allow cross-calibration on a small but invaluable subset of neutrino events detected by both the optical and radio methods. In addition to providing critical information in the identification of the source of UHECRs, such an observatory could also provide a unique probe of long baseline high energy neutrino interactions unattainable in any man-made neutrino beam.

Aims. We try to identify ultra-cool dwarfs from the seventh Data Release of the Sloan Digital Sky Survey (SDSS DR7) with SDSS i-z and r-z colors. We also obtain proper motion data from SDSS, 2MASS, and UKIDSS and improve spectral typing from SDSS and 2MASS photometric colors. Methods. We selected ultra-cool dwarf candidates from the SDSS DR7 with new photometric selection criteria, which are based on a parameterization study of known L and T dwarfs. The objects are then cross-identified with the Two Micron All Sky Survey and the Fourth Data Release of the UKIRT Infrared Deep Sky Survey (UKIDSS DR4). We derive proper motion constraints by combining SDSS, 2MASS, and UKIDSS positional information. In this way we are able to assess, to some extent, the credence of our sample using a multi epoch approach, which complements spectroscopic confirmation. Some of the proper motions are affected by short baselines, but, as a general tool, this method offers great potential to confirm faint L dwarfs as UKIDSS coverage increases. In addition we derive updated color-spectral type relations for L and T dwarfs with SDSS and 2MASS magnitudes. Results. We present 59 new nearby M and L dwarfs selected from the imaging catalog of the SDSS DR7, including proper motions and spectral types calculated from the updated color-spectral type relations. and obtain proper motions from SDSS, 2MASS, and UKIDSS for all of our objects.

On 2007 December 7, a small filament located in a small active region AR 10977 erupted and led to a B1.4 flare. An EIT wave associated with this eruption was observed both by SOHO/EIT and by EUVI on board STEREO. According to the observations from SOHO/LASCO and STEREO/COR A, we found that there was no CME associated with the EIT wave. This seems to challenge the argument that the cause of EIT waves is CME. However the data from STEREO/COR B indicated that there was a narrow CME associated with the EIT wave. This suggests that studying CMEs by investigating observations made in one direction alone may not be able to guarantee the reliability of the results.

The direct current (DC) electric field near the reconnection region has been proposed as an effective mechanism to accelerate protons and electrons in solar flares. A power-law energy spectrum was generally claimed in the simulations of electron acceleration by the reconnection electric field. However in most of the literature, the electric and magnetic fields were chosen independently. In this paper, we perform test-particle simulations of electron acceleration in a reconnecting magnetic field, where both the electric and magnetic fields are adopted from numerical simulations of the MHD equations. It is found that the accelerated electrons present a truncated power-law energy spectrum with an exponential tail at high energies, which is analogous to the case of diffusive shock acceleration. The influences of reconnection parameters on the spectral feature are also investigated, such as the longitudinal and transverse components of the magnetic field and the size of the current sheet. It is suggested that the DC electric field alone might not be able to reproduce the observed single or double power-law distributions.

We report initial results of the first flight of the Antarctic Impulsive Transient Antenna (ANITA-1) 2006-2007 Long Duration Balloon flight, which searched for evidence of a diffuse flux of cosmic neutrinos above energies of 3 EeV. ANITA-1 flew for 35 days looking for radio impulses due to the Askaryan effect in neutrino-induced electromagnetic showers within the Antarctic ice sheets. We report here on our initial analysis, which was performed as a blind search of the data. No neutrino candidates are seen, with no detected physics background. We set model-independent limits based on this result. Upper limits derived from our analysis rule out the highest cosmogenic neutrino models. In a background horizontal-polarization channel, we also detect six events consistent with radio impulses from ultra-high energy extensive air showers.

We have measured the dayside spectrum of HD 189733b between 1.5 and 2.5 microns using the NICMOS instrument on the Hubble Space Telescope. The emergent spectrum contains significant modulation, which we attribute to the presence of molecular bands seen in absorption. We find that water (H2O), carbon monoxide (CO), and carbon dioxide (CO2) are needed to explain the observations, and we are able to estimate the mixing ratios for these molecules. We also find temperature decreases with altitude in the ~0.01 < P < ~1 bar region of the dayside near-infrared photosphere and set an upper limit to the dayside abundance of methane (CH4) at these pressures.

We present a detailed report on the experimental details of the Antarctic Impulsive Transient Antenna (ANITA) long duration balloon payload, including the design philosophy and realization, physics simulations, performance of the instrument during its first Antarctic flight completed in January of 2007, and expectations for the limiting neutrino detection sensitivity. Neutrino physics results will be reported separately.

Radiowave detection of the Cherenkov radiation produced by neutrino-ice collisions requires an understanding of the radiofrequency (RF) response of cold polar ice. We herein report on a series of radioglaciological measurements performed approximately 10 km north of Taylor Dome Station, Antarctica from Dec. 6, 2006 - Dec. 16, 2006. Using RF signals broadcast from: a) an englacial discone, submerged to a depth of 100 meters and broadcasting to a surface dual polarization horn receiver, and b) a dual-polarization horn antenna on the surface transmitting signals which reflect off the underlying bed and back up to the surface receiver, we have made time-domain estimates of both the real (index-of-refraction) and imaginary (attenuation length) components of the complex ice dielectric constant. We have also measured the uniformity of ice response along two orthogonal axes in the horizontal plane. We observe a wavespeed asymmetry of order 0.1%, projected onto the vertical propagation axis, consistent with some previous measurements, but somewhat lower than others.

Magnetowave induced plasma wakefield acceleration (MPWA) in a relativistic astrophysical outflow has been proposed as a viable mechanism for the acceleration of cosmic particles to ultra high energies. Here we present simulation results that demonstrate the viability of this mechanism. We invoke the high frequency and high speed whistler mode for the driving pulse. The plasma wakefield so induced validates precisely the theoretical prediction. This mechanism is shown capable of accelerating charged particles to ZeV energies in Active Galactic Nuclei (AGN).

Giant flares from soft gamma-ray repeaters (SGRs) are one of the most violent phenomena in neutron stars. Quasi-periodic oscillations (QPOs) with frequencies ranging from 18 to 1840 Hz have been discovered in the tails of giant flares from two SGRs, and were ascribed to be seismic vibrations or torsional oscillations of magnetars. Here we propose an alternative explanation for the QPOs in terms of standing sausage mode oscillations of flux tubes in the magnetar coronae. We show that most of the QPOs observed in SGR giant flares could be well accounted for except for those with very high frequencies (625 and 1840 Hz).

It is well-known that the universe was in a plasma state both before decoupling and after reionization. However, the conventional wisdom has been that the plasma effects are largely Debye-shielded and can thus be safely ignored when considering large scale evolutions. Recently we showed that large scale structure formation in the universe may actually be suppressed by the plasma collective effect. Indeed, observational data indicate that the conventional theoretical formula tends to overestimate the matter power spectrum at scales $k>1h{\rm Mpc}^{-1}$. In this paper, we further develop our theory through a more thorough and general derivation of the Maxwell-Einstein-Boltzmann equation. In addition to baryon density perturbation post reionization, we apply this general formulation to investigate the possible plasma effect on CMB anisotropy. As expected, while the plasma effect does render an observable effect to the former, its impact on the latter is totally negligible.

Coronal mass ejections (CMEs), as a large-scale eruptive phenomenon, often reveal some precursors in the initiation phase, e.g., X-ray brightening, filament darkening, etc, which are useful for CME modeling and space weather forecast. With the SOHO/SUMER spectroscopic observations of the 2000 September 26 event, we propose another precursor for CME eruptions, namely, long-time prominence oscillations. The prominence oscillation-and-eruption event was observed by ground-based H$\alpha$ telescopes and space-borne white-light, EUV imaging and spectroscopic instruments. In particular, the SUMER slit was observing the prominence in a sit-and-stare mode. The observations indicate that a siphon flow was moving from the proximity of the prominence to a site at a projected distance of 270$''$, which was followed by repetitive H$\alpha$ surges and continual prominence oscillations. The oscillation lasted 4 hours before the prominence erupted as a blob-like CME. The analysis of the multiwavelength data indicates that the whole series of processes fits well into the emerging flux trigger mechanism for CMEs. In this mechanism, emerging magnetic flux drives a siphon flow due to increased gas pressure where the background polarity emerges. It also drives H$\alpha$ surges through magnetic reconnection where the opposite polarity emerges. The magnetic reconnection triggers the prominence oscillations, as well as its loss of equilibrium, which finally leads to the eruption of the prominence. It is also found that the reconnection between the emerging flux and the pre-existing magnetic loop proceeds in an intermittent, probably quasi-periodic, way.

This paper reviews recent progress in the research on the initiation and propagation of CMEs. In the initiation part, several trigger mechanisms are discussed; In the propagation part, the observations and modelings of EIT waves/dimmings, as the EUV counterparts of CMEs, are described.

There has been the suggestion that the cosmological constant as implied by the dark energy is related to the well-known hierarchy between the Planck scale, $M_{\rm Pl}$, and the Standard Model scale, $M_{\rm SM}$. Here we further propose that the same framework that addresses this hierarchy problem must also address the smallness problem of the cosmological constant. Specifically, we investigate the minimal supersymmetric (SUSY) extension of the Randall-Sundrum model where SUSY-breaking is induced on the TeV brane and transmitted into the bulk. We show that the Casimir energy density of the system indeed conforms with the observed dark energy scale.

The binary YY Gem shows many interesting properties, one of which is the periodicity in its flaring rate. The period, which is about $48 \pm 3$ min, was ever interpreted in terms of the oscillation of a filament. In this paper, we propose a new model to explain this phenomenon by means of 2.5-dimensional MHD numerical simulations. It is found that magnetic reconnection is induced as the coronal loops rooted on both stars inflate and approach each other, which is driven by the differential stellar rotation. The magnetic reconnection is modulated by fast-mode magnetoacoustic waves which are trapped between the surfaces of the two stars, so that the reconnection rate presents a periodic behaviour. With the typical parameters for the binary system, the observed period can be reproduced. We also derive an empirical formula to relate the period of the flaring rate to the coronal temperature and density, as well as the magnetic field.

Magnetowave induced plasma wakefield acceleration (MPWA) in a relativistic astrophysical outflow has been proposed as a viable mechanism for the acceleration of cosmic particles to ultra high energies. Here we present simulation results that clearly demonstrate the viability of this mechanism for the first time. We invoke the high frequency and high speed whistler mode for the driving pulse. The plasma wakefield so induced validates precisely the theoretical prediction. We show that under appropriate conditions, the plasma wakefield maintains very high coherence and can sustain high-gradient acceleration over a macroscopic distance. Invoking gamma ray burst (GRB) as the source, we show that MPWA production of ultra high energy cosmic rays (UHECR) beyond ZeV 10^21 eV is possible.

Measurements are reported of the yield and spectrum of fluorescence, excited by a 28.5 GeV electron beam, in air at a range of pressures of interest to ultra-high energy cosmic ray detectors. The wavelength range was 300 - 420 nm. System calibration has been performed using Rayleigh scattering of a nitrogen laser beam. In atmospheric pressure dry air at 304 K the yield is 20.8 +/- 1.6 photons per MeV.

We point out that during the reionization epoch of the cosmic history, the plasma collective effect among the ordinary matter would suppress the large scale structure formation. The imperfect Debye shielding at finite temperature would induce a residual long-range electrostatic potential which, working together with the baryon thermal pressure, would counter the gravitational collapse. As a result the effective Jean's length, $\tilde{\lambda}_J$, is increased by a factor, $\tilde{\lambda}_J/\lambda_J=\sqrt{8/5}$, relative to the conventional one. For scales smaller than the effective Jean's scale the plasma would oscillate at the ion-acoustic frequency. The modes that would be influenced by this effect depend on the starting time and the initial temperature of reionization, but roughly lie in the range $0.5 h{\rm Mpc}^{-1}< k$, which corresponds to the region of the Lyman-$\alpha$ forest from the inter-galactic medium. We predict that in the linear regime of density-contrast growth, the plasma suppression of the matter power spectrum would approach $1-(\Omega_{dm}/\Omega_m)^2\sim 1-(5/6)^2\sim 30%$.

We report results of a recently-completed pre-Formulation Phase study of SPIRIT, a candidate NASA Origins Probe mission. SPIRIT is a spatial and spectral interferometer with an operating wavelength range 25 - 400 microns. SPIRIT will provide sub-arcsecond resolution images and spectra with resolution R = 3000 in a 1 arcmin field of view to accomplish three primary scientific objectives: (1) Learn how planetary systems form from protostellar disks, and how they acquire their inhomogeneous composition; (2) characterize the family of extrasolar planetary systems by imaging the structure in debris disks to understand how and where planets of different types form; and (3) learn how high-redshift galaxies formed and merged to form the present-day population of galaxies. Observations with SPIRIT will be complementary to those of the James Webb Space Telescope and the ground-based Atacama Large Millimeter Array. All three observatories could be operational contemporaneously.

We investigate a possible new technique for microwave measurements of ultra-high energy cosmic ray (UHECR) extensive air showers which relies on detection of expected continuum radiation in the microwave range, caused by free-electron collisions with neutrals in the tenuous plasma left after the passage of the shower. We performed an initial experiment at the AWA (Argonne Wakefield Accelerator) laboratory in 2003 and measured broadband microwave emission from air ionized via high energy electrons and photons. A follow-up experiment at SLAC (Stanford Linear Accelerator Center) in summer of 2004 confirmed the major features of the previous AWA observations with better precision and made additional measurements relevant to the calorimetric capabilities of the method. Prompted by these results we built a prototype detector using satellite television technology, and have made measurements indicating possible detection of cosmic ray extensive air showers. The method, if confirmed by experiments now in progress, could provide a high-duty cycle complement to current nitrogen fluorescence observations of UHECR, which are limited to dark, clear nights. By contrast, decimeter microwave observations can be made both night and day, in clear or cloudy weather, or even in the presence of moderate precipitation.

We report on the first observations of the Askaryan effect in ice: coherent impulsive radio Cherenkov radiation from the charge asymmetry in an electromagnetic (EM) shower. Such radiation has been observed in silica sand and rock salt, but this is the first direct observation from an EM shower in ice. These measurements are important since the majority of experiments to date that rely on the effect for ultra-high energy neutrino detection are being performed using ice as the target medium. As part of the complete validation process for the Antarctic Impulsive Transient Antenna (ANITA) experiment, we performed an experiment at the Stanford Linear Accelerator Center (SLAC) in June 2006 using a 7.5 metric ton ice target, yielding results fully consistent with theoretical expectations.

The well-known hierarchy between the Planck scale (~10^{19}GeV) and the TeV scale, namely a ratio of ~10^{16} between the two, is coincidentally repeated in a inverted order between the TeV scale and the dark energy scale at \~10^{-3}eV implied by the observations. We argue that this is not a numerical coincidence. The same brane-world setups to address the first hierarchy problem may also in principle address this second hierarchy issue. Specifically, we consider supersymmetry in the bulk and its breaking on the brane and resort to the Casimir energy induced by the bulk graviton-gravitino mass-shift on the brane as the dark energy. For the ADD model we found that our notion is sensible only if the number of extra dimension n=2. We extend our study to the Randall-Sundrum model. Invoking the chirality-flip on the boundaries for SUSY-breaking, the zero-mode gravitino contribution to the Casimir energy does give rise to the double hierarchy. Unfortunately since the higher Kaluza-Klein modes acquire relative mass-shifts at the TeV level, the zero-mode contribution to Casimir energy is overshadowed.

We show that the cascade limit on ultra high energy cosmic neutrino (UHEC$\nu$) flux imposes a lower bound on the neutrino mass provided that super-GZK events of ultra high energy cosmic rays (UHECRs) are produced from Z-bursts. Based on the data from HiRes and AGASA, the obtained neutrino mass lower bound violates its existing cosmological upper bound. We conclude that the Z-bursts cannot be the dominant source for the observed super-GZK UHECR events. This is consistent with the recent ANITA-lite data.

Measurements are reported on the fluorescence of air as a function of depth in electromagnetic showers initiated by bunches of 28.5 GeV electrons. The light yield is compared with the expected and observed depth profiles of ionization in the showers. It validates the use of atmospheric fluorescence profiles in measuring ultra high energy cosmic rays.

We performed magnetohydrodynamic simulation of coronal mass ejections (CMEs) and associated giant arcade formations, and the results suggested new interpretations of observations of CMEs. We performed two cases of the simulation: with and without heat conduction. Comparing between the results of the two cases, we found that reconnection rate in the conductive case is a little higher than that in the adiabatic case and the temperature of the loop top is consistent with the theoretical value predicted by the Yokoyama-Shibata scaling law. The dynamical properties such as velocity and magnetic fields are similar in the two cases, whereas thermal properties such as temperature and density are very different.In both cases, slow shocks associated with magnetic reconnectionpropagate from the reconnection region along the magnetic field lines around the flux rope, and the shock fronts form spiral patterns. Just outside the slow shocks, the plasma density decreased a great deal. The soft X-ray images synthesized from the numerical results are compared with the soft X-ray images of a giant arcade observed with the Soft X-ray Telescope aboard {\it Yohkoh}, it is confirmed that the effect of heat conduction is significant for the detailed comparison between simulation and observation. The comparison between synthesized and observed soft X-ray images provides new interpretations of various features associated with CMEs and giant arcades.

In this paper we study various cosmological solutions for a D3/D7 system directly from M-theory with fluxes and M2-branes. In M-theory, these solutions exist only if we incorporate higher derivative corrections from the curvatures as well as G-fluxes. We take these corrections into account and study a number of toy cosmologies, including one with a novel background for the D3/D7 system whose supergravity solution can be completely determined. This new background preserves all the good properties of the original model and opens up avenues to investigate cosmological effects from wrapped branes and brane-antibrane annihilation, to name a few. We also discuss in some detail semilocal defects with higher global symmetries, for example exceptional ones, that could occur in a slightly different regime of our D3/D7 model. We show that the D3/D7 system does have the required ingredients to realise these configurations as non-topological solitons of the theory. These constructions also allow us to give a physical meaning to the existence of certain underlying homogeneous quaternionic Kahler manifolds.

The observed SEDs of a sample of 60 OH/IR stars are fitted using a radiative transfer model of a dusty envelope. Among the whole sample, 21 stars have reliable phase-lag distances while the others have less accurate distances. L*-P,Mlr-P and Mlr-L* relations have been plotted for these stars. It is found that type E (with emission feature at 10um and type A (with absorption feature at 10um) OH/IR stars have different L*-P and Mlr-L* relations while both of them follow a single Mlr-P relation. The type E stars are proven to be located in the area without large scale dense interstellar medium while the type A stars are located probably in dense interstellar medium. It is argued here that this may indicate the two types of OH/IR stars have different chemical composition or zero age main sequence mass and so evolve in different ways. This conclusion has reinforced the argument by Chen et al.(2001) who reached a similar conclusion from the galactic distribution of about 1000 OH/IR stars with the IRAS low-resolution spectra (LRS).

A new model for the origin of dark energy is proposed based on the Casimir effect in a supersymmetry-breaking brane-world. Supersymmetry is assumed to be preserved in the bulk while broken on a 3-brane. Due to the boundary conditions imposed on the compactified extra dimensions, there is an effective Casimir energy induced on the brane. The net Casimir energy contributed from the graviton and the gravitino modes as a result of supersymmetry-breaking on the brane is identified as the observed dark energy, which in our construction is a cosmological constant. We show that the smallness of the cosmological constant, which results from the huge contrast in the extra-dimensional volumes between that associated with the 3-brane and that of the bulk, is attainable under very relaxed conditions.

Most general relativity textbooks devote considerable space to the simplest example of a black hole containing a singularity, the Schwarzschild geometry. However only a few discuss the dynamical process of gravitational collapse, by which black holes and singularities form. We present here two types of analytic models for this process, which we believe are the simplest available; the first involves collapsing spherical shells of light, analyzed mainly in Eddington-Finkelstein coordinates; the second involves collapsing spheres filled with a perfect fluid, analyzed mainly in Painleve-Gullstrand coordinates. Our main goal is pedagogical simplicity and algebraic completeness, but we also present some results that we believe are new, such as the collapse of a light shell in Kruskal-Szekeres coordinates.

The 1612 MHz hydroxyl maser in circumstellar envelopes has long been thought to be pumped by 34.6um photons. Only recently, the Infrared Space Observatory has made possible spectroscopic observations which enable the direct confirmation of this pumping mechanism in a few cases. To look for the presence of this pumping line, we have searched the Infrared Space Observatory Data Archive and found 178 spectra with data around 34.6um for 87 galactic 1612MHz masers. The analysis performed showed that the noise level and the spectral resolution of the spectra are the most important factors affecting the detection of the 34.6um absorption line. Only 5 objects from the sample (3 red supergiants and 2 galactic center sources) are found to show clear 34.6um absorption (all of them already known) while two additional objects only tentatively show this line. The 3 supergiants show similar pump rates and their masers might be purely radiatively pumped. The pump rates of OH masers in late type stars are found to be about 0.05, only 1/5 of the theoretical value of 0.25 derived by Elitzur (1992). We have also found 16 maser sources which, according to the analysis assuming Elitzur's pump rate, should show the 34.6 $\mu$m absorption line but do not. These non-detections can be tentatively explained by far-infrared photon pumping, clumpy nature of the OH masing region or a limb-filling emission effect in the OH shell.

This paper analyzes the 53.3um line in the ISO LWS spectra towards a similar sample of OH/IR sources. We find 137 LWS spectra covering 53.3um and associated with 47 galactic OH/IR sources. Ten of these galactic OH/IR sources are found to show and another 5 ones tentatively show the 53.3um absorption while another 7 sources highly probably do not show this line. The source class is found to be correlated with the type of spectral profile: red supergiants (RSGs) and AGB stars tend to show strong blue-shifted filling emission in their 53.3um absorption line profiles while HII regions tend to show a weak red-shifted filling emission in the line profile. GC sources and megamasers do not show filling emission feature. It is argued that the filling emission might be the manifestation of an unresolved half emission half absorption profile of the 53.3um doublet. The 53.3 to 34.6um equivalent width (EW) ratio is close to unity for RSGs but much larger than unity for GC sources and megamasers while H II regions only show the 53.3um line. The pump rate defined as maser to IR photon flux ratio is approximately 5% for RSGs. The pump rates of GC sources are three order of magnitude smaller. Both the large 53.3 to 34.6um EW ratio and the small pump rate of the GC OH masers reflect that the two detected `pumping lines' in these sources are actually of interstellar origin. The pump rate of Arp 220 is 32%--much larger than that of RSGs, which indicates that the contribution of other pumping mechanisms to this megamaser is important.

While there exist various candidates, the identification of dark matter remains unresolved. Recently it was argued that the generalized uncertainty principle (GUP) may prevent a black hole from evaporating completely, and as a result there should exist a Planck-size BHR at the end of its evaporation. We speculate that the stability of BHR may be further protected by supersymmetry in the form of extremal black hole. If this is indeed the case and if a sufficient amount of small black holes can be produced in the early universe, then the resultant BHRs can be an interesting candidate for DM. We demonstrate that this is the case in the hybrid inflation model. By assuming BHR as DM, our notion imposes a constraint on the hybrid inflation potential. We show that such a constraint is not fine-tuned. Possible observational signatures are briefly discussed.

We perform a multi-wavelength study of a two-ribbon flare on 2002 September 29 and its associated filament eruption, observed simultaneously in the H-alpha line by a ground-based imaging spectrograph and in hard X-rays by RHESSI. The flare ribbons contain several H-alpha bright kernels that show different evolutional behaviors. In particular, we find two kernels that may be the footpoints of a loop. A single hard X-ray source appears to cover these two kernels and to move across the magnetic neutral line. We explain this as a result of the merging of two footpoint sources that show gradually asymmetric emission owing to an asymmetric magnetic topology of the newly reconnected loops. In one of the H-alpha kernels, we detect a continuum enhancement at the visible wavelength. By checking its spatial and temporal relationship with the hard X-ray emission, we ascribe it as being caused by electron beam precipitation. In addition, we derive the line-of-sight velocity of the filament plasma based on the Doppler shift of the filament-caused absorption in the H-alpha blue wing. The filament shows rapid acceleration during the impulsive phase. These observational features are in principal consistent with the general scenario of the canonical two-ribbon flare model.

There have been proposals that primordial black hole remnants (BHRs) are the dark matter, but the idea is somewhat vague. Recently we argued that the generalized uncertainty principle (GUP) may prevent black holes from evaporating completely, in a similar way that the standard uncertainty principle prevents the hydrogen atom from collapsing. We further noted that the hybrid inflation model provides a plausible mechanism for production of large numbers of small black holes. Combining these we suggested that the dark matter might be composed of Planck-size BHRs. In this paper we briefly review these arguments, and discuss the reheating temperature as a result of black hole evaporation.

There have been proposals that primordial black hole remnants (BHRs) are the dark matter, but the idea is somewhat vague. We argue here first that the generalized uncertainty principle (GUP) may prevent black holes from evaporating completely, in a similar way that the standard uncertainty principle prevents the hydrogen atom from collapsing. Secondly we note that the hybrid inflation model provides a plausible mechanism for production of large numbers of small black holes. Combining these we suggest that the dark matter might be composed of Planck-size BHRs and discuss the possible constraints and signatures associated with this notion.

Recent years have seen tremendous progress in our understanding of the extreme universe, which in turn points to even deeper questions to be further addressed. History has shown that the symbiosis between direct observations and laboratory investigations is instrumental in the progress of astrophysics. Current frontier astrophysical phenomena related to particle astrophysics and cosmology typically involve one or more of the following conditions: (1) extremely high energy events;(2) very high density, high temperature processes; (3) super strong field environments. Laboratory experiments using high intensity lasers and particle beams can calibrate astrophysical observation or detection processes, investigate the underlying dynamics of astrophysical phenomena, and probe fundamental physics in extreme limits. We give examples of possible laboratory experiments that investigate into the extreme universe.

The population inversion which leads to the 1612 MHz OH maser emission has long been thought to be radiatively pumped. Since OH rotational lines involved in this pumping scheme lie in the far-infrared they became observable only after the launch of the ISO satellite. With the aim to investigate the pumping conditions of the 1612 MHz OH maser emission in more details we have searched the ISO Archive for SWS observations around 34.6 um of 1024 OH/IR sources with IRAS LRS spectra from compilation of Chen et al. (2001). Surprisingly, among 81 OH/IR sources which have appropriate SWS data only already reported objects: VY CMa, IRC+10420 and the Galactic center, show clear 34.6 um absorption line. We discuss possible reasons for non-detection of this pumping line.

A cosmic acceleration mechanism is introduced which is based on the wakefields excited by the Alfven shocks in a relativistically flowing plasma, where the energy gain per distance of a test particle is Lorentz invariant. We show that there exists a threshold condition for transparency below which the accelerating particle is collision-free and suffers little energy loss in the plasma medium. The stochastic encounters of the random accelerating-decelerating phases results in a power-law energy spectrum: f(e) 1/e^2. The environment suitable for such plasma wakefield acceleration can be cosmically abundant. As an example, we discuss the possible production of super-GZK ultra high energy cosmic rays (UHECR) through this mechanism in the atmosphere of gamma ray bursts. We show that the acceleration gradient can be as high as G ~ 10^16 eV/cm. The estimated event rate in our model agrees with that from UHECR observations.

A sample of J-type carbon stars was searched for OH maser emission. The new detection of three OH lines towards two silicate carbon stars is reported. In V778 Cyg, previously known as the main-lines (1665 and 1667 MHz) maser source, the satellite 1612 MHz emission was discovered while in NSV 2814 the main OH lines were detected. The presence of OH maser lines confirms the former suggestion that oxygen-rich material is located in the vicinity ($\approx$ $10^{15-16}$ cm) of silicate carbon stars.

The Ultraviolet Imaging Telescope (UIT) was flown as part of the Astro observatory on the Space Shuttle Columbia in December 1990 and again on the Space Shuttle Endeavor in March 1995. Ultraviolet (1200-3300 Angstroms) images of a variety of astronomical objects, with a 40 arcmin field of view and a resolution of about 3 arcsec, were recorded on photographic film. The data recorded during the first flight are available to the astronomical community through the National Space Science Data Center (NSSDC); the data recorded during the second flight will soon be available as well. This paper discusses in detail the design, operation, data reduction, and calibration of UIT, providing the user of the data with information for understanding and using the data. It also provides guidelines for analyzing other astronomical imagery made with image intensifiers and photographic film.