Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica.

Machine learning has emerged as a popular and powerful approach for solving problems in astrophysics. We review applications of machine learning techniques for the analysis of ground-based gravitational-wave detector data. Examples include techniques for improving the sensitivity of Advanced LIGO and Advanced Virgo gravitational-wave searches, methods for fast measurements of the astrophysical parameters of gravitational-wave sources, and algorithms for reduction and characterization of non-astrophysical detector noise. These applications demonstrate how machine learning techniques may be harnessed to enhance the science that is possible with current and future gravitational-wave detectors.

We study surface effects of neutron $^{3}P_{2}$ superfluids in neutron stars. $^{3}P_{2}$ superfluids are in uniaxial nematic (UN), D$_{2}$ biaxial nematic (BN), or D$_{4}$ BN phase, depending on the strength of magnetic fields from small to large. We suppose a neutron $^{3}P_{2}$ superfluid in a ball with a spherical boundary. Adopting a suitable boundary condition for $^{3}P_{2}$ condensates, we solve the Ginzburg-Landau equation to find several surface properties for the neutron $^{3}P_{2}$ superfluid. First, the phase on the surface can be different from that of the bulk, and symmetry restoration or breaking occurs in general on the surface. Second, the distribution of the surface energy density has an anisotropy depending on the polar angle in the sphere, which may lead to the deformation of the geometrical shape of the surface. Third, the order parameter manifold induced on the surface, which is described by two-dimensional vector fields induced on the surface from the condensates, allows topological defects (vortices) on the surface, and there must exist such defects even in the ground state thanks to the Poincar\'{e}-Hopf theorem: although the numbers of the vortices and antivortices depend on the bulk phases, the difference between them is topologically invariant (the Euler number $\chi=2$) irrespective of the bulk phases. These vortices, which are not extended to the bulk, are called boojums in the context of liquid crystals and helium-3 superfluids. The surface properties of the neutron $^{3}P_{2}$ superfluid found in this paper may provide us useful information to study neutron stars.

In this paper, we report on the construction of a deep Artificial Neural Network (ANN) to localize simulated gravitational wave signals in the sky with high accuracy. We have modelled the sky as a sphere and have considered cases where the sphere is divided into 18, 50, 128, 1024, 2048 and 4096 sectors. The sky direction of the gravitational wave source is estimated by classifying the signal into one of these sectors based on it's right ascension and declination values for each of these cases. In order to do this, we have injected simulated binary black hole gravitational wave signals of component masses sampled uniformly between 30-80 solar mass into Gaussian noise and used the whitened strain values to obtain the input features for training our ANN. We input features such as the delays in arrival times, phase differences and amplitude ratios at each of the three detectors Hanford, Livingston and Virgo, from the raw time-domain strain values as well as from analytical versions of these signals, obtained through Hilbert transformation. We show that our model is able to classify gravitational wave samples, not used in the training process, into their correct sectors with very high accuracy (>90%) for coarse angular resolution using 18, 50 and 128 sectors. We also test our localization on test samples with injection parameters of the published LIGO binary black hole merger events GW150914, GW170818 and GW170823 for 1024, 2048 and 4096 sectors and compare the result with that from BAYESTAR and Parameter Estimation (PE). In addition, we report that the time taken by our model to localize one GW signal is around 0.018 secs on 14 Intel Xeon CPU cores.

The Ginzburg-Landau (GL) effective theory is a useful tool to study a superconductivity or superfluidity near the critical temperature, and usually the expansion up to the 4th order in terms of order parameters is sufficient for the description of the second-order phase transition. In this paper, we discuss the GL equation for the neutron $^{3}P_{2}$ superfluidity relevant for interior of neutron stars. We derive the GL expansion up to the 8th order in the condensates and find that this order is necessary for the system to have the unique ground state, unlike the ordinary cases. Starting from the $LS$ potential, which provides the dominant attraction between two neutrons at the high density, we derive the GL equation in the path-integral formalism, where the auxiliary field method and the Nambu-Gor'kov representation are used. We present the detailed description for the trace calculation necessary in the derivation of the GL equation. As numerical results, we show the phase diagram of the neutron $^{3}P_{2}$ superfluidity on the plane spanned by the temperature and magnetic field, and find that the 8th order terms lead to a first-order phase transition, whose existence was predicted in the Bogoliubov-de Gennes equation but has not been found thus far within the framework of the GL expansion up to the 6th order.The first-order phase transition will affect the interior structures inside the neutron stars.

We discuss the effect of a strong magnetic field on neutron $^{3}P_{2}$ superfluidity. Based on the attraction in the $^{3}P_{2}$ pair of two neutrons, we derive the Ginzburg-Landau equation in the path-integral formalism by adopting the bosonization technique and leave the next-to-leading order in the expansion of the magnetic field $B$. We determine the $(T,B)$ phase diagram with temperature $T$, comprising three phases: the uniaxial nematic (UN) phase for $B=0$, D$_{2}$-biaxial nematic (BN) and D$_{4}$-BN phases in finite $B$ and strong $B$ such as magnetars, respectively, where D$_{2}$ and D$_{4}$ are dihedral groups. We find that, compared with the leading order in the magnetic field known before, the region of the D$_{2}$-BN phase in the $(T,B)$ plane is extended by the effect of the next-to-leading-order terms of the magnetic field. We also present the thermodynamic properties, such as heat capacities and spin susceptibility, and find that the spin susceptibility exhibits anisotropies in the UN, D$_{2}$-BN, and D$_{4}$-BN phases. This information will be useful to understand the internal structures of magnetars.

