Heavy element contributions of rotating massive stars to interstellar medium

Employing the stellar evolution code Modules for Experiments in Stellar Astrophysics(MESA),we calculate yields of heavy elements from massive stars via stellar wind and core−collapse supernova(CCSN)ejecta to the interstellar medium(ISM).In our models,the initial masses(Mini)of massive stars are taken from 13 to 80 M⊙,their initial rotational velocities(V)are 0,300 and 500 km s^(−1),and their metallicities are[Fe/H]=−3,−2,−1 and 0.The yields of heavy elements coming from stellar winds are mainly affected by stellar rotation which changes the chemical abundances of stellar surfaces via chemically homogeneous evolution,and enhances mass-loss rate.We estimate that the stellar wind can produce heavy element yields of about 10^(−2)(for low metallicity models)to a mass of severalM⊙(for low metallicity and rapid rotation models).The yields of heavy elements produced by CCSN ejecta also depend on the large amount of remnantmass which is mainly determined by the mass of the CO-core.Our models calculate that the yields of heavy elements produced by CCSN ejecta can get up to several M⊙.Compared with stellar wind,CCSN ejecta has a greater contribution to the heavy elements in ISM.We also compare the 56Ni yields calculated in this work with the observational estimate.Our models only explain the 56Ni masses produced by faint SNe or normal SNe with progenitor mass lower than about 25 M⊙,and greatly underestimate the 56Ni masses produced by stars with masses higher than about 30M⊙.

An Asymmetric Galactic Stellar Disk Traced by OB-type Stars from LAMOST DR7

Using 9943 OB-type stars from LAMOST DR7 in the solar neighborhood,we fit the vertical stellar density profile with the model including a single exponential distribution at different positions(R,Φ).The distributions of the scale heights and scale length show that the young disk traced by the OB-type stars is not axisymmetric.The scale length decreases versus the azimuthal angleΦ,i.e.,from.■kpc withΦ=-3°to■kpc withΦ=9°.Meanwhile we find signal of non-symmetry in the distribution of the scale height of the north and south of the disk plane.The scale height in the north side shows signal of flaring of the disk,while that of the south disk stays almost constant around h_(s)=130 pc.The distribution of the displaceeent of the disk plane Z_(0)also shows variance versus the azimuthal angleΦ,which displays significant differences with the warp model constrained by the Cepheid stars.We also test different values for the position of the Sun,and the distance between the Sun and the Galactic center affects the scale heights and the displacement of the disk significantly,but that does not change our conclusion that the disk is not axisymmetric.

Spectroscopic Observations of Ten Galactic Wolf-Rayet Stars at Bosscha Observatory:Determination of Stellar Parameters and Mass-loss Rates

We present optical spectra of 10 Galactic Wolf-Rayet(WR)stars that consist of five WN and five WC stars.The optical observation was conducted using a low-resolution spectrograph NEO-R1000(λ/Δλ~1000)at GAO-ITB RTS(27.94 cm,F/10.0),Bosscha Observatory,Lembang.We implemented stellar atmosphere Postdam Wolf-Rayet(PoWR)grid modeling to derive stellar parameters.The normalized optical spectrum can be used to find the best model from the available PoWR grid,then we could derive stellar temperature and transformation radius.To derive luminosity,stellar radius and color excess,we conducted a Spectral Energy Distribution(SED)analysis with additional data on the near-ultraviolet spectrum from the International Ultraviolet Explorer(IUE)database,and UBV and 2MASS JHK broadband filter data.Additional analysis to derive asymptotic terminal wind velocity was conducted from the P-Cygni profile analysis of the high-resolution IUE ultraviolet spectrum.With previously derived parameters,we could determine the mass loss rate of the WR stars.Furthermore,we compared our results with previous work that used PoWR code and the differences are not more than 20%.We conclude that the PoWR spectral grid is sufficient to derive WR stellar parameters quickly and could provide more accurate initial parameter input to the PoWR program code.

