Pasta Structures Research Articles

Shear and interface modes in neutron stars with pasta structures

We carefully examine the shear and interface modes, which are excited due to the presence of crust elasticity, in neutron stars with pasta structures, adopting the relativistic Cowling approximation. We find that the shear modes are independent of the presence of the cylindrical-hole and spherical-hole nuclei at least up to a few kilohertz, while the interface modes strongly depend on the presence of the cylindrical-hole and spherical-hole nuclei. In addition, we find empirical relations for the interface mode frequencies multiplied by the stellar mass and for the shear mode frequencies multiplied by the stellar radius. These relations are expressed as a function of the stellar compactness almost independently of the stiffness in a higher-density region inside the neutron star, once one selects the crust equation of state. Thus, if one would simultaneously observe the shear and interface modes from a neutron star, one might extract the neutron star mass and radius with the help of the constraint on the crust stiffness obtained from terrestrial experiments.

Anisotropic electron transport in the nuclear pasta phase

ABSTRACTThe presence of nuclear pasta is expected to modify the transport properties in the mantle of neutron stars. The non-spherical geometry of the pasta nuclear clusters leads to anisotropies in the collision frequencies, impacting the thermal and electrical conductivity. We derive analytical expressions for the anisotropic collision frequencies using the Boltzmann equation in the relaxation time approximation. The average parallel, perpendicular, and Hall electrical conductivities are computed in the high-temperature regime above crustal melting, considering incoherent elastic electron-pasta scattering and randomly oriented pasta structures. Numerical values are obtained at different densities and temperatures by using the Indiana University–Florida State University (IUFSU) parametrization of the non-linear Walecka model to determine the crustal structure. We find that the anisotropy of the collision frequencies grows with the length of the pasta structures and, independently of the magnetic field, the presence of rod and slab phases decreases the conductivity by more than one order of magnitude. Our numerical results indicate that, even if the pasta structures might survive above the crustal melting point, no strong anisotropies are to be expected in the conduction properties in this temperature regime, even in the presence of a very high magnetic field.

Magnetized neutron star crust within effective relativistic mean-field model

Even though the crystallize nature of the neutron star crust plays a pivotal role in describing various fascinating astrophysical observations, its microscopic structure is not fully understood in the presence of a colossal magnetic field. In the present work, we study the crustal properties of a neutron star within an effective relativistic mean field framework in the presence of magnetic field strength $\sim 10^{17}$G. We calculate the equilibrium composition of the outer crust by minimizing the Gibbs free energy using the most recent atomic mass evaluations. The magnetic field significantly affects the equation of state (EoS) and the properties of the outer crust, such as neutron drip density, pressure, and melting temperature. For the inner crust, we use the compressible liquid drop model for the first time to study the crustal properties in a magnetic environment. The inner crust properties, such as mass and charge number distribution, isospin asymmetry, cluster density, etc., show typical quantum oscillations (De Haas-van Alphen effect) sensitive to the magnetic field's strength. The density-dependent symmetry energy influences the magnetic inner crust like the field-free case. We study the probable modifications in the pasta structures and it is observed that their mass and thickness changes by $\sim 10-15 \%$ depending upon the magnetic field strength. The fundamental torsional oscillation mode frequency is investigated for the magnetized crust in the context of quasiperiodic oscillations (QPO) in soft gamma repeaters. The magnetic field strengths considered in this work influences only the EoS of outer and shallow regions of the inner crust, which results in no significant change in global neutron star properties. However, the outer crust mass and its moment of inertia increase considerably with increase in magnetic field strength.

Nuclear pasta structures at high temperatures

