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

Andrew Cumming

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

ResearcherUnclaimed source record
astro-ph.HEastro-ph.SRastro-phastro-ph.EPphysics.plasm-phastro-ph.GAastro-ph.IMgr-qcnucl-th

Catalog footprint

What is connected

32works
9topics
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

32 published item(s)

preprint2020arXiv

The Challenge of Forming a Fuzzy Core in Jupiter

Recent structure models of Jupiter that match Juno gravity data suggest that the planet harbours an extended region in its deep interior that is enriched in heavy elements, a so-called dilute/fuzzy core. This finding raises the question of what possible formation pathways would lead to such a structure. We model Jupiter's formation and long-term evolution, starting at late-stage formation before runaway gas accretion. The formation scenarios we consider include both primordial composition gradients, as well as gradients that are built as proto-Jupiter rapidly acquires its gaseous envelope. We then follow Jupiter's evolution as it cools down and contracts, with a particular focus on the energy and material transport in the interior. We find that none of the scenarios we consider lead to a fuzzy core that is compatible with interior structure models. In all the cases, most of Jupiter's envelope becomes convective and fully mixed after at most a few million years. This is true even when we consider a case where the gas accretion leads to a cold planet, and large amounts of heavy elements are accreted. It is therefore concluded that it is very challenging to explain Jupiter's dilute core from standard formation models. We suggest that future work should consider more complex formation pathways as well as the modelling of additional physical processes that could lead to Jupiter's current-state internal structure.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The effect of diffusive nuclear burning in neutron star envelopes on cooling in accreting systems

Valuable information about the neutron star interior can be obtained by comparing observations of thermal radiation from a cooling neutron star crust with theoretical models. Nuclear burning of lighter elements that diffuse to deeper layers of the envelope can alter the relation between surface and interior temperatures and can change the chemical composition over time. We calculate new temperature relations and consider two effects of diffusive nuclear burning (DNB) for H-C envelopes. First, we consider the effect of a changing envelope composition and find that hydrogen is consumed on short timescales and our temperature evolution simulations correspond to those of a hydrogen-poor envelope within ~100 days. The transition from a hydrogen-rich to a hydrogen-poor envelope is potentially observable in accreting NS systems as an additional initial decline in surface temperature at early times after the outburst. Second, we find that DNB can produce a non-negligible heat flux, such that the total luminosity can be dominated by DNB in the envelope rather than heat from the deep interior. However, without continual accretion, heating by DNB in H-C envelopes is only relevant for <1-80 days after the end of an accretion outburst, as the amount of light elements is rapidly depleted. Comparison to crust cooling data shows that DNB does not remove the need for an additional shallow heating source. We conclude that solving the time-dependent equations of the burning region in the envelope self-consistently in thermal evolution models instead of using static temperature relations would be valuable in future cooling studies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Imprint of the Protoplanetary Disk in the Accretion of Super-Earth Envelopes

Super-Earths are by far the most dominant type of exoplanet, yet their formation is still not well understood. In particular, planet formation models predict that many of them should have accreted enough gas to become gas giants. Here we examine the role of the protoplanetary disk in the cooling and contraction of the protoplanetary envelope. In particular, we investigate the effects of 1) the thermal state of the disk as set by the relative size of heating by accretion or irradiation, and whether its energy is transported by radiation or convection, and 2) advection of entropy into the outer envelope by disk flows that penetrate the Hill sphere, as found in 3D global simulations. We find that, at 5 and 1 AU, this flow at the level reported in the non-isothermal simulations where it penetrates only to ~ 0.3 times the Hill radius has little effect on the cooling rate since most of the envelope mass is concentrated close to the core, and far from the flow. On the other hand, at 0.1 AU, the envelope quickly becomes fully-radiative, nearly isothermal, and thus cannot cool down, stalling gas accretion. This effect is significantly more pronounced in convective disks, leading to envelope mass orders of magnitude lower. Entropy advection at 0.1 AU in either radiative or convective disks could therefore explain why super-Earths failed to undergo runaway accretion. These results highlight the importance of the conditions and energy transport in the protoplanetary disk for the accretion of planetary envelopes.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

LRP2020: The cosmic origin and evolution of the elements

