scholarly journals Overshooting calibration and age determination from evolved binary systems

2018 ◽  
Vol 615 ◽  
pp. A62 ◽  
Author(s):  
G. Valle ◽  
M. Dell’Omodarme ◽  
P. G. Prada Moroni ◽  
S. Degl’Innocenti

Aims. The capability of grid-based techniques to estimate the age together with the convective core overshooting efficiency of stars in detached eclipsing binary systems for main sequence stars has previously been investigated. We have extended this investigation to later evolutionary stages and have evaluated the bias and variability on the recovered age and convective core overshooting parameter accounting for both observational and internal uncertainties. Methods. We considered synthetic binary systems, whose age and overshooting efficiency should be recovered by applying the SCEPtER pipeline to the same grid of models used to build the mock stars. We focus our attention on a binary system composed of a 2.50 M⊙ primary star coupled with a 2.38 M⊙ secondary. To explore different evolutionary scenarios, we performed the estimation at three different times: when the primary is at the end of the central helium burning, when it is at the bottom of the RGB, and when it is in the helium core burning phase. The Monte Carlo simulations have been carried out for two typical values of accuracy on the mass determination, that is, 1% and 0.1%. Results. Adopting typical observational uncertainties, we found that the recovered age and overshooting efficiency are biased towards low values in all three scenarios. For an uncertainty on the masses of 1%, the underestimation is particularly relevant for a primary in the central helium burning stage, reaching − 8.5% in age and − 0.04 (− 25% relative error) in the overshooting parameter β. In the other scenarios, an undervaluation of the age by about 4% occurs. A large variability in the fitted values between Monte Carlo simulations was found: for an individual system calibration, the value of the overshooting parameter can vary from β = 0.0 to β = 0.26. When adopting a 0.1% error on the masses, the biases remain nearly unchanged but the global variability is suppressed by a factor of about two. We also explored the effect of a systematic discrepancy between the artificial systems and the model grid by accounting for an offset in the effective temperature of the stars by ± 150 K. For a mass error of 1% the overshooting parameter is largely biased towards the edges of the explored range, while for the lower mass uncertainty it is basically unconstrained from 0.0 to 0.2. We also evaluate the possibility of individually recovering the β value for both binary stars. We found that this is impossible for a primary near to central hydrogen exhaustion owing to huge biases for the primary star of + 0.14 (90% relative error), while in the other cases the fitted β are consistent, but always biased by about − 0.04 (− 25% relative error). Finally, the possibility to distinguish between models computed with mild overshooting from models with no overshooting was evaluated, resulting in a reassuring power of distinction greater than 80%. However, the scenario with a primary in the central helium burning was a notable exception, showing a power of distinction lower than 5%.

2021 ◽  
Author(s):  
Masahide Sato

Abstract Performing isothermal-isochoric Monte Carlo simulations, I examine the types of clusters that dumbbell-like one–patch particles form in thin space between two parallel walls, assuming that each particle is synthesized through the merging of two particles, one non-attracting and the other attracting for which, for example, the inter-particle interaction is approximated by the DLVO model. The shape of these dumbbell-like particles is controlled by the ratio of the diameters q of the two spherical particles and by the dimensionless distance l between them. Using a modified Kern–Frenkel potential, I examine the dependence of the cluster shape on l and q. Large island-like clusters are created when q < 1. With increasing q, the clusters become chain-like. When q increases further, elongated clusters and regular polygonal clusters are created. In hte simulations, the cluster shape becomes three-dimensional with increasing l because the thickness of the thin system increases proportionally to l.


1976 ◽  
Vol 73 ◽  
pp. 27-34 ◽  
Author(s):  
C. De Loore ◽  
J. P. De Greve

