Study on period search methods of variable stars: Application to ASAS and CRTS databases

Study on variable stars is an important topic of modern astrophysics. After the invention of powerful telescopes and high resolving powered CCD’s, the variable star data is accumulating in the order of peta-bytes. The huge amount of data need lot of automated methods as well as human experts. This thesis is devoted to the data analysis on variable star’s astronomical time series data and hence belong to the inter-disciplinary topic, Astrostatistics. For an observer on earth, stars that have a change in apparent brightness over time are called variable stars. The variation in brightness may be regular (periodic), quasi periodic (semi-periodic) or irregular manner (aperiodic) and are caused by various reasons. In some cases, the variation is due to some internal thermo-nuclear processes, which are generally known as intrinsic vari- ables and in some other cases, it is due to some external processes, like eclipse or rotation, which are known as extrinsic variables. Intrinsic variables can be further grouped into pulsating variables, eruptive variables and flare stars. Extrinsic variables are grouped into eclipsing binary stars and chromospheri- cal stars. Pulsating variables can again classified into Cepheid, RR Lyrae, RV Tauri, Delta Scuti, Mira etc. The eruptive or cataclysmic variables are novae, supernovae, etc., which rarely occurs and are not periodic phenomena. Most of the other variations are periodic in nature. Variable stars can be observed through many ways such as photometry, spectrophotometry and spectroscopy. The sequence of photometric observa- xiv tions on variable stars produces time series data, which contains time, magni- tude and error. The plot between variable star’s apparent magnitude and time are known as light curve. If the time series data is folded on a period, the plot between apparent magnitude and phase is known as phased light curve. The unique shape of phased light curve is a characteristic of each type of variable star. One way to identify the type of variable star and to classify them is by visually looking at the phased light curve by an expert. For last several years, automated algorithms are used to classify a group of variable stars, with the help of computers. Research on variable stars can be divided into different stages like observa- tion, data reduction, data analysis, modeling and classification. The modeling on variable stars helps to determine the short-term and long-term behaviour and to construct theoretical models (for eg:- Wilson-Devinney model for eclips- ing binaries) and to derive stellar properties like mass, radius, luminosity, tem- perature, internal and external structure, chemical composition and evolution. The classification requires the determination of the basic parameters like pe- riod, amplitude and phase and also some other derived parameters. Out of these, period is the most important parameter since the wrong periods can lead to sparse light curves and misleading information. Time series analysis is a method of applying mathematical and statistical tests to data, to quantify the variation, understand the nature of time-varying phenomena, to gain physical understanding of the system and to predict future behavior of the system. Astronomical time series usually suffer from unevenly spaced time instants, varying error conditions and possibility of big gaps. This is due to daily varying daylight and the weather conditions for ground based observations and observations from space may suffer from the impact of cosmic ray particles. Many large scale astronomical surveys such as MACHO, OGLE, EROS, xv ROTSE, PLANET, Hipparcos, MISAO, NSVS, ASAS, Pan-STARRS, Ke- pler,ESA, Gaia, LSST, CRTS provide variable star’s time series data, even though their primary intention is not variable star observation. Center for Astrostatistics, Pennsylvania State University is established to help the astro- nomical community with the aid of statistical tools for harvesting and analysing archival data. Most of these surveys releases the data to the public for further analysis. There exist many period search algorithms through astronomical time se- ries analysis, which can be classified into parametric (assume some underlying distribution for data) and non-parametric (do not assume any statistical model like Gaussian etc.,) methods. Many of the parametric methods are based on variations of discrete Fourier transforms like Generalised Lomb-Scargle peri- odogram (GLSP) by Zechmeister(2009), Significant Spectrum (SigSpec) by Reegen(2007) etc. Non-parametric methods include Phase Dispersion Minimi- sation (PDM) by Stellingwerf(1978) and Cubic spline method by Akerlof(1994) etc. Even though most of the methods can be brought under automation, any of the method stated above could not fully recover the true periods. The wrong detection of period can be due to several reasons such as power leakage to other frequencies which is due to finite total interval, finite sampling interval and finite amount of data. Another problem is aliasing, which is due to the influence of regular sampling. Also spurious periods appear due to long gaps and power flow to harmonic frequencies is an inherent problem of Fourier methods. Hence obtaining the exact period of variable star from it’s time series data is still a difficult problem, in case of huge databases, when subjected to automation. As Matthew Templeton, AAVSO, states “Variable star data analysis is not always straightforward; large-scale, automated analysis design is non-trivial”. Derekas et al. 2007, Deb et.al. 2010 states “The processing of xvi huge amount of data in these databases is quite challenging, even when looking at seemingly small issues such as period determination and classification”. It will be beneficial for the variable star astronomical community, if basic parameters, such as period, amplitude and phase are obtained more accurately, when huge time series databases are subjected to automation. In the present thesis work, the theories of four popular period search methods are studied, the strength and weakness of these methods are evaluated by applying it on two survey databases and finally a modified form of cubic spline method is intro- duced to confirm the exact period of variable star. For the classification of new variable stars discovered and entering them in the “General Catalogue of Vari- able Stars” or other databases like “Variable Star Index“, the characteristics of the variability has to be quantified in term of variable star parameters.

