Synthetic Light Curves Research Articles

Robust Transiting Exoplanet Radii in the Presence of Starspots from Ingress and Egress Durations

Abstract We typically measure the radii of transiting exoplanets from the transit depth, which are given by the ratio of cross-sectional areas of the planet and star. However, if a star has dark starspots (or bright regions) distributed throughout the transit chord, the transit depth will be biased toward smaller (larger) values, and thus the inferred planet radius will be smaller (larger) if these are unaccounted for. We reparameterize the transit light curve to account for “self-contamination” by photospheric inhomogeneities by splitting the parameter R p /R ⋆ into two parameters: one for the radius ratio, which controls the duration of ingress and egress, and another which measures the possibly contaminated transit depth. We show that this is equivalent to the formulation for contamination by a second star (with positive or negative flux), and that it is sensitive to time-steady inhomogeneity of the stellar photosphere. We use synthetic light curves of spotted stars at high signal-to-noise to show that the radius recovered from measurement of the ingress/egress duration can recover the true radii of planets transiting spotted stars with axisymmetric spot distributions if the limb-darkening parameters are precisely known. We fit time-averaged high signal-to-noise transit light curves from Kepler and Spitzer of 10 planets to measure the planet radii and search for evidence of spot distributions. We find that this sample has a range of measured depths and ingress durations that are self-consistent, providing no strong evidence for contamination by spots. However, there is suggestive evidence for occultation of starspots on Kepler-17, and that relatively bright regions are occulted by the planets of Kepler-412 and HD 80606. Future observations with the James Webb Space Telescope may enable this technique to yield accurate planetary radii in the presence of stellar inhomogeneities.

On the use of wavelets to reveal oscillatory patterns in stellar flare emission.

Wavelet analysis is a powerful tool to investigate non-stationary signals such as amplitude modulated sinusoids or single events lasting for a small percentage of the observing time. Wavelet analysis can be used, for example, to reveal oscillations in the light curve of stars during coronal flares. A careful treatment of the background in the wavelet scalogram is necessary to determine robust confidence levels required to distinguish between patterns caused by actual oscillations and noise. This work describes the method using synthetic light curves and investigates the effect of background noise when determining confidence levels in the scalogram. The result of this analysis shows that the wavelet transform is able to reveal oscillatory patterns even when frequency-dependent noise is dominant. However, their significance in the wavelet scalogram may be reduced, depending on the assumed background spectrum. To show the power of wavelet analysis, the light curve of a well-known flaring star is analysed. It shows two oscillations overlapped. The lower-frequency oscillation is not mentioned in previous works in the literature. This result demonstrates the need for correctly characterizing the background noise of the signal.This article is part of the theme issue 'Redundancy rules: the continuous wavelet transform comes of age'.

The Influence of Accretion Disk Thickness on the Large-scale Magnetic Dynamo

Abstract The evolution of the magnetic field from the large-scale dynamo is considered a central feature of the accretion disk around a black hole. The resulting low-frequency oscillations introduced from the growth and decay of the field strength, along with the change in field orientation, play an integral role in the accretion disk behavior. Despite the importance of this process and how commonly it is invoked to explain variable features, it still remains poorly understood. We present a study of the dynamo using a suite of four global, high-resolution, MHD accretion disk simulations. We systematically vary the scale height ratio and find the large-scale dynamo fails to organize above a scale height ratio of h/r ≳ 0.2. Using spacetime diagrams of the azimuthal magnetic field, we show the large-scale dynamo is well ordered in the thinner accretion disk models, but fails to develop the characteristic “butterfly” pattern when the scale height ratio is increased, a feature which is also reflected in the power spectra. Additionally, we calculate the dynamo α-parameter and generate synthetic light curves. Using an emission proxy, we find the disks have markedly different characters as stochastic photometric fluctuations have a larger amplitude when the dynamo is unordered.

A GRB and Broad-lined Type Ic Supernova from a Single Central Engine