We discuss neutron $^{3}P_{2}$ phases in the core of neutron stars in strong magnetic field (magnetars). The neutron $^{3}P_{2}$ pairing provides a wide variety of condensates, such as the uniaxial nematic and (D$_{2}$ and D$_{4}$) biaxial nematic, with different symmetries stemming from the combinations of spin and momentum. Based on the spin-orbital angular momentum coupling and the spin-magnetic field coupling of the neutrons, we derive the Ginzburg-Landau equation containing higher order terms of the magnetic field. We investigate the phase diagram of the neutron $^{3}P_{2}$ superfluidity, and find that the D$_{2}$ biaxial nematic phase is extended by the higher order terms of the magnetic field. We also discuss the thermodynamic properties, the heat capacity and the spin susceptibility.

Axion is a promising candidate of dark matter. After the Peccei-Quinn symmetry breaking, axion strings are formed and attached by domain walls when the temperature of the universe becomes comparable to the QCD scale. Such objects can cause cosmological disasters if they are long-lived. As a solution for it, the Lazarides-Shafi mechanism is often discussed through introduction of a new non-Abelian (gauge) symmetry. We study this mechanism in detail and show configuration of strings and walls. Even if Abelian axion strings with a domain wall number greater than one are formed in the early universe, each of them is split into multiple Alice axion strings due to a repulsive force between the Alice strings even without domain wall. When domain walls are formed as the universe cools down, a single Alice string can be attached by a single wall because a vacuum is connected by a non-Abelian rotation without changing energy. Even if an Abelian axion string attached by domain walls are created due to the Kibble Zurek mechanism at the chiral phase transition, such strings are also similarly split into multiple Alice strings attached by walls in the presence of the domain wall tension. Such walls do not form stable networks since they collapse by the tension of the walls, emitting axions.

Quark-hadron continuity was proposed as crossover between hadronic matter and quark matter without a phase transition, based on the matching of the symmetry and excitations in both phases. In the limit of a light strange-quark mass, it connects hyperon matter and the color-flavor-locked (CFL) phase exhibiting color superconductivity. Recently, it was proposed that this conjecture could be generalized in the presence of superfluid vortices penetrating both phases (arXiv:1803.05115 [hep-ph]), and it was suggested that one hadronic superfluid vortex in hyperon matter could be connected to one non-Abelian vortex (color magnetic flux tube) in the CFL phase. Here, we argue that their proposal is consistent only at large distances; instead, we show that three hadronic superfluid vortices must combine with three non-Abelian vortices with different colors with the total color magnetic fluxes canceled out, where the junction is called a colorful boojum. We rigorously prove this in both a macroscopic theory based on the Ginzburg-Landau description in which symmetry and excitations match (including vortex cores), and a microscopic theory in which the Aharonov-Bohm phases of quarks around vortices match.

The quark-hadron crossover conjecture was proposed as a continuity between hadronic matter and quark matter with no phase transition. It is based on matching of symmetry and excitations in both the phases. It connects hyperon matter and color-flavor locked (CFL) phase of color superconductivity in the limit of light strange quark mass. We study generalization of this conjecture in the presence of topological vortices. We propose a picture where hadronic superfluid vortices in hyperon matter could be connected to non-Abelian vortices (color magnetic flux tubes) in the CFL phase during this crossover. We propose that three hadronic superfluid vortices must join together to three non-Abelian vortices with different color fluxes with the total color magnetic fluxes canceled out, where the junction is called a colorful boojum.

The LambdaCDM model of cosmology, though very successful at large scales, has some discrepancy with observations at the galactic and sub-galactic scales. These include the core-cusp problem, missing satellites problem etc. Spergel and Steingardt (2000) proposed that if dark matter undergoes feeble self interactions with each other, then such problems can be averted. In this thesis, a two-component Feebly Interacting Massive Particle (FIMP) dark matter model involving two singlet scalar fields capable of self-interactions has been proposed and its impact on large scale structure formation has been studied through cosmological simulations. The proposed model involves simple extensions of the Standard Model with two singlet scalar fields formed non-thermally through the decay of heavier particles in the very early universe. These particles acquire their relic abundance through a freeze-in mechanism. The coupled Boltzmann equation of the FIMP-FIMP model was solved and the relic densities for different values of the coupling parameters were obtained and matched with PLANCK results. The impact of dark matter self interactions was studied through cosmological simulations using a modified version of the parallel TreePM code GADGET-2 and the halo mass function and halo catalog for different dark matter self interaction cross sections were obtained. Lastly, the newly developed Effective Theory of Structure Formation (ETHOS) framework which is a new and innovative paradigm in the study of the cosmological effects of different dark matter models was studied and using the public code, ETHOS-CAMB the signatures of dark acoustic oscillations in the matter power spectrum for a particular dark matter model was obtained.

We discuss collective excitations (both fundamental and solitonic excitations) of quantized superfluid vortices in neutron $^3P_2$ superfluids, which likely exist in high density neutron matter such as neutron stars. Besides the well-known Kelvin modes (translational zero modes), we find a gapfull mode whose low-energy description takes the simple form of a double sine-Gordon model. The associated kink solution and its effects on spontaneous magnetization inside the vortex core are analyzed in detail.