The origin of many elements of the periodic table remains an unsolved problem. While many nucleosynthetic channels are broadly understood, significant uncertainties remain regarding certain groups of elements such as the intermediate and rapid neutron-capture processes, the p-process, or the origin of odd-Z elements in the most metal-poor stars. Canada has a long tradition of leadership in nuclear astrophysics, dating back to the work of Alastair Cameron in the 1950s. Recent faculty hires have further boosted activity in the field, including transient observation and theory, survey science on galactic nucleosynthesis, and nuclear experiments. This white paper contains a brief overview of recent activity in the community, highlighting strengths in each sub-field, and provides recommendations to improve interdisciplinary collaboration. Sustaining Canadian leadership in the next decade will require, on the observational side, access to transient and non-transient surveys like LSST, SKA, or MSE, support for target-of-opportunity observing in current and future Canadian telescopes, and participation in next-generation X-ray telescopes such as ATHENA. State-of-the-art theoretical predictions will require an ambitious succession plan for the Niagara supercomputer to support large parallel jobs. We propose a funding instrument for postdoctoral training that reflects the interdisciplinary nature of nuclear astrophysics research, and the creation of a national collaborative funding program that allows for joint projects and workshop organization.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Survey of Chemical Separation in Accreting Neutron Stars

The heavy element ashes of rp-process hydrogen and helium burning in accreting neutron stars are compressed to high density where they freeze, forming the outer crust of the star. We calculate the chemical separation on freezing for a number of different nuclear mixtures resulting from a range of burning conditions for the rp-process. We confirm the generic result that light nuclei are preferentially retained in the liquid and heavy nuclei in the solid. This is in agreement with the previous study of a 17-component mixture of rp-process ashes by Horowitz et al. (2007), but extends that result to a much larger range of compositions. We also find an alternate phase separation regime for the lightest ash mixtures which does not demonstrate this generic behaviour. With a few exceptions, we find that chemical separation reduces the expected $Q_{\rm imp}$ in the outer crust compared to the initial rp-process ash, where $Q_{\rm imp}$ measures the mean-square dispersion in atomic number $Z$ of the nuclei in the mixture. We find that the fractional spread of $Z$ plays a role in setting the amount of chemical separation and is strongly correlated to the divergence between the two/three-component approximations and the full component model. The contrast in $Y_e$ between the initial rp-process ashes and the equilibrium liquid composition is similar to that assumed in earlier two-component models of compositionally driven convection, except for very light compositions which produce nearly negligible convective driving. We discuss the implications of these results for observations of accreting neutron stars.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The thermal state of KS 1731-260 after 14.5 years in quiescence

Crustal cooling of accretion-heated neutron stars provides insight into the stellar interior of neutron stars. The neutron star X-ray transient, KS~1731$-$260, was in outburst for 12.5 years before returning to quiescence in 2001. We have monitored the cooling of this source since then through {\it Chandra} and {\it XMM-Newton} observations. Here, we present a 150 ks {\it Chandra} observation of KS~1731$-$260 taken in August 2015, about 14.5 years into quiescence, and 6 years after the previous observation. We find that the neutron star surface temperature is consistent with the previous observation, suggesting that crustal cooling has likely stopped and the crust has reached thermal equilibrium with the core. Using a theoretical crust thermal evolution code, we fit the observed cooling curves and constrain the core temperature (T$_c = 9.35\pm0.25\times10^7$ K), composition (Q$_{imp} = 4.4^{+2.2}_{-0.5}$) and level of extra shallow heating required (Q$_{sh} = 1.36\pm0.18$ MeV/nucleon). We find that the presence of a low thermal conductivity layer, as expected from nuclear pasta, is not required to fit the cooling curve well, but cannot be excluded either.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Urca cooling pairs in the neutron star ocean and their effect on superbursts

An accretion outburst onto a neutron star deposits hydrogen-rich and/or helium-rich material into the neutron star's envelope. Thermonuclear burning of accreted material robustly produces Urca pairs\textrm{---}pairs of nuclei that undergo cycles of $e^-$-capture and $β^-$-decay. The strong $T^5$ dependence of the Urca cooling neutrino luminosity means that Urca pairs in the neutron star interior potentially remove heat from accretion-driven nuclear reactions. In this study, we identify Urca pairs in the neutron star's ocean \textrm{---} a plasma of ions and electrons overlaying the neutron star crust \textrm{---} and demonstrate that Urca cooling occurs at all depths in the ocean. We find that Urca pairs in the ocean and crust lower the ocean's steady-state temperature during an accretion outburst and that unstable carbon ignition, which is thought to trigger superbursts, occurs deeper than it would otherwise. Cooling superburst light curves, however, are only marginally impacted by cooling from Urca pairs because the superburst peak radiative luminosity $L_{\rm peak}$ is always much greater than the Urca pair neutrino luminosity $L_ν$ in the hot post-superburst ocean.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A strong shallow heat source in the accreting neutron star MAXI J0556-332

