Variable Stars in the 50BiN Open Cluster Survey.Ⅲ.NGC 884

Open clusters are the basic building blocks that serve as a laboratory for the study of young stellar populations in the Milky Way.Variable stars in open clusters provide a unique way to accurately probe the internal structure,temporal and dynamical evolutionary stages of individual stars and the host cluster.The most powerful tool for such studies is time-domain photometric observations.This paper follows the route of our previous work,concentrating on a photometric search for variable stars in NGC 884.The target cluster is the companion of NGC869,forming the well-known double cluster system that is gravitationally bound.From the observation run in 2016 November,a total of 9247 B-band CCD images and 8218Ⅴ-band CCD images were obtained.We detected a total of 15 stars with variability in visual brightness,including five Be stars,three eclipsing binaries,and seven of unknown types.Two new variable stars were discovered in this work.We also compared the variable star content of NGC 884 with its companion NGC 869.

Revise Thermal Winds of Remnant Neutron Stars in Gamma-Ray Bursts

It seems that the wealth of information revealed by the multi-messenger observations of the binary neutron star(NS)merger event,GW170817/GRB 170817A/kilonova AT2017gfo,places irreconcilable constraints to models of the prompt emission of this gamma-ray burst(GRB).The observed time delay between the merger of the two NSs and the trigger of the GRB and the thermal tail of the prompt emission can hardly be reproduced by these models simultaneously.We argue that the merger remnant should be an NS(last for,at least,a large fraction of 1 s),and that the difficulty can be alleviated by the delayed formation of the accretion disk due to the absorption of high-energy neutrinos emitted by the NS and the delayed emergence of effective viscosity in the disk.Further,we extend the consideration of the effect of the energy deposition of neutrinos emitted from the NS.If the NS is the central object of a GRB with a distance and duration similar to that of GRB 170817A,thermal emission of the thermal bubble inflated by the NS after the termination of accretion may be detectable.If our scenario is verified,it would be of interest to investigate the cooling of nascent NSs.

以障碍点为指导的个性化综合训练对老年全髋关节置换伴肌少症患者运动能力及预后的影响

目的 探讨以障碍点为指导的个性化综合训练用于老年全髋关节置换伴肌少症围术期的康复效果。方法 选取2020年2月—2023年4月收治的84例老年全髋关节置换伴肌少症,根据康复训练方法不同分为对照组、观察组,各42例。对照组给予常规康复训练,观察组给予以障碍点为指导的个性化综合训练。比较2组术前及术后1个月Fugl-Meyer运动功能量表(FMA)、Berg平衡量表(BBS)、髋关节残疾和骨关节炎预后评分(HOOS)、肌力相关指标、肌肉相关康复指标及肺功能[用力肺活量(FVC)、第1秒用力呼气容积(FEV1)、FEV1/FVC、FEV1占预计值的百分比(FEV1%pred)]。结果 术后1个月,观察组FMA、BBS、HOOS评分较对照组升高,优势手握力、四肢骨骼肌质量、患肢骨骼肌质量和全身肌肉质量均较对照组升高,胫前屈肌厚度、肱二头肌厚度较对照组均降低(P<0.05);术后1个月,观察组FVC、FEV1、FEV1/FVC、FEV1%pred与对照组比较均升高(P<0.05)。结论 以障碍点为指导的个性化综合训练能提高老年全髋关节置换伴肌少症患者肌力与运动能力,利于改善肺功能及患者预后。

基于“三因制宜”理念的中医辨证护理在全髋关节置换术患者中的应用效果分析

目的:观察“三因制宜”理念的中医辨证护理在全髋关节置换术患者中的应用效果。方法:将88例全髋关节置换术(total hip arthroplasty,THA)患者按随机数字表法分为对照组和三因组各44例。对照组给予西医常规护理,三因组在对照组的基础上给予“三因制宜”理念的中医辨证护理。观察两组患者术前及术后6天视觉模拟量表(visual analogue scale,VAS)和汉密尔顿焦虑量表(hamilton anxiety scale,HAMA)评分、术后7天炎症指标[白细胞介素6(interleukin-6,IL-6)、肿瘤坏死因子α(tumor necrosis factorα,TNF-α)]、术后5个月Harris髋关节评分、并发症及护理满意度。结果:术后6天两组患者VAS和HAMA评分均低于术前(P<0.05),三因组低于对照组(P<0.05);术后7天两组患者血清IL-6、TNF-α水平均高于术前(P<0.05),三因组低于对照组(P<0.05);术后5个月,三因组患者Harris髋关节疼痛评分低于对照组(P<0.05),活动度、功能评分均高于对照组(P<0.05),畸形评分与对照组比较差异无统计学意义(P>0.05);三因组患者并发症总发生率为2.27%(1/44),低于对照组的13.64%(6/44)(P<0.05);三因组患者总满意度为95.45%(42/44),高于对照组的79.55%(35/44)(P<0.05)。结论:基于“三因制宜”理念的中医辨证护理可有效缓解THA患者疼痛,抑制炎症反应、改善髋关节功能,降低术后并发症,提高护理满意度。