Abstract Unusually high velocities (≳0.1c) and correspondingly high kinetic energies have been observed in a subset of Type Ic supernovae (so-called “broad-lined Ic” supernovae; SNe Ic-BL), prompting a search for a central engine model capable of generating such energetic explosions. A clue to the explosion mechanism may lie in the fact that all supernovae that accompany long-duration gamma-ray bursts (GRBs) belong to the SN Ic-BL class. Using a combination of two-dimensional relativistic hydrodynamics and radiation transport calculations, we demonstrate that the central engine responsible for long GRBs can also trigger an SN Ic-BL. We find that a reasonable GRB engine injected into a stripped Wolf–Rayet progenitor produces a relativistic jet with energy ∼1051 erg, as well as an SN whose synthetic light curves and spectra are fully consistent with observed SNe Ic-BL during the photospheric phase. As a result of the jet’s asymmetric energy injection, the SN spectra and light curves depend on viewing angle. The impact of viewing angle on the spectrum is particularly pronounced at early times, while the viewing-angle dependence for the light curves (∼10% variation in bolometric luminosity) persists throughout the photospheric phase.

Robust Period Estimation Using Mutual Information for Multiband Light Curves in the Synoptic Survey Era

The Large Synoptic Survey Telescope (LSST) will produce an unprecedented amount of light curves using six optical bands. Robust and efficient methods that can aggregate data from multidimensional sparsely-sampled time series are needed. In this paper we present a new method for light curve period estimation based on the quadratic mutual information (QMI). The proposed method does not assume a particular model for the light curve nor its underlying probability density and it is robust to non-Gaussian noise and outliers. By combining the QMI from several bands the true period can be estimated even when no single-band QMI yields the period. Period recovery performance as a function of average magnitude and sample size is measured using 30,000 synthetic multi-band light curves of RR Lyrae and Cepheid variables generated by the LSST Operations and Catalog simulators. The results show that aggregating information from several bands is highly beneficial in LSST sparsely-sampled time series, obtaining an absolute increase in period recovery rate up to 50%. We also show that the QMI is more robust to noise and light curve length (sample size) than the multiband generalizations of the Lomb Scargle and Analysis of Variance periodograms, recovering the true period in 10-30% more cases than its competitors. A python package containing efficient Cython implementations of the QMI and other methods is provided.

RAPTOR

Context. Observational efforts to image the immediate environment of a black hole at the scale of the event horizon benefit from the development of efficient imaging codes that are capable of producing synthetic data, which may be compared with observational data. Aims. We aim to present RAPTOR, a new public code that produces accurate images, animations, and spectra of relativistic plasmas in strong gravity by numerically integrating the equations of motion of light rays and performing time-dependent radiative transfer calculations along the rays. The code is compatible with any analytical or numerical spacetime. It is hardware-agnostic and may be compiled and run both on GPUs and CPUs. Methods. We describe the algorithms used in RAPTOR and test the code’s performance. We have performed a detailed comparison of RAPTOR output with that of other radiative-transfer codes and demonstrate convergence of the results. We then applied RAPTOR to study accretion models of supermassive black holes, performing time-dependent radiative transfer through general relativistic magneto-hydrodynamical (GRMHD) simulations and investigating the expected observational differences between the so-called fast-light and slow-light paradigms. Results. Using RAPTOR to produce synthetic images and light curves of a GRMHD model of an accreting black hole, we find that the relative difference between fast-light and slow-light light curves is less than 5%. Using two distinct radiative-transfer codes to process the same data, we find integrated flux densities with a relative difference less than 0.01%. Conclusions. For two-dimensional GRMHD models, such as those examined in this paper, the fast-light approximation suffices as long as errors of a few percent are acceptable. The convergence of the results of two different codes demonstrates that they are, at a minimum, consistent.

Spectroscopic and Photometric Analysis of the HW Vir Star PTF1 J011339.09+225739.1

Abstract HW Vir systems are rare eclipsing binary systems including a subdwarf B star (sdB) with a faint companion, mostly M-dwarfs. Up to now, 19 HW Vir systems have been published, three of them with substellar companions. We report the spectroscopic as well as photometric observation of the eclipsing sdB binary PTF1 J011339.09+225739.1 (PTF1 J0113) in a close (a=0.722 ± 0.023 R⊙), short period (P = 0.0933731(3)d) orbit. A quantitative spectral analysis of the sdB yields Te.=29280 ± 720 K, log(g)=5.77 ± 0.09 dex, and log(y)=−2.32 ± 0.12. The circular orbital velocity of the sdB of K1=74.2 ± 1.7 km s−1 is derived from the radial velocity curve. Except for the strong reflection effect, no other light contribution of the companion could be detected. The light curves - recorded with ULTRACAM - were analyzed using the Wilson-Devinney code. We find an inclination angle of i=79.88 ± 0.18∘. Because our first attempts to determine q failed, we calculated large grids of synthetic lightcurves for several mass ratios. Because of degeneracy, good solutions for different mass ratios were found - the one at q = 0.24 is consistent with the sdB’s canonical mass (MsdB = 0.47 M⊙). Accordingly, the mass of the companion is M2=0.112 ± 0.003 M⊙. The radii of the two components were also derived: RsdB=0.178 ± 0.006 R⊙ and R2 = 0.158 ± 0.009 R⊙. Thus, the results for the secondary are consistent with an M-dwarf as secondary