An accretion outburst in an X-ray transient deposits material onto the neutron star primary; this accumulation of matter induces reactions in the neutron star's crust. During the accretion outburst these reactions heat the crust out of thermal equilibrium with the core. When accretion halts, the crust cools to its long-term equilibrium temperature on observable timescales. Here we examine the accreting neutron star transient MAXI J0556-332, which is the hottest transient, at the start of quiescence, observed to date. Models of the quiescent light curve require a large deposition of heat in the shallow outer crust from an unknown source. The additional heat injected is $\approx 4\textrm{-}10\,\mathrm{MeV}$ per accreted nucleon; when the observed decline in accretion rate at the end of the outburst is accounted for, the required heating increases to $\approx 6\textrm{-}16\,\mathrm{MeV}$. This shallow heating is still required to fit the lightcurve even after taking into account a second accretion episode, uncertainties in distance, and different surface gravities. The amount of shallow heating is larger than that inferred for other neutron star transients and is larger than can be supplied by nuclear reactions or compositionally driven convection; but it is consistent with stored mechanical energy in the accretion disk. The high crust temperature ($T_b \gtrsim 10^{9} \, {\rm K}$) makes its cooling behavior in quiescence largely independent of the crust composition and envelope properties, so that future observations will probe the gravity of the source. Fits to the lightcurve disfavor the presence of Urca cooling pairs in the crust.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The link between coherent burst oscillations, burst spectral evolution and accretion state in 4U 1728-34

Coherent oscillations and the evolution of the X-ray spectrum during thermonuclear X-ray bursts in accreting neutron-star X-ray binaries have been studied intensively but separately. We analysed all the X-ray bursts of the source 4U 1728-34 with the Rossi X-ray Timing Explorer. We found that the presence of burst oscillations can be used to predict the behaviour of the blackbody radius during the cooling phase of the bursts. If a burst shows oscillations, during the cooling phase the blackbody radius remains more or less constant for ~2 - ~8s, whereas in bursts that do not show oscillations the blackbody radius either remains constant for more than ~2 - ~8s or it shows a rapid (faster than ~2s) decrease and increase. Both the presence of burst oscillations and the time-dependent spectral behaviour of the bursts are affected by accretion rate. We also found that the rise time and convexity of the bursts' light curve are different in bursts with and without oscillations in 4U 1728--34. Bursts with oscillations have a short rise time (~0.5s) and show both positive and negative convexity, whereas bursts without oscillations have a long rise time (~1s) and mostly positive convexity. This is consistent with the idea that burst oscillations are associated with off-equator ignition.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A signature of chemical separation in the cooling curves of transiently accreting neutron stars

We show that convection driven by chemical separation can significantly affect the cooling curves of accreting neutron stars after they go into quiescence. We calculate the thermal relaxation of the neutron star ocean and crust including the thermal and compositional fluxes due to convection. After the inward propagating cooling wave reaches the base of the neutron star ocean, the ocean begins to freeze, driving chemical separation. The resulting convection transports heat inward, giving much faster cooling of the surface layers than found assuming the ocean cools passively. The light curves including convection show a rapid drop in temperature weeks after outburst. Identifying this signature in observed cooling curves would constrain the temperature and composition of the ocean as well as offer a real time probe of the freezing of a classical multicomponent plasma.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Carbon Synthesis in Steady-State Hydrogen and Helium Burning On Accreting Neutron Stars

Superbursts from accreting neutron stars probe nuclear reactions at extreme densities ($ρ\approx 10^{9}~g\,cm^{-3}$) and temperatures ($T>10^9~K$). These bursts ($\sim$1000 times more energetic than type I X-ray bursts) are most likely triggered by unstable ignition of carbon in a sea of heavy nuclei made during the rp-process of regular type I X-ray bursts (where the accumulated hydrogen and helium are burned). An open question is the origin of sufficient amounts of carbon, which is largely destroyed during the rp-process in X-ray bursts. We explore carbon production in steady-state burning via the rp-process, which might occur together with unstable burning in systems showing superbursts. We find that for a wide range of accretion rates and accreted helium mass fractions large amounts of carbon are produced, even for systems that accrete solar composition. This makes stable hydrogen and helium burning a viable source of carbon to trigger superbursts. We also investigate the sensitivity of the results to nuclear reactions. We find that the $^{14}$O($α$,p)$^{17}$F reaction rate introduces by far the largest uncertainties in the $^{12}$C yield.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Constraining the initial entropy of directly-detected exoplanets