Analysis of the Impact of the Blazhko Effect Both on the Van Hoof Effect and Radial Velocity Amplitude in the Star RR Lyr

The Van Hoof effect is a phase shift existing between the radial velocity curves of hydrogen and metallic lines within the atmosphere of pulsating stars.In this article,we present a study of this phenomenon through the spectra of the brightest pulsating star RR Lyr of RR Lyrae stars recorded for 22 yr.We based ourselves,on the one hand,on 1268 spectra(41 nights of observation)recorded between the years 1994 and 1997 at the Observatory of Haute Provence(OHP,France)previously observed by Chadid and Gillet,and on the other hand on 1569 spectra(46nights of observation)recorded at our Oukaimeden Observatory(Morocco)between 2015 and 2016.Through this study,we have detected information on atmospheric dynamics that had not previously been detected.Indeed,the Van Hoof effect which results in a clear correlation between the radial velocities of hydrogen and those of the metallic lines has been observed and analyzed at different Blazhko phases.A correlation between the radial velocities of different metallic lines located in the lower atmosphere has been observed as well.For the first time,we were able to show that the amplitude of the radial velocity curves deduced from the lines of hydrogen and that of FeⅡ(λ4923.921?)increases toward the minimum of the Blazhko cycle and decreases toward the maximum of the same Blazhko cycle.Furthermore,we found that the Van Hoof effect is also modulated by the Blazhko effect.Thus,toward the minimum of the Blazhko cycle the Van Hoof effect is more visible and at the maximum of the Blazhko cycle,this effect is minimal.We also observed the temporal evolution of the amplitudes of the radial velocities of the lower and upper atmosphere.When observed over a long time,we can interpret it as a function of the Blazhko phases.

First Light Curve Analysis of NSVS 8294044,V1023 Her,and V1397 Her Contact Binary Systems

The first photometric light curve investigation of the NSVS 8294044,V1023 Her,and V1397 Her binary systems is presented.We used ground-based observations for the NSVS 8294044 system and Transiting Exoplanet Survey Satellite data for V1023 Her and V1397 Her.The primary and secondary times of minima were extracted from al the data,and,by collecting the literature,a new ephemeris was computed for each system.Linear fits for the O-C diagrams were conducted using the Markov Chain Monte Carlo (MCMC) method.Light curve solutions were performed using the PHysics Of Eclipsing BinariEs Python code and the MCMC approach.The systems were found to be contact binary stars based on the fillout factor and mass ratio.V1023 Her showed the O’Connell effect and a cold starspot on the secondary component was required for the light curve solution.The absolute parameters of the system were estimated based on an empirical relationship between orbital period and mass.We presented a new T–M equation based on a sample of 428 contact binary systems and found that our three target systems were in good agreement with the fit.The positions of the systems were also depicted on the M–L,M–R,q–L_(ratio),and M_(tot)–J_(0)diagrams in the logarithmic scales.

GRB 200612A:An Ultralong Gamma-Ray Burst Powered by Magnetar Spinning Down

GRB 200612A could be classified as an ultralong gamma-ray burst due to its prompt emission lasting up to~1020 s and the true timescale of the central engine activity t_(burst)≥4×10^(4) s.The late X-ray light curve with a decay index ofα=7.53 is steeper than the steepest possible decay from an external shock model.We propose that this X-ray afterglow can be driven by dipolar radiation from the magnetar spindown during its early stage,while the magnetar collapsed into the black hole before its spindown,resulting in a very steep decay of the late X-ray light curve.The optical data show that the light curve is still rising after 1.1 ks,suggesting a late onset.We show that GRB 200612A’s optical afterglow light curve is fitted with the forward shock model by Gaussian structured off-axis jet.This is a special case among GRBs,as it may be an ultralong gamma-ray burst powered by a magnetar in an off-axis observation scenario.

