Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: for the past several decades since the discovery of pulsars , many authors have investigated the properties of dense matter which exist inside neutron stars and supernova cores ( see , e.g. , refs . these objects have been shown that they consist of a variety of material phases whose physical properties reflect in many astrophysical phenomena of these objects . especially , the properties of nuclear matter under extreme conditions , which is one of the essential topics for understanding the mechanism of collapse - driven supernovae @xcite , the structure of neutron star crusts @xcite and its relating phenomena , have been studied actively . this subject is also interesting as one of the fundamental problems of the complex fluids of nucleons .    at subnuclear densities , nuclear matter exhibits the coexistence of a liquid phase with a gas phase due to the internucleon interaction which has an attractive part . at sufficiently low temperatures relevant to neutron star interiors , and sufficiently below the normal nuclear density , long - range coulomb interactions make the system divide periodically into gas and spherical liquid drops , adding a crystalline property to the liquid - gas coexistence .    in the density region where nuclei are about to melt into uniform matter , it is expected that the energetically favorable configuration of the mixed phase possesses interesting spatial structures such as rodlike and slablike nuclei and rodlike and spherical bubbles , etc . , which are referred to as nuclear `` pasta '' . this picture was originally proposed by ravenhall et al . @xcite and hashimoto et al . @xcite independently . their predictions both based on free energy calculations with liquid drop models assuming some specific nuclear shapes . these works clarify that the most energetically stable nuclear shape is determined by a subtle balance between the nuclear surface and coulomb energies . detailed aspects of equilibrium phase diagrams , such as a series of nuclear shapes which can be realized as the energetically most favorable and the density range corresponding to the phases with nonspherical nuclei , vary with nuclear models @xcite . however , the realization of the `` pasta '' phases as energy minimum states can be seen in a wide range of nuclear models and the phase diagrams possess an universal basic feature that , with increasing density , the shape of the nuclear matter region changes like  sphere @xmath4 cylinder @xmath4 slab @xmath4 cylindrical hole @xmath4 spherical hole @xmath4 uniform @xcite . this feature is also reproduced by the thomas - fermi calculations by several groups @xcite . the phases with these exotic nuclear structures , if they were realized in neutron star crusts or supernova cores , bring about many astrophysical consequences . as for those in neutron star phenomena , it is interesting to note the relevance of nonspherical nuclei to pulsar glitches and cooling of neutron stars . although the question whether the mechanism of pulsar glitches is depicted by vortex pinning model or star quake model has yet to be settled completely , the existence of nonspherical nuclei in neutron star matter ( nsm ) have significant effects in both cases . as for the former , while the force needed to pin vortices has yet to be clarified completely even for a bcc lattice of spherical nuclei mainly due to the uncertain properties of impurities and defects @xcite , the effect of spatial structure of normal nuclear matter on vortex dynamics can not be ignored . as for the latter , the existence of `` pasta '' phases with slablike and rodlike nuclei would change the elastic properties of inner crust matter from those of crystalline solid to those of liquid crystal as indicated by pethick and potekhin @xcite , which results in significant decrease of the maximum elastic energy that can be stored in the inner crust . the presence of nonspherical nuclei would also accelerate the cooling of the corresponding region of neutron stars by opening semileptonic weak processes which are unlikely to occur for spherical nuclei @xcite . `` pasta '' phases in supernova matter ( snm ) are expected to affect the neutrino transport and hydrodynamics in supernova cores . let us first note that the neutrino wavelengths , typically of order 20 fm , are comparable to or even greater than the internuclear spacing , leading to diffractive effects on the neutrino elastic scattering off such a periodic spatial structure of nuclear matter @xcite . these effects , induced by the internuclear coulombic correlations , would reduce the scattering rates and hence the lepton fraction @xmath5 . for the bcc lattice of spherical nuclei , such a reduction was examined by horowitz @xcite by calculating the associated static structure factor . it is also noteworthy that nonspherical nuclei and bubbles are elongated in specific direction . in such direction , the neutrino scattering processes are no longer coherent , in contrast to the case of roughly spherical nuclei whose finiteness in any direction yields constructive interference in the scattering . the final point to be mentioned is that the changes in the nuclear shape are accompanied by discontinuities in the adiabatic index , denoting how hard the equation of state of the material is . these discontinuities may influence the core hydrodynamics during the initial phase of the collapse @xcite . though the properties of `` pasta '' phases in equilibrium state have been investigated actively , the formation and the melting processes of these phases have not been discussed except for some limited cases which are based on perturbative approaches @xcite . it is important to adopt a microscopic and dynamical approach which allows arbitrary nuclear structures in order to understand these processes of nonspherical nuclei . at finite temperatures , it is considered that not only nuclear surface becomes obscure but also nuclei of various shapes may coexist . therefore , it is necessary to incorporate density fluctuations without any assumptions on nuclear shape to investigate the properties of `` pasta '' phases at finite temperatures . although the works done by williams and koonin @xcite and lassaut et al . @xcite do not assume nuclear structure , they can not incorporate fluctuations of nucleon distributions in a satisfying level because they are based on the thomas - fermi calculation , which is one - body approximation . in addition , only a single structure is contained in the simulation box in these works , there are thus possibilities that nuclear shape is strongly affected by boundary effect and some structures are prohibited implicitly .    in the present work , we study the structure of cold dense matter at subnuclear densities in the framework of quantum molecular dynamics ( qmd ) @xcite , which is one of the molecular dynamics ( md ) approaches for nucleon many - body systems ( see , e.g. , ref . @xcite for review ) . md for nucleons including qmd , which is a microscopic and dynamical method without any assumptions on nuclear structure , is suitable for incorporating fluctuations of particle distributions . previously , we have reported the first results of our study on nuclear `` pasta '' by qmd , which demonstrated that the `` pasta '' phases can be formed in a dynamical way for matter with proton fractions @xmath6 and 0.5 @xcite . in this full paper , we present new results for astrophysically interesting neutron - rich matter of @xmath7 in addition to the cases of @xmath6 and 0.5 reported before . the plan of this paper is as follows . in section ii , we describe the framework of the qmd model used in the present study . we then show the results of our simulations in section iii and analyze the structure of matter obtained by the simulations using two - point correlation functions and minkowski functionals in section iv . in section v , we investigate the properties of the effective nuclear interaction used in this work in order to examine the validity of our results in terms of nuclear forces . astrophysical discussions are given in section vi . summary and conclusions are presented in section vii . we have various types of molecular dynamics methods for nucleons including representative ones such as fermionic molecular dynamics ( fmd ) @xcite , antisymmetrized molecular dynamics ( amd ) @xcite and qmd , etc . in the present work , we choose qmd among them balancing between calculation cost and accuracy . the typical length scale @xmath8 of inter - structure is @xmath9 fm and the density region of interest is just below the normal nuclear density @xmath10 . the required nucleon number @xmath11 in order to reproduce @xmath12 unit structures in the simulation box is about @xmath13 ( for slabs ) . it is thus desirable that we prepare nucleons of order 10000 if we try to reduce boundary effects down to a satisfactory level by reproducing several several unit structures in the box . while it is a hard task to treat such a large system with , for example , fmd and amd whose calculation costs scale as @xmath14 , it is feasible to do it with qmd whose calculation costs scale as @xmath15 . this difference comes from summations in the slater determinants in the trial wave functions of the former models . in qmd , on the other hand , the total @xmath11-nucleon wave function @xmath16 is assumed to be a direct product of single - nucleon states @xmath17 : @xmath18 the single - nucleon state is represented by a gaussian wave packet : @xmath19}\ , \label{packet}\ ] ] where @xmath20 and @xmath21 are the centers of position and momentum of the packet @xmath22 , respectively , and @xmath23 is a parameter related to the extension of the wave packet in the coordinate space . it is also noted that we mainly focus on the macroscopic structures ; the exchange effect would not be so important for them . this can be seen by comparing the typical values of the exchange energy for the macroscopic scale and of the energy difference between two successive phases with nonspherical nuclei . suppose there are two identical nucleons , @xmath24 and 2 , bound in different nuclei each other . the exchange energy between these particles is calculated as an exchange integral : @xmath25 where @xmath26 is the potential energy . an asymptotic form of the wave function is given by @xmath27 with @xmath28 and @xmath29 , where @xmath30 is the binding energy and @xmath31 is the nucleon mass . the exchange integral reads @xmath32 } \sim 5 \times 10^{-6 } { \rm mev}\ ] ] for the internuclear distance @xmath33 fm and @xmath34 mev , which is extremely smaller than the typical energy difference per nucleon between the `` pasta '' phases of order 0.1 kev ( for nsm , see fig .  4 in ref . @xcite ) - 10 kev ( for snm , see fig .  4 in ref . therefore , it is expected that qmd , which is less elaborate in treating the exchange effect , is not bad approximation for investigating the nuclear `` pasta '' . consequently , qmd has the advantages over the other models in the present study . in the future , we will have to confirm the validity of the results obtained by qmd using other more elaborate model such as amd or fmd to treat the exchange effect more precicely . however , this problem is beyond the scope of the present work .      to simulate nuclear matter at subnuclear densities within the framework of qmd , we use a qmd model hamiltonian developed by maruyama et al . @xcite , which is constructed so as to reproduce bulk properties of nuclear matter and properties of finite nuclei . this model hamiltonian consists of the following six terms @xmath35 where @xmath36 is the kinetic energy , @xmath37 is the pauli potential introduced to reproduce the pauli principle effectively , @xmath38 is the skyrme potential which consists of an attractive two - body term and a repulsive three - body term , @xmath39 is the symmetry potential , @xmath40 is the momentum - dependent potential introduced as two fock terms of the yukawa interaction and @xmath41 is the coulomb potential . the expressions of these terms are given as    @xmath42 } \ \delta_{\tau_i \tau_j } \delta_{\sigma_i \sigma_j}\ , \label{pauli}\\    v_{\rm skyrme } & = &    { \alpha\over 2\rho_0}\sum_{i , j ( \neq i ) }    \rho_{ij }    +   { \beta\over ( 1+\tau)\ \rho_0^{\tau } }    \sum_i \left [ \sum_{j ( \neq i ) } \int { d^3{\bf r } \ \tilde{\rho_i}({\bf r } ) \                         \tilde{\rho_j}({\bf r } ) } \right]^{\tau}\ ,                     \label{skyrme}\\    v_{\rm sym } & = &    { c_{\rm s}\over 2\rho_0 } \sum_{i , j(\neq i ) } \ ,    ( 1 - 2 | c_i - c_j | ) \ \label{sym}\\    v_{\rm md }   & = &    v_{\rm md}^{(1 ) } + v_{\rm md}^{(2 ) } \nonumber \\    & = &           { c_{\rm ex}^{(1 ) } \over 2\rho_0 } \sum_{i , j(\neq i ) }         { 1 \over 1+\left[{{\bf p}_i-{\bf p}_j \over \hbar \mu_1}\right]^2 }         \ \rho_{ij }       +    { c_{\rm ex}^{(2 ) } \over 2\rho_0 } \sum_{i , j(\neq i ) }         { 1 \over 1+\left[{{\bf p}_i-{\bf p}_j \over \hbar \mu_2}\right]^2 }         \ \rho_{ij}\ , \label{md}\\    v_{\rm coulomb } & = &    { e^2 \over 2}\sum_{i , j(\neq i ) }    \left(\tau_{i}+\frac{1}{2}\right ) \ , \left(\tau_{j}+\frac{1}{2}\right )    \int\!\!\!\!\int d^3{\bf r}\,d^3{\bf r}^{\prime }     { 1 \over|{\bf r}-{\bf r}^{\prime}| } \ ,    \rho_i({\bf r})\rho_j({\bf r}^{\prime})\ , \label{coulomb}\end{aligned}\ ] ]    where @xmath43 means the overlap between the single - nucleon densities , @xmath44 and @xmath45 , for @xmath22-th and @xmath46-th nucleons given as @xmath47 @xmath48 is the nucleon spin and @xmath49 is the isospin ( @xmath50 for protons and @xmath51 for neutrons ) and @xmath52 and @xmath23 are model parameters determined to reproduce the properties of the ground states of the finite nuclei , especially heavier ones , and the saturation properties of nuclear matter @xcite . a parameter set used in this work is shown in table [ parameter ] . the single - nucleon densities @xmath44 and @xmath53 are given by @xmath54}\ , \quad \\ \tilde{\rho_i}({\bf r } ) & = &    \frac{1}{(2\pi \tilde{l})^{3/2}}\ \exp{\left [                  - \frac{({\bf r } - { \bf r}_i)^2}{2\tilde{l } } \right]}\ , \end{aligned}\ ] ] with @xmath55 the modified width @xmath56 in @xmath53 is introduced in the three - body term of skyrme interaction [ eq . ( [ skyrme ] ) ] to incorporate the effect of the repulsive density - dependent term . .effective interaction parameter set + ( incompressibility @xmath57=280 mev ; medium eos model in ref . @xcite ) [ cols="^,^,^,^ " , ]     [ parameter ]    we adopt qmd equations of motion with friction terms to simulate the dynamical relaxation : @xmath58 where the friction coefficients @xmath59 and @xmath60 are positive definite , which determine the relaxation time scale . the relaxation scheme given by eqs . ( [ qmdeom fric ] ) is referred to as the steepest descent method and it leads to the continuous decrease in @xmath61 as @xmath62 even though it is recognized that this method is not efficient , it is expected that the dynamics given by eqs . ( [ qmdeom fric ] ) with @xmath63 deviates slightly in a short period from the physically grounded dynamics given by qmd equations of motion without the friction terms [ equations without the second terms in the right - hand sides of eqs . ( [ qmdeom fric ] ) ] , which we would like to respect . we have performed qmd simulations of an infinite @xmath64 system with fixed proton fractions @xmath65 0.5 and 0.3 for various nucleon densities @xmath66 [ the density region is ( 0.05 - 1.0 ) @xmath67 . we set 2048 nucleons ( 1372 nucleons in some cases ) contained in a cubic box which is imposed by the periodic boundary condition . throughout this paper , the numbers of the protons ( neutrons ) with up - spin and with down - spin are equal . the relativistic degenerate electrons which ensure the charge neutrality are regarded as a uniform background and the coulomb interaction is calculated by the ewald method taking account of the gaussian charge distribution of each wave - packet ( see appendix [ ewald sum ] ) . this method enables us to efficiently sum up contributions of long - range interactions in a system with periodic boundary conditions . for nuclear interaction , we use the effective hamiltonian developed by maruyama et al . ( medium eos model ) @xcite whose expressions are given in the last section . we first prepare an uniform hot nucleon gas at @xmath68 mev as an initial condition equilibrated for @xmath69 fm/@xmath70 in advance . in order to realize the ground state of matter , we then cool it down slowly for @xmath71 fm/@xmath70 , keeping the nucleon density constant with the frictional relaxation method [ eqs . ( [ qmdeom fric ] ) ] , etc . @xcite until the temperature gets @xmath72 mev or less . note that no artificial fluctuations are given in the simulation . the qmd equations of motion with the friction terms given by eqs . ( [ qmdeom fric ] ) are solved using the fourth - order gear predictor - corrector method in conjunction with multiple time step algorithm @xcite . integration time steps @xmath73 are set to be adaptive in the range of @xmath74 fm/@xmath70 depending on the degree of convergence . at each step , the correcting operation is iterated until the error of position @xmath75 and the relative error of momentum @xmath76 become smaller than @xmath77 , where @xmath75 and @xmath76 are estimated as the maximum values of correction among all particles . we mainly use pcs(pentium iii ) equipped with mdgrape-2 , which accelerates calculations of momentum - independent forces including the long - range coulomb force .    shown in figs . [ fig pasta sym ] and [ fig pasta x0.3 ] are the resultant nucleon distributions of cold matter at @xmath78 0.5 and 0.3 , respectively . we can see from these figures that the phases with rodlike and slablike nuclei , cylindrical and spherical bubbles , in addition to the phase with spherical nuclei are reproduced in both the cases of @xmath78 0.5 and 0.3 . we here would like to mention the reasons of discrepancies between the present result and the result obtained by maruyama et al . which says `` the nuclear shape may not have these simple symmetries '' @xcite . one of the most crucial reasons seems to be the difference in treatment of the coulomb interaction . in the present simulation , we calculate the long range coulomb interaction in a consistent way using the ewald method . for the system of interest where the thomas - fermi screening length is comparable to or larger than the size of nuclei , this treatment is more adequate than that which inctroduces an artificial cutoff distance as in ref . the other crucial reason is the difference in the relaxation time scales @xmath79 fm/@xmath70 ; we set @xmath80 in the present work , but maruyama et al . set @xmath81 several@xmath82 fm/@xmath70 @xcite . in our simulation , we can reproduce the bubble - phases [ see ( d ) and ( e ) of figs . [ fig pasta sym ] and [ fig pasta x0.3 ] ] with @xmath83 fm/@xmath70 and the nucleus - phases [ see ( b ) and ( c ) of figs . [ fig pasta sym ] and [ fig pasta x0.3 ] ] with @xmath84 fm/@xmath70 . however , the matter in the density region corresponding to a nucleus - phase is quenched in an amorphous - like state when @xmath85 fm/@xmath70 . in the present work , we take @xmath79 much larger than typical time scale @xmath86 fm/@xmath70 for nucleons to thermally diffuse in the distance of @xmath9 fm at @xmath87 and @xmath88 mev . this temperature is lower than the typical value of the liquid - gas phase transition temperature in the density region of interest , it is thus considered that our results are thermally relaxed in a satisfying level . phase diagrams of matter in the ground state are shown in figs . [ phase diagram 1](a ) and ( b ) for @xmath78 0.5 and 0.3 , respectively . as can be seen from these figures , the obtained phase diagrams basically reproduce the sequence of the energetically favored nuclear shapes predicted by simple discussions @xcite which only take account of the coulomb and surface effects ; this prediction is that the nuclear shape changes like  sphere @xmath4 cylinder @xmath4 slab @xmath4 cylindrical hole @xmath4 spherical hole @xmath4 uniform , with increasing density . comparing figs .  3 ( a ) and ( b ) , we can see that the phase diagram shifts towards the lower density side with decreasing @xmath89 , which is due to the tendency that the saturation density is lowered as the neutron excess increases . it is remarkable that the density dependence of the nuclear shape , except for spherical nuclei and bubbles , is quite sensitive and phases with intermediate nuclear shapes which are not simple as shown in figs . 1 and 2 are observed in two density regions : one is between the cylinder phase and the slab phase , the other is between the slab phase and the cylindrical hole phase . we note that these phases are different from coexistence phases with nuclei of simple shapes , which will be referred to as `` intermediate phases '' . we have also performed qmd simulations of matter with proton fraction @xmath7 as a more realistic condition for the neutron star matter . in this case , we have to deal with a larger system than in the cases of @xmath0 and 0.3 because enough number of protons for reproducing several nuclei in a simulation box are required to obtain significant results ; protons play an important role in generating the long - range order due to their electric charge . we have investigated the neutron - rich matter at @xmath7 with 10976 nucleons , in which 1098 protons and 9878 neutrons are contained . following basically the same procedure that used for the cases of @xmath0 and 0.3 ( see section [ sim p - rich ] for detail ) , we tried to obtain the ground - state matter . however , in the present case , we quickly relax from the initial state at @xmath90 mev to the state at @xmath91 mev , at which matter is still uniform , with a nos - hoover - like thermostat @xcite which will be discussed in another paper @xcite . after the relaxation at @xmath91 mev for @xmath92  7000 fm/@xmath70 , we then cool down the system with a relaxation time scale @xmath84 fm/@xmath70 using the qmd equations of motion with friction terms ( [ qmdeom fric ] ) . these simulations are performed by fujitsu vpp 5000 equipped with mdgrape-2 . some resultant nucleon distributions are shown in figs . [ fig sphere x0.1 ] and [ fig cylinder x0.1 ] , which correspond to the sphere phase and the cylinder phase , respectively . as can be seen in fig . [ fig sphere x0.1 ] , dripped neutrons spread over the whole region in the simulation box , which lead to smaller density contrast compared with that for the cases of @xmath0 and 0.3 depicted in figs . [ fig pasta sym ] and [ fig pasta x0.3 ] , respectively . the results obtained for @xmath7 are summarized in the phase diagram shown in fig . [ phase diagram 2 ] . a striking feature is that the wide density region from @xmath93 to @xmath94 is occupied by an intermediate phase . the structure of matter seems to change rather continuously from that consists of branching rodlike nuclei connected to each other [ obtained in the lower density region of the intermediate phase denoted by ( c , u ) ] to that consists of branching bubbles connected to each other [ higher density region of the intermediate phase ( c , u ) ] . however , in the present neutron - rich case , `` pasta '' phase with slablike nuclei can not be obtained as far as we have investigated , which will be discussed at the end of the next section . it is also noted that the density at which matter turns into uniform is lower than those in the cases of @xmath0 and 0.3 , which is consistent with the tendency that matter becomes the more neutron - rich , the saturation density becomes the lower . to analyze the spatial distribution of nucleons , we calculate two - point correlation function @xmath95 for nucleon density field @xmath96 @xmath97 where @xmath11 stands for nucleons@xmath98 . @xmath95 is here defined as @xmath99 where @xmath100 denotes an average over the position @xmath101 and the direction of @xmath102 , and @xmath103 is the fluctuation of the density field @xmath104 given by @xmath105 with @xmath106    we construct the nucleon density distribution @xmath104 from a data set of the centers of position of the nucleons by the following procedure . we first set @xmath107 ( for @xmath0 and 0.3 ) or @xmath108 ( for @xmath7 ) grid points in the simulation box and then distribute particle numbers on each grid point using the cloud - in - cell method ( see , e.g. , ref . next , we carry out the smoothing procedure in the discrete fourier space with a gaussian smoothing function corresponding to the distribution of the wave packet given by eq . ( [ packet ] ) . the density distributions @xmath109 , where @xmath110 or 1 projects on particle type @xmath22 , in the discrete real space can be obtained by the inverse fourier transformation . the fourier transformations are performed using the fft algorithm .    the resultant two - point correlation functions @xmath111 , @xmath112 and @xmath113 at various densities below @xmath114 at which matter becomes uniform at zero temperature for @xmath0 , 0.3 and 0.1 are plotted in figs . [ corr x0.5 ] , [ corr x0.3 ] and [ corr x0.1 ] , respectively . we can see the general tendency , which is common for the different values of @xmath89 , that the amplitude of @xmath115 decreases with increasing the density . it is noted that even though the change in the amplitude of @xmath115 is quite noticeable , the smallest zero - point @xmath116 of @xmath115 takes similar values at various densities especially for @xmath0 and 0.3 . this feature means that the typical length scales of the nuclear structures , i.e. , the internuclear distance and the nuclear radius , keep comparable at subnuclear densities from @xmath117 to @xmath114 , which is consistent with the results obtained by the previous works ( see , e.g. , refs . this behavior just below @xmath114 will be discussed further concerning a problem about the properties of the transition to uniform matter . we also note that a strong attractive force acting between a proton and a neutron leads to the good agreement of the zero - points of @xmath118 and @xmath119 even for @xmath6 and 0.1 as well as for @xmath0 although the zero - point @xmath120 of @xmath119 is @xmath121 fm ( @xmath122 fm ) larger than that of @xmath118 for @xmath6 ( 0.1 ) at each density . this shows that the phase of the density fluctuations of protons and neutrons correlate so strongly with each other at zero temperature that they almost coincide . as can be seen by comparing @xmath111 for different values of the proton fraction @xmath89 , the amplitude decreases and the value of the smallest zero - point @xmath120 increases with decreasing @xmath89 . this behavior means that , as matter becomes more neutron - rich , not only the nucleon density distribution gets smoother but also the spatial structure becomes larger . let us then examine @xmath95 for each value of the proton fraction . for symmetric matter ( @xmath0 ) , protons and neutrons are equivalent except for the mass difference and the coulomb interaction . therefore , the two - point correlation functions @xmath123 , @xmath118 and @xmath119 are almost the same at larger values of @xmath124 . at smaller values of @xmath125 , @xmath118 is slightly smaller than @xmath119 because the repulsive coulomb interaction among protons tends to reduce the proton density inhomogeneity especially in the smaller scale . for asymmetric matter ( @xmath6 and 0.1 ) , on the contrary , the amplitude of @xmath119 is much smaller than that of @xmath118 due to the dripped neutrons which distribute rather uniformly outside the nuclei . let us here examine the properties of the transition from the phase with spherical bubbles to uniform matter for @xmath0 . for this purpose , two - point correlation function of the nucleon density fluctuation is useful . in fig . [ corr2 melt ] , we thus plot the two - point correlation function of the nucleon density fluctuation @xmath111 for several densities around the melting density @xmath114 . to compute @xmath111 , we use a 1372-nucleon system cooled down until the temperature gets @xmath126 0.05 mev by qmd equations of motion ( [ qmdeom fric ] ) with friction terms .    for uniform phase , @xmath111 should be zero except for the contribution of short - range correlation . the behavior of @xmath111 shows that @xmath114 lies between 0.7@xmath127 and @xmath128 at @xmath0 [ see also fig . [ phase diagram 1](a ) ] above which long - range correlation disappears . it is noted that the smallest value of @xmath116 at which @xmath129 keeps around 8 fm even at densities just below @xmath114 . this means that , for @xmath0 , the phase with spherical bubbles whose radii are around @xmath120 suddenly disappear rather than they shrink gradually and the system turns into uniform with increasing the density because the quantity @xmath120 nearly corresponds to the half wave - length of the inhomogeneous density profile . the discontinuous change in the density profile indicates that the transition between the phase with spherical bubbles and the uniform phase is of first order . this conclusion is also obtained in the previous calculations for which the spatial structure of the nuclear matter region and/or the shape of the density profile are assumed @xcite . in the present work , we have confirmed the first order nature by qmd simulations without these assumptions as several authors have done so without the assumptions by the thomas - fermi approximation @xcite .    for the cases of @xmath6 and 0.1 , we could not see the significant sign of the first order nature of the transition between the mixed phase and the uniform phase because the amplitude of @xmath111 is quite small just below @xmath114 . further study is necessary to determine the properties of the transition for theses cases of asymmetric matter .      to extract the morphological characteristics of the nuclear shape changes and the intermediate phases , we introduce the minkowski functionals ( see , e.g. , ref . @xcite and references therein ; a concise review is provided by ref . @xcite ) as geometrical and topological measures of the nuclear surface . let us consider a homogeneous body @xmath130 in the @xmath131-dimensional eucledian space , where @xmath132 is the class of such bodies . morphological measures are defined as functionals @xmath133 which satisfy the following three general properties :    * _ motion invariance _ : the functional is independent of the position and the direction of the body , i.e. , @xmath134 where @xmath135 denotes any translations and rotations . * _ additivity _ : the functional of the union of two bodies should behave like a volume . the contribution of the overlapping region should be subtracted , i.e. , @xmath136 where @xmath137 . * _ continuity _ : if the body is approximated with pixels , the functional of the approximate body converges to that of the original body when the pixels get smaller , i.e. , @xmath138 where @xmath57 is a convex body and @xmath139 is a sequence of convex bodies . hadwiger s theorem in integral geometry states that there are just @xmath140 independent functionals which satisfy the above properties ; they are known as minkowski functionals . in three dimensional space , four minkowski functionals are related to the volume , the surface area , the integral mean curvature and the euler characteristic .    in classifying nuclear shapes including those of the intermediate phases obtained in our simulations , the integral mean curvature and the euler characteristic are useful , which will be discussed later . both are described by surface integrals of the following local quantities : the mean curvature @xmath141 and the gaussian curvature @xmath142 , i.e. , @xmath143 and @xmath144 , where @xmath145 and @xmath146 are the principal curvatures and @xmath147 is the area element of the surface of the body @xmath57 . the euler characteristic @xmath148 is a purely topological quantity and is given by @xmath149 thus @xmath150 for the sphere and the spherical hole phases and the coexistence phase of spheres and cylinders , and @xmath151 for the other ideal `` pasta '' phases , i.e. the cylinder , the slab and the cylindrical hole phases which consist of infinitely long rods , infinitely extending slabs and infinitely long cylindrical holes , respectively . we introduce the area - averaged mean curvature , @xmath152 , and the euler characteristic density , @xmath153 , as normalized quantities , where @xmath154 is the volume of the whole space . we calculate the normalized minkowski functionals , i.e. , the volume fraction @xmath155 , the surface area density @xmath156 , the area - averaged mean curvature @xmath157 and the euler characteristic density @xmath158 for @xmath0 and 0.3 by the following procedure . as described in subsection [ 2pcf ] , we first construct proton and nucleon density distributions @xmath159 and @xmath160 , where @xmath161 or 1 is the isospin projection on the proton state . we set a threshold proton density @xmath162 and then calculate @xmath163 , where @xmath164 and @xmath165 are the volume and the surface area of the regions in which @xmath166 . we find out the value @xmath167 where @xmath168 and define the regions in which @xmath169 as nuclear regions . for spherical nuclei , for example , @xmath170 corresponds to a point of inflection of a radial density distribution . in the most phase - separating region , the values of @xmath171 distribute in the range of about @xmath172 @xmath173 in the both cases of @xmath174 and 0.3 , where @xmath171 is the threshold nucleon density corresponds to @xmath170 . we then calculate @xmath155 , @xmath175 , @xmath176 and @xmath148 for the identified nuclear surface . we evaluate @xmath155 by counting the number of pixels at which @xmath177 is higher than @xmath170 , @xmath175 by the triangle decomposition method , @xmath176 by the algorithm shown in ref . @xcite in conjunction with a calibration by correction of surface area , and @xmath148 by the algorithm of ref . @xcite and by that of counting deficit angles @xcite , which are confirmed that both of them give the same results . we have plotted the resultant @xmath66 dependence of @xmath155 , @xmath156 , @xmath157 and @xmath153 for the isodensity surface of @xmath178 in figs . [ fig minko 1 ] and [ fig minko 2 ] . in addition to the values of @xmath155 , @xmath156 and @xmath157 for the isodensity surface of @xmath178 , we have also investigated those for the isodensity surfaces of @xmath179 to examine the extent of the uncertainties of these quantities which stem from the arbitrariness in the definition of the nuclear surface . as shown in fig . [ fig minko 1 ] , these uncertainties are at most @xmath180 for @xmath155 and @xmath181 @xmath182 for @xmath156 . for @xmath157 , we could not observe remarkable differences from the values for @xmath178 ( they were smaller than 0.015 @xmath182 ) . we could not see these kinds of uncertainties in @xmath153 , except for the densities near below @xmath114 .    as shown in fig . [ fig minko 1 ] , the volume fraction @xmath155 of the nuclear regions increases almost monotonically under @xmath114 in both cases of @xmath0 and 0.3 . this feature reflects the incompressible nature of nuclear matter . it is interesting to see the density dependence of the nuclear surface density @xmath156 because this quantity is directly related to the surface energy density , which is one of the key factors in determining the nuclear shape . figs . [ fig minko 1](c ) and [ fig minko 1](d ) show that , as the nucleon density increases , @xmath156 increases at a nearly constant rate until @xmath183 , and then its increasing rate becomes rather smaller around the density region of the slab phase , and finally it begins to decrease in the density region of the cylindrical hole phase or the spherical hole phase . this general behavior can be understood from the density dependence of the surface energy density obtained by simple arguments , which only allow for the nuclear surface and the coulomb effects ( see , e.g. , refs . @xcite ) . we also plot @xmath155 and @xmath156 for the nucleon density distribution @xmath184 as functions of the threshold density @xmath185 evaluated at various values of @xmath66 . the results for @xmath0 and 0.3 are shown in figs . [ fig minko thres x0.5 ] and [ fig minko thres x0.3 ] , respectively . features of nuclear shape changes can be seen in the behavior of the curves of @xmath156 . peaks in the higher @xmath185 region are attributed to nucleons in the nuclear matter regions and broad bumps in the lower @xmath185 region around 0.05  0.1@xmath127 observed for @xmath6 are due to the dripped neutrons outside nuclei . @xmath156 in the intermediate @xmath185 region in which its slope is nearly constant mainly comes from contribution of nuclear surfaces . as the nucleon density @xmath66 increases , the higher @xmath185 peak becomes more clear and the position of the center of the peak finally coincides with @xmath66 in the uniform phase . this feature shows that the nuclear matter regions become more uniform with increasing the density . when the dispersion of the internucleon distance is small , large surface area caused by the gaussian density distribution of each nucleon can be picked up with a single value of @xmath185 . as can be seen in fig . [ fig minko thres x0.3 ] , the lower @xmath185 bump , in turn , disappears with increasing @xmath66 . this is because , as @xmath66 increases , the dripped neutron gas becomes more inhomogeneous and tends to distribute close to the nuclear surface leading to a lower proton fraction in the nuclear matter regions . it is also noted that , as the nucleon density increases , the slope in the intermediate @xmath185 region changes from negative to positive at the density corresponding to the phase of slablike nuclei ( 0.4@xmath127 for @xmath0 and 0.35@xmath127 for @xmath6 ) , which is consistent with what is expected from the sign of @xmath157 for the nuclear surface [ see figs .  [ fig minko 2](a ) and [ fig minko 2](b ) ] . let us then focus on @xmath157 and @xmath158 to classify the nuclear shape . the behavior of @xmath157 shows that it decreases almost monotonically from positive to negative with increasing @xmath66 until the matter turns into uniform . the densities corresponding to @xmath186 are about 0.4 and 0.35@xmath127 for @xmath0 and 0.3 , respectively ; these values are consistent with the density regions of the phase with slablike nuclei ( see fig . [ phase diagram 1 ] ) . as mentioned previously , @xmath153 is actually positive in the density regions corresponding to the phases with spherical nuclei , coexistence of spherical and cylindrical nuclei , and spherical holes because of the existence of isolated regions . as for those corresponding to the phases with cylindrical nuclei , planar nuclei and cylindrical holes , @xmath187 . the fact that the values of @xmath153 are not exactly zero for nucleon distributions shown as the slab phase in figs . [ fig pasta sym ] and [ fig pasta x0.3 ] reflects the imperfection of these `` slabs '' , which is due to the small nuclear parts connecting the neighboring slabs . however , we can say that the behaviors of @xmath153 plotted in figs . [ fig minko 2](c ) and [ fig minko 2](d ) show that @xmath153 is negative in the density region of the intermediate phases , even if we take into account the imperfection of the obtained nuclear shapes and the uncertainties of the definition of the nuclear surface . this means that the intermediate phases consist of nuclear surfaces which are saddle - like at each point on average and they consist of each highly connected nuclear and gas regions due to a lot of tunnels [ see eq . ( [ euler ] ) ] . using the quantities @xmath157 and @xmath158 , the sequence of the nuclear shapes with increasing the density can be described as follows : @xmath188 @xmath4 @xmath189 @xmath4 @xmath190 @xmath4 @xmath191 @xmath4 @xmath192 @xmath4 @xmath193 @xmath4 @xmath194 @xmath4 uniform . let us now consider the discrepancy from the results of previous works which do not assume nuclear structure @xcite ; the intermediate phases can not be seen in these works . we can give following two reasons for the discrepancy . * these previous calculations are based on the thomas - fermi approximation which can not sufficiently incorporate fluctuations of nucleon distributions . this shortcoming may result in favoring nuclei of smoothed simple shapes than in the real situation . * there is a strong possibility that some highly connected structures which have two or more substructures in a period are neglected in these works because only one structure is contained in a simulation box . it is not unnatural that the phases with highly connected nuclear and bubble regions are realized as the most energetically stable state @xcite . it is considered that , for example , a phase with perforated slablike nuclei , which has negative @xmath153 , could be more energetically stable than that with extremely thin slablike nuclei . the thin planar nucleus costs surface - surface energy which stems from the fact that nucleons bound in the nucleus feel its surfaces of both sides . the surface - surface energy brings about an extra energy increase in addition to the contribution of the surface energy . we have to examine the existence of the intermediate phases by more extensive simulations with larger nucleon numbers and with longer relaxation time scales in the future .      in the case of @xmath7 , the criterion for identification of the isodensity surface corresponding to the nuclear surface using the second derivative of @xmath195 does not work at higher densities . we thus use another method to calculate the normalized minkowski functionals of the nuclear surface for @xmath7 .    in fig . [ fig euler typical ] , we have plotted the @xmath185 dependence of the euler characteristic density @xmath158 at @xmath196 as an example . we can see that this curve consists of three components : the peaks of the lower @xmath185 region , the plateau region and the peaks of the higher @xmath185 region , which are due to dripped neutrons ( thus these peaks can not be observed for @xmath0 ) , nuclear surfaces and nucleons in nuclei , respectively . these components can also be seen at the other values of @xmath66 lower than @xmath114 . however , we have to mention that the higher the density becomes , the smaller the plateau region gets , which means that the density contrast between the dripped neutron gas region and the nuclear matter region becomes obscure . here , we take the mean values of the normalized minkowski functionals in the plateau region as those for the nuclear surface , which are plotted as crosses in fig . [ fig minko x0.1 ] . the error bars shown in this figure are the standard deviations of these quantities in the plateau region . consistency between this method and the one using the second derivative of @xmath195 has been confirmed for @xmath6 . [ fig minko x0.1 ] shows the resultant normalized minkowski functionals for the nuclear surface at various values of @xmath66 below @xmath114 . the qualitative behaviors of @xmath155 and @xmath157 for @xmath7 are the same as those for @xmath0 and 0.3 ; as @xmath66 increases , @xmath155 increases and @xmath157 decreases ( from positive to negative ) almost monotonically in the density region of 0.1@xmath197 . however , in the behaviors of @xmath156 and @xmath158 , qualitative differences can be observed between the present case and the cases of @xmath0 and 0.3 . as can be seen in fig . [ fig minko x0.1](b ) , @xmath156 increases almost linearly until just below @xmath114 .    the absence of the phases with cylindrical bubbles and with spherical bubbles in the phase diagram of @xmath7 ( fig . [ phase diagram 2 ] ) is well characterized by the behavior of @xmath158 shown in fig . [ fig minko x0.1](d ) . in the cases of @xmath0 and 0.3 , @xmath158 increases from negative to positive with increasing density in the density region higher than that of the slab phase . however , for @xmath7 , we can not observe the tendency that @xmath158 starts increasing even at just below @xmath114 . as a result , @xmath158 keeps negative until matter becomes uniform . ( further discussion will be given in section [ int ] . ) let us then consider the phase with slablike nuclei , which have not been obtained in the simulations for @xmath7 . if it were realized by using a longer relaxation time scale , it is expected to be obtained at @xmath198  0.34 @xmath127 according to the behaviors of the @xmath157 and @xmath158 . however , we can not see any signs from fig . [ fig minko x0.1](b ) that @xmath156 stops increasing in this density region unlike the behaviors of @xmath156 in the cases of @xmath0 and 0.3 . here , we would like to mention that , according to the landau - peierls argument , thermal fluctuations are effective at destroying the long - range order of one - dimensional layered lattice of slablike nuclei rather than that of triangular lattice of rodlike nuclei and of bcc lattice of spherical ones . thus , the melting temperature of the planar phase would be lower than the other phases , which leads to a longer time scale for formation of the slablike nuclei by the thermal diffusion . therefore , a further investigation with a longer relaxation time scale is necessary to determine whether or not the phase with slablilke nuclei is really prohibited in such neutron - rich matter in the present model .    in fig . [ fig minko thres x0.1 ] , we have also plotted @xmath155 and @xmath156 for the density distribution @xmath184 as functions of @xmath185 like figs . [ fig minko thres x0.5 ] and [ fig minko thres x0.3 ] . in comparison with fig . [ fig minko thres x0.3 ] , the contribution of the dripped neutrons is shown more clearly in this case . we can see that the peak in the lower @xmath185 region due to the dripped neutrons combines to the peak in the higher @xmath185 region . this behavior stems from the fact that a part of the dripped neutrons at lower densities are absorbed into nuclear matter region with increasing the density at fixed @xmath89 ; finally , all the neutrons are contained there in the uniform phase . we can also expect from the @xmath185 dependence of @xmath156 that the phase with slablike nuclei might be obtained around @xmath199 , where the slope of the plateau region is close to zero . let us here examine the effective nuclear interaction used in this work . structure of matter at subnuclear densities is affected by the properties of neutron - rich nuclei and of the pure neutron gas resulted from the nuclear interaction . key quantities are the energy per nucleon @xmath1 of the pure neutron matter , the proton chemical potential @xmath2 in the pure neutron matter , and the nuclear surface tension @xmath3 . there is a tendency , especially in the case of neutron star matter , that the higher @xmath1 , the density @xmath114 at which matter becomes uniform is lowered . this is because larger @xmath1 tends to favor uniform nuclear matter without dripped neutron gas regions than mixed phases with dripped neutron gas regions . in the neutron star matter , there is also a tendency that the lower @xmath2 , the smaller @xmath114 . this is because @xmath200 represents the degree to which the neutron gas outside the nuclei favors the presence of protons in itself . the quantity @xmath3 controls the size of the nuclei and bubbles , and hence the sum of the coulomb and surface energies . with increasing @xmath3 and so this energy sum , @xmath114 gets lowered . it is important to check whether the effective nuclear force given by eqs . ( [ hamiltonian])([coulomb ] ) yields unrealistic values of these quantities or not . if @xmath1 , @xmath201 and @xmath3 for the present model are unrealistically small in comparison with those for the other models , our results which have reproduced the `` pasta '' phases might be quite limited for the present model hamiltonian .    in order to evaluate @xmath1 , we set 1372 neutrons in a simulation box imposed of the periodic boundary condition . this system is cooled down by the qmd equations of motion with friction terms [ see eqs . ( [ qmdeom fric ] ) ] until the temperature becomes @xmath126 1 kev . the resultant values of @xmath1 are plotted in fig . [ fig e_n ] . we note that our results for @xmath2020.2 , 0.6 and 1.0 @xmath127 ( the result for @xmath203 is not plotted in fig . [ fig e_n ] ) coincide with the results for zero - proton ratio plotted in fig .  9 of ref . @xcite . the values of @xmath1 for the present model behave like those for the skm skyrme interaction especially in the density region of @xmath204 @xmath173 ; they are close to the result of the variational chain summation obtained by akmal , pandharipande and ravenhall @xcite at @xmath205 . the steep rise in @xmath1 in the higher neutron density region ( @xmath206 @xmath173 ) compared to those obtained from the hartree - fock theory using various skyrme interactions would help neutron - rich matter , which have larger dripped neutron density of @xmath206 @xmath173 , to be uniform . therefore , we can say that this behavior of @xmath1 for the present qmd model hamiltonian suppresses the density region in which the `` pasta '' phases are the most energetically favorable in neutron star matter and in the case of @xmath7 in the present study . we also note that @xmath1 at lower neutron densities of @xmath207 @xmath173 is relatively small . this , in turn , would lead to increase the density at which matter with lower dripped neutron density ( e.g. , the case of @xmath6 in the present study ) turns into uniform . next , we calculate the proton chemical potential @xmath2 in the pure neutron matter . we use the cold neutron matter prepared for the above calculation of @xmath1 as an initial condition . we insert a proton into this pure neutron matter , and then minimize the total energy by the frictional relaxation method with fixing the positions and momenta of the other neutrons . the position of the inserted proton is chosen randomly in the simulation box , and its momentum is chosen randomly from @xmath208 mev@xmath209 . we evaluate @xmath2 as the difference in the total energy between that before the insertion of the proton and that after the optimization of the position and the momentum of the proton .    in fig . [ fig mu_p0 ] , we plot @xmath2 for the present model hamiltonian . as can be seen from this figure , the result for the present model hamiltonian generally reproduce the data of the other results obtained from the hartree - fock theory using the various skyrme interactions at densities @xmath210 @xmath173 . at lower densities of @xmath211 @xmath173 , errors are quite large and data scatter significantly . this is because density fluctuations in pure neutron matter obtained by qmd would be unrealistically large at such low densities due to the fixed width of the wave packets in this model . however , it is noted that even in such a density region , our data are generally consistent with the other results mentioned above . finally , we turn to the surface tension , which affects energetically favorable nuclear shape most directly among the three quantities discussed here . we have calculated @xmath3 from the energy change induced by the change in the area of planar nuclear matter . we use nucleon gas with proton fractions @xmath7 , 0.2 , 0.3 and 0.5 composed of 1372 nucleons . in the calculation of @xmath3 , the coulomb interaction is excluded . in order to prepare slablike nuclear matter , we first cool down the above nucleon gas from @xmath90 mev to @xmath212 mev using eqs . ( [ qmdeom fric ] ) in a shallow trapping harmonic potential @xmath213 where @xmath214 mev @xmath215 . the simulation box here is imposed of the periodic boundary condition in the @xmath89- and @xmath216-directions , and is imposed of the open boundary condition in the @xmath217-direction . the box size @xmath218 in the @xmath89- and @xmath216-directions is set 20.26 fm . when the temperature reaches @xmath212 mev , we remove the trapping potential and , except for the case with @xmath0 , we change the boundary condition in the @xmath217-direction from open to periodic . the box size in the @xmath217-direction @xmath219 is chosen so as to at least all nucleons dripped outside the nuclear matter region can be contained in the box : @xmath220 fm ( @xmath7 ) , 79.12 fm ( @xmath221 ) , and 82.85 fm ( @xmath6 ) . after we relax the system for @xmath126 7000 fm/@xmath70 , we prepare three kinds of samples : one is nothing changed ( sample 1 ) and the others have the area of the @xmath222 side of the simulation box increased ( decreased ) by 1 @xmath223 with the total volume of the box kept constant [ sample 2 ( sample 3 ) ] . we further cool them down until @xmath224 kev . the resultant nucleon density profiles for the sample 1 of @xmath7 projected on the @xmath217-axis is shown in fig . [ fig profile ] . as can be seen from this figure , @xmath225 fm is much larger than the thickness of the slab @xmath226 fm , and thus the volume of the dripped neutron gas region is almost the same among the three kinds of samples because the volume change in the nuclear matter region is negligible . it is also noted that @xmath131 is much larger than the surface thickness @xmath227 fm , which ensures that the surfaces at both sides of the slab are separated well . therefore , we can say that the energy difference between the sample 2 and 3 is just due to the difference in the surface area of the planar nuclear matter . following the spirit of ravenhall , bennett and pethick ( hereafter rbp ) @xcite and by using the sample 1 , we define the proton fraction @xmath228 in the nuclear matter region as an averaged value for the region of the width of @xmath229 fm in the central part of the slab , where the proton and neutron density profiles ripple around constant values . we extract @xmath3 from the total energy @xmath230 of the sample 2 and @xmath231 of the sample 3 as follows : @xmath232 where @xmath233 is the area of the @xmath222 side of the sample 2 ( sample 3 ) given by @xmath234 @xmath235 ( @xmath236 @xmath235 ) . the factor 2 in the denominator represents , of course , the contribution of the two sides of the planar nuclear matter . as shown in fig . [ fig esurf ] , all the results plotted in this figure almost coincide with each other at @xmath237 , where uncertainty in the nuclear surface tension is rather small . these results deviate significantly at lower @xmath228 , and the values of @xmath3 of the present calculation lie between those obtained by baym , bethe and pethick ( hereafter bbp ) @xcite and by rbp @xcite , which is based on the hartree - fock calculation using a skyrme interaction , at @xmath238 . thus we can say that , in comparison with the result by rbp taken to be standard here , the contribution of @xmath3 of the present qmd model tends to favor uniform nuclear matter rather than the inhomogeneous `` pasta '' phases for lower @xmath228 of @xmath239 . now let us discuss the fact that the intermediate phases with a negative value of @xmath158 are obtained instead of the bubble phases in a wide density region for @xmath7 . in considering this problem , we should note the general tendency that @xmath201 increases and @xmath3 decreases as matter becomes neutron - rich ( see figs . [ fig mu_p0 ] and [ fig esurf ] ) . the density dependence of @xmath156 monotonically increases until @xmath114 as shown in fig . [ fig minko x0.1 ] would partly stem from the small value of @xmath3 . the apparently lower melting density @xmath114 in this case than in the cases of @xmath0 and 0.3 , even though @xmath3 is small , is due to a large @xmath201 , which increases typically of order 10 mev as @xmath240 increases . the small @xmath3 and the large @xmath201 in neutron - rich matter would help nuclear matter regions and neutron gas regions mix each other at the cost of small surface energy below @xmath114 . as a result , the structures with a negative @xmath158 could be favored . according to fig . [ fig mu_p0 ] , the quantity @xmath2 obtained for the present model hamiltonian is consistent with those for the other skyrme - hartree - fock calculations in the relevant region of @xmath241 @xmath173 . it is thus possible that a result which shows the structure of matter changes from negative @xmath158 to uniform without undergoing `` swiss cheese '' structure will be obtained by another calculation for neutron - rich matter using some framework without assuming nuclear shape .    in closing this subsection , we summarize the consequences of the resultant @xmath1 , @xmath2 and @xmath3 for the present qmd hamiltonian .    * for symmetric matter ( @xmath242 ) + according to @xmath3 at @xmath237 , the present model is consistent with the other results , and is an appropriate effective interaction for the study of the `` pasta '' phases at @xmath0 . * for neutron - rich cases + for neutron - rich cases such as @xmath243 , in which the dripped neutron density @xmath240 grows @xmath244 @xmath173 just before matter turns into uniform , the present qmd model can be taken as a conservative one in reproducing the `` pasta '' phases . its @xmath1 , @xmath2 and @xmath3 act to suppress the density region of the `` pasta '' phases compared to other skyrme - hartree - fock results . * for intermediate cases + at intermediate proton fraction of @xmath245 , @xmath1 of the present model acts to favor the inhomogeneous `` pasta''phases rather than the uniform phase and @xmath3 acts in the opposite way in comparison with other skyrme - hartree - fock results . here we would like to discuss astrophysical consequences of our results . pethick and potekhin have pointed out that elastic properties of `` pasta '' phases with rodlike and slablike nuclei are similar to those of liquid crystals , which stems from the similarity in the geometrical structures @xcite . it can also be said that the intermediate phases observed in the present work are `` spongelike '' ( or `` rubberlike '' for @xmath246 ) phases because they have both highly connected nuclear and bubble regions shown as @xmath247 . the elastic properties of the spongelike intermediate phases are qualitatively different from those of the liquid - crystal - like `` pasta '' phases because the former ones do not have any directions in which the restoring force does not act ; while the latter ones have . our results imply that the intermediate phases occupy a significant fraction of the density region in which nonspherical nuclei can be seen ( see figs . [ phase diagram 1 ] and [ phase diagram 2 ] ) . according to figs . [ phase diagram 1 ] and [ phase diagram 2 ] , we expect that the maximum elastic energy that can be stored in the neutron star crust and supernova inner core is higher than that in the case where all nonspherical nuclei have simple `` pasta '' structures . besides , the cylinder and the slab phases , which are liquid - crystal - like , lie between the spongelike intermediate phases or the crystalline solidlike phase , and the releasing of the strain energy would , in consequence , concentrate in the domain of these liquid - crystal - like phases . the above effects of the intermediate phases should be taken into account in considering the crust dynamics of starquakes and hydrodynamics of the core collapse , etc . if these phases exist in neutron star matter and supernova matter . in the context of pulsar glitch phenomena , the effects of the spongelike nuclei on the pinning rate and the creep velocity of superfluid neutron vortices also have yet to be investigated .    for neutrino cooling of neutron stars , some version of the direct urca process which is suggested by lorenz et al . @xcite , that this might be allowed in the `` pasta '' phases , would be suppressed in the intermediate phases . this is due to the fact that the proton spectrum at the fermi surface is no longer continuous in the spongelike nuclei . an important topic which we would like to mention is about the effects of the intermediate phases on neutrino trapping in supernova cores . the nuclear parts connect over a wide region which is much larger than that characterized by the typical neutrino wave length @xmath248 fm . thus the neutrino scattering processes are no longer coherent in contrast to the case of the spherical nuclei , and this may , in consequence , reduce the diffusion time scale of neutrinos as in the case of `` pasta '' phases with simple structures . this reduction softens the supernova matter and would thus act to enhance the amount of the released gravitational energy . it would be interesting to estimate the neutrino opacity of the spongelike phases and the `` pasta '' phases . finally , we would like to mention the thermal fluctuations with long wavelengths leading to displacements of `` pasta '' nuclei , which can not be incorporated into the simulations using a finite - size box . even if we succeed in reproducing the phase with slablike nuclei in neutron - rich matter in the future study , we should remind the above effect of thermal fluctuations to consider the real situation of matter in inner crusts of neutron stars . following the discussions in refs . @xcite by a liquid - drop model , it is likely that the extension of slablike nuclei is limited to a finite length scale of @xmath249 fm in the temperature regions typical for neutron star crusts and supernova cores . we have performed qmd simulations for matter with fixed proton fractions @xmath0 , 0.3 and 0.1 at various densities below the normal nuclear density . our calculations without any assumptions on the nuclear shape demonstrate that the `` pasta '' phases with rodlike nuclei , with slablike nuclei , with cylindrical bubbles and with spherical bubbles can be formed dynamically from hot uniform matter within the time scale of @xmath80 fm/@xmath70 in the proton - rich cases of @xmath0 and 0.3 . we also demonstrate that the `` pasta '' phase with cylindrical nuclei can be formed dynamically within the time scale of @xmath84 fm/@xmath70 for the neutron - rich case of @xmath7 . our results imply the existence of at least the phase with cylindrical nuclei in neutron star crusts because they cool down keeping the local thermal equilibrium after proto - neutron stars are formed and their cooling time scale , which is macroscopic one , is much larger than the relaxation time scale of our simulations . in addition to these `` pasta '' phases with simple structures , our results obtained here also suggest the existence of intermediate phases which have complicated nuclear shapes . we have systematically analyzed the structure of matter with two - point correlation functions and with morphological measures `` minkowski functionals '' , and have demonstrated how structure changes with increasing density . making use of a topological quantity called euler characteristic , which is one of the minkowski functionals , it has been found that the intermediate phases can be characterized as those with negative euler characteristic . this means that the intermediate phases have `` spongelike '' ( or `` rubberlike '' for @xmath246 ) structures which have both highly connected nuclear and bubble regions . the elastic properties of the spongelike intermediate phases are qualitatively different from those of the liquid - crystal - like `` pasta '' phases . we have also investigated the properties of the effective qmd interaction used in the present work in order to examine the validity of our results . important quantities which affect the structure of matter are the energy per nucleon @xmath1 of the pure neutron matter , the proton chemical potential @xmath2 in pure neutron matter and the nuclear surface tension @xmath3 . the resultant these quantities show that the present qmd interaction has generally reasonable properties at subnuclear densities among other nuclear interactions . it is thus concluded that our results are not exceptional ones in terms of nuclear forces . our results which suggest the existence of the highly connected intermediate phases as well as the simple `` pasta '' phases provide a vivid picture that matter in neutron star inner crusts and supernova inner cores has a variety of material phases . the stellar region which we have tried to understand throughout this paper is relatively tiny , but there has quite rich properties which stem from the fancy structures of dense matter . g. w. is grateful to t. maruyama , k. iida , a. tohsaki , h. horiuchi , k. oyamatsu , h. takemoto , k. niita , s. chikazumi , i. kayo , c. hikage , t. buchert , n. itoh and k. kotake for helpful discussions and comments . g. w. also appreciate j. lee for checking the manuscript . this work was supported in part by the junior research associate program in riken through research grant no . j130026 and by grants - in - aid for scientific research provided by the ministry of education , culture , sports , science and technology through research grant ( s ) no . 14102004 , no . 14079202 and no . 14 - 7939 . in calculating a long - range interaction , such as the coulomb interaction , it is necessary to sum up all contributions of particles at a sufficiently far distance . the ewald sum is a familiar technique for efficiently computing the long - range contributions in a system with the periodic boundary condition ( see , e.g. , refs . @xcite ; recent mathematically careful discussion relating to the conditionally convergence of the coulomb energy can be seen , e.g. , in ref . the basic idea of the ewald sum is that the contributions of particles in a long distance in real space can be calculated as contributions in the neighborhood in fourier space : the contributions of particles in a short distance is summed up in real space and those of particles in a long distance is summed up in fourier space . let us consider a system consisting of charged particles which have a gaussian charge distribution and a uniform background charge which cancels the total charge of charged particles . these @xmath11 particles are assumed to be in a cubic simulation box with volume @xmath250 which is imposed by the periodic boundary condition .    if every particle @xmath22 with total charge @xmath251 is surrounded by a gaussian charge distribution with total charge @xmath252 , the electrostatic interaction of particle @xmath22 turns into a screened short - range interaction . thus the total coulomb energy @xmath253 of this system can be decomposed into as follows ( see fig . [ fig ewald ] ) : @xmath254 where @xmath255 is the sum of the coulomb energy between an each unscreened charged particle @xmath22 and the other charged particles with screening charges , @xmath256 is that between an each unscreened charged particle @xmath22 and compensating charges which cancels the screening gaussian charges , and @xmath257 is the sum of spurious self interactions between an each charged particle @xmath22 and its compensating charge @xmath22 . charge densities of the real charged particles @xmath184 and of the screening charges @xmath258 can be written as @xmath259 where @xmath260 and @xmath261 are the reciprocals of the widths of charge distributions for a charged particle and a screening charge , respectively , and * n * denotes a position vector of a periodic image normalized by @xmath262 . thus a distribution of screened charges @xmath263 is @xmath264 the electrostatic potential @xmath265 due to @xmath266 can be obtained as      because a solution of the poisson equation @xmath268\ ] ] is @xmath269 where erf(@xmath89 ) and erfc(@xmath89 ) are the error function and the complementary error function , respectively , and they are defined as @xmath270 and @xmath271 . here , we determine the constant @xmath272 so that the average value of @xmath273 in the simulation box @xmath154 be zero : @xmath274 thus , @xmath265 leads to              next , we calculate @xmath256 , which is the sum of the coulomb energy between an unscreened charged particle @xmath22 and a charge density @xmath281 which consists of the compensating charges and the background charge : @xmath282 fourier transforming the charge distribution @xmath281 yields @xmath283    using @xmath284 and the poisson equation @xmath285 in the fourier form @xmath286 we can at once obtain the electrostatic potential @xmath287 due to the charge density @xmath288 : @xmath289 the term with @xmath290 is canceled due to charge neutrality . thus the coulomb energy @xmath256 is given by @xmath291 we have to subtract a sum of spurious self interactions @xmath292 between a charged particle @xmath22 and its compensating charge from @xmath256 . according to eq . ( [ solution ] ) , an electrostatic potential @xmath293 due to a gaussian compensating charge is @xmath294 , thus the self interaction of particle @xmath22 reads to @xmath295 and hence , @xmath296 finally , the total coulomb energy can be calculated by eqs . ( [ coul ] ) , ( [ u_short ] ) , ( [ u1 ] ) and ( [ u_self ] ) . the positive background charge does not appear explicitly because the average value of @xmath273 within the simulation box is set to be zero .    the total coulomb energy per particle @xmath297 and the @xmath89 component of the coulomb force @xmath298 acting on a particle for various values of @xmath261 are plotted in fig . [ fig ewald ] . in this calculation , we use 1024 positive charged particles ( protons ) distributed randomly in a simulation box of @xmath299 fm ( i.e. @xmath300 ) , which is imposed of the periodic boundary condition . figs . [ fig ewald ] ( a ) and [ fig ewald ] ( b ) show the results for point charges and for gaussian charge distributions , respectively , which are calculated for the same configuration of the particle positions @xmath301 . the width of the gaussian charge distributions is set to be @xmath302 with @xmath303 @xmath235 , which corresponds to the width of the wave packets in the qmd model used in this work . we note that there are plateau regions of @xmath297 and @xmath298 whose values do not depend on @xmath261 . these constant values give results to be obtained . we note that the range of @xmath304 of the plateau regions become larger with increasing @xmath305 , where @xmath305 is the cutoff radius in the unit of @xmath306 for the summation in fourier space . as can be seen from fig . [ fig ewald ] , @xmath261 dependences of @xmath297 and @xmath298 for the present qmd model with a finite width of the gaussian charge distributions are weaker than those for the point charges . these features are also confirmed for different proton number densities of 0.2 and 0.3 @xmath127 . in our simulations , @xmath261 is set 13 or 14 , which are considered to be large enough to calculate the total coulomb energy per particle in accuracy less than @xmath307 kev , which is the typical value of the energy difference between successive `` pasta '' phases in neutron star matter obtained by previous works . bethe , g. borner and k. sato , astron .  astrophys . * 7 * , 279 ( 1970 ) . g. baym , h. a. bethe and c. j. pethick , nucl . a175 * , 225 ( 1971 ) . k. sato , prog . * 54 * , 1325 ( 1975 ) . h. a. bethe , rev . phys . * 62 * , 801 ( 1990 ) . c. j. pethick and d. g. ravenhall , annu . nucl .  part . * 45 * , 429 ( 1995 ) . d. g. ravenhall , c. j. pethick and j. r. wilson , phys . lett . * 50 * , 2066 ( 1983 ) . m. hashimoto , h. seki and m. yamada , prog . 71 * , 320 ( 1984 ) . lorenz , d.g . ravenhall and c.j . pethick , phys . lett . * 70 * , 379 ( 1993 ) . g. watanabe , k. iida and k. sato , nucl . * a676 * , 455 ( 2000 ) ; erratum , nucl .  phys . a ( in press ) . g. watanabe , k. iida and k. sato , nucl . * a687 * , 512 ( 2001 ) ; erratum , nucl . a ( in press ) . r. d. williams and s. e. koonin , nucl . a435 * , 844 ( 1985 ) . m. lassaut , h. flocard , p. bonche , p.h . heenen and e. suraud , astron .  astrophys . * 183 * , l3 ( 1987 ) . k. oyamatsu , nucl .  phys . * a561 * , 431 ( 1993 ) . b. jones , mon .  not . .  soc . * 296 * , 217 ( 1998 ) ; mon .  not . .  soc . * 306 * , 327 ( 1998 ) ; phys . lett . * 83 * , 3589 ( 1999 ) . c. j. pethick and a.y . potekhin , phys . b427 * , 7 ( 1998 ) . c. j. horowitz , phys . d * 55 * ( 1997 ) 4577 . k. iida , g. watanabe and k. sato , prog . * 106 * , 551 ( 2001 ) . j. aichelin and h. stcker , phys .  lett . * b176 * , 14 ( 1986 ) . h. feldmeier and j. schnack , rev . phys . * 72 * , 655 ( 2000 ) . g. watanabe , k. sato , k. yasuoka and t. ebisuzaki , phys . c * 66 * , 012801(r ) ( 2002 ) . h. feldmeier , nucl .  phys . * a515 * , 147 ( 1990 ) . a. ono , h. horiuchi , t. maruyama and a. ohnishi , prog .  theor .  phys . * 87 * , 1185 ( 1992 ) ; phys .  rev .  lett . * 68 * , 2898 ( 1992 ) . t. maruyama , k. niita , k. oyamatsu , t. maruyama , s. chiba and a. iwamoto , phys . c * 57 * , 655 ( 1998 ) . in some cases , the nos - hoover - like thermostat which will be described in another paper @xcite is also used . in the present study , employment of this thermostat is almost limited in a higher temperature region where matter is completely uniform . g. watanabe , k. sato , k. yasuoka and t. ebisuzaki , in preparation . m. p. allen and d. j. tildesley , _ computer simulation of liquids _ ( clarendon , oxford , 1987 ) . t. maruyama , private communication . s. nos , j.  chem . 81 * , 511 ( 1984 ) . w. g. hoover , phys . rev .  a * 31 * , 1695 ( 1985 ) . r. w. hockney and j. w. eastwood , _ computer simulation using particles _ ( iop publishing , bristol , 1988 ) . r. ogasawara and k. sato , prog . * 68 * , 222 ( 1982 ) . k. michielsen and h. de raedt , phys .  rep . * 347 * , 461 ( 2001 ) . j. schmalzing , s. engineer , t. buchert , v. sahni and s. shandarin , ( 1996 ) unpublished . j. r. gott iii , a. l. melott and m. dickinson , astrophys . j.  * 306 * , 341 ( 1986 ) ; d. h. weinberg , publ . pacific * 100 * , 1373 ( 1988 ) . g. watanabe , phd thesis , univ . of tokyo ( 2003 ) . existence of some phase with a complicated structure is also suggested by p. magierski and p. h. heenen , phys . c * 65 * , 045804 ( 2002 ) . they focus on shell effects of dripped neutrons . non - integral values of the `` dimensionality '' @xmath131 in the early work done by ravenhall et al . @xcite also suggest the possibilty of complicated nuclear shapes . akmal , v. r. pandharipande and d. g. ravenhall , phys . c * 58 * , 1804 ( 1998 ) . b. friedman and v. r. pandharipande , nucl . a361 * , 501 ( 1981 ) . c. j. pethick , d. g. ravenhall and c. p. lorenz , nucl .  phys . * a584 * , 675 ( 1995 ) . f. douchin and p. haensel , phys .  lett . * b485 * , 107 ( 2000 ) . o. sjberg , nucl . a222 * , 161 ( 1974 ) . j. siemens and v. r. pandharipande , nucl . a173 * , 561 ( 1971 ) . d. g. ravenhall , c. d. bennett and c. j. pethick , phys . rev .  lett . * 28 * , 978 ( 1972 ) . m. p. tosi , in : f. seitz and d. turnbull ( eds . ) , _ solid state physics : advances in research and applications _ ( academic press , new york , 1964 ) p. 1 . h. takemoto , t. ohyama and a. tohsaki , prog . .  phys . * 109 * , 563 ( 2003 ) .
Abstract:

Answer: structure of cold nuclear matter at subnuclear densities for the proton fraction @xmath0 , 0.3 and 0.1 is investigated by quantum molecular dynamics ( qmd ) simulations . we demonstrate that the phases with slablike and rodlike nuclei , etc . can be formed dynamically from hot uniform nuclear matter without any assumptions on nuclear shape , and also systematically analyze the structure of cold matter using two - point correlation functions and minkowski functionals . in our simulations , we also observe intermediate phases , which has complicated nuclear shapes . it is found out that these phases can be characterized as those with negative euler characteristic . our result implies the existence of these kinds of phases in addition to the simple `` pasta '' phases in neutron star crusts and supernova inner cores . in addition , we investigate the properties of the effective qmd interaction used in the present work to examine the validity of our results . the resultant energy per nucleon @xmath1 of the pure neutron matter , the proton chemical @xmath2 in pure neutron matter and the nuclear surface tension @xmath3 are generally reasonable in comparison with other nuclear interactions .