It is well known that the outcome of case B evolution of the primaries of massive close binary systems (M1 ≥ 9 M⊙) depends on the initial primary mass. The most massive primaries finally ignite carbon, form iron cores and presumably end in a supernova explosion, whereas the lighter ones presumably end as white dwarfs, without carbon ignition. This paper derives an estimate of the mass boundary separating these two kinds of evolution.As an example of the first case, the evolution of a 20 M⊙ + 14 M⊙ system was computed; after the mass exchange, the primary star (with M = 5.43 M⊙) evolves through the helium-burning (Wolf-Rayet) stage towards a supernova explosion; finally the system evolves into an X-ray binary (BWRX-evolution).As a representative for the second case the evolution of a 10 M⊙ + 8 M⊙ system was examined. After the first stage of mass exchange, the primary (with a mass of 1.66 M⊙) approaches the helium main sequence; during later phases of helium burning the radius increases again, and a second stage of mass transfer starts; after this the star (with a mass of 1.14 M⊙) again evolves towards the left in the Hertzsprung-Russell diagram and ends as a white dwarf (BSWD-evolution). A system of 15 M⊙ + 8 M⊙ is found to evolve very similar to the 20 M⊙ + 14 M⊙ system. The mass Mu, separating the two types of evolution, must therefore be situated between 10 and 15 solar masses. An initial chemical composition X = 0.70, Z = 0.03 was used for all systems.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Masahide Sato

AbstractPerforming isothermal-isochoric Monte Carlo simulations, I examine the types of clusters that dumbbell-like one–patch particles form in thin space between two parallel walls, assuming that each particle is synthesized through the merging of two particles, one non-attracting and the other attracting for which, for example, the inter-particle interaction is approximated by the DLVO model . The shape of these dumbbell-like particles is controlled by the ratio of the diameters q of the two spherical particles and by the dimensionless distance l between these centers. Using a modified Kern–Frenkel potential, I examine the dependence of the cluster shape on l and q. Large island-like clusters are created when $$q<1$$ q < 1 . With increasing q, the clusters become chain-like . When q increases further, elongated clusters and regular polygonal clusters are created. In the simulations, the cluster shape becomes three-dimensional with increasing l because the thickness of the thin system increases proportionally to l.


Author(s):  
Binyamin V Naiman ◽  
Efrat Sabach ◽  
Avishai Gilkis ◽  
Noam Soker

Abstract We simulate the evolution of binary systems with a massive primary star of 15M⊙ where we introduce an enhanced mass loss due to jets that the secondary star might launch, and find that in many cases the enhanced mass loss brings the binary system to experience the grazing envelope evolution (GEE) and form a progenitor of Type IIb supernova (SN IIb). The jets, the Roche lobe overflow (RLOF), and a final stellar wind remove most of the hydrogen-rich envelope, leaving a blue-compact SN IIb progenitor. In many cases without this jet-driven mass loss the system enters a common envelope evolution (CEE) and does not form a SN IIb progenitor. We use the stellar evolutionary code MESA binary and mimic the jet-driven mass loss with a simple prescription and some free parameters. Our results show that the jet-driven mass loss, that some systems have during the GEE, increases the parameter space for stellar binary systems to form SN IIb progenitors. We estimate that the binary evolution channel with GEE contributes about a quarter of all SNe IIb, about equal to the contribution of each of the other three channels, binary evolution without a GEE, fatal CEE (where the secondary star merges with the core of the giant primary star), and the single star channel.


2011 ◽  
Vol 172-174 ◽  
pp. 658-663 ◽  
Author(s):  
Mohamed Briki ◽  
Jérôme Creuze ◽  
Fabienne Berthier ◽  
Bernard Legrand

In order to build the phase diagram of Cu-Ag nanoalloys, we study a 405-atom nanoparticle by means of Monte Carlo simulations with relaxations usingN-body interatomic potentials. We focus on a range of nominal concentrations for which the cluster core remains Cu-pure and the (001) facets of the outer shell exhibit two original phenomena. Within the (N,mAg-mCu,P,T) ensemble, a structural and chemical bistability is observed, which affects all the (001) facets together. For a nanoparticle assembly, this will result in a bimodal distribution of clusters, some of them having their (001) facets Cu-rich with the usual square shape, the other ones having their (001) facets Ag-rich with a diamond shape. This bistability is replaced in the (NAg,NCu,P,T) ensemble by a continuous evolution of both the structure and the concentration of the (001) facets from Cu-rich square-shaped to Ag-rich diamond-shaped facets as the number of Ag atoms increases. For a nanoparticle assembly, this will result in an unimodal distribution of the cluster population concerning the properties of the (001) facets. This comparison between pseudo grand canonical and isothermal-isobaric results shows that the distribution of a population of bimetallic nanoparticles depends strongly on the conditions under it is elaborated.


1998 ◽  
Vol 14 (2) ◽  
pp. 200-221 ◽  
Author(s):  
Jörg Breitung ◽  
Philip Hans Franses

