BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

D. K. Berry

D. K. Berry appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
nucl-thastro-ph.SRastro-ph.HEcond-mat.stat-mechphysics.plasm-phastro-phBiological Physicscond-mat.othercond-mat.softhep-ph

Catalog footprint

What is connected

14works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2016arXiv

Effect of topological defects on "nuclear pasta" observables

[Background] The "pasta" phase of nuclear matter may play an important role in the structure and evolution of neutron stars. Recent works suggest nuclear pasta has a high resistivity which could be explained by the presence of long lived topological defects. The defects act as impurities that decrease thermal and electrical conductivity of the pasta. [Purpose] To quantify how topological defects affect transport properties of nuclear pasta and estimate this effect using an impurity parameter $Q_{\text{imp}}$. [Methods] Contrast molecular dynamics simulations of up to 409\,600 nucleons arranged in parallel nuclear pasta slabs (perfect pasta) with simulations of pasta slabs connected by topological defects (impure pasta). From these simulations compare the viscosity and heat conductivity of perfect and impure pasta to obtain an effective impurity parameter $Q_{\text{imp}}$ due to the presence of defects. [Results] Both the viscosity and thermal conductivity calculated for both perfect and impure pasta are anisotropic, peaking along directions perpendicular to the slabs and reaching a minimum close to zero parallel to them. In our 409\,600 nucleon simulation topological defects connecting slabs of pasta reduce both the thermal conductivity and viscosity on average by about 37\%. We estimate an effective impurity parameter due to the defects of order $Q_{\text{imp}}\sim30$. [Conclusions] Topological defects in the pasta phase of nuclear matter have an effect similar to impurities in a crystal lattice. The irregularities introduced by the defects reduce the thermal and electrical conductivities and the viscosity of the system. This effect can be parameterized by a large impurity parameter $Q_{\text{imp}}\sim30$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nuclear pasta and supernova neutrinos at late times

Nuclear pasta, with nucleons arranged into tubes, sheets, or other complex shapes, is expected in core collapse supernovae (SNe) at just below nuclear density. We calculate the additional opacity from neutrino-pasta coherent scattering using molecular dynamics simulations. We approximately include this opacity in simulations of SNe. We find that pasta slows neutrino diffusion and greatly increases the neutrino signal at late times of 10 or more seconds after stellar core collapse. This signal, for a galactic SN, should be clearly visible in large detectors such as Super-Kamiokande.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Parking-garage structures in astrophysics and biophysics

A striking shape was recently observed for the cellular organelle endoplasmic reticulum consisting of stacked sheets connected by helical ramps. This shape is interesting both for its biological function, to synthesize proteins using an increased surface area for ribosome factories, and its geometric properties that may be insensitive to details of the microscopic interactions. In the present work, we find very similar shapes in our molecular dynamics simulations of the nuclear pasta phases of dense nuclear matter that are expected deep in the crust of neutron stars. There are dramatic differences between nuclear pasta and terrestrial cell biology. Nuclear pasta is 14 orders of magnitude denser than the aqueous environs of the cell nucleus and involves strong interactions between protons and neutrons, while cellular scale biology is dominated by the entropy of water and complex assemblies of biomolecules. Nonetheless the very similar geometry suggests both systems may have similar coarse-grained dynamics and that the shapes are indeed determined by geometrical considerations, independent of microscopic details. Many of our simulations self-assemble into flat sheets connected by helical ramps. These ramps may impact the thermal and electrical conductivities, viscosity, shear modulus, and breaking strain of neutron star crust. The interaction we use, with Coulomb frustration, may provide a simple model system that reproduces many biologically important shapes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Disordered nuclear pasta, magnetic field decay, and crust cooling in neutron stars

Nuclear pasta, with non-spherical shapes, is expected near the base of the crust in neutron stars. Large scale molecular dynamics simulations of pasta show long lived topological defects that could increase electron scattering and reduce both the thermal and electrical conductivities. We model a possible low conductivity pasta layer by increasing an impurity parameter Q_{imp}. Predictions of light curves for the low mass X-ray binary MXB 1659-29, assuming a large Q_{imp}, find continued late time cooling that is consistent with Chandra observations. The electrical and thermal conductivities are likely related. Therefore observations of late time crust cooling can provide insight on the electrical conductivity and the possible decay of neutron star magnetic fields (assuming these are supported by currents in the crust).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nuclear Waffles

