This is a community draft white paper for submission to the Decadal Survey in Planetary Science and Astrobiology, reflecting the views of the NASA Astrobiology Program's Research Coordination Network for Ocean Worlds (NOW). We recommend the establishment of a dedicated Ocean Worlds Exploration Program within NASA, providing sustained funding to support the science, engineering, research, development, and mission planning needed to implement a multi-decadal, multi-mission program to explore Ocean Worlds for life and understand the conditions for habitability. The two new critical flagship missions within this program would 1) land on Europa or Enceladus in the decade 2023-2032 to investigate geophysical and geochemical environments while searching for biosignatures, and 2) access a planetary ocean to directly search for life in the decade 2033-2042. The technological solutions for a landed mission are already in-hand, evidenced by the successful delta-Mission Concept Review of the Europa Lander pre-flight project in the fall of 2018. Following an initial landed mission, an ocean access mission will require substantial research, development, and analog testing this decade enable the initiation of a pre-flight project at the start of the following decade.

The first results of a new three-dimensional, finite temperature Skyrme-Hartree-Fock+BCS study of the properties of inhomogeneous nuclear matter at densities and temperatures leading to the transition to uniform nuclear matter are presented. Calculations are carried out in a cubic box representing a unit cell of the locally periodic structure of the matter. A constraint is placed on the two independent components of the quadrupole moment of the neutron density in order to investigate the dependence of the total energy-density of matter on the geometry of the nuclear structure in the unit cell. This approach allows self-consistent modeling of effects such as (i) neutron drip, resulting in a neutron gas external to the nuclear structure, (ii) shell effects of bound and unbound nucleons, (iii) the variety of exotic nuclear shapes that emerge, collectively termed `nuclear pasta' and (iv) the dissolution of these structures into uniform nuclear matter as density and/or temperature increase. In part I of this work the calculation of the properties of inhomogeneous nuclear matter in the core collapse of massive stars is reported. Calculations are performed at baryon number densities of $n_{\rm b}$ = 0.04 - 0.12 fm$^{\rm -3}$, a proton fraction of $y_{\rm p}=0.3$ and temperatures in the range 0 - 7.5 MeV. A wide variety of nuclear shapes are shown to emerge. It is suggested that thermodynamical properties change smoothly in the pasta regime up to the transition to uniform matter; at that transition, thermodynamic properties of the matter vary discontinuously.