The simplest model of mirror sector dark matter maintains exact mirror symmetry, but has a baryon abundance $\Omega_{b'} = \beta \Omega_b$ and a suppressed temperature $T' = x T$ in the mirror sector; hence it depends only on two parameters, $\beta,x$. For sufficiently small $x$, early cosmological observables may not constrain mirror baryons from constituting all of the dark matter despite their strong self-interactions, depending on the unknown details of structure formation in the hidden sector. Here we close this loophole by simulating mirror structure formation, mapping out the allowed regions of parameter space using cosmological and astronomical data. We find that the Milky Way disk surface density and bulge mass constrain $\Omega_{b'}\lesssim 0.3 \Omega_{b}$ at the highest $T'$ allowed by BBN and CMB ($T'=0.5 T$), or $\Omega_{b'}\lesssim 0.8 \Omega_{b}$ at lower values of $T'$. We also briefly discuss the realization of the necessary temperature asymmetry between the SM and the mirror sector in our model with unbroken mirror symmetry.

Small-scale magnetic flux ropes, in the solar wind, have been studied for decades via the approach of both simulation and observation. Statistical analysis utilizing various in-situ spacecraft measurements is the main observational approach. In this study, we extend the automated detection of small-scale flux ropes based on the Grad-Shafranov reconstruction to the complete dataset of \emph{Ulysses} spacecraft in-situ measurements. We first discuss the temporal variation of the bulk properties of 22,719 flux ropes found through our approach, namely, the average magnetic field and plasma parameters, etc., as functions of the heliographical latitudes and heliocentric radial distances. We then categorize all identified events into three groups based on event distributions in different latitudes separated at 30$^{\circ}$, at different radial distances, and under different solar activities, respectively. By the detailed statistical analysis, we conclude as follows. (1) The properties of flux ropes, such as the duration, scale size, etc., follow the power-law distributions, but with different slope indices, especially for distributions at different radial distances. (2) Also, they are affected by the solar wind speed which has different distributions under the different solar activities, which is manifested as the latitudinal effect. (3) The main difference in flux rope properties between the low and high latitudes is attributed to possible Alfv\'enic structures or waves. (4) Flux ropes with longer duration and larger scale sizes occur more often at larger radial distances. (5) With more strict Wal\'en slope threshold, more events are excluded at higher latitudes. The entire database is published online at \url{http://www.fluxrope.info}.

We have developed a new automated small-scale magnetic flux rope (SSMFR) detection algorithm based on the Grad-Shafranov (GS) reconstruction technique. We have applied this detection algorithm to the Wind spacecraft in-situ measurements during 1996 - 2016, covering two solar cycles, and successfully detected a total number of 74,241 small-scale magnetic flux rope events with duration from 9 to 361 minutes. This large number of small-scale magnetic flux ropes has not been discovered by any other previous studies through this unique approach. We perform statistical analysis of the small-scale magnetic flux rope events based on our newly developed database, and summarize the main findings as follows. (1) The occurrence of small-scale flux ropes has strong solar cycle dependency with a rate of a few hundreds per month on average. (2) The small-scale magnetic flux ropes in the ecliptic plane tend to align along the Parker spiral. (3) In low speed ($<$ 400 km/s) solar wind, the flux ropes tend to have lower proton temperature and higher proton number density, while in high speed ($\ge$ 400 km/s) solar wind, they tend to have higher proton temperature and lower proton number density. (4) Both the duration and scale size distributions of the small-scale magnetic flux ropes obey a power law. (5) The waiting time distribution of small-scale magnetic flux ropes can be fitted by an exponential function (for shorter waiting times) and a power law function (for longer waiting times). (6) The wall-to-wall time distribution obeys double power laws with the break point at 60 minutes (corresponding to the correlation length). (7) The small-scale magnetic flux ropes tend to accumulate near the heliospheric current sheet (HCS). The entire database is available at \url{http://fluxrope.info} and in machine readable format in this article.

This paper summarizes the results obtained by the team "Heliosheath Processes and the Structure of the Heliopause: Modeling Energetic Particles, Cosmic Rays, and Magnetic Fields" supported by the International Space Science Institute in Bern, Switzerland. We focus on the physical processes occurring in the outer heliosphere, especially at its boundary called the heliopause (HP), and in the LISM. The importance of magnetic field, charge exchange between atoms and ions, and solar cycle on the heliopause topology and observed heliocentric distances to different heliospheric discontinuities are discussed. It is shown that time-dependent boundary conditions are necessary to describe the heliospheric asymmetries detected by the Voyager spacecraft. We also discuss the structure of the HP, especially due to its instability and magnetic reconnection. It is demonstrated that the Rayleigh-Taylor instability of the nose of the HP creates consecutive layers of the interstellar and heliospheric plasma which are magnetically connected to different sources. This may be a possible explanation of abrupt changes in the galactic and anomalous cosmic ray fluxes observed by Voyager 1 when it was crossing the HP structure for a period of about one month in the summer of 2012. This paper also discusses the plausibility of fitting simulation results to a number of observational data sets obtained by in situ and remote measurements. The distribution of magnetic field in the vicinity of the HP is discussed in the context of Voyager measurements. We discuss the transport of energetic particles in the inner and outer heliosheath, concentrating on the anisotropic spatial diffusion diffusion tensor and the pitch-angle dependence of perpendicular diffusion and demonstrate that the latter can explain the observed pitch-angle anisotropies of both the anomalous and galactic cosmic rays in the outer heliosheath.

The transport of charged energetic particles in the presence of strong intermittent heliospheric turbulence is computationally analyzed based on known properties of the interplanetary magnetic field and solar wind plasma at 1 Astronomical Unit (AU). The turbulence is assumed to be static, composite, and quasi-three-dimensional with a varying energy distribution between a one-dimensional Alfv\'enic (slab) and a structured two-dimensional component. The spatial fluctuations of the turbulent magnetic field are modeled either as homogeneous with a Gaussian probability distribution function (PDF), or as intermittent on large and small scales with a q-Gaussian PDF. Simulations showed that energetic particle diffusion coefficients both parallel and perpendicular to the background magnetic field are significantly affected by intermittency in the turbulence. This effect is especially strong for parallel transport where for large-scale intermittency results show an extended phase of subdiffusive parallel transport during which cross-field transport diffusion dominates. The effects of intermittency are found to depend on particle rigidity and the fraction of slab energy in the turbulence, yielding a perpendicular to parallel mean free path ratio close to 1 for large-scale intermittency. Investigation of higher order transport moments (kurtosis) indicates that non-Gaussian statistical properties of the intermittent turbulent magnetic field are present in the parallel transport, especially for low rigidity particles at all times.

We investigate heliosheath energetic neutral atom (ENA) fluxes at keV energies, by assuming that the heliosheath proton distribution can be approximated by a kappa-distribution. The choice of the kappa parameter derives from observational data of the solar wind (SW). This has direct applications to the upcoming IBEX mission. We will look at all-sky ENA maps within the IBEX energy range (10 eV to 6 keV), as well as ENA energy spectra in several directions. We find that the use of kappa, as opposed to a Maxwellian, gives rise to greatly increased ENA fluxes above 1 keV, while medium energy fluxes are somewhat reduced. We show how IBEX data can be used to estimate the spectral slope in the heliosheath, and that the use of kappa reduces the differences between ENA maps at different energies. We also investigate the effect introducing a kappa-distribution has on the global interaction between the SW and the local interstellar medium (LISM), and find that there is generally an increase in energy transport from the heliosphere into the LISM, due to the modified profile of ENA's energies. This results in a termination shock that moves out by 4 AU, a heliopause that moves in by 9 AU and a bow shock 25 AU farther out, in the nose direction.