In this review we highlight a few physical properties of neutron stars and their theoretical treatment inasmuch as they can be useful for nuclear and particle physicists concerned with matter at finite density (and newly, temperature). Conversely, we lay out some of the hadron physics necessary to test General Relativity with binary mergers including at least one neutron star, in view of the event GW170817: neutron stars and their mergers reach the highest matter densities known, offering access to the matter side of Einstein's equations. In addition to minimum introductory material for those interested in starting research in the field of neutron stars, we dedicate quite some effort to a discussion of the Equation of State of hadron matter in view of gravitational wave developments; we address phase transitions and how the new data may help; we show why transport is expected to be dominated by turbulence instead of diffusion through most if not all of the star, in view of the transport coefficients that have been calculated from microscopic hadron physics; and we relate many of the interesting physics topics in neutron stars to the radius and tidal deformability.

We contribute a publicly available set of tables and code to provide Equations of State (EoS) for matter at neutron star densities. Our EoSes are constrained only by input from hadron physics and fundamental principles, without feedback from neutron star observations, and so without relying on General Relativity. They can therefore be used to test General Relativity itself, as well as modified gravity theories, with neutron star observables, without logical circularity. We have adapted state of the art results from NN chiral potentials for the low--density limit, pQCD results for the asymptotically high-density EoS, and use monotonicity and causality as the only restrictions for intermediate densities, for the EoS sets to remain as model-independent as is feasible today.

In the context of f(R) gravity theories, we show that the apparent mass of a neutron star as seen from an observer at infinity is numerically calculable but requires careful matching, first at the star's edge, between interior and exterior solutions, none of them being totally Schwarzschild-like but presenting instead small oscillations of the curvature scalar R; and second at large radii, where the Newtonian potential is used to identify the mass of the neutron star. We find that for the same equation of state, this mass definition is always larger than its general relativistic counterpart. We exemplify this with quadratic $R^2$ and Hu-Sawicki-like modifications of the standard General Relativity action. Therefore, the finding of two-solar mass neutron stars basically imposes no constraint on stable f(R) theories. However, star radii are in general smaller than in General Relativity, which can give an observational handle on such classes of models at the astrophysical level. Both larger masses and smaller matter radii are due to much of the apparent effective energy residing in the outer metric for scalar-tensor theories. Finally, because the f(R) neutron star masses can be much larger than General Relativity counterparts, the total energy available for radiating gravitational waves could be of order several solar masses, and thus a merger of these stars constitutes an interesting wave source.

We point out that in the early universe, for temperatures in the approximate interval 175-80 MeV (after the quark-gluon plasma), pions carried a large share of the entropy and supported the largest inhomogeneities. Thus, we examine the production of entropy in a pion gas, particularizing to inhomogeneities of the temperature, for which we benefit from the known thermal conductivity. We finally put that entropy produced in relaxing such thermal inhomogeneities in the broad context of this relatively unexplored phase of early-universe cosmology.

We hypothesize that in a non-metallic crystalline structure under extreme pressures, atomic wavefunctions deform to adopt a reduced rotational symmetry consistent with minimizing interstitial space in the crystal. We exemplify with a simple numeric variational calculation that yields the energy cost of this deformation for Helium to 25%. Balancing this with the free energy gained by tighter packing we obtain the pressures required to effect such deformation. The consequent modification of the structure suggests a decrease in the resistance to tangential stress, and an associated decrease of the crystal's shear modulus. The atomic form factor is also modified. We also compare with neutron matter in the interior of compact stars.

We highlight the progress, current status, and open challenges of QCD-driven physics, in theory and in experiment. We discuss how the strong interaction is intimately connected to a broad sweep of physical problems, in settings ranging from astrophysics and cosmology to strongly-coupled, complex systems in particle and condensed-matter physics, as well as to searches for physics beyond the Standard Model. We also discuss how success in describing the strong interaction impacts other fields, and, in turn, how such subjects can impact studies of the strong interaction. In the course of the work we offer a perspective on the many research streams which flow into and out of QCD, as well as a vision for future developments.

The neutron is largely spherical and incompressible in atomic nuclei. These two properties are however challenged in the extreme pressure environment of a neutron star. Our variational computation within the Cornell model of Coulomb gauge QCD shows that the neutron (and also the Delta-3/2 baryon) can adopt cubic symmetry at an energy cost of about 150 MeV. Balancing this with the free energy gained by tighter neutron packing, we expose the possible softening of the equation of state of neutron matter.

If cooled-down neutron stars have a thin atomic crystalline-iron crust, they must diffract X-rays of appropriate wavelength. If the diffracted beam is to be visible from Earth, the illuminating source must be very intense and near the reflecting star. An example is a binary system composed of two neutron stars in close orbit, one of them inert, the other an X-ray pulsar (perhaps an "anomalous" X-ray pulsar or magnetar, not powered by gas absorption from the companion or surrounding space, would be the cleanest example). The observable to be searched for is a secondary peak added (quasi-) periodically to the main X-ray pulse. The distinguishing feature of this secondary peak is that it appears at wavelengths related by simple integer numbers, lambda, lambda/2, lambda/3... lambda/n because of Bragg's diffraction law.