We investigate nuclear pasta structures at high temperatures in the framework of relativistic mean field model with Thomas-Fermi approximation. Typical pasta structures (droplet, rod, slab, tube, and bubble) are obtained, which form various crystalline configurations. The properties of those nuclear pastas are examined in a three-dimensional geometry with reflection symmetry, where the optimum lattice constants are fixed by reproducing the droplet/bubble density that minimizes the free energy adopting spherical or cylindrical approximations for Wigner-Seitz cells. It is found that different crystalline structures can evolve into each other via volume conserving deformations. For fixed densities and temperatures, the differences of the free energies per baryon of nuclear pasta in various shapes and lattice structures are typically on the order of tens of keV, suggesting the possible coexistence of those structures. As temperature increases, the thermodynamic fluctuations are expected to disrupt the long-range ordering in nuclear pasta structures. We then estimate the critical conditions for nuclear pasta to become disordered and behave like liquid, which are found to be sensitive to the densities, temperatures, proton fractions, and nuclear shapes. If we further increase temperature, eventually the nonuniform structures of nuclear pasta become unstable and are converted into uniform nuclear matter. The phase diagrams of nuclear matter are then estimated, which should be useful for understanding the evolutions of neutron stars, supernova dynamics, and binary neutron star mergers.

Pasta properties of the neutron star within effective relativistic mean-field model

We study the properties of pasta structures and their influence on the neutron star observables employing the effective relativistic mean-field (E-RMF) model. The compressible liquid drop model is used to incorporate the finite size effects, considering the possibility of nonspherical structures in the inner crust. The unified equation of states are constructed for several E-RMF parameters to study various properties such as pasta mass and thickness in the neutron star's crust. The majority of the pasta properties are sensitive to the symmetry energy in the subsaturation density region. Using the results from Monte Carlo simulations, we estimate the shear modulus of the crust in the context of quasiperiodic oscillations from soft gamma-ray repeaters and calculate the frequency of fundamental torsional oscillation mode in the inner crust. Global properties of the neutron star such as mass-radius profile, the moment of inertia, crustal mass, crustal thickness, and fractional crustal moment of inertia are worked out. The results are consistent with various observational and theoretical constraints.

Pasta phases in neutron stars under strong magnetic fields

In the present work, we consider nuclear matter in the innermost crust of neutron stars under the presence of a strong magnetic field within the framework of a relativistic mean-field description. Two models with a different slope of the symmetry energy are considered in order to discuss the density dependence of the equation of state on the crust structure. The nonhomogeneous matter in $\ensuremath{\beta}$ equilibrium is described within the coexisting phases method, and the effect of including the anomalous magnetic moment is discussed. Five different geometries for the pasta structures are considered. It is shown that strong magnetic fields cause an extension of the inner crust of the neutron stars, with the occurrence of a series of disconnected nonhomogeneous matter regions above the one existing for a null magnetic field. Moreover, we observed that in these disconnected regions, for some values of the magnetic field, all five different cluster geometrical shapes occur, and the gas density is close to the cluster density. Also, the pressure at the neutron star crust-core transition is much larger than the pressure obtained for a zero magnetic field. Another noticeable effect of the presence of strong magnetic fields is the increase of the proton fraction, favoring the appearance of protons in the gas background.

Nuclear pasta structures and symmetry energy

In the framework of the relativistic mean field model with Thomas-Fermi approximation, we study the structures of low density nuclear matter in a three-dimensional geometry with reflection symmetry. The numerical accuracy and efficiency are improved by expanding the mean fields according to fast cosine transformation and considering only one octant of the unit cell. The effect of finite cell size is treated carefully by searching for the optimum cell size. Typical pasta structures (droplet, rod, slab, tube, and bubble) arranged in various crystalline configurations are obtained for both fixed proton fractions and $\beta$-equilibration. It is found that the properties of droplets/bubbles are similar in body-centered cubic (BCC) and face-centered cubic (FCC) lattices, where the FCC lattice generally becomes more stable than BCC lattice as density increases. For the rod/tube phases, the honeycomb lattice is always more stable than the simple one. By introducing an $\omega$-$\rho$ cross coupling term, we further examine the pasta structures with a smaller slope of symmetry energy $L = 41.34$ MeV, which predicts larger onset densities for core-crust transition and non-spherical nuclei. Such a variation due to the reduction of $L$ is expected to have impacts on various properties in neutron stars, supernova dynamics, and binary neutron star mergers.

Nuclear pasta and symmetry energy in the relativistic point-coupling model