The dense neutron-rich matter found in supernovae and neutron stars is expected to form complex nonuniform phases referred to as nuclear pasta. The pasta shapes depend on density, temperature and proton fraction and determine many transport properties in supernovae and neutron star crusts. We use two recently developed hybrid CPU/GPU codes to perform large scale molecular dynamics (MD) simulations with $51200$ and $409600$ nucleons of nuclear pasta. From the output of the MD simulations we characterize the topology and compute two observables, the radial distribution function $g(r)$ and the structure factor $S(q)$, for systems with proton fractions $Y_p=0.10, 0.20, 0.30$ and $0.40$ at about one third of nuclear saturation density and temperatures near $1.0$ MeV. We observe that the two lowest proton fraction systems simulated, $Y_p=0.10$ and $0.20$, equilibrate quickly and form liquid-like structures. Meanwhile, the two higher proton fraction systems, $Y_p=0.30$ and $0.40$, take a longer time to equilibrate and organize themselves in solid-like periodic structures. Furthermore, the $Y_p=0.40$ system is made up of slabs, lasagna phase, interconnected by defects while the $Y_p=0.30$ systems consist of a stack of perforated plates, the nuclear waffle phase. The periodic configurations observed in our MD simulations for proton fractions $Y_p\ge0.30$ have important consequences for the structure factors $S(q)$ of protons and neutrons, which relate to many transport properties of supernovae and neutron star crust. A detailed study of the waffle phase and how its structure depends on temperature, size of the simulation and the screening length showed that finite-size effects appear to be under control and, also, that the plates in the waffle phase merge at temperatures slightly above $1.0$ MeV and the holes in the plates form an hexagonal lattice at temperatures slightly lower than $1.0$ MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Pasta Nucleosynthesis: Molecular dynamics simulations of nuclear statistical equilibrium

Background: Exotic non-spherical nuclear pasta shapes are expected in nuclear matter at just below saturation density because of competition between short range nuclear attraction and long range Coulomb repulsion. Purpose: We explore the impact of nuclear pasta on nucleosynthesis, during neutron star mergers, as cold dense nuclear matter is ejected and decompressed. Methods: We perform classical molecular dynamics simulations with 51200 and 409600 nucleons, that are run on GPUs. We expand our simulation region to decompress systems from an initial density of 0.080 fm^{-3} down to 0.00125 fm^{-3}. We study proton fractions of Y_P=0.05, 0.10, 0.20, 0.30, and 0.40 at T =0.5, 0.75, and 1.0 MeV. We calculate the composition of the resulting systems using a cluster algorithm. Results: We find final compositions that are in good agreement with nuclear statistical equilibrium models for temperatures of 0.75 and 1 MeV. However, for proton fractions greater than Y_P=0.2 at a temperature of T = 0.5 MeV, the MD simulations produce non-equilibrium results with large rod-like nuclei. Conclusions: Our MD model is valid at higher densities than simple nuclear statistical equilibrium models and may help determine the initial temperatures and proton fractions of matter ejected in mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nuclear Pasta Formation

The formation of complex nonuniform phases of nuclear matter, known as nuclear pasta, is studied with molecular dynamics simulations containing 51200 nucleons. A phenomenological nuclear interaction is used that reproduces the saturation binding energy and density of nuclear matter. Systems are prepared at an initial density of 0.10fm$^{-3}$ and then the density is decreased by expanding the simulation volume at different rates to densities of 0.01 fm$^{-3}$ or less. An originally uniform system of nuclear matter is observed to form spherical bubbles ("swiss cheese"), hollow tubes, flat plates ("lasagna"), thin rods ("spaghetti") and, finally, nearly spherical nuclei with decreasing density. We explicitly observe nucleation mechanisms, with decreasing density, for these different pasta phase transitions. Topological quantities known as Minkowski functionals are obtained to characterize the pasta shapes. Different pasta shapes are observed depending on the expansion rate. This indicates non equilibrium effects. We use this to determine the best ways to obtain lower energy states of the pasta system from MD simulations and to place constrains on the equilibration time of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Direct MD simulation of liquid-solid phase equilibria for three-component plasma

The neutron rich isotope 22Ne may be a significant impurity in carbon and oxygen white dwarfs and could impact how the stars freeze. We perform molecular dynamics simulations to determine the influence of 22Ne in carbon-oxygen-neon systems on liquid-solid phase equilibria. Both liquid and solid phases are present simultaneously in our simulation volumes. We identify liquid, solid, and interface regions in our simulations using a bond angle metric. In general we find good agreement for the composition of liquid and solid phases between our MD simulations and the semi analytic model of Medin and Cumming. The trace presence of a third component, neon, does not appear to strongly impact the chemical separation found previously for two component carbon and oxygen systems. This suggests that small amounts of 22Ne may not qualitatively change how the material in white dwarf stars freezes. However, we do find systematically lower melting temperatures (higher Gamma) in our MD simulations compared to the semi analytic model. This difference seems to grow with impurity parameter Q_imp and suggests a problem with simple corrections to the linear mixing rule for the free energy of multicomponent solid mixtures that is used in the semi analytic model.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Direct MD simulation of liquid-solid phase equilibria for two-component plasmas

