No Arabic abstract
An accurate knowledge of the neutron capture cross sections of 62,63Ni is crucial since both isotopes take key positions which affect the whole reaction flow in the weak s process up to A=90. No experimental value for the 63Ni(n,gamma) cross section exists so far, and until recently the experimental values for 62Ni(n,gamma) at stellar temperatures (kT=30 keV) ranged between 12 and 37 mb. This latter discrepancy could now be solved by two activations with following AMS using the GAMS setup at the Munich tandem accelerator which are also in perfect agreement with a recent time-of-flight measurement. The resulting (preliminary) Maxwellian cross section at kT=30 keV was determined to be <sigma>30keV = 23.4 +/- 4.6 mb. Additionally, we have measured the 64Ni(gamma,n)63Ni cross section close to threshold. Photoactivations at 13.5 MeV, 11.4 MeV and 10.3 MeV were carried out with the ELBE accelerator at Forschungszentrum Dresden-Rossendorf. A first AMS measurement of the sample activated at 13.5 MeV revealed a cross section smaller by more than a factor of 2 compared to NON-SMOKER predictions.
The stellar (n,gamma) cross section of 40Ca at kT=25 keV has been measured with a combination of the activation technique and accelerator mass spectrometry (AMS). This combination is required when direct off-line counting of the produced activity is compromised by the long half-life and/or missing gamma-ray transitions. The neutron activations were performed at the Karlsruhe Van de Graaff accelerator using the quasistellar neutron spectrum of kT=25 keV produced by the 7Li(p,n)7Be reaction. The subsequent AMS measurements were carried out at the Vienna Environmental Research Accelerator (VERA) with a 3 MV tandem accelerator. The doubly magic 40Ca is a bottle-neck isotope in incomplete silicon burning, and its neutron capture cross section determines the amount of leakage, thus impacting on the eventual production of iron group elements. Because of its high abundance, 40Ca can also play a secondary role as neutron poison for the s-process. Previous determinations of this value at stellar energies were based on time-of-flight measurements. Our method uses an independent approach, and yields for the Maxwellian-averaged cross section at kT=30 keV a value of <sigma>30 keV= 5.73+/-0.34 mb.
Thanks to missions like Kepler and TESS, we now have access to tens of thousands of high precision, fast cadence, and long baseline stellar photometric observations. In principle, these light curves encode a vast amount of information about stellar variability and, in particular, about the distribution of starspots and other features on their surfaces. Unfortunately, the problem of inferring stellar surface properties from a rotational light curve is famously ill-posed, as it often does not admit a unique solution. Inference about the number, size, contrast, and location of spots can therefore depend very strongly on the assumptions of the model, the regularization scheme, or the prior. The goal of this paper is twofold: (1) to explore the various degeneracies affecting the stellar light curve inversion problem and their effect on what can and cannot be learned from a stellar surface given unresolved photometric measurements; and (2) to motivate ensemble analyses of the light curves of many stars at once as a powerful data-driven alternative to common priors adopted in the literature. We further derive novel results on the dependence of the null space on stellar inclination and limb darkening and show that single-band photometric measurements cannot uniquely constrain quantities like the total spot coverage without the use of strong priors. This is the first in a series of papers devoted to the development of novel algorithms and tools for the analysis of stellar light curves and spectral time series, with the explicit goal of enabling statistically robust inference about their surface properties.
We present new methods to solve the Riemann problem both exactly and approximately for general equations of state (EoS) to facilitate realistic modeling and understanding of astrophysical flows. The existence and uniqueness of the new exact general EoS Riemann solution can be guaranteed if the EoS is monotone regardless of the physical validity of the EoS. We confirm that: (1) the solution of the new exact general EoS Riemann solver and the solution of the original exact Riemann solver match when calculating perfect gas Euler equations; (2) the solution of the new Harten-Lax-van Leer-Contact (HLLC) general EoS Riemann solver and the solution of the original HLLC Riemann solver match when working with perfect gas EoS; and (3) the solution of the new HLLC general EoS Riemann solver approaches the new exact solution. We solve the EoS with two methods, one is to interpolate 2D EoS tables by the bi-linear interpolation method, and the other is to analytically calculate thermodynamic variables at run-time. The interpolation method is more general as it can work with other monotone and realistic EoS while the analytic EoS solver introduced here works with a relatively idealized EoS. Numerical results confirm that the accuracy of the two EoS solvers is similar. We study the efficiency of these two methods with the HLLC general EoS Riemann solver and find that analytic EoS solver is faster in the test problems. However, we point out that a combination of the two EoS solvers may become favorable in some specific problems. Throughout this research, we assume local thermal equilibrium.
Occurring in protoplanetary discs composed of dust and gas, streaming instabilities are a favoured mechanism to drive the formation of planetesimals. The Polydispserse Streaming Instability is a generalisation of the Streaming Instability to a continuum of dust sizes. This second paper in the series provides a more in-depth derivation of the governing equations and presents novel numerical methods for solving the associated linear stability problem. In addition to the direct discretisation of the eigenproblem at second order introduced in the previous paper, a new technique based on numerically reducing the system of integral equations to a complex polynomial combined with root finding is found to yield accurate results at much lower computational cost. A related method for counting roots of the dispersion relation inside a contour without locating those roots is also demonstrated. Applications of these methods show they can reproduce and exceed the accuracy of previous results in the literature, and new benchmark results are provided. Implementations of the methods described are made available in an accompanying Python package psitools.
The 62Ni(n,gamma)63Ni(t_1/2=100+-2 yrs) reaction plays an important role in the control of the flow path of the slow neutron-capture (s-) nucleosynthesis process. We have measured for the first time the total cross section of this reaction for a quasi-Maxwellian (kT = 25 keV) neutron flux. The measurement was performed by fast-neutron activation, combined with accelerator mass spectrometry to detect directly the 63Ni product nuclei. The experimental value of 28.4+-2.8 mb, fairly consistent with a recent theoretical estimate, affects the calculated net yield of 62Ni itself and the whole distribution of nuclei with 62<A <90 produced by the weak s-process in massive stars.