Nonuniform structure of low-density nuclear matter, known as nuclear pasta, is expected to appear not only in the inner crust of neutron stars but also in core-collapse supernova explosions and neutron-star mergers. We perform fully three-dimensional calculations of inhomogeneous nuclear matter and neutron-star matter in the low-density region using the Thomas-Fermi approximation. The nuclear interaction is described in the relativistic mean-field approach with the point-coupling interaction, where the meson exchange in each channel is replaced by the contact interaction between nucleons. We investigate the influence of nuclear symmetry energy and its density dependence on pasta structures by introducing a coupling term between the isoscalar-vector and isovector-vector interactions. It is found that the properties of pasta phases in the neutron-rich matter are strongly dependent on the symmetry energy and its slope. In addition to typical shapes like droplets, rods, slabs, tubes, and bubbles, some intermediate pasta structures are also observed in cold stellar matter with a relatively large proton fraction. We find that nonspherical shapes are unlikely to be formed in neutron-star crusts, since the proton fraction obtained in $\beta$ equilibrium is rather small. The inner crust properties may lead to a visible difference in the neutron-star radius.

Properties of nuclear pastas

In this review we study the nuclear pastas as they are expected to be formed in neutron star crusts. We start with a study of the pastas formed in nuclear matter (composed of protons and neutrons), we follow with the role of the electron gas on the formation of pastas, and we then investigate the pastas in neutron star matter (nuclear matter embedded in an electron gas).Nuclear matter (NM) at intermediate temperatures (1 MeV ≲ T ≲ 15 MeV), at saturation and sub-saturation densities, and with proton content ranging from 30% to 50% was found to have liquid, gaseous and liquid-gas mixed phases. The isospin-dependent phase diagram was obtained along with the critical points, and the symmetry energy was calculated and compared to experimental data and other theories. At low temperatures (T ≲ 1 MeV) NM produces crystal-like structures around saturation densities, and pasta-like structures at sub-saturation densities. Properties of the pasta structures were studied with cluster-recognition algorithms, caloric curve, the radial distribution function, the Lindemann coefficient, Kolmogorov statistics, Minkowski functionals; the symmetry energy of the pasta showed a connection with its morphology.Neutron star matter (NSM) is nuclear matter embedded in an electron gas. The electron gas is included in the calculation by the inclusion of an screened Coulomb potential. To connect the NM pastas with those in neutron star matter (NSM), the role the strength and screening length of the Coulomb interaction have on the formation of the pastas in NM was investigated. Pasta was found to exist even without the presence of the electron gas, but the effect of the Coulomb interaction is to form more defined pasta structures, among other effects. Likewise, it was determined that there is a minimal screening length for the developed structures to be independent of the cell size.Neutron star matter was found to have similar phases as NM, phase transitions, symmetry energy, structure function and thermal conductivity. Like in NM, pasta forms at around T ≈ 1.5 MeV, and liquid-to-solid phase changes were detected at T ≈ 0.5 MeV. The structure function and the symmetry energy were also found to depend on the pasta structures.

A mixing interpolation method to mimic pasta phases in compact star matter

We present a new method to interpolate between two matter phases that allows for a description of mixed phases and can be used, e.g., for mimicking transitions between pasta structures occurring in the crust as well as in the inner core of compact stars. This interpolation method is based on assuming switch functions that are used to define a mixture of subphases while fulfilling constraints of thermodynamic stability. The width of the transition depends on a free parameter, the pressure increment relative to the critical pressure of a Maxwell construction. As an example we present a trigonometric function ansatz for the switch function together with a pressure increment during the transition. We note that the resulting mixed phase equation of state bears similarities with the appearance of substitutional compounds in neutron star crusts and with the sequence of transitions between different pasta phases in the hadron-to-quark matter transition. We apply this method to the case of a hadron-to-quark matter transition and test the robustness of the compact star mass twin phenomenon against the appearance of pasta phases modelled in this way.

Chiral pasta: Mixed phases at the chiral phase transition