Impact of ejecta morphology and composition on the electromagnetic signatures of neutron star mergers

The electromagnetic transients accompanying compact binary mergers ($\gamma$-ray bursts, afterglows and 'macronovae') are crucial to pinpoint the sky location of gravitational wave sources. Macronovae are caused by the radioactivity from freshly synthesised heavy elements, e.g. from dynamic ejecta and various types of winds. We study macronova signatures by using multi-dimensional radiative transfer calculations. We employ the radiative transfer code SuperNu and state-of-the-art LTE opacities for a few representative elements from the wind and dynamical ejecta (Cr, Pd, Se, Te, Br, Zr, Sm, Ce, Nd, U) to calculate synthetic light curves and spectra for a range of ejecta morphologies. The radioactive power of the resulting macronova is calculated with the detailed input of decay products. We assess the detection prospects for our most complex models, based on the portion of viewing angles that are sufficiently bright, at different cosmological redshifts ($z$). The brighter emission from the wind is unobscured by the lanthanides (or actinides) in some of the models, permitting non-zero detection probabilities for redshifts up to $z=0.07$. We also find the nuclear mass model and the resulting radioactive heating rate are crucial for the detectability. While for the most pessimistic heating rate (from the FRDM model) no reasonable increase in the ejecta mass or velocity, or wind mass or velocity, can possibly make the light curves agree with the observed nIR excess after GRB130603B, a more optimistic heating rate (from the Duflo-Zuker model) leads to good agreement. We conclude that future reliable macronova observations would constrain nuclear heating rates, and consequently help constrain nuclear mass models.

Interaction of a Supernova with a Circumstellar Disk

Abstract Interaction between a supernova (SN) ejecta and a dense circumstellar medium (CSM) can power a luminous light curve and create narrow emission lines in the spectra. While theoretical studies of interaction often assume a spherically symmetric CSM, there are observational indications that the gas surrounding some SNe has a disk-like geometry. Here, we use moving-mesh hydrodynamics simulations to study the interaction of an SN with a disk and determine how the dynamics and observable signatures may depend on the disk mass, thickness, and radial extent. We find that simple modifications to standard spherically-symmetric scaling laws can be used to describe the propagation and heating rate of the interaction shock. We use the resulting shock heating rates to derive approximate bolometric light curves, and provide analytic formulas that can be used to generate simple synthetic light curves for general supernova-disk interactions. For certain disk parameters and explosion energies, we are able to produce luminosities akin to those seen in super-luminous SNe. Because the SN ejecta can flow around and engulf the CSM disk, the interaction region may become embedded and, from certain viewing angles, the narrow emission lines indicative of interaction may be hidden.

The life cycles of Be viscous decretion discs: fundamental disc parameters of 54 SMC Be stars

Be stars are main-sequence massive stars with emission features in their spectrum, which originates in circumstellar gaseous discs. Even though the viscous decretion disc (VDD) model can satisfactorily explain most observations, two important physical ingredients, namely the magnitude of the viscosity ($\alpha$) and the disk mass injection rate, remain poorly constrained. The light curves of Be stars that undergo events of disc formation and dissipation offer an opportunity to constrain these quantities. A pipeline was developed to model these events that uses a grid of synthetic light curves, computed from coupled hydrodynamic and radiative transfer calculations. A sample of 54 Be stars from the OGLE survey of the Small Magellanic Cloud (SMC) was selected for this study. Because of the way our sample was selected (bright stars with clear disc events), it likely represents the densest discs in the SMC. Like their siblings in the Galaxy, the mass of the disc in the SMC increases with the stellar mass. The typical mass and angular momentum loss rates associated with the disk events are of the order of $\sim$$10^{-10}\, M_\odot\,\mathrm{yr^{-1}}$ and $\sim$$5\times 10^{36}\, \mathrm{g\, cm^{2}\, s^{-2}}$, respectively. The values of $\alpha$ found in this work are typically of a few tenths, consistent with recent results in the literature and with the ones found in dwarf novae, but larger than current theory predicts. Considering the sample as a whole, the viscosity parameter is roughly two times larger at build-up ($\left\langle\alpha_\mathrm{bu}\right\rangle = 0.63$) than at dissipation ($\left\langle\alpha_\mathrm{d}\right\rangle = 0.26$). Further work is necessary to verify whether this trend is real or a result of some of the model assumptions.