The post-formation, initial entropy S_i of a gas giant planet is a key witness to its mass-assembly history and a crucial quantity for its early evolution. However, formation models are not yet able to predict reliably S_i, making unjustified the use solely of traditional, 'hot-start' cooling tracks to interpret direct imaging results and calling for an observational determination of initial entropies to guide formation scenarios. Using a grid of models in mass and entropy, we show how to place joint constraints on the mass and initial entropy of an object from its observed luminosity and age. This generalises the usual estimate of only a lower bound on the real mass through hot-start tracks. Moreover, we demonstrate that with mass information, e.g. from dynamical stability analyses or radial velocity, tighter bounds can be set on the initial entropy. We apply this procedure to 2M1207 b and find that its initial entropy is at least 9.2 kB/baryon, assuming that it does not burn deuterium. For the planets of the HR 8799 system, we infer that they must have formed with S_i > 9.2 kB/baryon, independent of uncertainties about the age of the star. Finally, a similar analysis for beta Pic b reveals that it must have formed with S_i > 10.5 kB/baryon, using the radial-velocity mass upper limit. These initial entropy values are respectively ca. 0.7, 0.5, and 1.5 kB/baryon higher than the ones obtained from core accretion models by Marley et al., thereby quantitatively ruling out the coldest starts for these objects and constraining warm starts, especially for beta Pic b.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hall Attractor in Axially Symmetric Magnetic Fields in Neutron Star Crusts

We have found an attractor for an axially symmetric magnetic field evolving under the Hall effect and subdominant ohmic dissipation, resolving the question of the long term fate of the magnetic field in neutron star crusts. The electron fluid is in isorotation, analogous to Ferraro's law, with its angular velocity being approximately proportional to the poloidal magnetic flux, $Ω\propto Ψ$. This equilibrium is the long term configuration of a magnetic field evolving because of the Hall effect and ohmic dissipation. For an initial dipole dominated field the attractor consists mainly of a dipole and an octupole component accompanied by an energetically negligible quadrupole toroidal field. The field dissipates in a self-similar way: although higher multipoles should have been decaying faster, the toroidal field mediates transfer of energy into them from the lower ones, leading to an advection diffusion equilibrium and keeping the ratio of the poloidal multipoles almost constant. This has implications for the structure of the intermediate age neutron stars, suggesting that their poloidal field should consist of a dipole and a octupole component accompanied by a very weak toroidal quadrupole. For initial conditions that have a higher multipole $\ell$ structure the attractor consists mainly of $\ell$ and $\ell+2$ poloidal components.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Shedding Light on the Eccentricity Valley: Gap Heating and Eccentricity Excitation of Giant Planets in Protoplanetary Disks

We show that the first order (non co-orbital) corotation torques are significantly modified by entropy gradients in a non-barotropic protoplanetary disk. Such non-barotropic torques can dramatically alter the balance that, for barotropic cases, results in the net eccentricity damping for giant gap-clearing planets embedded in the disk. We demonstrate that stellar illumination can heat the gap enough for the planet's orbital eccentricity to instead be excited. We also discuss the "Eccentricity Valley" noted in the known exoplanet population, where low-metallicity stars have a deficit of eccentric planets between $\sim 0.1$ and $\sim 1$ AU compared to metal-rich systems (Dawson & Murray-Clay 2013). We show that this feature in the planet distribution may be due to the self-shadowing of the disk by a rim located at the dust sublimation radius $\sim 0.1$ AU, which is known to exist for several T Tauri systems. In the shadowed region between $\sim 0.1$ and $\sim 1$ AU lack of gap insolation allows disk interactions to damp eccentricity. Outside such shadowed regions stellar illumination can heat the planetary gaps and drive eccentricity growth for giant planets. We suggest that the self-shadowing does not arise at higher metallicity due to the increased optical depth of the gas interior to the dust sublimation radius.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Thermal Stability of Helium Burning on Accreting Neutron Stars

Thermonuclear burning on the surface of accreting neutron stars is observed to stabilize at accretion rates almost an order of magnitude lower than theoretical models predict. One way to resolve this discrepancy is by including a base heating flux that can stabilize the layer. We focus our attention on pure helium accretion, for which we calculate the effect of a base heating flux on the critical accretion rate at which thermonuclear burning stabilizes. We use the MESA stellar evolution code to calculate $\dot m_{\rm crit}$ as a function of the base flux, and derive analytic fitting formulae for $\dot m_{\rm crit}$ and the burning temperature at that critical accretion rate, based on a one-zone model. We also investigate whether the critical accretion rate can be determined by examining steady-state models only, without time-dependent simulations. We examine the argument that the stability boundary coincides with the turning point $dy_{\rm burn}/d\dot m=0$ in the steady-state models, and find that it does not hold outside of the one-zone, zero base flux case. A linear stability analysis of a large suite of steady-state models is also carried out, which yields critical accretion rates a factor of $\sim3$ larger than the MESA result, but with a similar dependence on base flux. Lastly, we discuss the implications of our results for the ultracompact X-ray binary 4U~1820-30.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Time-dependent, compositionally driven convection in the oceans of accreting neutron stars