A Note on the Mixing Length Theory and Massive Star Evolution

In this paper, we investigate the problem of supersonic convection caused by the application of the usual Mixing Length Theory (MLT), and give a modification to the original expression of the MLT. In the case of the usual MLT, a zone of supersonic convection exists when the stellar model goes into yellow-red supergiant phase. When the modification is applied, the models of late type supergiants possess shallower convective envelopes for a given temperature compared with the normal formalism of the MLT. Therefore the stellar models made with our new formalism move to lower effective temperature by 100-400 K depending on the luminosity of the star. Such a modification does not affect the convective envelope structure of lower luminosity late type stars, as both expressions of the MLT give the same result. In these stars.

New Massive Contact Twin Binary in a Radio-quiet HⅡRegion Associated with the M17 Complex

Early-B stars,much less energetic than O stars,may create an HⅡregion that appears as radio-quiet.We report the identification of new early-B stars associated with the radio-quiet HⅡregion G014.645–00.606 in the M17 complex.The radio-quiet HⅡregion G014.645–00.606 is adjacent to three radio-quiet WISE HⅡregion candidates.The ionizing sources of the radio-quiet HⅡregions are expected to be later than B1 V,given the sensitivity about 1–2 m Jy of the MAGPIS 20 cm survey.The stars were first selected if their parallaxes of Gaia EDR3 match that of the 22 GHz H2 O maser source within the same region.We used the color–magnitude diagram made from the Zwicky Transient Facility photometric catalog to select the candidates for massive stars because the intrinsic g-r colors of massive stars change little from B-type to O-type stars.Five stars lie in the areas of the color–magnitude diagram where either reddened massive stars or evolved post-main sequence stars of lower masses are commonly found.Three of the five stars,sources 1,2,and 3,are located at the cavities of the three IR bubbles,and extended Hαemission is detected around the three IR bubbles.We suggest that sources 1,2,and 3 are candidates for early-B stars associated with the radio-quiet region G014.645–00.606.Particularly,source 1 is an EW type eclipsing binary with a short period of 0.825 day,while source 2 is an EA type eclipsing binary with a short period of 0.919 day.The physical parameters of the two binary systems have been derived through the PHOEBE model.Source 1 is a twin binary of two stars with Teff≈23,500 K,and source 2 contains a hotter component(T_(eff)≈20,100 K)and a cooler one(T_(eff)≈15,500 K).The O-C values of source 1 show a trend of decline,implying that the period of the source is deceasing.Source 1 is likely a contact early-B twin binary,for which mass transfer might cause its orbit to shrink.

Searching for the brightest stars in galaxies outside the Local Group

This paper introduces a technique for searching for bright massive stars in galaxies beyond the Local Group.To search for massive stars,we processed the results of stellar photometry from the Hubble Space Telescope(HST)images using the DAOPHOT and DOLPHOT packages.The results of such searches are demonstrated with examples of the galaxies DDO 68,M94 and NGC 1672.In the galaxy DDO 68,the LBV star changes its brightness,and massive stars in M94 can be identified by excess in the Ha band.For the galaxy NGC 1672,we measure the distance for the first time by the TRGB method,which enabled determining the luminosities of the brightest stars,likely hypergiants,in the young star formation region.So far,we have performed stellar photometry on HST images of 320 northern sky galaxies located at a distance less than 12 Mpc.This allowed us to identify 53 galaxies with probable hypergiants.Further photometric and spectral observations of these galaxies are planned to search for massive stars.

The Depletion of the Red Supergiant Envelope Radiative Zone During Common Envelope Evolution