Sub-Chandrasekhar-mass White Dwarf Detonations Revisited

Abstract The detonation of a sub-Chandrasekhar-mass white dwarf (WD) has emerged as one of the most promising Type Ia supernova (SN Ia) progenitor scenarios. Recent studies have suggested that the rapid transfer of a very small amount of helium from one WD to another is sufficient to ignite a helium shell detonation that subsequently triggers a carbon core detonation, yielding a “dynamically driven double-degenerate double-detonation” SN Ia. Because the helium shell that surrounds the core explosion is so minimal, this scenario approaches the limiting case of a bare C/O WD detonation. Motivated by discrepancies in previous literature and by a recent need for detailed nucleosynthetic data, we revisit simulations of naked C/O WD detonations in this paper. We disagree to some extent with the nucleosynthetic results of previous work on sub-Chandrasekhar-mass bare C/O WD detonations; for example, we find that a median-brightness SN Ia is produced by the detonation of a 1.0 WD instead of a more massive and rarer 1.1 WD. The neutron-rich nucleosynthesis in our simulations agrees broadly with some observational constraints, although tensions remain with others. There are also discrepancies related to the velocities of the outer ejecta and light curve shapes, but overall our synthetic light curves and spectra are roughly consistent with observations. We are hopeful that future multidimensional simulations will resolve these issues and further bolster the dynamically driven double-degenerate double-detonation scenario’s potential to explain most SNe Ia.

PyTranSpot: A tool for multiband light curve modeling of planetary transits and stellar spots

Several studies have shown that stellar activity features, such as occulted and non-occulted starspots, can affect the measurement of transit parameters biasing studies of transit timing variations and transmission spectra. We present PyTranSpot, which we designed to model multiband transit light curves showing starspot anomalies, inferring both transit and spot parameters. The code follows a pixellation approach to model the star with its corresponding limb darkening, spots, and transiting planet on a two dimensional Cartesian coordinate grid. We combine PyTranSpot with a Markov chain Monte Carlo framework to study and derive exoplanet transmission spectra, which provides statistically robust values for the physical properties and uncertainties of a transiting star-planet system. We validate PyTranSpot’s performance by analyzing eleven synthetic light curves of four different star-planet systems and 20 transit light curves of the well-studied WASP-41b system. We also investigate the impact of starspots on transit parameters and derive wavelength dependent transit depth values for WASP-41b covering a range of 6200−9200 Å, indicating a flat transmission spectrum.

Lebedev acceleration and comparison of different photometric models in the inversion of lightcurves for asteroids

In the lightcurve inversion process where asteroid's physical parameters such as rotational period, pole orientation and overall shape are searched, the numerical calculations of the synthetic photometric brightness based on different shape models are frequently implemented. Lebedev quadrature is an efficient method to numerically calculate the surface integral on the unit sphere. By transforming the surface integral on the Cellinoid shape model to that on the unit sphere, the lightcurve inversion process based on the Cellinoid shape model can be remarkably accelerated. Furthermore, Matlab codes of the lightcurve inversion process based on the Cellinoid shape model are available on Github for free downloading. The photometric models, i.e., the scattering laws, also play an important role in the lightcurve inversion process, although the shape variations of asteroids dominate the morphologies of the lightcurves. Derived from the radiative transfer theory, the Hapke model can describe the light reflectance behaviors from the viewpoint of physics, while there are also many empirical models in numerical applications. Numerical simulations are implemented for the comparison of the Hapke model with the other three numerical models, including the Lommel-Seeliger, Minnaert, and Kaasalainen models. The results show that the numerical models with simple function expressions can fit well with the synthetic lightcurves generated based on the Hapke model; this good fit implies that they can be adopted in the lightcurve inversion process for asteroids to improve the numerical efficiency and derive similar results to those of the Hapke model.