One-sided jet at milliarcsecond scales in LS I +61˚ 303

We present Very Long Baseline Interferometry (VLBI) observations of the high mass X-ray binary LS I +61˚303, carried out with the European VLBI Network (EVN). Over the 11 hour observing run, performed ~10 days after a radio outburst, the radio source showed a constant flux density, which allowed sensitive imaging of the emission distribution. The structure in the map shows a clear extension to the southeast. Comparing our data with previous VLBI observations we interpret the extension as a collimated radio jet as found in several other X-ray binaries. Assuming that the structure is the result of an expansion that started at the onset of the outburst, we derive an apparent expansion velocity of 0:003 c, which, in the context of Doppler boosting, corresponds to an intrinsic velocity of at least 0:4 c for an ejection close to the line of sight. From the apparent velocity in all available epochs we are able to establish variations in the ejection angle which imply a precessing accretion disk. Finally we point out that LS I +61˚303, like SS 433 and Cygnus X-1, shows evidence for an emission region almost orthogonal to the relativistic jet

SIMULATIONS OF WINDS OF WEAK-LINED T TAURI STARS. II. THE EFFECTS OF A TILTED MAGNETOSPHERE AND PLANETARY INTERACTIONS

Based on our previous work, we investigate here the effects on the wind and magnetospheric structures of weak-lined T Tauri stars due to a misalignment between the axis of rotation of the star and its magnetic dipole moment vector. In such a configuration, the system loses the axisymmetry presented in the aligned case, requiring a fully three-dimensional (3D) approach. We perform 3D numerical magnetohydrodynamic simulations of stellar winds and study the effects caused by different model parameters, namely the misalignment angle theta(t), the stellar period of rotation, the plasma-beta, and the heating index.. Our simulations take into account the interplay between the wind and the stellar magnetic field during the time evolution. The system reaches a periodic behavior with the same rotational period of the star. We show that the magnetic field lines present an oscillatory pattern. Furthermore, we obtain that by increasing theta(t), the wind velocity increases, especially in the case of strong magnetic field and relatively rapid stellar rotation. Our 3D, time-dependent wind models allow us to study the interaction of a magnetized wind with a magnetized extrasolar planet. Such interaction gives rise to reconnection, generating electrons that propagate along the planet`s magnetic field lines and produce electron cyclotron radiation at radio wavelengths. The power released in the interaction depends on the planet`s magnetic field intensity, its orbital radius, and on the stellar wind local characteristics. We find that a close-in Jupiter-like planet orbiting at 0.05 AU presents a radio power that is similar to 5 orders of magnitude larger than the one observed in Jupiter, which suggests that the stellar wind from a young star has the potential to generate strong planetary radio emission that could be detected in the near future with LOFAR. This radio power varies according to the phase of rotation of the star. For three selected simulations, we find a variation of the radio power of a factor 1.3-3.7, depending on theta(t). Moreover, we extend the investigation done in Vidotto et al. and analyze whether winds from misaligned stellar magnetospheres could cause a significant effect on planetary migration. Compared to the aligned case, we show that the timescale tau(w) for an appreciable radial motion of the planet is shorter for larger misalignment angles. While for the aligned case tau(w) similar or equal to 100 Myr, for a stellar magnetosphere tilted by theta(t) = 30 degrees, tau(w) ranges from similar to 40 to 70 Myr for a planet located at a radius of 0.05 AU. Further reduction on tau(w) might occur for even larger misalignment angles and/or different wind parameters.