V0405 Dra:A New Deep and Low Mass Ratio Contact Binary with Extremely Fast Decrease in the Orbital Period

V0405 Dra is a W UMa-type binary star.Based on the TESS data,we have conducted an orbital period study and performed a light curve analysis for the system.The orbital period study reveals that the O-C curve for V0405 Dra exhibits secular decrease at an extremely high rate of d P/dt=-2.71×10^(-6)day year^(-1),along with periodic variations characterized by an amplitude of A_(3)=0.0032 day and a period of P_(3)=1.413 years.The orbital periodic change is possibly due to the light-travel time effect resulting from an additional third body in the system,for which we estimate a minimum mass of M_(3)=0.77M_(⊙).By employing the 2013 version of the Wilson-Devinney(W-D)method to synthesize a light curve,we derived photometric solutions indicating that V0405 Dra is a new deep(f=68.7%)and low-mass ratio(q=0.175)contact binary.The fast decrease in its orbital period is likely caused by mass transfer from the more massive primary star to the less massive secondary star,or due to angular momentum loss.With further mass transfer and loss of angular momentum,the binary will gradually evolve into a tighter contact configuration,eventually leading to a merger into a single star,following the evolutionary paths suggested for such deep and low mass ratio contact binaries.

A Study of 10 Rotating Radio Transients Using Parkes Radio Telescope

Rotating Radio Transients(RRATs)are a relatively new subclass of pulsars that emit detectable radio bursts sporadically.We analyzed 10 RRATs observed using the Parkes telescope,with eight of these observed via the ultra-wide-bandwidth low-frequency(UWL)receiver.We measured the burst rate and produced integrated profiles spanning multiple frequency bands for three RRATs.We also conducted a spectral analysis on both integrated pulses and individual pulses of three RRATs.All of their integrated pulses follow a simple power law,consistent with the known range of pulsar spectral indices.Their average spectral indices of single pulses are-0.9,-1.2,and-1.0 respectively,which are within the known range of pulsar spectral indices.Additionally,we find that the spreads of single-pulse spectral indices for these RRATs(ranging from-3.5 to+0.5)are narrower compared to what has been observed in other RRATs.Notably,the average spectral index and scatter of single pulses are both relatively small.For the remaining five RRATs observed at the UWL receiver,we also provide the upper limits on fluence and flux density.In addition,we obtain the timing solution of PSR J1709-43.Our analysis shows that PSRs J1919+1745,J1709-43,and J1649-4653 are potentially nulling pulsars or weak pulsars with sparse strong pulses.

BSN:Photometric Light Curve Analysis of Two Contact Binary Systems LS Del and V997 Cyg

The light curve analyses and orbital period variations for two contact binary stars,LS Del and V997 Cyg,are presented in this work which was conducted in the frame of the Binary Systems of South and North project Ground-based photometric observations were performed at two observatories in France.We used the Transiting Exoplanet Survey Satellite(TESS)data for extracting times of minima and light curve analysis of the targe systems.The O-C diagram for both systems displays a parabolic trend.LS Del and V997 Cyg’s orbital periods are increasing at rates of dP/dt=7.20×10^(-08)days yr^(-1)and dP/dt=2.54×10^(-08)days yr^(-1),respectively Therefore,it can be concluded that mass is being transferred from the less massive star to the more massive component with a rate of dM/dt=-1.96×10^(-7)M_(⊙)yr^(-1)for the LS Del system,and dM/dt=-3.83×10^(-7)M_(⊙)yr-1for V997 Cyg.The parameters of a third possible object in the system are also considered.The PHysics Of Eclipsing BinariEs Python code was used to analyze the light curves.The light curve solutions needed a cold starspot due to the asymmetry in the LS Del system’s light curve maxima.The mass ratio fill-out factor,and star temperature all indicate that both systems are contact binary types in this investigation.Two methods were applied to estimate the absolute parameters of the systems:one method relied on the parallax of Gaia DR3,and the other used a P-M relationship.The positions of the systems are also depicted on the M-L,M-R q-L_(ratio),and logM_(tot)-logJ_(0)diagrams.We recommend that further observations and investigations be done on the existence of a fourth body in this system.