The super-orbital modulation of supergiant high-mass X-ray binaries

The long-term X-ray lightcurves of classical supergiant X-ray binaries and supergiant fast X-ray transients show relatively similar super-orbital modulations, which are still lacking a sound interpretation. We propose that these modulations are related to the presence of corotating interaction regions (CIRs) known to thread the winds of OB supergiants. To test this hypothesis, we couple the outcomes of 3-D hydrodynamic models for the formation of CIRs in stellar winds with a simplified recipe for the accretion onto a neutron star. The results show that the synthetic X-ray light curves are indeed modulated by the presence of the CIRs. The exact period and amplitude of these modulations depend on a number of parameters governing the hydrodynamic wind models and on the binary orbital configuration. To compare our model predictions with the observations, we apply the 3-D wind structure previously shown to well explain the appearance of discrete absorption components in the UV time series of a prototypical B0.5I-type supergiant. Using the orbital parameters of IGRJ16493-4348 which has the same B0.5I donor spectral type, the period and modulations in the simulated X-ray light-curve are similar to the observed ones, thus providing support to our scenario. We propose, that the presence of CIRs in donor star winds should be considered in future theoretical and simulation efforts of wind-fed X-ray binaries.

Pair-instability Supernova Simulations: Progenitor Evolution, Explosion, and Light Curves

Abstract In recent years, the viability of the pair-instability supernova (PISN) scenario for explaining superluminous supernovae has all but disappeared except for a few slowly-evolving examples. However, PISNe are not predicted to be superluminous throughout the bulk of their mass range. In fact, it is more likely that the first PISN we see (if we have not seen one already) will not be superluminous. Here, we present hydrodynamic simulations of PISNe for four stellar models with unique envelope properties spanning the PISN mass range. In addition, we compute synthetic light curves (LCs) for comparison with current and future observations. We also investigate, in the context of our most massive model, the prospect of mixing in the supernova ejecta, alleviating discrepancies between current PISN models and the remaining superluminous candidate events. To this end, we present the first published 3D hydrodynamic simulations of PISNe. After achieving convergence between 1D, 2D, and 3D simulations, we examine mixing in the supernova ejecta and its affect on the bolometric LC. We observe slight deviations from spherical symmetry, which increase with the number of dimensions. We find no significant effects on the bolometric LC; however, we conclude that mixing between the silicon and oxygen rich layers caused by the Rayleigh–Taylor instability may affect spectra.

Scattering Law Analysis Based on Hapke and Lommel-Seeliger Models for Asteroidal Taxonomy

In deriving the physical properties of asteroids from their photometric data, the scattering law plays an important role, although the shape variations of asteroids result in the main variations in lightcurves. By following the physical behaviors of light reflections, Hapke et al. deduced complex functions to represent the scattering process, however, it is very hard to accurately simulate the surface scattering law in reality. For simplicity, other numerical scattering models are presented for efficiently calculating the physical properties of asteroids, such as the Lommel-Seeliger (LS) model. In this article, these two models are compared numerically. It is found that in some numerical applications the LS model in simple form with four parameters can be exploited to replace the Hapke model in complex form with five parameters. Furthermore, the generated synthetic lightcurves by the Cellinoid shape model also show that the LS model can perform as well as the Hapke model in the inversion process. Finally, by applying the Principal Component Analysis (PCA) technique to the parameters of the LS model, we present an efficient method to classify C and S type asteroids, instead of the conventional method using the parameters of the Hapke model.

Early light curves for Type Ia supernova explosion models

Upcoming high-cadence transient survey programmes will produce a wealth of observational data for Type Ia supernovae. These data sets will contain numerous events detected very early in their evolution, shortly after explosion. Here, we present synthetic light curves, calculated with the radiation hydrodynamical approach Stella for a number of different explosion models, specifically focusing on these first few days after explosion. We show that overall the early light curve evolution is similar for most of the investigated models. Characteristic imprints are induced by radioactive material located close to the surface. However, these are very similar to the signatures expected from ejecta-CSM or ejecta-companion interaction. Apart from the pure deflagration explosion models, none of our synthetic light curves exhibit the commonly assumed power-law rise. We demonstrate that this can lead to substantial errors in the determination of the time of explosion. In summary, we illustrate with our calculations that even with very early data an identification of specific explosion scenarios is challenging, if only photometric observations are available.