ON THE NATURE OF THE PROTOTYPE LUMINOUS BLUE VARIABLE AG CARINAE. II. WITNESSING A MASSIVE STAR EVOLVING CLOSE TO THE EDDINGTON AND BISTABILITY LIMITS

We show that the significantly different effective temperatures (T(eff)) achieved by the luminous blue variable AG Carinae during the consecutive visual minima of 1985-1990 (T(eff) similar or equal to 22,800 K) and 2000-2001 (T(eff) similar or equal to 17,000 K) place the star on different sides of the bistability limit, which occurs in line-driven stellar winds around T(eff) similar to 21,000 K. Decisive evidence is provided by huge changes in the optical depth of the Lyman continuum in the inner wind as T(eff) changes during the S Dor cycle. These changes cause different Fe ionization structures in the inner wind. The bistability mechanism is also related to the different wind parameters during visual minima: the wind terminal velocity was 2-3 times higher and the mass-loss rate roughly two times smaller in 1985-1990 than in 2000-2003. We obtain a projected rotational velocity of 220 +/- 50 km s(-1) during 1985-1990 which, combined with the high luminosity (L(star) = 1.5 x 10(6) L(circle dot)), puts AG Car extremely close to the Eddington limit modified by rotation (Omega Gamma limit): for an inclination angle of 90 degrees, Gamma(Omega) greater than or similar to 1.0 for M(circle dot) less than or similar to 60. Based on evolutionary models and mass budget, we obtain an initial mass of similar to 100 M(circle dot) and a current mass of similar to 60-70 M(circle dot) for AG Car. Therefore, AG Car is close to, if not at, the Omega Gamma limit during visual minimum. Assuming M = 70 M(circle dot), we find that Gamma(Omega) decreases from 0.93 to 0.72 as AG Car expands toward visual maximum, suggesting that the star is not above the Eddington limit during maximum phases.

A possible stellar metallic enhancement in post-T Tauri stars by a planetesimal bombardment

The photospheres of stars hosting planets have larger metallicity than stars lacking planets. This could be the result of a metallic star contamination produced by the bombarding of hydrogen-deficient solid bodies. In the present work we study the possibility of an earlier metal enrichment of the photospheres by means of impacting planetesimals during the first 20-30 Myr. Here we explore this contamination process by simulating the interactions of an inward migrating planet with a disc of planetesimal interior to its orbit. The results show the percentage of planetesimals that fall on the star. We identified the dependence of the planet's eccentricity (e(p)) and time-scale of migration (tau) on the rate of infalling planetesimals. For very fast migrations (tau= 10(2) and 10(3) yr) there is no capture in mean motion resonances, independently of the value of e(p). Then, due to the planet's migration the planetesimals suffer close approaches with the planet and more than 80 per cent of them are ejected from the system. For slow migrations (tau= 10(5)and 10(6) yr) the percentage of collisions with the planet decreases with the increase of the planet's eccentricity. For e(p) = 0 and 0.1 most of the planetesimals were captured in the 2:1 resonance and more than 65 per cent of them collided with the star. Whereas migration of a Jupiter mass planet to very short pericentric distances requires unrealistic high disc masses, these requirements are much smaller for smaller migrating planets. Our simulations for a slowly migrating 0.1 M-Jupiter planet, even demanding a possible primitive disc three times more massive than a primitive solar nebula, produces maximum [Fe/H] enrichments of the order of 0.18 dex. These calculations open possibilities to explain hot Jupiter exoplanet metallicities.

Planet formation in a triple stellar system: implications of the third star's orbital inclination

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

NGC 2782: a merger remnant with young stars in its gaseous tidal tail