We discuss the effect of chemical separation as matter freezes at the base of the ocean of an accreting neutron star, and the subsequent enrichment of the ocean in light elements and inward transport of heat through convective mixing. We extend the steady-state results of Medin & Cumming 2011 to transiently accreting neutron stars, by considering the time-dependent cases of heating during accretion outbursts and cooling during quiescence. Convective mixing is extremely efficient, flattening the composition profile in about one convective turnover time (weeks to months at the base of the ocean). During accretion outbursts, inward heat transport has only a small effect on the temperature profile in the outer layers until the ocean is strongly enriched in light elements, a process that takes hundreds of years to complete. During quiescence, however, inward heat transport rapidly cools the outer layers of the ocean while keeping the inner layers hot. We find that this leads to a sharp drop in surface emission at around a week followed by a gradual recovery as cooling becomes dominated by the crust. Such a dip should be observable in the light curves of these neutron star transients, if enough data is taken at a few days to a month after the end of accretion. If such a dip is definitively observed, it will provide strong constraints on the chemical composition of the ocean and outer crust.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Chandra observations of SGR 1627-41 near quiescence

We report on an observation of SGR 1627-41 made with the Chandra X-ray Observatory on 2011 June 16. Approximately three years after its outburst activity in 2008, the source's flux has been declining, as it approaches its quiescent state. For an assumed power-law spectrum, we find that the absorbed 2--10 keV flux for the source is $1.0^{+0.3}_{-0.2} \times 10^{-13} erg cm^{-2} s^{-1}$ with a photon index of $2.9 \pm 0.8$ ($N_H=1.0\times10^{23}$ cm^{-2}). This flux is approximately consistent with that measured at the same time after the source's outburst in 1998. With measurements spanning 3 years after the 2008 outburst, we analyze the long-term flux and spectral evolution of the source. The flux evolution is well described by a double exponential with decay times of 0.5 $\pm$ 0.1 and 59 $\pm$ 6 days, and a thermal cooling model fit suggests that SGR 1627-41 may have a hot core ($T_c ~ 2\times 10^8$ K). We find no clear correlation between flux and spectral hardness as found in other magnetars. We consider the quiescent X-ray luminosities of magnetars and the subset of rotation-powered pulsars with high magnetic fields ($B >~ 10^{13}$ G) in relation to their spin-inferred surface magnetic-field strength, and find a possible trend between the two quantities.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraints on neutron star mass and radius in GS 1826-24 from sub-Eddington X-ray bursts

We investigate the constraints on neutron star mass and radius in GS 1826-24 from models of lightcurves and spectral evolution of type I X-ray bursts. This source shows remarkable agreement with theoretical calculations of burst energies, recurrence times, and lightcurves. We first exploit this agreement to set the overall luminosity scale of the observed bursts. When combined with a measured blackbody normalization, this leads to a distance and anisotropy independent measurement of the ratio between the redshift 1+z and color correction factor f_c. We find 1+z=1.19-1.28 for f_c=1.4-1.5. We then compare the evolution of the blackbody normalization with flux in the cooling tail of bursts with predictions from spectral models of Suleimanov et al. (2011b). The observations are well described by the models at luminosities greater than about one third of the peak luminosity, with deviations emerging at luminosities below that. We show that this comparison leads to distance independent upper limits on R_\infty and neutron star mass of R_\infty\lesssim 9.0-13.2 km and M<1.2-1.7 M_\odot, respectively, for solar abundance of hydrogen at the photosphere and a range of metallicity and surface gravity. The radius limits are low in comparison to previous measurements. This may be indicative of a subsolar hydrogen fraction in the GS 1826-24 photosphere, or of larger color corrections than that predicted by spetral models. Our analysis also gives an upper limit on the distance to GS 1826-24 of d<4.0-5.5 kpc ξ_b^{-1/2}, where ξ_b is the degree of anisotropy of the burst emission.

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

Ohmic Dissipation in the Interiors of Hot Jupiters

We present models of ohmic heating in the interiors of hot jupiters in which we decouple the interior and the wind zone by replacing the wind zone with a boundary temperature Tiso and magnetic field Bphi0. Ohmic heating influences the contraction of gas giants in two ways: by direct heating within the convection zone, and by heating outside the convection zone which increases the effective insulation of the interior. We calculate these effects, and show that internal ohmic heating is only able to slow the contraction rate of a cooling gas giant once the planet reaches a critical value of internal entropy. We determine the age of the gas giant when ohmic heating becomes important as a function of mass, Tiso and induced Bphi0. With this survey of parameter space complete, we then adopt the wind zone scalings of Menou (2012) and calculate the expected evolution of gas giants with different levels of irradiation. We find that,with this prescription of magnetic drag, it is difficult to inflate massive planets or those with strong irradiation using ohmic heating, meaning that we are unable to account for many of the observed hot jupiter radii. This is in contrast to previous evolutionary models that assumed that a constant fraction of the irradiation is transformed into ohmic power.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Compositionally-driven convection in the oceans of accreting neutron stars