We conduct one-dimensional stellar evolution simulations of red supergiant(RSG)stars that mimic common envelope evolution(CEE)and find that the inner boundary of the envelope convective zone moves into the initial envelope radiative zone.The envelope convection practically disappears only when the RSG radius decreases by about an order of magnitude or more.The implication is that one cannot split the CEE into one stage during which the companion spirals-in inside the envelope convective zone and removes it,and a second slower phase when the companion orbits the initial envelope radiative zone and a stable mass transfer takes place.At best,this might take place when the orbital separation is about several solar radii.However,by that time other processes become important.We conclude that as of yet,the commonly used alpha-formalism that is based on energy considerations is the best phenomenological formalism.

Classifying Core Collapse Supernova Remnants by Their Morphology as Shaped by the Last Exploding Jets

Under the assumption that jets explode all core collapse supernovae(CCSNe),I classify 14 CCSN remnants(CCSNRs)into five groups according to their morphology as shaped by jets,and attribute the classes to the specific angular momentum of the pre-collapse core.Point-symmetry(one CCSNR):According to the jittering jets explosion mechanism(JJEM)when the pre-collapse core rotates very slowly,the newly born neutron star(NS)launches tens of jet-pairs in all directions.The last several jet-pairs might leave an imprint of several pairs of“ears,”i.e.,a point-symmetric morphology.One pair of ears(eight CCSNRs):More rapidly rotating cores might force the last pair of jets to be long-lived and shape one pair of jet-inflated ears that dominates the morphology.S-shaped(one CCSNR):The accretion disk might precess,leading to an S-shaped morphology.Barrel-shaped(three CCSNRs):Even more rapidly rotating pre-collapse cores might result in a final energetic pair of jets that clear the region along the axis of the pre-collapse core rotation and form a barrel-shaped morphology.Elongated(one CCSNR):A very rapidly rotating pre-collapse core forces all jets to be along the same axis such that the jets are inefficient in expelling mass from the equatorial plane and the long-lasting accretion process turns the NS into a black hole.The two new results of this study are the classification of CCSNRs into five classes based on jet-shaped morphological features,and the attribution of the morphological classes mainly to the pre-collapse core rotation in the frame of the JJEM.

Probing the Progenitor of High-z Short-duration GRB 201221D and its Possible Bulk Acceleration in Prompt Emission

The growing observed evidence shows that the long-and short-duration gamma-ray bursts(GRBs) originate from massive star core-collapse and the merger of compact stars,respectively.GRB 201221 D is a short-duration GRB lasting~0.1 s without extended emission at high redshift z=1.046.By analyzing data observed with the Swift/BAT and Fermi/GBM,we find that a cutoff power-law model can adequately fit the spectrum with a soft E=113keV,and isotropic energy E=1.36× 10erg.In order to reveal the possible physical origin of GRB 201221 D,we adopted multi-wavelength criteria(e.g.,Amati relation,ε-parameter,amplitude parameter,local event rate density,luminosity function,and properties of the host galaxy),and find that most of the observations of GRB 201221 D favor a compact star merger origin.Moreover,we find that α is larger than 2+βin the prompt emission phase which suggests that the emission region is possibly undergoing acceleration during the prompt emission phase with a Poynting-flux-dominated jet.

Influence of Entropy on Composition and Structure of Massive Protoneutron Stars

Adjusting the suitable coupling constants in relativistic mean field(RMF) theory and focusing on thermal effect of an entropy per baryon(S) from 0 to 3, we investigate the composition and structure of massive protoneutron stars corresponding PSR J1614-2230 and PSR J0348+0432. It is found that massive protoneutron stars(PNSs) have more hyperons than cold neutron stars. The entropy per baryon will stiffen the equation of state, and the influence on the pressure is more obvious at low density than high density, while the influence on the energy density is more obvious at high density than low density. It is found that higher entropy will give higher maximum mass, higher central temperature and lower central density. The entropy per baryon changes from 0 to 3, the radius of a PNS corresponding PSR J0348+0432 will increase from 12.86 km to 19.31 km and PSR J1612-2230 will increase from 13.03 km to 19.93 km.The entropy per baryon will raise the central temperature of massive PNSs in higher entropy per baryon, but the central temperature of massive PNSs maybe keep unchanged in lower entropy per baryon. The entropy per baryon will increase the moment of inertia of a massive protoneutron star, while decrease gravitational redshift of a massive neutron star.