We have searched for young star-forming regions around the merger remnant NGC 2782. By using Galaxy Evolution Explorer far-ultraviolet and near-ultraviolet imaging and H i data we found seven ultraviolet sources, located at distances greater than 26 kpc from the centre of NGC 2782, and coinciding with its western H i tidal tail. These regions were resolved in several smaller systems when Gemini/Gemini multi-object spectrograph (GMOS) r-band images were used. We compared the observed colours to stellar population synthesis models and found that these objects have ages of similar to 1 to 11 Myr and masses ranging from 103.9 to 104.6 M circle dot. By using Gemini/GMOS spectroscopic data we confirm memberships and derive high metallicities for three of the young regions in the tail (12+log(O/H) = 8.74 +/- 0.20, 8.81 +/- 0.20 and 8.78 +/- 0.20). These metallicities are similar to the value presented by the nuclear region of NGC 2782 and also similar to the value presented for an object located close to the main body of NGC 2782. The high metallicities measured for the star-forming regions in the gaseous tidal tail of NGC 2782 could be explained if they were formed out of highly enriched gas which was once expelled from the centre of the merging galaxies when the system collided. An additional possibility is that the tail has been a nursery of a few generations of young stellar systems which ultimately polluted this medium with metals, further enriching the already pre-enriched gas ejected to the tail when the galaxies collided.

EVOLUTIONARY TRAJECTORIES OF ULTRACOMPACT "BLACK WIDOW" PULSARS WITH VERY LOW MASS COMPANIONS

The existence of millisecond pulsars with planet-mass companions in close orbits is challenging from the stellar evolution point of view. We calculate in detail the evolution of binary systems self-consistently, including mass transfer, evaporation, and irradiation of the donor by X-ray feedback, demonstrating the existence of a new evolutionary path leading to short periods and compact donors as required by the observations of PSR J1719-1438. We also point out the alternative of an exotic nature of the companion planet-mass star.

THE INTERSTELLAR MAGNETIC FIELD CLOSE TO THE SUN. II.

The magnetic field in the local interstellar medium (ISM) provides a key indicator of the galactic environment of the Sun and influences the shape of the heliosphere. We have studied the interstellar magnetic field (ISMF) in the solar vicinity using polarized starlight for stars within 40 pc of the Sun and 90 degrees of the heliosphere nose. In Frisch et al. (Paper I), we developed a method for determining the local ISMF direction by finding the best match to a group of interstellar polarization position angles obtained toward nearby stars, based on the assumption that the polarization is parallel to the ISMF. In this paper, we extend the analysis by utilizing weighted fits to the position angles and by including new observations acquired for this study. We find that the local ISMF is pointed toward the galactic coordinates l, b = 47 degrees +/- 20 degrees, 25 degrees +/- 20 degrees. This direction is close to the direction of the ISMF that shapes the heliosphere, l, b = 33 degrees +/- 4 degrees, 55 degrees +/- 4 degrees, as traced by the center of the "Ribbon" of energetic neutral atoms discovered by the Interstellar Boundary Explorer (IBEX) mission. Both the magnetic field direction and the kinematics of the local ISM are consistent with a scenario where the local ISM is a fragment of the Loop I superbubble. A nearby ordered component of the local ISMF has been identified in the region l approximate to 0 degrees -> 80 degrees and b approximate to 0 degrees -> 30 degrees, where PlanetPol data show a distance-dependent increase of polarization strength. The ordered component extends to within 8 pc of the Sun and implies a weak curvature in the nearby ISMF of +/- 0 degrees.25 pc(-1). This conclusion is conditioned on the small sample of stars available for defining this rotation. Variations from the ordered component suggest a turbulent component of +/- 23 degrees. The ordered component and standard relations between polarization, color excess, and H-o column density predict a reasonable increase of N(H) with distance in the local ISM. The similarity of the ISMF directions traced by the polarizations, the IBEX Ribbon, and pulsars inside the Local Bubble in the third galactic quadrant suggest that the ISMF is relatively uniform over spatial scales of 8-200 pc and is more similar to interarm than spiral-arm magnetic fields. The ISMF direction from the polarization data is also consistent with small-scale spatial asymmetries detected in GeV-TeV cosmic rays with a galactic origin. The peculiar geometrical relation found earlier between the cosmic microwave background dipole moment, the heliosphere nose, and the ISMF direction is supported by this study. The interstellar radiation field at +/- 975 angstrom does not appear to play a role in grain alignment for the low-density ISM studied here.

DYNAMICAL EVOLUTION OF VISCOUS DISKS AROUND Be STARS. I. PHOTOMETRY