Spin Evolution of the Magnetar SGR J1935+2154

Fast radio bursts(FRBs)are short pulses observed in radio frequencies usually originating from cosmological distances.The discovery of FRB 200428 and its X-ray counterpart from the Galactic magnetar SGR J1935+2154suggests that at least some FRBs can be generated by magnetars.However,the majority of X-ray bursts from magnetars are not associated with radio emission.The fact that only in rare cases can an FRB be generated raises the question regarding the special triggering mechanism of FRBs.Here we report long time spin evolution of SGR J1935+2154 until the end of 2022.According to v and v,the spin evolution of SGR J1935+2154 could be divided into two stages.The first stage evolves relatively steady evolution until 2020 April 27.After the burst activity in2020,the spin of SGR J1935+2154 shows strong variations,especially for v.After the burst activity in 2022October,a new spin-down glitch with△v/v=(-7.2±0.6)×10^(-6)is detected around MJD 59876,which is the second event in SGR J1935+2154.At the end,spin frequency and pulse profile do not show variations around the time of FRB 200428 and radio bursts 221014 and 221021,which supply strong clues to constrain the trigger mechanism of FRBs or radio bursts.

The Bright Single Pulse Emission from PSR B1133+16

We have conducted a comprehensive investigation into the bright single pulse emission from PSR B1133+16using the Giant Metrewave Radio Telescope.High time resolution data(61μs)were obtained at a center frequency of 322 MHz with a bandwidth of 32 MHz over a continuous observation period of 7.45 hr.A total of 1082 bright pulses were sporadically detected with peak flux densities ranging from 10 to 23 times stronger than the average pulse profile.However,no giant pulse-like emission with a relative pulse energy larger than 10 and extremely short duration was detected,indicating that these bright pulses cannot be categorized as giant pulse emission.The majority of these bright pulses are concentrated in pulse phases at both the leading and trailing windows of the average pulse profile,with an occurrence ratio of approximately 2.74.The pulse energy distribution for all individual pulses can be described by a combination of two Gaussian components and a cutoff power-law with an index of α=-3.2.An updated nulling fraction of 15.35%±0.45% was determined from the energy distribution.The emission of individual pulses follows a log-normal distribution in peak flux density ratio.It is imperative that regular phase drifting in bright pulse sequence is identified in both the leading and trailing components for the first time.Possible physical mechanisms are discussed in detail to provide insights into these observations.

Study of Complex Nitrogen and Oxygen-bearing Molecules toward the High-mass Protostar IRAS 18089–1732

The observation of oxygen(O)-and nitrogen(N)-bearing molecules gives an idea about the complex prebiotic chemistry in the interstellar medium.Recent millimeter and submillimeter wavelength observations have shown the presence of complex O-and N-bearing molecules in the star formation regions.So,the investigation of those molecules is crucial to understanding the chemical complexity in the star-forming regions.In this article,we present the identification of the rotational emission lines of N-bearing molecules ethyl cyanide(C_(2)H_(5)CN)and cyanoacetylene(HC_(3)N),and O-bearing molecule methyl formate(CH_(3)OCHO)toward high-mass protostar IRAS18089–1732 using the Atacama Compact Array.We also detected the emission lines of both the N-and O-bearing molecule formamide(NH_(2)CHO)in the envelope of IRAS 18089–1732.We have detected the v=0 and 1 state rotational emission lines of CH_(3)OCHO.We also detected the two vibrationally excited states of HC_(3)N(v7=1 and v7=2).The estimated fractional abundances of C_(2)H_(5)CN,HC_(3)N(v7=1),HC_(3)N(v7=2),and NH_(2)CHO toward IRAS 18089–1732 are(1.40±0.5)×10^(-10),(7.5±0.7)×10^(-11),(3.1±0.4)×10^(-11),and(6.25±0.82)×10^(-11)respectively.Similarly,the estimated fractional abundances of CH_(3)OCHO(v=0)and CH_(3)OCHO(v=1)are(1.90±0.9)×10^(-9)and(8.90±0.8)×10^(-10),respectively.We also created the integrated emission maps of the detected molecules,and the observed molecules may have originated from the extended envelope of the protostar.We show that C_(2)H_(5)CNand HC_(3)N are most probably formed via the subsequential hydrogenation of the CH_(2)CHCNand the reaction between C_(2)H_(2)and CN on the grain surface of IRAS 18089–1732.We found that NH_(2)CHO is probably produced due to the reaction between NH_(2)and H_(2)CO in the gas phase.Similarly,CH_(3)OCHO is possibly created via the reaction between radical CH_(3)O and radical HCO on the grain surface of IRAS 18089–1732.