Influence of the weakly interacting light U boson on the properties of massive protoneutron stars

Considering the octet baryons in relativistic mean field theory and selecting entropy per baryon S=1, we calculate and discuss the influence of U bosons on the equation of state, mass-radius, moment of inertia and gravitational redshift of massive protoneutron stars （PNSs）. The effective coupling constant gu of U bosons and nucleons is selected from 0 to 70 GeV-2. The results indicate that U bosons will stiffen the equation of state （EOS）. The influence of U bosons on the pressure is more obvious at low density than high density, while the influence of U bosons on the energy density is more obvious at high density than low density. The U bosons play a significant role in increasing the maximum mass and radius of PNS. When the value of gu changes from 0 to 70 GeV-2, the maximum mass of a massive PNS increases from 2.11M to 2.58Me, and the radius of a PNS corresponding to PSR J0348＋0432 increases from 13.71 km to 24.35 km. The U bosons will increase the moment of inertia and decrease the gravitational redshift of a PNS. For the PNS of the massive PSR J0348＋0432, the radius and moment of inertia vary directly with gu, and the gravitational redshift varies approximately inversely with gu.

Massive neutron stars and A-hypernuclei in relativistic mean field models

Based on relativistic mean field（RMF） models, we study finite A-hypernuclei and massive neutron stars. The effective N-N interactions PK1 and TM1 are adopted, while the N-A interactions are constrained by reproducing the binding energy of A-hyperon at 1 s orbit of Λ^40Ca. It is found that the A-meson couplings follow a simple relation, indicating a fixed A potential well for symmetric nuclear matter at saturation densities, i.e., around VΛ=-29.786 MeV. With those interactions, a large mass range of Λ-hypernuclei can be described well. Furthermore,the masses of PSR J1614-2230 and PSR J0348＋0432 can be attained adopting the Λ-meson couplings gσΛ/gσN≥0.73,gωΛ/gωN≥0.80 for PK1 and gσΛ/gσN≥0.81,gωΛ/gωN≥0.90 for TM1, respectively. This resolves the hyperon puzzle without introducing any additional degrees of freedom.

Radiative Diffusion in a Time-dependent Outflow:a Model for Fast Blue Optical Transients

Fast Blue Optical Transients(FBOTs)are luminous transients with fast evolving(typically t;<12 days)light curve and blue color(usually-0.2>g-r>-0.3)that cannot be explained by a supernova-like explosion.We propose a radiative diffusion in a time-dependent outflow model to interpret such special transients.In this model,we assume that a central engine ejects continuous outflow during a few days.We consider the ejection of the outflow to be time-dependent.The outflow is optically thick initially and photons are frozen in it.As the outflow expands over time,photons gradually escape,and our work is to model such an evolution.Numerical and analytical calculations are considered separately,and the results are consistent.We apply the model to three typical FBOTs:PS1-10 bjp,ZTF18 abukavn,and ATLAS 19 dqr.The modeling finds the total mass of the outflow(~1-5 M;),and the total time of the ejection(~a few days)for them,leading us to speculate that they may be the result of the collapse of massive stars.

Pollution tracks of r-process material in[Sr/Ba]vs.[Ba/Fe]space for the early Galaxy