We discuss the effect of chemical separation as matter freezes at the base of the ocean of an accreting neutron star, and argue that the retention of light elements in the liquid acts as a source of buoyancy that drives a slow but continual mixing of the ocean, enriching it substantially in light elements, and leading to a relatively uniform composition with depth. We first consider the timescales associated with different processes that can redistribute elements in the ocean, including convection, sedimentation, crystallization, and diffusion. We then calculate the steady state structure of the ocean of a neutron star for an illustrative model in which the accreted hydrogen and helium burns to produce a mixture of O and Se. Even though the H/He burning produces only 2% oxygen by mass, the steady state ocean has an oxygen abundance more than ten times larger, almost 40% by mass. Furthermore, we show that the convective motions transport heat inwards, with a flux of ~ 0.2 MeV per nucleon for an O-Se ocean, heating the ocean and steepening the outwards temperature gradient. The enrichment of light elements and heating of the ocean due to compositionally-driven convection likely have important implications for carbon ignition models of superbursts.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Continued cooling of the crust in the neutron star low-mass X-ray binary KS 1731-260

Some neutron star low-mass X-ray binaries have very long outbursts (lasting several years) which can generate a significant amount of heat in the neutron star crust. After the system has returned to quiescence, the crust then thermally relaxes. This provides a rare opportunity to study the thermal properties of neutron star crusts, putting constraints on the thermal conductivity and hence the structure and composition of the crust. KS 1731-260 is one of only four systems where this crustal cooling has been observed. Here, we present a new Chandra observation of this source approximately 8 years after the end of the last outburst, and 4 years since the last observation. We find that the source has continued to cool, with the cooling curve displaying a simple power-law decay. This suggests that the crust has not fully thermally relaxed yet, and may continue to cool further. A simple power law decay is in contrast to theoretical cooling models of the crust, which predict that the crust should now have cooled to the same temperature as the neutron star core.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Crystallization of classical multi-component plasmas

We develop a method for calculating the equilibrium properties of the liquid-solid phase transition in a classical, ideal, multi-component plasma. Our method is a semi-analytic calculation that relies on extending the accurate fitting formulae available for the one-, two-, and three-component plasmas to the case of a plasma with an arbitrary number of components. We compare our results to those of Horowitz, Berry, & Brown (Phys. Rev. E, 75, 066101, 2007), who use a molecular dynamics simulation to study the chemical properties of a 17-species mixture relevant to the ocean-crust boundary of an accreting neutron star, at the point where half the mixture has solidified. Given the same initial composition as Horowitz et al., we are able to reproduce to good accuracy both the liquid and solid compositions at the half-freezing point; we find abundances for most species within 10% of the simulation values. Our method allows the phase diagram of complex mixtures to be explored more thoroughly than possible with numerical simulations. We briefly discuss the implications for the nature of the liquid-solid boundary in accreting neutron stars.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Radiative Hydrodynamic Simulations of HD209458b: Temporal Variability

We present a new approach for simulating the atmospheric dynamics of the close-in giant planet HD209458b that allows for the decoupling of radiative and thermal energies, direct stellar heating of the interior, and the solution of the full 3D Navier Stokes equations. Simulations reveal two distinct temperature inversions (increasing temperature with decreasing pressure) at the sub-stellar point due to the combined effects of opacity and dynamical flow structure and exhibit instabilities leading to changing velocities and temperatures on the nightside for a range of viscosities. Imposed on the quasi-static background, temperature variations of up to 15% are seen near the terminators and the location of the coldest spot is seen to vary by more than 20 degrees, occasionally appearing west of the anti-solar point. Our new approach introduces four major improvements to our previous methods including simultaneously solving both the thermal energy and radiative equations in both the optical and infrared, incorporating updated opacities, including a more accurate treatment of stellar energy deposition that incorporates the opacity relevant for higher energy stellar photons, and the addition of explicit turbulent viscosity.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

An Integrated Analysis of Radial Velocities in Planet Searches