A Possible γ-Ray Pulsation from PSR J1740-5340B in the Globular Cluster NGC 6397

Recently, a new radio millisecond pulsar(MSP) J1740-5340B, hosted in the globular cluster(GC) NGC 6397,was reported with a 5.78 ms spin period in an eclipsing binary system with a 1.97 days orbital period. Based on a modified radio ephemeris updated by tool tempo2, we analyze the ~15 yr γ-ray data obtained from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope and detect PSR J1740-5340B's γ-ray pulsation at a confidence level of ~4σ with a weighted H-test value of ~26. By performing a phase-resolved analysis, the γ-ray luminosity in on-pulse interval of PSR J1740-5340B is L_(γ)~ 3.8 × 10^(33) erg s^(-1) using NGC 6397's distance of 2.48 kpc. And γ-rays from the on-pulse part of PSR J1740-5340B contribute ~90% of the total observed γ-ray emissions from NGC 6397. No significant γ-ray pulsation of another MSP J1740-5340A in the GC is detected.Considering that the previous four cases of MSPs in GCs, more data in γ-ray, X-ray, and radio are encouraged to finally confirm the γ-ray emissions from MSP J1740-5340B, especially starving for a precise ephemeris.

The RFI Fast Mitigation Algorithm Based on Block LMS Filter

The radio telescope possesses high sensitivity and strong signal collection capabilities.While receiving celestial radiation signals,it also captures Radio Frequency Interferences(RFIs)introduced by human activities.RFI,as signals originating from sources other than the astronomical targets,significantly impacts the quality of astronomical data.This paper presents an RFI fast mitigation algorithm based on block Least Mean Square(LMS)algorithm.It enhances the traditional adaptive LMS filter by grouping L adjacent time-sampled points into one block and applying the same filter coefficients for filtering within each block.This transformation reduces multiplication calculations and enhances algorithm efficiency by leveraging the time-domain convolution theorem.The algorithm is tested using baseband data from the Parkes 64 m radio telescope's pulsar observations and simulated data.The results confirm the algorithm's effectiveness,as the pulsar profile after RFI mitigation closely matches the original pulsar profile.

Isochrone Fitting of Galactic Globular Clusters—Ⅵ.High-latitude Clusters NGC 5024 (M53),NGC 5053,NGC 5272 (M3),NGC 5466,and NGC7099 (M30)

We fit various color–magnitude diagrams(CMDs) of the high-latitude Galactic globular clusters NGC 5024(M53),NGC 5053,NGC 5272(M3),NGC 5466,and NGC 7099(M30) by isochrones from the Dartmouth Stellar Evolution Database and Bag of Stellar Tracks and Isochrones for α–enrichment [α/Fe] = +0.4.For the CMDs,we use data sets from Hubble Space Telescope,Gaia,and other sources utilizing,at least,25 photometric filters for each cluster.We obtain the following characteristics with their statistical uncertainties for NGC 5024,NGC 5053,NGC 5272,NGC 5466,and NGC 7099,respectively:metallicities [Fe/H] =-1.93 ± 0.02,-2.08 ± 0.03,-1.60 ± 0.02,-1.95 ± 0.02,and-2.07 ± 0.04 dex with their systematic uncertainty 0.1 dex;ages 13.00 ± 0.11,12.70 ± 0.11,11.63 ± 0.07,12.15 ± 0.11,and 12.80 ± 0.17 Gyr with their systematic uncertainty 0.8 Gyr;distances(systematic uncertainty added) 18.22 ± 0.06 ± 0.60,16.99 ± 0.06 ± 0.56,10.08 ± 0.04 ± 0.33,15.59 ±0.03 ± 0.51,and 8.29 ± 0.03 ± 0.27 kpc;reddenings E(B-V) = 0.023 ± 0.004,0.017 ± 0.004,0.023 ± 0.004,0.023 ± 0.003,and 0.045 ± 0.002 mag with their systematic uncertainty 0.01 mag;extinctions AV= 0.08 ± 0.01,0.06 ± 0.01,0.08 ± 0.01,0.08 ± 0.01,and 0.16 ± 0.01 mag with their systematic uncertainty 0.03 mag,which suggest the total Galactic extinction AV= 0.08 across the whole Galactic dust to extragalactic objects at the North Galactic Pole.The horizontal branch morphology difference of these clusters is explained by their different metallicity,age,mass-loss efficiency,and loss of low-mass members in the evolution of the core-collapse cluster NGC 7099 and loose clusters NGC 5053 and NGC 5466.