Be stars possess gaseous circumstellar disks that modify in many ways the spectrum of the central B star. Furthermore, they exhibit variability at several timescales and for a large number of observables. Putting the pieces together of this dynamical behavior is not an easy task and requires a detailed understanding of the physical processes that control the temporal evolution of the observables. There is an increasing body of evidence that suggests that Be disks are well described by standard alpha-disk theory. This paper is the first of a series that aims at studying the possibility of inferring several disk and stellar parameters through the follow-up of various observables. Here we study the temporal evolution of the disk density for different dynamical scenarios, including the disk buildup as a result of a long and steady mass injection from the star, the disk dissipation that occurs after mass injection is turned off, as well as scenarios in which active periods are followed by periods of quiescence. For those scenarios, we investigate the temporal evolution of continuum photometric observables using a three-dimensional non-LTE radiative transfer code. We show that light curves for different wavelengths are specific of a mass loss history, inclination angle, and alpha viscosity parameter. The diagnostic potential of those light curves is also discussed.

THE GEMINI NICI PLANET-FINDING CAMPAIGN: DISCOVERY OF A MULTIPLE SYSTEM ORBITING THE YOUNG A STAR HD 1160

We report the discovery of two low-mass companions to the young A0V star HD 1160 at projected separations of 81 +/- 5 AU (HD 1160 B) and 533 +/- 25 AU (HD 1160 C) by the Gemini NICI Planet-Finding Campaign. Very Large Telescope images of the system taken over a decade for the purpose of using HD 1160 A as a photometric calibrator confirm that both companions are physically associated. By comparing the system to members of young moving groups and open clusters with well-established ages, we estimate an age of 50(-40)(+50) Myr for HD 1160 ABC. While the UVW motion of the system does not match any known moving group, the small magnitude of the space velocity is consistent with youth. Near-IR spectroscopy shows HD 1160 C to be an M3.5 +/- 0.5 star with an estimated mass of 0.22(-0.04)(+0.03) M-circle dot, while NIR photometry of HD 1160 B suggests a brown dwarf with a mass of 33(-9)(+12) M-Jup. The very small mass ratio (0.014) between the A and B components of the system is rare for A star binaries, and would represent a planetary-mass companion were HD 1160 A to be slightly less massive than the Sun.

The Gemini planet-finding campaign: the frequency of giant planets aroud debris disk stars

We have completed a high-contrast direct imaging survey for giant planets around 57 debris disk stars as part of the Gemini NICI Planet-Finding Campaign. We achieved median H-band contrasts of 12.4 mag at 0.''5 and 14.1 mag at 1'' separation. Follow-up observations of the 66 candidates with projected separation <500 AU show that all of them are background objects. To establish statistical constraints on the underlying giant planet population based on our imaging data, we have developed a new Bayesian formalism that incorporates (1) non-detections, (2) single-epoch candidates, (3) astrometric and (4) photometric information, and (5) the possibility of multiple planets per star to constrain the planet population. Our formalism allows us to include in our analysis the previously known β Pictoris and the HR 8799 planets. Our results show at 95% confidence that <13% of debris disk stars have a ≥5 M Jup planet beyond 80 AU, and <21% of debris disk stars have a ≥3 M Jup planet outside of 40 AU, based on hot-start evolutionary models. We model the population of directly imaged planets as d 2 N/dMdavpropm α a β, where m is planet mass and a is orbital semi-major axis (with a maximum value of a max). We find that β < –0.8 and/or α > 1.7. Likewise, we find that β < –0.8 and/or a max < 200 AU. For the case where the planet frequency rises sharply with mass (α > 1.7), this occurs because all the planets detected to date have masses above 5 M Jup, but planets of lower mass could easily have been detected by our search. If we ignore the β Pic and HR 8799 planets (should they belong to a rare and distinct group), we find that <20% of debris disk stars have a ≥3 M Jup planet beyond 10 AU, and β < –0.8 and/or α < –1.5. Likewise, β < –0.8 and/or a max < 125 AU. Our Bayesian constraints are not strong enough to reveal any dependence of the planet frequency on stellar host mass. Studies of transition disks have suggested that about 20% of stars are undergoing planet formation; our non-detections at large separations show that planets with orbital separation >40 AU and planet masses >3 M Jup do not carve the central holes in these disks.

Laser guide stars wavefront sensors for the EELT

