• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    A Common Envelope Jets Supernova (CEJSN) Impostor Scenario for Fast Blue Optical Transients

    2022-05-24 06:34:08NoamSoker

    Noam Soker

    Department of Physics,Technion,Haifa,3200003,Israel;soker@physics.technion.ac.il

    Abstract I propose a new scenario,the polar common envelope jets supernova (CEJSN) impostor scenario,to account for AT2018cow-like fast blue optical transients (FBOTs).The polar CEJSN impostor scenario evolves through four main phases.(1)A red supergiant(RSG)star expands to tidally interact with a neutron star(NS)companion(or a black hole).The interaction increases the RSG mass loss rate to form a circumstellar matter (CSM) halo to r ?0.1 pc.(2)Shortly before the onset of a common envelope evolution(CEE)and about a year before explosion the NS accretes mass from the RSG envelope and launches jets that inflate two opposite lobes in the CSM within ≈100 au.(3)The NS-RSG system enters a CEE phase during which the system ejects most of the envelope mass in a dense equatorial outflow.(4)At the termination of the CEE the leftover envelope forms a circumbinary disk around the NS-core system.The NS accretes mass from the circumbinary disk and launches energetic jets that,when colliding with the fronts of the CSM lobes,power an FBOT event.The low mass of the jets-lobes interaction zones and their large distance,of about 100 au,from the center account for the fast transient.In the future the core collapses to form a second NS.In the far future the two NSs might merge.I suggest that FBOTs and similar fast transients are CEJSN impostors which compose a large fraction of the progenitors of NS-NS merger binaries.

    Key words:neutron stars– supernovae– stellar jets– common envelope binary stars– transient sources

    1.Introduction

    In the common envelope jets supernova (CEJSN) scenario the source of explosion energy is the accretion of gas from a red supergiant (RSG) star onto a neutron star (NS) or a black hole (BH),hereafter NS/BH,that spirals-in inside the RSG,first inside the envelope and then inside the RSG core (e.g.,Papish et al.2015;Soker &Gilkis 2018;Gilkis et al.2019;Grichener &Soker 2019a;López-Cámara et al.2019,2020;Soker 2021).The accretion process proceeds via an accretion disk that launches energetic jets that explode the star.

    In case the NS/BH spirals-in inside the envelope of the giant star but does not enter(does not merge with)the core,the event is a CEJSN impostor (Gilkis et al.2019).

    The jets that the NS/BH launches at late times interact with the remaining RSG envelope and with the circumstellar matter(CSM) that the earlier common envelope evolution (CEE)ejected.The interaction excites two shock waves,the forward shock that expands into the ambient gas (the envelope and/or the CSM)and the reverse shock that shocks the outflowing jets’gas.The cocoon is the hot post-shock regions of the two shocks,which have a contact discontinuity between them.The hot cocoon radiates some of its energy (e.g.,Schreier et al.2021),leading to a bright transient event that mimics a core collapse supernova(CCSN).I do not always explicitly mention the cocoon.However,it should be taken for granted that any jet interaction with the envelope and/or with the CSM leads to the formation of a cocoon (e.g.,López-Cámara et al.2019;Schreier et al.2021).For further discussion of the accretion rate in the CEJSN scenario see Grichener et al.(2021) and Hillel et al.(2022).

    Because they mimic CCSNe and have a rich variety of properties,e.g.,the companion can be an NS or a BH with a range of masses,the RSG envelope mass and radius can vary by factors of about an order of magnitude,the core mass and composition can vary,and the eccentricity of the orbit before the strong interaction can be from zero to a high value (e.g.,Soker et al.2019).CEJSN and CEJSN impostor events are very likely to be behind several types of transient events that CCSNe cannot account for.Th?ne et al.(2011) explain the unusual gamma-ray burst GRB 101225A by an NS that merged with a helium star.The enigmatic supernovae iPTF14hls(Arcavi et al.2017) and SN 2020faa (Yang et al.2021) might be CEJSN events (Soker &Gilkis 2018).Adopting the CEJSN scenario and processes from these earlier studies,Dong et al.(2021)proposed that the luminous radio transient VT J121001+495647 was a CEJSN event.CEJSNe might account also for some processes that require extreme conditions,like some fraction of the r-process nucleosynthesis (Grichener &Soker 2019a,2019b;Grichener et al.2022) and,for a CEJSN impostor with a BH companion,the formation of a fraction of the very high-energy neutrinos (Grichener &Soker 2021).

    In this study I concentrate on AT2018cow-like fast blue optical transients (FBOTs).AT2018cow-like FBOTs are fastrising,only a few days (e.g.,Prentice et al.2018;Perley et al.2019),bright transients,i.e.,they might be brighter than superluminous CCSNe (e.g.,Margutti et al.2019).They display high velocities of0.1c with a total kinetic energy of?1051–1052erg (e.g.,Coppejans et al.2020),with hydrogen lines (e.g.,Margutti et al.2019).They tend to occur in starforming galaxies (e.g.,Prentice et al.2018).They might display rapid X-ray variability (e.g.,Pasham et al.2021;Yao et al.2022) and have dense CSM (e.g.,Bright et al.2021;Nayana &Chandra 2021).In Section 4 I return to discuss in more detail the properties of FBOTs,including their uncertain rate.

    In the CEJSN scenario an NS/BH that was formed in an earlier CCSN of the initially more massive star in a binary system enters the envelope of an RSG that evolved from the initially less massive star of the binary.The main power of the event comes from the very energetic jets that the NS/BH launches as it spirals-in inside the envelope and then as it spirals-in inside the core of the RSG star.The very powerful jets are the main difference from the many types of binary interactions with and without CEE(e.g.,Han et al.2020),from cases where a main sequence companion launches jets in a CEE (e.g.,Shiber et al.2016),and from the more closely related NS/BH scenarios of CEE without jets (e.g.,Fryer &Woosley 1998;Zhang &Fryer 2001;Barkov &Komissarov 2011;Th?ne et al.2011;Chevalier 2012;Schr?der et al.2020).

    In the polar CEJSN channel of the CEJSN scenario that Soker et al.(2019) constructed,the early jets remove most of the envelope gas along the polar directions.Late jets then expand almost freely to large distances such that they might account for the fast rise and decline of FBOTs,as well as for the very high outflow velocities of up to >0.1c (e.g.,Margutti et al.2019;Perley et al.2019;Coppejans et al.2020).Jet-CSM interaction must take place to convert kinetic energy to thermal energy of the post-shock jet and CSM gas (the cocoon).This hot gas radiates a large fraction of its energy to power a bright event.

    The enigmatic AT2018cow event (Prentice et al.2018) and similar FBOTs have promoted the development of several scenarios for AT2018cow-like FBOTs (e.g.,Liu et al.2018;Fox&Smith 2019;Kuin et al.2019;Lyutikov&Toonen 2019;Margutti et al.2019;Quataert et al.2019;Yu et al.2019;Leung et al.2020;Mohan et al.2020;Piro&Lu 2020;Uno&Maeda 2020;Kremer et al.2021;Xiang et al.2021;Chen &Shen 2022;Gottlieb et al.2022).In the present study I do not compare the different scenarios with each other but rather present the polar-CEJSN impostor channel (Section 2) and its new ingredients (Section 3).While Soker et al.(2019)presented a polar CEJSN scenario,I here present a polar CEJSN impostor scenario,i.e.,the NS/BH does not enter (nor destroy) the RSG core.The scenario I propose here does not replace the one Soker et al.(2019) proposed,but rather adds another channel to the CEJSN and CEJSN impostor scenarios.In Section 4 I discuss how the polar CEJSN impostor scenario might account for the observational properties of AT2018cowlike FBOTs.I summarize in Section 5.

    2.The Polar-CEJSN Impostor Scenario

    I describe the general polar CEJSN impostor scenario here.I will present the details of some ingredients in Section 3.I emphasize that the main difference between the polar CEJSN impostor that I propose here for FBOTs and the polar CEJSN scenario for FBOTs that Soker et al.(2019)proposed is that,as the“impostor”in the name implies,in the present scenario the NS does not enter (does not merge with) the core of the RSG star.

    The system will experience the CEJSN impostor scenario if the system manages to remove the entire RSG envelope before the NS reaches the core.This in turns depends on the radius and mass of the RSG and on the eccentricity of the orbit when the strong binary interaction of the NS and the RSG starts.A lower envelope mass,a higher RSG radius and a higher eccentricity favor envelope removal,hence CEJSN impostor evolution.However,there is yet no quantitative study of these effects.

    The main evolutionary phases after the formation of a binary system of an NS and an RSG star are as follows.

    1.Pre-CEE.The NS perturbs the RSG envelope by tidal forces,that both deform the envelope and spin it up.As a result of that the mass loss rate increases by possibly up to about an order of magnitude.This wind forms the CSM halo (Section 3.1).

    2.The onset of the CEE.As the NS enters the envelope(the onset of the CEE)it starts to accrete mass via an accretion disk and launches jets (see López-Cámara et al.2020 for simulations and Hillel et al.2022 for discussion of the accretion process).During this phase the jets interact with the dense wind and shape two opposite lobes along the bipolar directions,as observed in some planetary nebulae(Section 3.1).Since the NS is just about to enter the envelope and accretion is similar to a Roche Lobe overflow,the jets manage to expand and accelerate gas along the two opposite polar directions.This takes place about a year before the explosion,which is the duration of the CEE phase τCEE(about the dynamical time on the surface of the RSG star,e.g.,Lau et al.2022),and might lead to a faint precursor at tPreC??τCEE≈?1 yr relative to the FBOT event itself.

    3.The main CEE phase.During the main CEE phase the NS spirals-in inside the RSG envelope,and when it removes the entire envelope it ends at a final orbital radius from the RSG core that I mark as aNC.The jets do not penetrate the envelope and most mass is ejected near the equatorial plane and up to mid-latitudes (e.g.,Schreier et al.2021).The average jets’ power is ≈1043erg s?1.The outflow does not destroy the polar lobes,but rather mainly adds mass to the dense equatorial outflow.

    4.Accretion from a circumbinary disk as the energy source of the explosion.Following Kashi &Soker (2011),I assume that at the end of the main CEE phase most of the leftover envelope mass forms a circumbinary disk around the NS-core binary system(Section 3.2).I further assume that there are two opposite openings in the envelope along the polar directions (funnels;e.g.,Soker 1992;Zou et al.2020).Therefore,at the end of the CEE the jets that the NS launches as it accretes hydrogen-rich gas from the circumbinary disk manage to escape the star.The “polar” in the name of the polar CEJSN scenario that Soker et al.(2019)proposed and in the present polar CEJSN impostor scenario comes from the (almost) empty polar directions inside the RSG envelope that allow the jets to propagate to large distances.The difference is that in the polar CEJSN the NS launches jets as it accretes mass from the destroyed massive core of the RSG,leading to (1) a very energetic event with a kinetic energy of1052erg,much above the kinetic energy of typical CCSNe (≈few×1050?few×1051erg),and (2) hydrogen-poor jets.The typical explosion energy of the FBOTs that I consider here is that of typical CCSNe and the fast outflow is hydrogen rich (Section 4),which lead me to propose that the NS does not enter the RSG core,hence the“impostor”in the name of the scenario.

    5.Late accretion phase from the circumbinary disk.The accretion process depletes the circumbinary disk,which is the reservoir for the accreted mass,and therefore the accretion rate decreases.The accretion rate decreases over a timescale of weeks,namely,it might continue for weeks(Section 3.2).This explains the late X-ray rapid variability (Section 4).

    6.Far-future events.This scenario predicts that in the future(up to about hundreds of thousands of years in the future)the RSG core experiences a striped-CCSN (type Ib or type Ic CCSN)event,leaving behind a second NS.At an even later time the two NSs might merge as in the channel-I CEE scenario that Vigna-Gómez et al.(2018)study for NS-NS merger events.

    3.New Ingredients

    Here I describe the new ingredients of the polar-CEJSN impostor scenario with respect to the polar-CEJSN scenario that Soker et al.(2019)introduced.The main difference,as the“impostor” in the name implies,is that the NS does not enter(does not merge with)the core of the RSG.Note,however,that the bipolar CSM structure that I describe next might be also the CSM structure in the polar CEJSN scenario.

    Several emission sources and properties as well as the dense equatorial CSM of the polar-CEJSN impostor scenario are similar to the general picture that Margutti et al.(2019)describe.These similar properties include a bipolar structure of the interaction,a central engine,late X-ray variability that results from changes in the power of the central engine,a fast polar outflow that explains the absorption lines and an equatorial outflow that accounts for the ?several×1000 km s?1broad hydrogen emission lines at later times.Another process that Soker et al.(2019)adopted from Margutti et al.(2019) is the receding of the photosphere from the fast polar ejecta to the slower equatorial ejecta,which might explain the transition from X-ray to ultraviolet/visible/infrared emission.Here I further note that some processes and properties that I describe here are also similar to some processes in the FBOT scenario that Gottlieb et al.(2022) develop.There are,however,some basic differences between the scenarios,in particular that in the polar-CEJSN impostor scenario the NS is an old one,rather than a newly born NS or BH as in the scenario of Gottlieb et al.(2022),and that in the polar CEJSN impostor scenarios the CSM contains two opposite lobes along the polar directions.The strong interaction of the jets that yields the early optical emission is with the CSM,rather than with the stellar envelope that Gottlieb et al.(2022) consider.At late times,i.e.,after a few weeks,some emission properties,like radio emission,are similar as in both scenarios the interaction of the jets is with the CSM.

    3.1.Pre-explosion Bipolar CSM

    I take the pre-explosion CSM structure just after the NS reaches close to the RSG core to be similar to that of some bipolar planetary nebulae,e.g.,the Owl Nebula (NGC 3587;e.g.,García-Díaz et al.2018).The structure of this type of planetary nebula contains two opposite low-density lobes (the white ellipse in the upper half of Figure 1)inside a spherical or an elliptical shell (black dots with beige background in Figure 1).Jets inflate such lobes.In the present scenario the NS launches jets just before it enters the envelope of the RSG star.At the same time its interaction with the envelope increases the mass loss rate in the RSG wind by about three orders of magnitude to be ≈0.01 M⊙yr?1.These jets that are active for several months inflate the lobes as they interact with the intensive wind from the RSG star.The shell is optically thick.

    The inflation of these two lobes compresses the equatorial gas (e.g.,Akashi &Soker 2008),adding to the action of the gravity of the companion,here an NS,and jets at the onset of the CEE (e.g.,Shiber et al.2019).I schematically draw this CSM structure in Figure 1.Note that there is a mirror symmetry about the equatorial plane,but in Figure 1 the upper loberepresents the case where the two lobes are closed by the optically-thick shell (the walls of the lobes) while the lower lobe corresponds to the case where the two lobes are open.The ratio of the energy of the pre-CEE jets that inflate the bubbles to the density of the dense wind is the main factor that determines the shape of the lobes(other parameters include the opening angle of the jets and the wind velocity).Very energetic jets will break out from the dense wind and form open lobes.

    Figure 1.A schematic drawing (not to scale) of the CSM just before the explosion,which is also very shortly after the system exits the CEE.The plane of the figure is the meridional plane that momentarily contains the RSG core(red filled-circle) and the NS (green dot).There is an axial-symmetry about a vertical axis in the figure through the center of mass of the core-NS binary system(the symmetry axis coincides with the RL arrow in the figure).There is also a mirror symmetry about the equatorial plane.However,in this figure the upper lobe represents the case of the two lobes being closed by the opticallythick shell,while the lower lobe represents the case where the two lobes are open.The main structural components are the low-density lobes(in white),the optically thick shell that encloses the two lobes from all or most sides (dottedbeige),the dense equatorial CSM (red),the circumbinary disk (the two beige zones in the inner part)and the CSM halo that extends beyond the boundary of the figure(dotted-white area).RL ≈1015 cm is the length of one lobe along the polar direction,aNC ?2 R⊙is the post-CEE orbital separation of the NS and the core,and RCBD ?several×aNC is the radius of the post-CEE circumbinary disk.

    The duration of the CEE phase,τCEE,from the onset of the CEE until the compact companion spirals-in deeply into the giant envelope is about the orbital (Keplerian) time on the surface of the giant (e.g.,Glanz &Perets 2021).In the present study the orbital time on the surface of an RSG is about a year.During that time the CSM lobes that the jets expel reach a distance of

    where for the expansion velocity of the lobe,vL,I take a velocity larger by about a factor of three than the escape velocity from the RSG star,vL?3vesc(RRSG),where RRSG?few×AU is the radius of the RSG.The equatorial outflow is slower and closer to the escape velocity from the RSG.

    I emphasize the following properties.

    1.The formation of the two opposite polar lobes of the CSM can take place before the CEE.Examples are the Red Rectangle Nebula around a post-asymptotic giant branch binary system and the Homunculus Nebula around the massive binary system Eta Carinae.Both binary systems have bipolar CSM around them,e.g.,Cohen et al.(1975)and Smith (2006),respectively,but did not enter a full CEE.I suggest the same in the case of some FBOTs.

    2.Likely,there is a more spherical CSM and of lower density at much larger distances,as observed in many planetary nebulae(e.g.,Corradi et al.2003),e.g.,the Owl Nebula.This is the CSM halo (black dots with a white background in Figure 1),which is optically thin.

    3.It is during the CEE that the binary system ejects the massive and dense equatorial CSM.

    4.During most of the CEE the jets do not break out from the envelope (e.g.,Grichener &Soker 2021;Hillel et al.2022).They can break out only at the beginning of the CEE just before the NS enters the very dense parts of the envelope,or at the end of the CEE when the CEE cleared the polar directions (e.g.,Soker 2019).

    5.Because the companion star spins-up the envelope during the CEE,in some cases the fast envelope rotation leads to the opening of a funnel along the two polar directions(e.g.,Soker 1992;Zou et al.2020).This ensures that the jets expand almost freely before they encounter the dense parts of the lobes.The two opposite openings(funnels)in the polar directions give this scenario the name “polar CEJSN” (Soker et al.2019) or polar CEJSN impostor(present study).The funnels in the polar directions are another property that distinguishes the polar CEJSN scenario from the scenario that Gottlieb et al.(2022)propose in which the jets strongly interact with the envelope.

    6.It is possible that the fronts of the two lobes are open.I schematically present this possibility in the lower lobe in Figure 1.

    3.2.Post CEE Accretion from a Circumbinary Disk

    Margutti et al.(2019) find that the central power source(“engine”) of AT2018cow should release a total energy of?1050–1051.5erg over a characteristic timescale of ≈103–105s.They also deduce that the ejecta mass is Mej?0.1–1 M⊙a(bǔ)nd that it contains hydrogen and helium,but a limited56Ni mass of MNi<0.04 M⊙.The ejecta velocity spans a large range,from vej?0.01 km s?1to vej?0.2c.I raise the possibility here that the NS launches these outflows as it accretes mass from a circumbinary disk at the end of the CEE inside the RSG envelope,i.e.,the RSG core is still intact.

    Kashi &Soker (2011) argue that ηCBD?0.01–0.1 of the common envelope might remain bound in a circumbinary disk after the compact companion ends the CEE inside the giant envelope (later it might enter the core or destroy the core).A dynamically stable circumbinary disk extends from RCBD?2.5aNCto RCBD?several×aNC,where aNCis the orbital separation of the NS and the core (assuming a circular orbit).Here I do not require the circumbinary disk to be dynamically stable,and its inner boundary will be much closer to the NS-core binary system because the NS-core system has just emerged from the CEE.I assume that the NS launches fast jets at velocity vjshortly after the CEE phase and that the jets carry a fraction of η2j?0.05–0.1 of the circumbinary mass.For an RSG envelope mass Menvthe mass in the two jets and their energy are then

    respectively.This energy is about equal to the jets’ energy in the scenario of Gottlieb et al.(2022),but here the jets are mainly baryonic while Gottlieb et al.(2022) consider relativistic jets.

    The duration of the accretion process is about the viscosity time of the circumbinary disk,which is ≈10–100 times the Keplerian orbital time of the disk around the binary.The accretion time period of the entire post-CEE circumbinary disk,from its inner boundary to its outer boundary,is

    whereMcoreis the mass of the RSG core and I took the inner and outer boundaries of the circumbinary disk to beRCBD,in?aNCand RCBD,out?5aNC,respectively.I expect a high accretion rate in the first few days,˙Macc≈that slowly decreases with timescales of weeks to months tobut not necessarily monotonically.The decline in the mass accretion rate might at best launch weak jets that cannot propagate to large distances.This might be related to the finding of Bietenholz et al.(2020) that there are no observed long-lived relativistic jets in AT2018cow.

    The phase of mass accretion from the circumbinary disk is a short phase that precedes the BB mass transfer phase in the channel-I CEE scenario that Vigna-Gómez et al.(2018)describe (their Figure 5) for the formation of double NSs that later merge by gravitational wave emission.Namely,I suggest here that some of the AT1018cow-like FBOTs are progenitors of binary NS systems that much later merge by gravitational waves.In the channel-I scenario that Vigna-Gómez et al.(2018) describe,the NS ends at aNC?1 ?few×R⊙from the helium core.The helium core masses in most cases they consider for channel-I areThey consider a phase of stable mass transfer from the core to the NS,which implies aNC1.5 R⊙(e.g.,Tauris et al.2015).Since I do not require here a stable mass transfer at a phase after the explosion,I allow for smaller orbital separations of even aNC?1 R⊙.

    3.3.The Interaction of the Jets with the Lobes

    As the jets expand they cool adiabatically.To channel a large fraction of the jets’kinetic energy to radiation they must collide with an ambient gas.In the polar CEJSN impostor scenario the jets collide with the bipolar CSM as I schematically present in Figure 2 (only for one lobe).Due to the orbital motion of the NS and the interaction of the jets with some tenuous gas in the lobes’interior(the lobes are not completely empty)I expect the jets not to be collimated,and they might even be wide.Specifically,interaction of jets with a dense close (to the launching point) gas can collimate the jets.Here the polar directions are almost empty,and the tenuous gas along the polar directions will be accelerated and entrained by the jets to a wider polar outflow(wide jets).In addition,the orbital motion prevents the build-up of a dense gas very close to the NS.As well,when they collide with the dense front of the lobe the jets might be slower and more massive than at their origin near the NS,as they entrain the tenuous gas along the polar directions.

    Consider then that a jet collides with the front of the lobe over an angle αI(Figure 2)such that the solid-angle of the two interaction regions (one in each lobe) isThe jets shock the front of the close lobe,and when the forward shock breaks out from the lobe there is a hot cocoon of width ΔRI.At that time the photosphere starts to move inward with respect to the mass of the cocoon,and a weaker shockcontinues to propagate into the optically-thin halo,as I schematically illustrate in Figure 3.

    Figure 2.A schematic drawing(not to scale)of the interaction of the post-CEE hydrogen-rich jets with the lobes before shock break-out.The NS launches the jets as it orbits the RSG core and accretes from the circumbinary disk.Interaction takes place in both lobes,but the drawing is only in the upper one.The elliptical CSM shell(beige background)is optically thick,while the CSM halo is optically thin.

    The combined volume of the two hot cocoons that break out from the lobe in the jet-lobe interaction regions isThe mass inside each of the two hot cocoons includes the mass of the jet,the mass that the jet drags with it as it expands through the lobe’s interior and the mass of the dense shell of the lobe within the interaction region.This mass will be several times the original mass of the jet.I scale this mass with MI?0.1 M⊙.The photon diffusion time from this region is

    where τ=ρIkΔRIis the optical depth of the interaction region,ρIis the density of the cocoon and k is the opacity.Equation(5)gives the timescale for the variation of the luminosity and photosphere size at rise and early decline.

    Due to this geometry a spherically-symmetric model of the photosphere will yield a smaller distance from the center,depending on the viewing angle of the observer.Approximately,the inferred radius of a spherical model will beFor example,the earliest photosphere radius that Perley et al.(2019) deduce for AT2018cow is Rph,sph=8×1014cm.For the scaling I use here of ΩI=π the lobe radius is RL?1.6×1015cm.Perley et al.(2019)find that the radius of the photosphere that they calculate decreases from the first time they calculate this radius.I attribute this monotonic decrease of the photosphere radius to the structure of the lobes,as I draw schematically in Figure 4.At the decline phase the low density regions near the symmetry axis are already optically thin and the jets clean these regions.The shocks that run into the walls of the lobes at lower latitudes form hot cocoons there that continue to radiate,but from smaller and smaller photosphere areas.

    Figure 3.Similar to Figure 2 but after shock break-out and before maximum optical luminosity.The dashed blue line marks the photosphere on one side(another exists on the opposite lobe but is not drawn).The hatched blue volume outside the photosphere that once was a hot cocoon is already cold due to radiative cooling.

    With the simple scenario that I describe here I cannot follow to late times when the luminosity and photosphere area have dropped by a large factor.The reason is that the photosphere area becomes small,and therefore it is sensitive to the initial conditions and the jet-lobe interaction that I cannot follow with analytical means at late times.Perley et al.(2019) find that at t ?40 days the luminosity and photosphere area of AT2018cow have dropped by about three and two orders of magnitude,respectively.At these late times the optical emission properties depend on small-scale interaction of the jets with left-over RSG envelope gas near the polar directions.As well,there might be a slow disk-wind from the circumbinary disk.

    At late time an observer that is not at too low latitudes can see the other side of the lobes’ walls,and then the central region,namely,the core-NS binary system and the circumbinary disk (that has its mass decreasing due to accretion and disk-wind).

    I end this section by mentioning the case of open lobes.In this case the jets interact with less mass of the lobes,but still interact with a mass.I expect that in that case the fast outflow(jets) will be more pronounced and that the channelling of the kinetic energy of the jets to radiation will be less efficient.Therefore,the FBOT will be fainter.The exact properties depend strongly on the density of the gas inside the lobes and the opening angle of the lobes,as well as the jets’properties of course.

    4.Accounting for Observational Properties

    In this section I discuss some of the properties of AT2018-like FBOTs alongside the processes in the polar CEJSN impostor scenario that might account for these properties.

    In comparing the theoretical expectations with observations one should note the following.First,there are several processes that require quantitative studies to determine more accurate values and to check the suggestions that I raise here.In the present study I propose the scenario but there are hydrodynamical simulations,radiative transfer calculations and population synthesis studies to conduct in the future.

    Second,the parameter space of the polar CEJSN impostor scenario is large.Coppejans et al.(2020)point out the diversity of FBOTs.Indeed,the polar CEJSN impostor scenario has properties that can change from one FBOT to another,including the CSM structure (compare the two lobes in Figure 1),the mass in the circumbinary disk,the companion type (an NS or a BH),and the evolutionary phase of the RSG and its mass at the onset of the CEE.

    Figure 4.Similar to Figure 3 but during the decline phase of the light curve as the photosphere recedes into lower latitudes of the walls of the lobes.The jets are weaker now and have already cleaned the polar directions.Note that now an observer at a high enough latitude can see the other side of the lobe,as well as the central region.(As before,the lower lobe shows the structure before the interaction with jets for comparison,although interaction takes place there also).

    In what follows I consider the compact object that spirals-in inside the RSG star to be an NS.However,in some cases it might be a BH.

    4.1.Star-forming Galaxies

    AT2018cow-like FBOTs tend to occur in star-forming galaxies,namely,coming from massive stars (e.g.,Prentice et al.2018;Perley et al.2019;Lyman et al.2020;Perley et al.2021).Soker et al.(2019) noted that from their population synthesis study Mapelli et al.(2018) find that in the local universe NS-NS mergers tend to occur shortly after star formation.Since many CEJSN events lead to NS-NS close binary systems,Soker et al.(2019) argued that we do expect CEJSN events to take place in star-forming galaxies.I here add that the envelope of the RSG star must be massive enough,crudely M2,env10 M⊙,to force the NS companion to spiral-in down to final core-NS orbital separation of about aNS?1–3 R⊙.This implies that not only the progenitor of the NS was a massive star,but the initially less massive star,the progenitor of the RSG star,should also be a massive star,i.e.,with an initial mass of8 M⊙.However,in the present scenario the RSG envelope cannot be too massive as to force the NS to spiral all the way to the core.I return to my suggestion that some FBOTs are progenitors of binary NSs that later merge in Section 4.10.

    4.2.Hydrogen in the Fast Ejecta

    There are two sources of hydrogen in the fast ejecta.(1)One of the new ingredients that I add here is that the NS accretes mass from a circumbinary disk when it orbits the core at a very small radius of aNS?1–3 R⊙(Section 3.2).This circumbinary disk is the leftover of the RSG envelope,hence it is hydrogenrich.Its inner and outer radii are RCBD,in?aNCand RCBD,out?5–30 R⊙,respectively,while its mass is crudely MCBD?0.1–1 M⊙.(2) The jets that the NS launches at the explosion as it accretes mass from the circumbinary disk(Section 3.2;it might accrete some mass from the core of the RSG) sweep hydrogen-rich CSM that was ejected during the CEE of the NS inside the RSG envelope.

    4.3.High Velocities of vej >0.1c

    The outflow velocities in the FBOTs AT2018cow,ZTF18abvkwla and CRTSCSS161010 J045834-081803 (CSS16 hereafter)are ?0.1c(e.g.,Margutti et al.2019),?0.3c(e.g.,Ho et al.2020) and ?0.5c (Coppejans et al.2020),respectively.I attribute these high velocities(relative to CCSNe)to the clean polar directions that allow the jets to expand while interacting with a relatively low envelope and CSM mass,at least during part of the event.

    Perley et al.(2019) find that in AT2018cow the highvelocity absorption lines disappear after about two weeks.I attribute this to the fact that the fast material has expanded to large distances and the dominant absorbing gas is the slower gas at large angle with respect to the symmetry axis,namely,a gas that comes mainly from the slower walls of the polar lobes(Figure 4).This suggestion requires further study by hydrodynamical simulations.

    4.4.A Small Mass of Ejecta at vej >0.1c

    Coppejans et al.(2020),for example,estimate the fast ejecta mass in the FBOT CSS16 to be ?0.01–0.1 M⊙a(bǔ)nd the kinetic energy ?1051–1052erg.The circumbinary disk at the end of the CEE inside the envelope is expected to contain a small fraction of the original envelope mass,and therefore the mass that the jets carry is small,as Equation (2) affirms.The jets expand almost freely,but not totally so that the lobes cannot be completely empty.There is a gas inside the lobes before explosion,with lower densities and higher temperatures than those of the walls of the lobes (the optically-thick shell).As such,the jets entrain gas and slow down.For example,in AT2018cow I expect the jets to interact with gas in the lobes that is a few times more massive than the jets,so the jets slow down from vj?105km s?1at launching to a few times slower,?0.1c.As well,parts of the jets,in particular at large angles with respect to the symmetry axis,might be chocked by the wider walls of the lobes(the shell)at low latitudes.Overall,the mass in the fast ejecta might crudely be

    where in the second equality I inserted Equation (2).This is compatible with the fast ejecta in the FBOT CSS16(Coppejans et al.2020),?0.01–0.1 M⊙,as the value of η2jηCBDmight be larger by up to a factor of few than the scaling of Equations(2)and (6).

    4.5.Total Event Energy of ≈1050-1052 erg

    I attribute this energy of <1052erg to the accretion from a low mass circumbinary disk,as Equation(3)shows.I note that in accretion disks of collapsars,i.e.,the collapse of the core in a CCSN that forms a BH with an accretion disk around it,nuclear burning with helium might take place (e.g.,Zenati et al.2020).Here the accretion disk is hydrogen-rich and has lower densities.Future numerical simulations should examine whether nuclear reactions take place in the accretion disk.

    4.6.A Fast Rise (a Day to a Few Days)

    The rise time of AT2018cow,as an example,was less than three days,and at a time of Δt <1.3 days its magnitude raised by 4.2 mag(Prentice et al.2018;Perley et al.2019).I estimate this typical timescale in Equation (5) as the photon diffusion time.With that comes the observations that AT2018cow-like transients are the most luminous and fast type of FBOTs (e.g.,Ho et al.2021).This requires an efficient channelling of kinetic energy to thermal energy and then radiation.The low densities at large distances account for that (Equation (5)).

    The shock breakout through the walls of the lobes at large distances accounts for the properties of AT2018cow as Perley et al.(2019) deduced.They concluded that at shock breakout the photosphere should be unbound and results from a preexplosion dense wind or shell ejection.Perley et al.(2019)concluded also that the CSM shell should be localized in extent.The fronts of the lobes (Figures 1–4) have these properties.

    4.7.Decreasing Photosphere in AT2018cow

    In AT2018cow the photosphere decreases during the first several weeks of observations (Perley et al.2019).In Section 3.3 I attributed this structure to the pre-collapse lobes and the nature of the interaction with the jets that removes the fronts of the lobes such that the area of photosphere decreases(schematically in the transition from Figures 3 to 4).The complicated structure of the photosphere might account for the increase in photospheric temperature after several weeks in AT2018cow (Perley et al.2019).Future hydrodynamical numerical simulations will determine the exact behavior.The large parameter space of CSM shapes and densities and the jets’ properties will require an intensive study.

    4.8.Rapidly Variable X-ray Source

    The X-ray emission of AT2018cow-like FBOTs varies from timescales of days(e.g.,Yao et al.2022 for AT2020mrf)down to a fraction of a second (e.g.,Pasham et al.2021 for AT2018cow).Some X-rays might come from the jets as they pass through shocks in optically-thin regions for X-ray.This might account for hours to days timescales of variability.This is a subject of future hydrodynamical simulations.

    The rapid variability can come from the accretion disk.The X-ray emission of FBOTs requires a central compact source(engine),e.g.,AT2020xnd (Ho et al.2022).The rapid X-ray variability even weeks after explosion points to a central energy course (engine;e.g.,Pasham et al.2021).The central engine most likely involves jets as in the polar CEJSN scenario(Soker et al.2019) and in the scenario where the inner parts of a star collapse to form a BH that launches jets (e.g.,Perley et al.2021;Gottlieb et al.2022).This is also the case with the polar CEJSN impostor scenario that I propose here where a circumbinary disk feeds the accretion disk around the NS (or BH in some cases) that launches the jets.Another central engine is a magnetar (e.g.,Mohan et al.2020).However,the formation of an energetic magnetar most likely is accompanied by more energetic jets (e.g.,Soker &Gilkis 2017).Note that the tidal disruption scenario(e.g.,Kuin et al.2019;Perley et al.2019) seems unable to explain FBOTs because FBOTs come from massive stars and have dense CSMs (e.g.,Huang et al.2019;Yao et al.2022).

    4.9.Possible Dense CSM at r ?1017 cm

    AT2018cow (e.g.,Nayana&Chandra 2021)and AT2020xnd(e.g.,Bright et al.2021) have a dense CSM that extends up to r ≈1017cm.Coppejans et al.(2020) conclude that the ejecta of the FBOT CSS16 interacts with a dense wind at r ≈1017cm.From the deceleration of the outflow they conclude that the CSM mass that the outflow(blast wave)sweeps is comparable or larger than the mass of the fast material (ejecta).In the polar CEJSN impostor scenario the extended CSM is the CSM halo(Figure 1)that I discussed in Section 3.1,which extends to Rhalo≈1017cm and which,although containing a mass larger than a regular RSG wind,i.e.,corresponding to a mass loss rate of ˙Mhalo≈is optically thin.The close CSM,at r <RL≈1015cm,is much more massive,MCSM10 M⊙,and is concentrated in the equatorial plane.

    4.10.The Rate of AT2018cow-like FBOTs

    Ho et al.(2022) estimate the rate of AT2018cow-like events to be ?0.01%–0.1% of the rate of CCSNe (see also Ho et al.2021).In Soker et al.(2019)we estimate the rate of all polar-CEJSN and impostor events to be 0.2%–0.5% of all CCSNe.The polar CEJSN impostors are a fraction of these,and so the expected rate is compatible with the estimate of Ho et al.(2022).The outcome of the polar CEJSN impostor scenario after the core explodes is a binary system of two NSs that might merge at a later time.The NS-NS merger rate is≈1% of all CCSNe (e.g.,see discussion by Mapelli et al.2018).If these rates hold,then AT2018cow-like events account for ≈10% of the progenitors of NS-NS mergers.However,the rate of all types of FBOTs is ≈1% of all CCSNe (e.g.,Coppejans et al.2020).As well,some CEJSN impostor channels that are similar,but not identical,to those I study here might account for other fast transients,e.g.,AT2018lqh that Ofek et al.(2021) observed.Tsuna et al.(2021) propose a scenario for AT2018lqh where a rotating blue supergiant collapses to form a BH of ?30 M⊙a(bǔ)nd blows a disk wind of?0.8 M⊙.I instead propose that AT2018lqh is a type of CEJSN impostor similar to the polar CEJSN impostors that explain AT2018cow-like transients.Ofek et al.(2021)estimate the rate of such fast transients to be of the order of magnitude of the rate of NS-NS mergers.

    I raise therefore the possibility that FBOTs and similar fast transients are CEJSN impostors which compose a large fraction of the progenitors of NS-NS merger binaries.

    4.11.Optical and Radio Emission

    I will neither study here the observed radio properties(e.g.,Huang et al.2019;Ho et al.2022) nor the exact optical emission,as these require more detailed radiative transfer calculations following hydrodynamical simulations.I referred to the light curve in Sections 3.3 and 4.6,and to the photosphere in Section 4.7.I also pointed out (Section 3) that I expect the emission properties to be similar to some of those in the studies of Margutti et al.(2019) and Gottlieb et al.(2022).

    The common properties of the polar CEJSN impostor scenario and of the geometrical model of Margutti et al.(2019) are the presence of a dense equatorial outflow,which can explain the late outflow velocity of ≈1000-several×1000 km s?1in AT2018cow,the possibility that the ejecta shocks this dense equatorial outflow to contribute to the X-ray emission,and the opening view along the polar directions that allows an observer to view the central engine.The common properties with the scenario that Gottlieb et al.(2022)propose are the hot cocoons along the two polar directions that explain the early optical emission and the jets along the polar directions that might contribute to the X-ray emission even at late times.

    To summarize this section I list the properties and processes that I discussed above in Table 1.I note that Soker et al.(2019)already discussed some of these properties and processes in the frame of the polar CEJSN scenario,i.e.,for which the NS enters the RSG core.

    Table 1 The Explanation of the Polar-CEJSN Impostor Scenario for the Different Properties of AT2018cow-like FBOTs

    Further hydrodynamical simulations together with radiative transfer calculations are required to determine the optical and radio properties of the polar CEJSN impostor scenario.

    5.Summary

    The puzzling properties of FBOTs,like fast rise and decline,X-ray variability,and fast outflows,have lead different researchers to propose several different theoretical scenarios to account for FBOTs,in particular AT2018cow-like FBOTs(see references in Section 1).In the present study I added another scenario,the polar CEJSN impostor scenario.I described the main evolutionary phases of this scenario in Section 2.

    The main new ingredients of the polar CEJSN impostor scenario are a pre-explosion bipolar CSM (Section 3.1,Figure 1),the post-CEE accretion from a hydrogen-rich circumbinary disk that feeds the accretion disk around the NS that launches the jets (Section 3.2),and the interaction of these jets with the lobes (Section 3.3;Figures 2–4).The main energy source of FBOTs in the scenario I propose is the gravitational energy of the accretion process onto a pre-existing NS that takes place immediately after the termination of the CEE.The NS accretes mass from an accretion disk that launches jets.The post-CEE circumbinary disk feeds the accretion disk for a time period of weeks(Section 2;Figure 1).The collision of the jets with the CSM gives rise to thermal energy and then radiation (Figures 2–4).

    In Section 4 I listed different properties of AT2018cow-like FBOTs and discussed the way by which the polar CEJSN impostor scenario might account for these properties(Table 1).

    At the end of the FBOT there is a bare system of an NS and the core of the RSG.The core later explodes,leading to the formation of a binary NS system that might be bound.Such a binary might much later experience NS-NS merger.In Section 4.10 I crudely estimated,based on the proposed scenario,that AT2018cow-like events are progenitors of ≈10%of the NS-NS mergers.From that I raised––the possibility that FBOTs and similar fast transients are CEJSN impostors that are progenitors of a large fraction of NS-NS merger binaries.

    Acknowledgments

    I thank Ore Gottlieb and Aldana Grichener for useful comments,and an anonymous referee for detailed comments that substantially improved the manuscript.This research was supported by the Amnon Pazy Research Foundation and the Asher Space Research Fund at the Technion.

    欧美日韩福利视频一区二区| 狠狠精品人妻久久久久久综合| 国产一区二区 视频在线| 中文字幕色久视频| 精品一区二区免费观看| 精品国产一区二区久久| av福利片在线| 91aial.com中文字幕在线观看| 十八禁人妻一区二区| 欧美日韩成人在线一区二区| 亚洲天堂av无毛| 国产精品 国内视频| 18禁裸乳无遮挡动漫免费视频| av卡一久久| 在线亚洲精品国产二区图片欧美| 国产一区二区激情短视频 | 国产毛片在线视频| 国产精品无大码| 国产亚洲最大av| 成年人午夜在线观看视频| 伊人久久国产一区二区| 国产伦人伦偷精品视频| 最近中文字幕高清免费大全6| 国产在视频线精品| 国产免费又黄又爽又色| 电影成人av| 日本猛色少妇xxxxx猛交久久| 2021少妇久久久久久久久久久| 久久人人97超碰香蕉20202| 国产成人一区二区在线| 久久性视频一级片| 国产亚洲欧美精品永久| 国产精品免费视频内射| 国产成人一区二区在线| 日本一区二区免费在线视频| 欧美乱码精品一区二区三区| 精品一区二区三卡| 国产黄色免费在线视频| 午夜免费鲁丝| 999久久久国产精品视频| 亚洲在久久综合| 制服诱惑二区| 国产免费视频播放在线视频| 亚洲少妇的诱惑av| 永久免费av网站大全| 成年女人毛片免费观看观看9 | 丝袜脚勾引网站| 国产又爽黄色视频| 亚洲精品视频女| 丰满乱子伦码专区| 亚洲av综合色区一区| 午夜福利,免费看| 一区二区三区精品91| 搡老乐熟女国产| 亚洲欧美日韩另类电影网站| 亚洲成人国产一区在线观看 | 青春草亚洲视频在线观看| 在线免费观看不下载黄p国产| 精品一品国产午夜福利视频| 国产黄频视频在线观看| 一二三四中文在线观看免费高清| 日韩制服丝袜自拍偷拍| 欧美xxⅹ黑人| 日韩制服骚丝袜av| 毛片一级片免费看久久久久| 日本91视频免费播放| 亚洲成人av在线免费| 精品一区二区三区四区五区乱码 | 久久久久精品性色| 免费在线观看完整版高清| 久久99精品国语久久久| 精品一区二区三区av网在线观看 | 中文天堂在线官网| 一级毛片电影观看| 如何舔出高潮| 精品国产一区二区三区久久久樱花| 亚洲伊人久久精品综合| 色婷婷久久久亚洲欧美| 国产一卡二卡三卡精品 | 老司机影院成人| 国产又爽黄色视频| 在线观看免费日韩欧美大片| 国产免费视频播放在线视频| tube8黄色片| 免费女性裸体啪啪无遮挡网站| 国产一区亚洲一区在线观看| 亚洲欧美日韩另类电影网站| 九草在线视频观看| 日日撸夜夜添| 国产成人a∨麻豆精品| 久久精品久久久久久久性| 久久久久久久大尺度免费视频| 日本一区二区免费在线视频| 国产精品一国产av| 又大又黄又爽视频免费| 激情五月婷婷亚洲| 黑人猛操日本美女一级片| 在线观看免费高清a一片| 久久久久网色| 悠悠久久av| videos熟女内射| 考比视频在线观看| 人人妻人人爽人人添夜夜欢视频| 国产熟女欧美一区二区| 大香蕉久久成人网| 精品午夜福利在线看| 久久综合国产亚洲精品| 99久久99久久久精品蜜桃| 老司机在亚洲福利影院| 人人妻,人人澡人人爽秒播 | av电影中文网址| 亚洲av综合色区一区| 国产野战对白在线观看| 亚洲国产毛片av蜜桃av| 国产又色又爽无遮挡免| 午夜福利免费观看在线| 国产成人精品久久二区二区91 | 日韩中文字幕视频在线看片| 国产成人欧美在线观看 | av福利片在线| 桃花免费在线播放| 男女高潮啪啪啪动态图| 亚洲国产欧美网| 9热在线视频观看99| 在线亚洲精品国产二区图片欧美| 亚洲国产毛片av蜜桃av| 99re6热这里在线精品视频| 精品福利永久在线观看| 操出白浆在线播放| 日韩欧美精品免费久久| 精品视频人人做人人爽| 观看av在线不卡| 精品国产国语对白av| 国产又爽黄色视频| 国产精品国产三级专区第一集| av国产精品久久久久影院| 国产 一区精品| 天美传媒精品一区二区| 欧美av亚洲av综合av国产av | 999久久久国产精品视频| 国产熟女欧美一区二区| 男女高潮啪啪啪动态图| 亚洲国产精品国产精品| 国产极品天堂在线| 亚洲国产最新在线播放| 亚洲精华国产精华液的使用体验| 国产一区二区在线观看av| 久久精品亚洲熟妇少妇任你| 亚洲,欧美精品.| 99热国产这里只有精品6| 日韩免费高清中文字幕av| 久久婷婷青草| 日本wwww免费看| 久久久久网色| 久久久久久久久久久免费av| 伦理电影免费视频| 亚洲少妇的诱惑av| 久久精品亚洲av国产电影网| 色播在线永久视频| 久久久欧美国产精品| 美女主播在线视频| 欧美日韩国产mv在线观看视频| 男女无遮挡免费网站观看| 黄片无遮挡物在线观看| 亚洲精品日本国产第一区| 宅男免费午夜| 久久久久视频综合| 欧美成人精品欧美一级黄| 叶爱在线成人免费视频播放| 国产一区二区 视频在线| 黑人巨大精品欧美一区二区蜜桃| 亚洲av中文av极速乱| 亚洲精品av麻豆狂野| tube8黄色片| 这个男人来自地球电影免费观看 | 在线看a的网站| 丰满迷人的少妇在线观看| 伊人久久大香线蕉亚洲五| 亚洲欧美精品自产自拍| 亚洲精品av麻豆狂野| 日日啪夜夜爽| 久久天躁狠狠躁夜夜2o2o | 男人舔女人的私密视频| 男女边摸边吃奶| 91老司机精品| 国产淫语在线视频| 国产麻豆69| av女优亚洲男人天堂| 高清视频免费观看一区二区| 国产精品免费大片| 亚洲av成人不卡在线观看播放网 | 男女之事视频高清在线观看 | 大片免费播放器 马上看| 精品一品国产午夜福利视频| 极品人妻少妇av视频| 欧美日韩精品网址| 欧美日韩亚洲综合一区二区三区_| 欧美国产精品一级二级三级| 在线观看人妻少妇| 免费在线观看黄色视频的| www日本在线高清视频| 一本一本久久a久久精品综合妖精| 操美女的视频在线观看| 韩国高清视频一区二区三区| 亚洲七黄色美女视频| 欧美日韩亚洲国产一区二区在线观看 | 人妻一区二区av| 天堂中文最新版在线下载| 国产精品秋霞免费鲁丝片| av天堂久久9| 久热爱精品视频在线9| 国产日韩欧美亚洲二区| 成人免费观看视频高清| 91aial.com中文字幕在线观看| 制服人妻中文乱码| 久久精品久久久久久久性| a级片在线免费高清观看视频| 免费高清在线观看视频在线观看| 国产精品 欧美亚洲| 国产亚洲av高清不卡| 亚洲视频免费观看视频| 日韩一本色道免费dvd| 久久ye,这里只有精品| 国产乱来视频区| 人妻一区二区av| 亚洲一卡2卡3卡4卡5卡精品中文| 日本爱情动作片www.在线观看| 晚上一个人看的免费电影| 国产一区二区三区综合在线观看| 美女主播在线视频| 日日爽夜夜爽网站| 欧美日韩亚洲高清精品| avwww免费| 成年美女黄网站色视频大全免费| 亚洲美女视频黄频| 无限看片的www在线观看| 赤兔流量卡办理| 少妇人妻 视频| 热99国产精品久久久久久7| 在线观看一区二区三区激情| av国产精品久久久久影院| 亚洲国产欧美网| 久久久久久久大尺度免费视频| 在线精品无人区一区二区三| 亚洲色图 男人天堂 中文字幕| 最近中文字幕高清免费大全6| 国产一级毛片在线| 黄色一级大片看看| 观看av在线不卡| 亚洲国产欧美网| 免费黄频网站在线观看国产| 51午夜福利影视在线观看| 一边亲一边摸免费视频| 欧美日韩av久久| 国产一卡二卡三卡精品 | 又粗又硬又长又爽又黄的视频| 国产免费一区二区三区四区乱码| videos熟女内射| 99精国产麻豆久久婷婷| 777米奇影视久久| 97人妻天天添夜夜摸| 人妻一区二区av| 亚洲伊人久久精品综合| 狠狠精品人妻久久久久久综合| 婷婷色麻豆天堂久久| 国产在视频线精品| 美女中出高潮动态图| 永久免费av网站大全| 男女高潮啪啪啪动态图| 亚洲av电影在线观看一区二区三区| 精品国产国语对白av| 丝袜喷水一区| 人成视频在线观看免费观看| av一本久久久久| 欧美亚洲日本最大视频资源| 伦理电影大哥的女人| 一本色道久久久久久精品综合| 2021少妇久久久久久久久久久| h视频一区二区三区| 亚洲精品久久久久久婷婷小说| 精品少妇久久久久久888优播| 18禁观看日本| 亚洲国产精品一区二区三区在线| 搡老乐熟女国产| 国产一级毛片在线| 大陆偷拍与自拍| 2018国产大陆天天弄谢| 欧美变态另类bdsm刘玥| 免费人妻精品一区二区三区视频| 最近中文字幕高清免费大全6| 老熟女久久久| av网站免费在线观看视频| 国产一区二区 视频在线| 老司机靠b影院| 欧美xxⅹ黑人| 人人妻人人添人人爽欧美一区卜| 久久国产精品大桥未久av| 波多野结衣一区麻豆| 亚洲精品国产一区二区精华液| 咕卡用的链子| 国产在线一区二区三区精| 精品免费久久久久久久清纯 | 国产精品香港三级国产av潘金莲 | 男人操女人黄网站| 卡戴珊不雅视频在线播放| 午夜福利视频在线观看免费| av网站免费在线观看视频| 久久影院123| 国产精品久久久久久久久免| 看十八女毛片水多多多| 亚洲国产欧美在线一区| 青春草国产在线视频| 色94色欧美一区二区| 色婷婷av一区二区三区视频| 中文乱码字字幕精品一区二区三区| 久久久精品免费免费高清| 亚洲天堂av无毛| 高清欧美精品videossex| av视频免费观看在线观看| 操出白浆在线播放| 免费日韩欧美在线观看| 精品人妻熟女毛片av久久网站| 中文精品一卡2卡3卡4更新| 国产亚洲av片在线观看秒播厂| 亚洲一卡2卡3卡4卡5卡精品中文| 女的被弄到高潮叫床怎么办| 国产欧美日韩综合在线一区二区| 亚洲国产成人一精品久久久| bbb黄色大片| 国产亚洲av高清不卡| 精品国产乱码久久久久久男人| 大片免费播放器 马上看| 亚洲国产精品999| 国产成人啪精品午夜网站| 天天影视国产精品| 十分钟在线观看高清视频www| 如何舔出高潮| 成年美女黄网站色视频大全免费| 成人国语在线视频| 国产熟女欧美一区二区| 黑人欧美特级aaaaaa片| 一区二区三区激情视频| 亚洲精品美女久久久久99蜜臀 | 亚洲精品aⅴ在线观看| 久久热在线av| 悠悠久久av| 51午夜福利影视在线观看| 我的亚洲天堂| 日本91视频免费播放| 一级,二级,三级黄色视频| 丝袜美足系列| 男人添女人高潮全过程视频| 免费人妻精品一区二区三区视频| 久久影院123| 亚洲成人免费av在线播放| 大话2 男鬼变身卡| 久久久久久久精品精品| 无遮挡黄片免费观看| 免费女性裸体啪啪无遮挡网站| 中文字幕最新亚洲高清| 亚洲精品国产区一区二| 成人毛片60女人毛片免费| 91精品伊人久久大香线蕉| 精品少妇黑人巨大在线播放| 久久久久人妻精品一区果冻| 天天躁狠狠躁夜夜躁狠狠躁| 久久性视频一级片| 日本wwww免费看| 男女国产视频网站| 久久精品久久久久久久性| 岛国毛片在线播放| 日韩中文字幕视频在线看片| 亚洲成人一二三区av| 女人久久www免费人成看片| 中文精品一卡2卡3卡4更新| 亚洲成人一二三区av| 成年美女黄网站色视频大全免费| 久久久久精品国产欧美久久久 | 亚洲一码二码三码区别大吗| 青春草亚洲视频在线观看| 在线观看人妻少妇| 午夜av观看不卡| 亚洲精品国产av成人精品| 一本一本久久a久久精品综合妖精| 一级毛片 在线播放| 麻豆精品久久久久久蜜桃| 最近中文字幕2019免费版| 综合色丁香网| 青春草视频在线免费观看| 日本wwww免费看| 国产精品一二三区在线看| 天堂中文最新版在线下载| xxxhd国产人妻xxx| 秋霞伦理黄片| e午夜精品久久久久久久| 亚洲av成人不卡在线观看播放网 | 999久久久国产精品视频| 久久精品国产综合久久久| 久久狼人影院| 国产精品秋霞免费鲁丝片| 日韩中文字幕欧美一区二区 | 亚洲精品美女久久av网站| 亚洲一级一片aⅴ在线观看| 婷婷成人精品国产| 免费女性裸体啪啪无遮挡网站| 久久精品久久精品一区二区三区| 日韩精品免费视频一区二区三区| 国产一区二区三区av在线| 99久久99久久久精品蜜桃| 免费少妇av软件| 波多野结衣一区麻豆| 亚洲精品日本国产第一区| 日本黄色日本黄色录像| 久久久国产一区二区| 日日啪夜夜爽| 亚洲av中文av极速乱| 一区二区三区乱码不卡18| 亚洲视频免费观看视频| 一本色道久久久久久精品综合| www.自偷自拍.com| 操美女的视频在线观看| 亚洲国产毛片av蜜桃av| 熟女少妇亚洲综合色aaa.| 精品一品国产午夜福利视频| 精品一区在线观看国产| 亚洲视频免费观看视频| 天堂8中文在线网| 91精品三级在线观看| 日本爱情动作片www.在线观看| 国产在线免费精品| 国语对白做爰xxxⅹ性视频网站| 亚洲七黄色美女视频| 热99国产精品久久久久久7| 最近最新中文字幕大全免费视频 | 在线观看人妻少妇| xxx大片免费视频| 精品一区二区三卡| 黑人猛操日本美女一级片| 国产午夜精品一二区理论片| 丰满迷人的少妇在线观看| 精品人妻熟女毛片av久久网站| 丁香六月欧美| 欧美激情高清一区二区三区 | 久热这里只有精品99| 2018国产大陆天天弄谢| 国产精品久久久人人做人人爽| 多毛熟女@视频| 免费看不卡的av| 久久热在线av| 亚洲av日韩精品久久久久久密 | 免费观看性生交大片5| 亚洲av日韩精品久久久久久密 | 妹子高潮喷水视频| 少妇人妻久久综合中文| 欧美97在线视频| 国产日韩一区二区三区精品不卡| 久久精品人人爽人人爽视色| 国产精品熟女久久久久浪| 又大又爽又粗| 九九爱精品视频在线观看| 成年女人毛片免费观看观看9 | 欧美最新免费一区二区三区| 下体分泌物呈黄色| 午夜av观看不卡| 97在线人人人人妻| 黑人猛操日本美女一级片| 十八禁网站网址无遮挡| 婷婷色综合www| 麻豆av在线久日| 黄网站色视频无遮挡免费观看| 超色免费av| 一边摸一边做爽爽视频免费| 丝袜在线中文字幕| 久久国产精品大桥未久av| 精品一区二区三区av网在线观看 | 国产欧美日韩一区二区三区在线| 青草久久国产| 一边摸一边做爽爽视频免费| 亚洲精品日韩在线中文字幕| 国产精品久久久久久精品电影小说| 成年美女黄网站色视频大全免费| 一级毛片我不卡| 人体艺术视频欧美日本| 亚洲国产最新在线播放| 老司机影院毛片| 嫩草影院入口| 美女高潮到喷水免费观看| 一区二区三区四区激情视频| 看免费成人av毛片| 亚洲精品av麻豆狂野| 日韩,欧美,国产一区二区三区| 精品少妇久久久久久888优播| 中文字幕人妻丝袜一区二区 | 中文字幕人妻丝袜一区二区 | 精品少妇黑人巨大在线播放| 亚洲精品久久午夜乱码| 亚洲精品美女久久av网站| av在线播放精品| 黄色一级大片看看| 啦啦啦视频在线资源免费观看| 19禁男女啪啪无遮挡网站| 啦啦啦中文免费视频观看日本| 国产一区二区三区综合在线观看| 欧美日韩视频高清一区二区三区二| a级片在线免费高清观看视频| 色吧在线观看| 亚洲精品国产区一区二| 91aial.com中文字幕在线观看| 天堂俺去俺来也www色官网| 国产精品av久久久久免费| 久久国产亚洲av麻豆专区| 亚洲成人国产一区在线观看 | 制服诱惑二区| 丝瓜视频免费看黄片| 爱豆传媒免费全集在线观看| 久热这里只有精品99| 黑人猛操日本美女一级片| 亚洲久久久国产精品| 日本av手机在线免费观看| 在线精品无人区一区二区三| 在线观看免费高清a一片| 在线观看免费视频网站a站| 在线观看国产h片| 国产不卡av网站在线观看| 成人漫画全彩无遮挡| 久久精品熟女亚洲av麻豆精品| 夫妻午夜视频| 亚洲第一青青草原| 校园人妻丝袜中文字幕| 色94色欧美一区二区| 2021少妇久久久久久久久久久| 久久久久久久久久久久大奶| 在现免费观看毛片| 日韩av不卡免费在线播放| 日韩一区二区三区影片| 欧美亚洲 丝袜 人妻 在线| 亚洲国产成人一精品久久久| 亚洲免费av在线视频| 日日啪夜夜爽| 高清av免费在线| 国产又色又爽无遮挡免| 校园人妻丝袜中文字幕| 老司机影院毛片| 美女国产高潮福利片在线看| 欧美亚洲 丝袜 人妻 在线| 国产日韩欧美视频二区| 麻豆av在线久日| 青春草国产在线视频| 亚洲av男天堂| 午夜福利视频在线观看免费| 在线观看三级黄色| 热99久久久久精品小说推荐| 色视频在线一区二区三区| 日韩熟女老妇一区二区性免费视频| 日韩欧美精品免费久久| 久久99一区二区三区| 在现免费观看毛片| 一级,二级,三级黄色视频| 亚洲精品在线美女| 又大又黄又爽视频免费| 国产深夜福利视频在线观看| 亚洲成国产人片在线观看| h视频一区二区三区| 亚洲情色 制服丝袜| 亚洲av日韩精品久久久久久密 | 9色porny在线观看| 黄色一级大片看看| 美国免费a级毛片| 伊人亚洲综合成人网| 久久精品久久久久久久性| 操出白浆在线播放| 日本91视频免费播放| 夫妻午夜视频| 高清视频免费观看一区二区| 中文字幕av电影在线播放| 成人亚洲精品一区在线观看| 我要看黄色一级片免费的| 国产黄色免费在线视频| 国产欧美日韩一区二区三区在线| 国产国语露脸激情在线看| 国产又爽黄色视频| 精品第一国产精品| 亚洲成av片中文字幕在线观看| 久久精品国产综合久久久| 亚洲一码二码三码区别大吗| 亚洲成人一二三区av| 色婷婷av一区二区三区视频| 久久狼人影院| 日韩成人av中文字幕在线观看| 高清视频免费观看一区二区| 精品国产国语对白av| 国产精品 国内视频| 色婷婷久久久亚洲欧美| 欧美精品一区二区免费开放| 女人被躁到高潮嗷嗷叫费观| 欧美亚洲 丝袜 人妻 在线| 丰满饥渴人妻一区二区三| 麻豆av在线久日| 麻豆乱淫一区二区| 又大又爽又粗| 999久久久国产精品视频| 精品午夜福利在线看| 亚洲第一青青草原| 观看美女的网站| 最近中文字幕2019免费版| 最近2019中文字幕mv第一页| 日本av免费视频播放| 新久久久久国产一级毛片| 久久久久精品国产欧美久久久 | 一级毛片电影观看| 大香蕉久久网| 欧美日韩国产mv在线观看视频| 欧美激情高清一区二区三区 | 18在线观看网站| av线在线观看网站|