In the early Galaxy,elemental abundances of the extremely metal-poor(EMP)stars contain abundant information about the neutron-capture nucleosynthesis and the chemical enrichment history.In this work,we study the abundance characteristics of Sr and Ba for the EMP stars in the[Sr/Ba]vs.[Ba/Fe]space.We find that there are three boundaries for the distribution region of the EMP stars.The weak rprocess star CS 22897-008 lies on the upper end and the main r-process stars lie on the right end of the region.Near the right boundary of the distribution region,there is an Fe-normal belt.For the EMP stars in the belt,element Fe dominantly originates from the normal massive stars.The low-Sr stars([Sr/Fe]6−0.3)distribute in the region of the lower left of the Fe-normal belt and their Fe should originate partly from the prompt inventory.We find that the formation of the lower boundary of the distribution region is due to the pollution of the main r-process material and the formation of the right boundary could be explained by the combination of the weak r-and main r-process material.Furthermore,the formation of the left boundary is due to the pollution of the weak r-process material.Although the[Sr/Ba]ratios are related to the relative importance of the weak r-process material,the scatter of[Sr/Ba]ratios for the EMP stars mainly depends on the abundance ratio of the weak r-process.

Search for age pattern across spiral arms of the Milky Way

The age pattern across spiral arms is one of the key observational features utilised to study the dynamic nature of the Galaxy’s spiral structure.With the most updated samples of high-mass star formation region(HMSFR)masers,O stars and open clusters,we investigated their distributions and kinematic properties in the vicinity of the Sun.We found that the Sagittarius-Carina Arm traced by HMSFRs,O stars((?)10 Myr)and young open clusters(<30 Myr)seem to deviate gradually towards the Galactic Anticenter(GAC)direction.The Local Arm traced by HMSFRs,O stars,young clusters and also mediumyoung clusters(30-100 Myr)are inclined to gradually deviate toward the Galactic Center(GC)direction.The properties for the Local Arm are supported by a simplified simulation of cluster motions in the Galaxy.Indications of systematic motions in the circular and radial velocities are noticed for the old open clusters(>200 Myr).These results are consistent with the idea that star formation can be triggered by spiral shocks of density waves,and indicate that the corotation radius of the Galaxy is located between the SagittariusCarina Arm and the Local Arm,close to the Solar circle.

Photometric and Spectroscopic Analysis of LBV Candidate J004341.84+411112.0 in M31

We study Luminous Blue Variable(LBV)candidate J004341.84+411112.0 in the Andromeda galaxy.We present optical spectra of the object obtained with the 6 m telescope of BTA SAO RAS.The candidate shows typical LBV features in its spectra:broad and strong hydrogen lines and the He I lines with P Cygni profiles.Its remarkable spectral resemblance to the well known LBV P Cygni suggests a common nature of the objects and supports LBV classification of J004341.84+411112.0.We estimate the temperature,reddening,radius and luminosity of the star using its spectral energy distribution.Obtained bolometric luminosity of the candidate(M_(bol)=-10.41±0.12 mag)is quite similar to those of known LBV stars in the Andromeda galaxy.We analyzed a ten year light curve of the object in R filter.The candidate demonstrates photometric variations of the order of 0.4 mag,with an overall brightness increasing trendΔR>0.1 mag.Therewith,the corresponding color variations of the object are fully consistent with LBV behavior when a star become cooler and brighter in the optical spectral range with a nearly constant bolometric luminosity.LBV-type variability of the object,similarity of its spectrum and estimated luminosity to those of known LBVs allow us to classify J004341.84+411112.0 as an LBV.

CSO CO(2–1) and Spitzer IRAC observations of a bipolar outflow in high-mass star-forming region IRAS 22506+5944

We present Caltech Submillimeter Observatory CO （2-1） and Spitzer IRAC observations toward IRAS 22506＋5944, which is a 104 Lo massive star-forming region. The CO （2-1） maps show an east-west bipolar molecular outflow originating from the 3 mm dust continuum peak. The Spitzer IRAC color-composite image reveals a pair of bow-shaped tips which are prominent in excess 4.5 p.m emission and are located at the leading fronts of the bipolar outflow, providing compelling evidence for the existence of bow-shocks as the driving agents of the molecular outflow. By comparing our CO （2- 1） observations with previously published CO （1-0） data, we find that the CO （2-1）/（1-0） line ratio increases from low （-5 km s- 1） to moderate （- 8-12 km s- 1） velocities, and then decreases at higher velocities. This is qualitatively consistent with the scenario that the molecular outflow is driven by multiple bow-shocks. We also revisit the position-velocity diagram of the CO （1-0） data, and find two spur structures along the outflow axis, which are further evidence for the presence of multiple jet bow- shocks. Finally, power-law fittings to the mass spectrum of the outflow gives power law indexes more consistent with the jet bow-shock model than the wide-angle wind model.