We discuss a Bayesian approach to the analysis of radial velocities in planet searches. We use a combination of exact and approximate analytic and numerical techniques to efficiently evaluate chi-squared for multiple values of orbital parameters, and to carry out the marginalization integrals for a single planet including the possibility of a long term trend. The result is a robust algorithm that is rapid enough for use in real time analysis that outputs constraints on orbital parameters and false alarm probabilities for the planet and long term trend. The constraints on parameters and odds ratio that we derive compare well with previous calculations based on Markov Chain Monte Carlo methods, and we compare our results with other techniques for estimating false alarm probabilities and errors in derived orbital parameters. False alarm probabilities from the Bayesian analysis are systematically higher than frequentist false alarm probabilities, due to the different accounting of the number of trials. We show that upper limits on the velocity amplitude derived for circular orbits are a good estimate of the upper limit on the amplitude of eccentric orbits for eccentricities less than about 0.5.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Mapping crustal heating with the cooling lightcurves of quasi-persistent transients

The monitoring of quiescent emission from neutron star transients with accretion outbursts long enough to significantly heat the neutron star crust has opened a new vista onto the physics of dense matter. In this paper we construct models of the thermal relaxation of the neutron star crust following the end of a protracted accretion outburst. We confirm the finding of Shternin et al., that the thermal conductivity of the neutron star crust is high, consistent with a low impurity parameter. We describe the basic physics that sets the broken power-law form of the cooling lightcurve. The initial power law decay gives a direct measure of the temperature profile, and hence the thermal flux during outburst, in the outer crust. The time of the break, at hundreds of days post-outburst, corresponds to the thermal time where the solid transitions from a classical to quantum crystal, close to neutron drip. We calculate in detail the constraints on the crust parameters of both KS 1731-260 and MXB 1659-29 from fitting their cooling lightcurves. Our fits to the lightcurves require that the neutrons do not contribute significantly to the heat capacity in the inner crust, and provide evidence in favor of the existence of a neutron superfluid throughout the inner crust. Our fits to both sources indicate an impurity parameter of order unity in the inner crust.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

An Accurate Determination of the Optical Periodic Modulation in the X-Ray Binary SAX J1808.4-3658

We report on optical imaging of the X-ray binary SAX J1808.4-3658 with the 8-m Gemini South Telescope. The binary, containing an accretion-powered millisecond pulsar, appears to have a large periodic modulation in its quiescent optical emission. In order to clarify the origin of this modulation, we obtained three time-resolved $r'$-band light curves (LCs) of the source in five days. The LCs can be described by a sinusoid, and the long time-span between them allows us to determine optical period P=7251.9 s and phase 0.671 at MJD 54599.0 (TDB; phase 0.0 corresponds to the ascending node of the pulsar orbit), with uncertainties of 2.8 s and 0.008 (90 % confidence), respectively. This periodicity is highly consistent with the X-ray orbital ephemeris. By considering this consistency and the sinusoidal shape of the LCs, we rule out the possibility of the modulation arising from the accretion disk. Our study supports the previous suggestion that the X-ray pulsar becomes rotationally powered in quiescence, with its energy output irradiating the companion star, causing the optical modulation. While it has also been suggested that the accretion disk would be evaporated by the pulsar, we argue that the disk exists and gives rise to the persistent optical emission. The existence of the disk can be verified by long-term, multi-wavelength optical monitoring of the source in quiescence, as an increasing flux and spectral changes from the source would be expected based on the standard disk instability model.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The Keck Planet Search: Detectability and the Minimum Mass and Orbital Period Distribution of Extrasolar Planets

We analyze 8 years of precise radial velocity measurements from the Keck Planet Search, characterizing the detection threshold, selection effects, and completeness of the survey. We carry out a systematic search for planets by assessing the false alarm probability associated with Keplerian orbit fits to the data. This allows us to understand the detection threshold for each star in terms of the number and time baseline of the observations, and size of measurement errors and stellar jitter. We show that all planets with orbital periods <2000 days, velocity amplitudes >20 m/s, and eccentricities <0.6 have been announced, and summarize the candidates at lower amplitudes and longer orbital periods. For the remaining stars, we calculate upper limits on the velocity amplitude of a companion, typically 10 m/s, and use the non-detections to derive completeness corrections at low amplitudes and long orbital periods. We give the fraction of stars with a planet as a function of planet mass and orbital period, and extrapolate to long period orbits and low planet masses. A power law fit for planet masses >0.3 Jupiter masses and periods <2000 days gives a mass-period distribution dN=C M^αP^βdlnM dlnP with α=-0.31 \pm 0.2, β=0.26\pm 0.1, and the normalization constant C such that 10.5% of solar type stars have a planet with mass in the range 0.3-10 Jupiter masses and orbital period 2-2000 days. The orbital period distribution shows an increase in the planet fraction by a factor of 5 for orbital periods beyond 300 days. Extrapolation gives 17-20% of stars having gas giant planets within 20 AU. Finally, taking into account differences in detectability, we find that M dwarfs are 3 to 10 times less likely to harbor a Jupiter mass planet than solar type stars.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Helium-rich thermonuclear bursts and the distance to the accretion-powered millisecond pulsar SAX J1808.4-3658