In this paper we consider a semiparametric version of the test for seasonal unit roots suggested by Hylleberg, Engle, Granger, and Yoo (1990, Journal of Econometrics 44, 215–238). The asymptotic theory is based on the analysis of a simple regression problem, and the results apply to tests at any given frequency in the range (0,π]. Monte Carlo simulations suggest that the test may have more power than the parametric test of Hylleberg et al. (1990). On the other hand, the semiparametric version suffers from severe size distortions in some situations.


2006 ◽  
Vol 17 (06) ◽  
pp. 785-793 ◽  
Author(s):  
F. W. S. LIMA

On directed Barabási–Albert networks with two and seven neighbours selected by each added site, the Ising model does not seem to show a spontaneous magnetisation. Instead, the decay time for flipping of the magnetisation follows an Arrhenius law for Metropolis and Glauber algorithms, but for Wolff cluster flipping the magnetisation decays exponentially with time. On these networks the magnetisation behaviour of the Ising model, with Glauber, HeatBath, Metropolis, Wolf or Swendsen–Wang algorithm competing against Kawasaki dynamics, is studied by Monte Carlo simulations. We show that the model exhibits the phenomenon of self-organisation (= stationary equilibrium) defined in Ref. 8 when Kawasaki dynamics is not dominant in its competition with Glauber, HeatBath and Swendsen–Wang algorithms. Only for Wolff cluster flipping the magnetisation, this phenomenon occurs after an exponentially decay of magnetisation with time. The Metropolis results are independent of competition. We also study the same process of competition described above but with Kawasaki dynamics at the same temperature as the other algorithms. The obtained results are similar for Wolff cluster flipping, Metropolis and Swendsen–Wang algorithms but different for HeatBath.


2021 ◽  
Vol 922 (2) ◽  
pp. 105
Author(s):  
Yong Yang ◽  
Jingkun Zhao ◽  
Jiajun Zhang ◽  
Xianhao Ye ◽  
Gang Zhao

Abstract We revisit the moving groups (MGs) in the solar neighborhood with a sample of 91,969 nearby stars constructed from LAMOST DR7. Using the wavelet technique and Monte Carlo simulations, five MGs together with a new candidate located at V≃−130 km s−1 are detected simultaneously in V − U 2 + 2 V 2 space. Taking into account the other known MGs, we conclude that MGs in the Galactic disk are spaced by approximately 15–25 km s−1 along V velocity. The origin of detected MGs is analyzed through the distributions of [Fe/H]−[Mg/Fe] and ages. Our results support attributing the origin to the continuous resonant mechanisms probably induced by the bar or spiral arms of the Milky Way.


2021 ◽  
Author(s):  
Elyes Balti

Reed-Solomon and related codes have recently become very important for erasure correction in large disk arrays used in data centers. In this paper, we will implement a 3-error correcting Reed-Solomon encoder and decoder over the field GF(9) generated by the primitive polynomial D^2 + D + 2 over GF(3) and the decoding is carried out by the Berlekamp. We simulate the encoder and decoder using Monte-Carlo simulations over the 9-ary symmetric channel that outputs the correct symbol with probability (1-p), and outputs one of the other 8 possible incorrect symbols with probability p/8. Then, we compare the simulated probability of symbol error P(E) of out code with the union upper bound.


2013 ◽  
Vol 23 ◽  
pp. 157-160
Author(s):  
ZHENG CHENG ◽  
CHENGMIN ZHANG ◽  
ALI TAANI

We perform a Monte Carlo simulation to investigate the distribution and evolution of 66 mass measured pulsars. We get the best fits to the distribution at 1.35 ± 0.27M⨀ (1σ confidence level). In addition, we notice bimodal distributions in 1.34 ± 0.15M⨀ and 1.48 ± 0.53M⨀, this can be led to the idea that radio pulsars in binary systems have recycled. Thus we divide the data according to the characteristic spin period into two groups, millisecond pulsars (MSPs), P Spin ≤ 20 ms and less recycled pulsars P Spin ≥ 20 ms , respectively. We show that the distributions of MSPs at 1.42 ± 0.36M⨀, and 1.32 ± 0.18M⨀ for less recycled pulsars. As such, the mass of MSPs are heavier than those in less recycled pulsars by ~ 0.1M⨀, since they accreting material from their companions. On the other hand, the formation of heavier pulsars from the accretion induced collapse of accreting white dwarfs, must be invoked.


Sign in / Sign up

Export Citation Format

Share Document