Interiors of neutron stars are ultra-dense and may contain a core of deconfined quark matter. Such a core connects to the outer layers smoothly or through a sharp microscopic interface or through an intermediate macroscopic layer of inhomogeneous mixed phases, which is globally neutral but locally contains electrically charged domains. Here I employ a nucleon-meson model under neutron star conditions that shows a first-order chiral phase transition at large densities. In the vicinity of this chiral transition I calculate the free energies of various mixed phases - different 'pasta structures' - in the Wigner-Seitz approximation. Crucially, chirally broken nuclear matter and the approximately chirally symmetric phase (loosely interpreted as quark matter) are treated on the same footing. This allows me to compute the interface profiles of bubbles, rods, and slabs fully consistently, taking into account electromagnetic screening effects and without needing the surface tension as an input. I find that the full numerical results tend to disfavor mixed phases compared to a simple step-like approximation used frequently in the literature and that the predominantly favored pasta structure consists of slabs with a surface tension $\Sigma \simeq 6\, {\rm MeV}/{\rm fm}^2$.

Crustal torsional oscillations inside the deeper pasta structures

AbstractThe quasiperiodic oscillations (QPOs) have been discovered in the afterglow of the giant flares observed in the soft‐gamma repeaters. As QPOs are strongly associated with the neutron star oscillations, one could extract the neutron star properties via identification of QPOs as specific neutron star oscillations. In this study, we take into account the torsional oscillations inside the phase composed of cylindrical holes (tubes) and spherical holes (bubbles), which correspond to the inner layer in the neutron star crust. We remark that such oscillations are independent of the oscillations inside the phase composed of spherical and cylindrical nuclei, which correspond to the outer layer. Then, we find that all the observed QPOs can be explained with torsional oscillations excited in either layer, which is consistent with the severe constraint on the slope parameter, L, of nuclear symmetry energy as L ≃ 58 − 73 MeV obtained in our previous study.

Domains and defects in nuclear pasta

Nuclear pasta topology is an essential ingredient to determine transport properties in the inner crust of neutron stars. We perform semi-classical molecular dynamics simulations of nuclear pasta for proton fractions $Y_p=0.30$ and $Y_p=0.40$ near one third of nuclear saturation density, $n=0.05\,\mathrm{fm}^{-3}$, at a temperature $T=1.0\,\mathrm{MeV}$. Our simulations are, to our knowledge, the largest nuclear pasta simulations to date and contain up to $3\,276\,800$ nucleons in the $Y_p=0.30$ and $819\,200$ nucleons in the $Y_p=0.40$ case. An algorithm to determine which nucleons are part of a given sub-domain in the system is presented. By comparing runs of different sizes we study finite size effects, equilibration time, the formation of multiple domains and defects in the pasta structures, as well as the structure factor dependence on simulation size. Although we find qualitative agreement between the topological structure and the structure factors of runs with $51\,200$ nucleons and those with $819\,200$ nucleons or more, we show that simulations with hundreds of thousands of nucleons may be necessary to accurately predict pasta transport properties.

Phase transitions and symmetry energy in nuclear pasta

Cold and isospin-symmetric nuclear matter at sub-saturation densities is known to form the so-called pasta structures, which, in turn, are known to undergo peculiar phase transitions. Here we investigate if such pastas and their phase changes survive in isospin asymmetric nuclear matter, and whether the symmetry energy of such pasta configurations is connected to the isospin content, the morphology of the pasta and to the phase transitions. We find that indeed pastas are formed in isospin asymmetric systems with proton to neutron ratios of x=0.3, 0.4 and 0.5, densities in the range of 0.05 fm−3<ρ<0.08 fm−3, and temperatures T<2 MeV. Using tools (such as the caloric curve, Lindemann coefficient, radial distribution function, Kolmogorov statistic, and Euler functional) on the composition of the pasta, determined the existence of homogeneous structures, tunnels, empty regions, cavities and transitions among these regions. The symmetry energy was observed to attain different values in the different phases showing its dependence on the morphology of the nuclear matter structure.

"Do-it-in-classroom" fabrication of microfluidic systems by replica moulding of pasta structures.

Here, we describe a novel method for fabrication of microfluidic structures in classroom environments. This method is based on replica moulding of pasta structures in polydimethylsiloxane. Placing pasta structures on a petroleum jelly base layer enables templating round-shaped structures with controllable cross-sectional profiles. The pasta structures can be easily deformed and combined to create more complex 3D microfluidic structures. Proof-of-concept experiments indicate the capability of this method for studying the mixing of neighbouring flows, generation of droplets, lateral migration of particles, as well as culturing, shear stress stimulation, and imaging of cells. Our "do-it-in-classroom" method bridges the gap between the classroom and the laboratory.