We analysed Rossi X-ray Timing Explorer observations of the accretion-powered 401 Hz pulsar SAX J1808.4-3658, in order to precisely determine the source distance. While the fluences for the five transient outbursts observed from 1996 were constant to within the uncertainties, the outburst interval varied signficantly, so that the time-averaged flux (and accretion rate) decreased by around 40%. By equating the time-averaged X-ray flux with the expected mass transfer rate from gravitational radiation, we derived a lower limit on the distance of 3.4 kpc. Combined with an upper limit from assuming that the four radius-expansion thermonuclear bursts observed during the 2002 October outburst reached at most the Eddington limit for a pure He atmosphere, we found that the probable distance range for the source is 3.4-3.6 kpc. The implied inclination, based on the optical/IR properties of the counterpart, is i<~30 degrees. We compared the properties of the bursts with an ignition model. The time between bursts was long enough for hot CNO burning to significantly deplete the accreted hydrogen, so that ignition occurred in a pure helium layer underlying a stable hydrogen burning shell. This is the first time that this burning regime has been securely observationally identified. The observed energetics of the bursts give a mean hydrogen fraction at ignition of <X> approx. 0.1, and require that the accreted hydrogen fraction X_0 and the CNO metallicity Z_CNO are related by Z_CNO approx. 0.03(X_0/0.7)^2. We show that in this burning regime, a measurement of the burst recurrence time and energetics allows the local accretion rate onto the star to be determined independently of the accreted composition, giving a new method for estimating the source distance which is in good agreement with our other estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Millihertz Quasi-Periodic Oscillations from Marginally Stable Nuclear Burning on an Accreting Neutron Star

We investigate marginally stable nuclear burning on the surface of accreting neutron stars as an explanation for the mHz quasi-periodic oscillations (QPOs) observed from three low mass X-ray binaries. At the boundary between unstable and stable burning, the temperature dependence of the nuclear heating rate and cooling rate almost cancel. The result is an oscillatory mode of burning, with an oscillation period close to the geometric mean of the thermal and accretion timescales for the burning layer. We describe a simple one-zone model which illustrates this basic physics, and then present detailed multizone hydrodynamical calculations of nuclear burning close to the stability boundary using the KEPLER code. Our models naturally explain the characteristic 2 minute period of the mHz QPOs, and why they are seen only in a very narrow range of X-ray luminosities. The oscillation period is sensitive to the accreted hydrogen fraction and the surface gravity, suggesting a new way to probe these parameters. A major puzzle is that the accretion rate at which the oscillations appear in the theoretical models is an order of magnitude larger than the rate implied by the X-ray luminosity when the mHz QPOs are seen. We discuss the implications for our general understanding of nuclear burning on accreting neutron stars. One possibility is that the accreted material covers only part of the neutron star surface at luminosities Lx > ~1E37 erg/s.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

The Thermal Evolution following a Superburst on an Accreting Neutron Star

Superbursts are very energetic Type I X-ray bursts discovered in recent years by long term monitoring of X-ray bursters, and believed to be due to unstable ignition of carbon in the deep ocean of the neutron star. In this Letter, we follow the thermal evolution of the surface layers as they cool following the burst. The resulting lightcurves agree very well with observations for layer masses and energy releases in the range expected from ignition calculations. At late times, the cooling flux from the layer decays as a power law in time, giving timescales for quenching of normal Type I bursting of weeks, in good agreement with observational limits. We show that simultaneous modelling of superburst lightcurves and quenching times promises to constrain both the thickness of the fuel layer and the energy deposited.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

The Lick Planet Search: Detectability and Mass Thresholds

We analyse 11 years of precise radial velocities for 76 solar type stars from the Lick survey. Eight stars in this sample have previously reported planetary-mass companions, all with mass (m sin i) less than 8 Jupiter masses (MJ). For the stars without a detected companion, we place upper limits on possible companion mass. For most stars, we can exclude companions with m sin i > 0.7 MJ (a/AU)^1/2 for orbital radii a < 5 AU. We use our results to interpret the observed masses and orbital radii of planetary-mass companions. For example, we show that the finite duration of the observations makes detection of Jupiter mass companions more and more difficult for orbital radii beyond 3 AU. Thus it is possible that the majority of solar type stars harbor Jupiter-mass companions much like our own, and if so these companions should be detectable in a few years. To search for periodicities, we adopt a "floating-mean" periodogram, which improves on the traditional Lomb-Scargle periodogram by accounting for statistical fluctuations in the mean of a sampled sinusoid. We discuss in detail the normalization of the periodogram, an issue which has been of some debate in the literature.

0 citations0 reviewsNo code linkedOpen access