Effects of relativistic parameter sets on tidal deformabilities and f-mode oscillations of neutron stars

The implications of relativistic parameter sets established at saturation density on the tidal deformabilities and f-mode oscillations of neutron stars(NSs)are examined using constraints from the gravitational wave(GW)event GW170817 and NICER.According to our findings,the isovector saturation parameters have a greater impact on the radii and tidal deformabilities of NSs than the isoscalar saturation parameters.Our analysis also examines the impact of saturation properties on f-mode frequencies and finds that f-mode frequencies with 1.4 M⊙(solar mass)are roughly between 1.95 and 2.15 kHz.These findings could be confirmed by future advanced GW detectors.A good linear parameter-independent correlation between f-mode frequencies inferred from saturation parameters in the entire region is also observed,and we attempt to fit an updated version of this universal relationship.Furthermore,we used chiral effective theory(χEFT)together with the multi-messenger astronomy constraints to further reinforce the rationality of the conclusions we have reached.

Impact of the new ^(12)C+^(12)C reaction rate on presupernova nucleosynthesis

The ^(12)C+^(12)C reaction rate plays an essential role in stellar evolution and nucleosynthesis.Nevertheless,the uncertainties of this reaction rate are still large.We calculate a series of stellar evolution models with the near solar abundance from the zero-age main-sequence through presupernova stages for initial masses of 20 M_(⊙) to 40 M_(⊙).The ^(12)C+^(12)C reaction rates from two different studies are used in our investigation.One is the rate obtained using the Trojan Horse Method(THM)by Tumino et al.[Nature 557(7707),687(2018)],and the other was obtained by Mukhamedzhanov et al.[Physical Review C 99(6),064618(2019)](Muk19).Then,comparisons of the nucleosynthesis and presupernova isotopic abundances are conducted.In particular,we find that in the C burning shell,models with the THM produce a smaller amount of ^(23)Na and some neutron-rich isotopes than Muk19.The difference in the abundance ratios of Na/Mg,S/Mg,Ar/Mg,and K/Mg between the two models are apparent.We compare Na/Mg obtained from our theoretical presupernovae models with Na/Mg in stellar atmospheres observed with high-resolution spectra as well as from the latest galactic chemical evolution model.Although Na/Mg obtained using the THM is within 2σ of the observed stellar ratio,the theoretical uncertainty on Na/Mg introduced by the uncertainty of the ^(12)C+^(12)C reaction rate is almost equivalent to the standard deviation of astronomical observations.Therefore,a more accurate ^(12)C+^(12)C reaction rate is crucial.

Short-range correlation effects in neutron star's radial and non-radial oscillations

In this study,we determine the influence of the nucleon-nucleon short range correlation(SRC)on static spherically symmetric neutron stars(NSs)from the perspectives of radial and nonradial oscillations for the first time.We revise the equation of state and coupling parameters in the relativistic mean field theory after considering the SRC effect,and select the hyperon coupling parameters as the SU(3)model.For the non-radial oscillations,the SRC effect decreases the f-mode frequency by 0.2~0.3 kHz.For the radial oscillations,it decreases the fundamental radial frequency f_(1) by 0.75~0.85 kHz.Additionally,we refit the linear relationship between the average density and f-mode frequency for SRC.Combining the characteristics of the radial and non-radial frequencies,we provide a view of inferring the maximum mass of NSs.Owing to the characteristics of the SRC influence on the radial frequency,we expect that the SRC can be tested by future observation and can also be used as a probe for the structure inside NSs.