The Adaptive Optics is the measurement and correction in real time of the wavefront aberration of the star light caused by the atmospheric turbulence, that limits the angular resolution of ground based telescopes and thus their capabilities to deep explore faint and crowded astronomical objects. The lack of natural stars enough bright to be used as reference sources for the Adaptive Optics, over a relevant fraction of the sky, led to the introduction of artificial reference stars. The so-called Laser Guide Stars are produced by exciting the Sodium atoms in a layer laying at 90km of altitude, by a powerful laser beam projected toward the sky. The possibility to turn on a reference star close to the scientific targets of interest has the drawback in an increased difficulty in the wavefront measuring, mainly due to the time instability of the Sodium layer density. These issues are increased with the telescope diameter. In view of the construction of the 42m diameter European Extremely Large Telescope a detailed investigation of the achievable performances of Adaptive Optics becomes mandatory to exploit its unique angular resolution . The goal of this Thesis was to present a complete description of a laboratory Prototype development simulating a Shack-Hartmann wavefront sensor using Laser Guide Stars as references, in the expected conditions for a 42m telescope. From the conceptual design, through the opto-mechanical design, to the Assembly, Integration and Test, all the phases of the Prototype construction are explained. The tests carried out shown the reliability of the images produced by the Prototype that agreed with the numerical simulations. For this reason some possible upgrades regarding the opto-mechanical design are presented, to extend the system functionalities and let the Prototype become a more complete test bench to simulate the performances and drive the future Adaptive Optics modules design.

Brown Dwarf Binaries

Nearly 500 brown dwarfs have been discovered in recent years. The majority of these brown dwarfs exist in the solar neighborhood, yet determining their fundamental properties (mass, age, temperature & metallicity) has proved to be quite difficult, with current estimates relying heavily on theoretical models. Binary brown dwarfs provide a unique opportunity to empirically determine fundamental properties, which can then be used to test model predictions. In addition, the observed binary fractions, separations, mass ratios, & orbital eccentricities can provide insight into the formation mechanism of these low-mass objects. I will review the results of various brown dwarf multiplicity studies, and will discuss what we have learned about the formation and evolution of brown dwarfs by examining their binary properties as a function of age and mass.

A VLT/FLAMES survey for massive binaries in Westerlund 1. III. The WC9d binary W239 and implications for massive stellar evolution

Context. There is growing evidence that a treatment of binarity amongst OB stars is essential for a full theory of stellar evolution. However the binary properties of massive stars – frequency, mass ratio & orbital separation – are still poorly constrained. Aims. In order to address this shortcoming we have undertaken a multiepoch spectroscopic study of the stellar population of the young massive cluster Westerlund 1. In this paper we present an investigation into the nature of the dusty Wolf-Rayet star and candidate binary W239. Methods. To accomplish this we have utilised our spectroscopic data in conjunction with multi-year optical and near-IR photometric observations in order to search for binary signatures. Comparison of these data to synthetic non-LTE model atmosphere spectra were used to derive the fundamental properties of the WC9 primary. Results. We found W239 to have an orbital period of only ~5.05 days, making it one of the most compact WC binaries yet identified. Analysis of the long term near-IR lightcurve reveals a significant flare between 2004-6. We interpret this as evidence for a third massive stellar component in the system in a long period (>6 yr), eccentric orbit, with dust production occuring at periastron leading to the flare. The presence of a near-IR excess characteristic of hot (~1300 K) dust at every epoch is consistent with the expectation that the subset of persistent dust forming WC stars are short (<1 yr) period binaries, although confirmation will require further observations. Non-LTE model atmosphere analysis of the spectrum reveals the physical properties of the WC9 component to be fully consistent with other Galactic examples. Conclusions. The simultaneous presence of both short period Wolf-Rayet binaries and cool hypergiants within Wd 1 provides compelling evidence for a bifurcation in the post-Main Sequence evolution of massive stars due to binarity. Short period O+OB binaries will evolve directly to the Wolf-Rayet phase, either due to an episode of binary mediated mass loss – likely via case A mass transfer or a contact configuration – or via chemically homogenous evolution. Conversely, long period binaries and single stars will instead undergo a red loop across the HR diagram via a cool hypergiant phase. Future analysis of the full spectroscopic dataset for Wd 1 will constrain the proportion of massive stars experiencing each pathway; hence quantifying the importance of binarity in massive stellar evolution up to and beyond supernova and the resultant production of relativistic remnants.