Trabecular Plate Loss and Deteriorating Elastic Modulus of Femoral Trabecular Bone in Intertrochanteric Hip Fractures

Osteoporotic hip fracture is associated with significant trabecular bone loss, which is typically characterized as low bone density by dual-energy X-ray absorptiometry （DXA） and altered microstructure by micro-computed tomography （pCT）. Emerging morphological analysis techniques, e.g. individual trabecula segmentation （ITS）, can provide additional insights into changes in plate-like and rod-like trabeculae, two major micro- structural types serving different roles in determining bone strength. Using ITS, we evaluated trabecular microstructure of intertrochanteric bone cores obtained from 23 patients undergoing hip replacement surgery for intertrochanteric fracture and 22 cadaveric controls. Micro-finite element （~FE） analyses were performed to further understand how the abnormalities seen by ITS might translate into effects on bone strength. ITS analyses revealed that, near fracture site, plate-like trabeculae were seriously depleted in fracture patients, but trabecular rod volume was maintained. Besides, decreased plate area and rod length were observed in fracture patients. Fracture patients also showed decreased elastic moduli and shear moduli of trabecular bone. These results provided evidence that in intertrochanteric hip fracture, preferential loss of plate-like trabeculae led to more rod-like microstructure and deteriorated mechanical competence adjacent to the fracture site, which increased our understanding of the biomechanical pathogenesis of hip fracture in osteoporosis.

Physical and geometrical parameters of VCBS XIII: HIP 105947

The best physical and geometrical parameters of the main sequence close visual binary system(CVBS), HIP 105947, are presented. These parameters have been constructed conclusively using Al-Wardat’s complex method for analyzing CVBSs, which is a method for constructing a synthetic spectral energy distribution(SED) for the entire binary system using individual SEDs for each component star. The model atmospheres are in its turn built using the Kurucz(ATLAS9) line-blanketed plane-parallel models. At the same time, the orbital parameters for the system are calculated using Tokovinin’s dynamical method for constructing the best orbits of an interferometric binary system. Moreover, the mass-sum of the components, as well as the ?θ and ?ρ residuals for the system, is introduced. The combination of Al-Wardat’s and Tokovinin’s methods yields the best estimations of the physical and geometrical parameters. The positions of the components in the system on the evolutionary tracks and isochrones are plotted and the formation and evolution of the system are discussed.

Physical and orbital properties of the stellar system HIP 43766

In this paper,we present the analysis of the stellar binary system HIP 43766 to determine its properties.We rely on dynamical modeling and atmospheric modeling with recent data to determine the orbital solution and the physical properties of the system.There is a consistency between observed and synthetic photometry obtained using atmospheric modeling.The calculated dynamical mass sum of the system ranged between 1.691 and 2.609 solar masses,while it ranges between 2.0 and 2.1 as estimated utilizing atmospheric modeling.This could be due to inaccuracy in estimating the orbit,which could be modified with future observations with more relative positional measurements.The parameters of the system and the position of the components on the evolutionary tracks show that the system consists of F5 and G5 subgiant stars,mostly formed by fragmentation.A dynamical mass sum is predicted for the system.

The stellar system HIP 101227:is it a binary,a triple or a quadruple system?

In this paper,we present the analysis of the stellar system HIP 101227 to determine the actual number of components in the system,and their properties.We use dynamical modeling and complex spectrophotometric(involving atmospheric modeling)techniques with recent data,to determine the physical properties and orbital solution for the system,respectively,with better accuracy than past studies.Based on our analysis,we found that the system is more consistent with being a quadruple rather than a binary or a triple system as suggested by previous studies.The total mass of the system determined from our SED analysis is 3.42±0.20 M,which is distributed almost equally between the four stars.The stars are found to be zero-age main sequence stars;i.e.,at the last stage of pre-main sequence,with age less than 200 Myr and spectral types K0.All four stars have very similar physical characteristics,suggesting that the fragmentation process is the most likely theory for the formation and evolution of the system.