Robustness of third family solutions for hybrid stars against mixed phase effects

We investigate the robustness of third family solutions for hybrid compact stars with a quark matter core that correspond to the occurrence of high-mass twin stars against a softening of the phase transition by means of a construction that mimics the effects of pasta structures in the mixed phase. We consider a class of hybrid equations of state that exploits a relativistic mean-field model for the hadronic as well as for the quark matter phase. We present parametrizations that correspond to branches of high-mass twin star pairs with maximum masses between $2.05~M_\odot$ and $1.48~M_\odot$ having radius differences between 3.2 km and 1.5 km, respectively. When compared to a Maxwell construction with a fixed value of critical pressure $P_c$ the effect of the mixed phase construction consists in the occurrence of a region of pressures around $P_c$ belonging to the coexistence of hadronic and quark matter phases between the onset pressure at $P_H$ and the end of the transition at $P_Q$. The maximum broadening which would still allow mass twin compact stars is found to be $(P_Q-P_H)_{\rm max} \approx P_c$ for all parametrizations within the present class of models. At least the heavier of the neutron stars of the binary merger GW170817 could have been a member of the third family of hybrid stars.

Inhomogeneous Chiral and Coulomb Crystal in Neutron Stars

AbstractWe study the possibility of two types of inhomogeneous phases in core of neutron stars: one is the Coulomb crystal, which is known as quark-hadron pasta structures, and another one is chiral crystal. In the Coulomb crystal, the inhomogeneous phase appears as the result of the balance between the surface tension and the Coulomb interaction. In chiral crystal, we study the inhomogeneous chiral condensate, which has spatial modulation. In the simple model in 1+1 dimensions, this condensate has the same feature with the FFLO state, which is well known in the condensed matter physics.

Rotational stellar structures based on the Lagrangian variational principle

A new method for multi-dimensional stellar structures is proposed in this study. As for stellar evolution calculations, the Heney method is the defacto standard now, but basically assumed to be spherical symmetric. It is one of the difficulties for deformed stellar-evolution calculations to trace the potentially complex movements of each fluid element. On the other hand, our new method is very suitable to follow such movements, since it is based on the Lagrange coordinate. This scheme is also based on the variational principle, which is adopted to the studies for the pasta structures inside of neutron stars. Our scheme could be a major break through for evolution calculations of any types of deformed stars: proto-planets, proto-stars, and proto-neutron stars, etc.

Mixed phase effects on high-mass twin stars

Recently it has been found that a certain class of hybrid star equations of state with a large latent heat (strong first order phase transition obtained by a Maxwell construction) between stiff hadronic hadronic and stiff quark matter phases allows for the appearance of a third family of compact stars (including "twins") at high mass of $\sim 2 M_\odot$. We investigate how robust this high-mass twin phenomenon is against a smoothing of the transition which would occur, e.g., due to pasta structures in the mixed phase. To this end we propose a simple construction of a pasta-like equation of state with a parameter that quantifies the degree of smoothing of the transition and could eventually be related to the surface tension of the pasta structures. It is interesting to note that the range of energy densities for the transition as well as the pressure at the onset of the transition of this class of hybrid star matter at zero temperature corresponds well to values of the same quantities found in finite temperature lattice QCD simulations for the 1 $\sigma$ region at the pseudocritical temperature $T_c=154 \pm 9$ MeV. The pattern of the speed of sound as a function of energy density is very different.

Appearance of the single gyroid network phase in “nuclear pasta” matter

Nuclear matter under the conditions of a supernova explosion unfolds into a rich variety of spatially structured phases, called nuclear pasta. We investigate the role of periodic network-like structures with negatively curved interfaces in nuclear pasta structures, by static and dynamic Hartree-Fock simulations in periodic lattices. As the most prominent result, we identify for the first time the {\it single gyroid} network structure of cubic chiral $I4_123$ symmetry, a well known configuration in nanostructured soft-matter systems, both as a dynamical state and as a cooled static solution. Single gyroid structures form spontaneously in the course of the dynamical simulations. Most of them are isomeric states. The very small energy differences to the ground state indicate its relevance for structures in nuclear pasta.