Limb Darkening and Planetary Transits: Testing Center-to-limb Intensity Variations and Limb-darkening Directly from Model Stellar Atmospheres

Abstract The transit method, employed by Microvariability and Oscillation of Stars (MOST), Kepler, and various ground-based surveys has enabled the characterization of extrasolar planets to unprecedented precision. These results are precise enough to begin to measure planet atmosphere composition, planetary oblateness, starspots, and other phenomena at the level of a few hundred parts per million. However, these results depend on our understanding of stellar limb darkening, that is, the intensity distribution across the stellar disk that is sequentially blocked as the planet transits. Typically, stellar limb darkening is assumed to be a simple parameterization with two coefficients that are derived from stellar atmosphere models or fit directly. In this work, we revisit this assumption and compute synthetic planetary-transit light curves directly from model stellar atmosphere center-to-limb intensity variations (CLIVs) using the plane-parallel Atlas and spherically symmetric SAtlas codes. We compare these light curves to those constructed using best-fit limb-darkening parameterizations. We find that adopting parametric stellar limb-darkening laws leads to systematic differences from the more geometrically realistic model stellar atmosphere CLIV of about 50–100 ppm at the transit center and up to 300 ppm at ingress/egress. While these errors are small, they are systematic, and they appear to limit the precision necessary to measure secondary effects. Our results may also have a significant impact on transit spectra.

Modeling the Infrared Reverberation Response of the Circumnuclear Dusty Torus in AGNs: The Effects of Cloud Orientation and Anisotropic Illumination

Abstract The obscuring circumnuclear torus of dusty molecular gas is one of the major components of active galactic nuclei (AGN). The torus can be studied by analyzing the time response of its infrared (IR) dust emission to variations in the AGN continuum luminosity, a technique known as reverberation mapping. The IR response is the convolution of the AGN ultraviolet/optical light curve with a transfer function that contains information about the size, geometry, and structure of the torus. Here, we describe a new computer model that simulates the reverberation response of a clumpy torus. Given an input optical light curve, the code computes the emission of a 3D ensemble of dust clouds as a function of time at selected IR wavelengths, taking into account light travel delays. We present simulated dust emission responses at 3.6, 4.5, and 30 μm that explore the effects of various geometrical and structural properties, dust cloud orientation, and anisotropy of the illuminating radiation field. We also briefly explore the effects of cloud shadowing (clouds are shielded from the AGN continuum source). Example synthetic light curves have also been generated, using the observed optical light curve of the Seyfert 1 galaxy NGC 6418 as input. The torus response is strongly wavelength-dependent, due to the gradient in cloud surface temperature within the torus, and because the cloud emission is strongly anisotropic at shorter wavelengths. Anisotropic illumination of the torus also significantly modifies the torus response, reducing the lag between the IR and optical variations.

The K spectral type contact binary NSVS 1557555: Photometric solution and preliminary elements

We present the first CCD sets of complete light curves for the W Ursae Majoris system NSVS 1557555.The observations were performed in the B, V and Ic bands using the 0.25 m telescope of the Stazione Astronomica Betelgeuse Northern Italy, during 8 nights in October and November 2016.Based on our new eleven Time of Minima (ToM), and two recent ones found in bibliography, the short orbital period of the system is confirmed and revised to P = 0.2725163 days .A reasonable fit of the synthetic light curves of the data indicate that NSVS 1557555 is a late-type (K1+K3) shallow contact binary system of W-Subtype of the W Ursae Majoris systems, with a mass ratio of q = 1.8, a degree of contact factor f = 12.5%, a temperature difference between the components of 240K and inclination i = 85°.The light curves show asymmetries at the maxima with the maximum at phase 0.75 higher the other one (inverse O’Connell effect).To explain the light asymmetries we used a model that involves an hot spotted region on the surface of the cooler star.The definitive solution is only possible with a large amount of third light (L3 = 0.58 in B Filter). It may come from a hot tertiary component.The absolute dimensions of the system are estimated. From the logM-logL diagram it is seen that both components of NSVS 1557555 follow the general pattern of the W subtype W Ursae Majoris systems.The orbital angular momentum is compared with those of other W UMa type binaries and is normal.