We determine the liquid-solid phase diagram for carbon-oxygen and oxygen-selenium plasma mixtures using two-phase MD simulations. We identified liquid, solid, and interface regions using a bond angle metric. To study finite size effects, we perform 27648 and 55296 ion simulations. To help monitor non-equilibrium effects, we calculate diffusion constants $D_i$. For the carbon-oxygen system we find that $D_O$ for oxygen ions in the solid is much smaller than $D_C$ for carbon ions and that both diffusion constants are 80 or more times smaller than diffusion constants in the liquid phase. There is excellent agreement between our carbon-oxygen phase diagram and that predicted by Medin and Cumming. This suggests that errors from finite size and non-equilibrium effects are small and that the carbon-oxygen phase diagram is now accurately known. The oxygen-selenium system is a simple two-component model for more complex rapid proton capture nucleosynthesis ash compositions for an accreting neutron star. Diffusion of oxygen, in a predominately selenium crystal, is remarkably fast, comparable to diffusion in the liquid phase. We find a somewhat lower melting temperature for the oxygen-selenium system than that predicted by Medin and Cumming. This is probably because of electron screening effects.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Diffusion in Coulomb Crystals

Diffusion in coulomb crystals can be important for the structure of neutron star crusts. We determine diffusion constants $D$ from molecular dynamics simulations. We find that $D$ for coulomb crystals with relatively soft-core $1/r$ interactions may be larger than $D$ for Lennard-Jones or other solids with harder-core interactions. Diffusion, for simulations of nearly perfect body-centered-cubic lattices, involves the exchange of ions in ring-like configurations. Here ions "hop" in unison without the formation of long lived vacancies. Diffusion, for imperfect crystals, involves the motion of defects. Finally, we find that diffusion, for an amorphous system rapidly quenched from coulomb parameter $Γ=175$ to coulomb parameters up to $Γ=1750$, is fast enough so that the system starts to crystallize during long simulation runs. These results strongly suggest that coulomb solids in cold white dwarf stars, and the crust of neutron stars, will be crystalline and not amorphous.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Neutron Star Crust and Molecular Dynamics Simulation

In this book chapter we review plasma crystals in the laboratory, in the interior of white dwarf stars, and in the crust of neutron stars. We describe a molecular dynamics formalism and show results for many neutron star crust properties including phase separation upon freezing, diffusion, breaking strain, shear viscosity and dynamics response of nuclear pasta. We end with a summary and discuss open questions and challenges for the future.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Crystallization of Carbon Oxygen Mixtures in White Dwarf Stars

We determine the phase diagram for dense carbon/ oxygen mixtures in White Dwarf (WD) star interiors using molecular dynamics simulations involving liquid and solid phases. Our phase diagram agrees well with predictions from Ogata et al. and Medin and Cumming and gives lower melting temperatures than Segretain et al. Observations of WD crystallization in the globular cluster NGC 6397 by Winget et al. suggest that the melting temperature of WD cores is close to that for pure carbon. If this is true, our phase diagram implies that the central oxygen abundance in these stars is less than about 60%. This constraint, along with assumptions about convection in stellar evolution models, limits the effective S factor for the $^{12}$C($α,γ$)$^{16}$O reaction to S_{300} <= 170 keV barns.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Diffusion of Neon in White Dwarf Stars

Sedimentation of the neutron rich isotope $^{22}$Ne may be an important source of gravitational energy during the cooling of white dwarf stars. This depends on the diffusion constant for $^{22}$Ne in strongly coupled plasma mixtures. We calculate self-diffusion constants $D_i$ from molecular dynamics simulations of carbon, oxygen, and neon mixtures. We find that $D_i$ in a mixture does not differ greatly from earlier one component plasma results. For strong coupling (coulomb parameter $Γ>$ few), $D_i$ has a modest dependence on the charge $Z_i$ of the ion species, $D_i \propto Z_i^{-2/3}$. However $D_i$ depends more strongly on $Z_i$ for weak coupling (smaller $Γ$). We conclude that the self-diffusion constant $D_{\rm Ne}$ for $^{22}$Ne in carbon, oxygen, and neon plasma mixtures is accurately known so that uncertainties in $D_{\rm Ne}$ should be unimportant for simulations of white dwarf cooling.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Dynamical response of the nuclear pasta in neutron star crusts

The nuclear pasta -- a novel state of matter having nucleons arranged in a variety of complex shapes -- is expected to be found in the crust of neutron stars and in core-collapse supernovae at subnuclear densities of about $10^{14}$ g/cm$^3$. Due to frustration, a phenomenon that emerges from the competition between short-range nuclear attraction and long-range Coulomb repulsion, the nuclear pasta displays a preponderance of unique low-energy excitations. These excitations could have a strong impact on many transport properties, such as neutrino propagation through stellar environments. The excitation spectrum of the nuclear pasta is computed via a molecular-dynamics simulation involving up to 100,000 nucleons. The dynamic response of the pasta displays a classical plasma oscillation in the 1-2 MeV region. In addition, substantial strength is found at low energies. Yet this low-energy strength is missing from a simple ion model containing a single-representative heavy nucleus. The low-energy strength observed in the dynamic response of the pasta is likely to be a density wave involving the internal degrees of freedom of the clusters.

0 citations0 reviewsNo code linkedOpen access