Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STAR MODELS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Towards a self-consistent model of the convective core boundary in upper main sequence stars: I. 2.5D and 3D simulations

There is strong observational evidence that the convective cores of intermediate-mass and massive main sequence stars are substantially larger than those predicted by standard stellar-evolution models. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-timescale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number S EVEN -L EAGUE H YDRO code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a 15 M ⊙ zero-age main sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor, ranging from 10 3 to 10 6 . After reaching thermal equilibrium, all of our models show a clear penetration layer; its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. The penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification in simulations at the lower end of our luminosity range. This structure corresponds to the adiabatic ‘step overshoot’ model often employed in stellar-evolution calculations. The simulations with the highest and lowest luminosity differ by less than a factor of two in the penetration distance. The high computational cost of 3D simulations makes our current 3D data set rather sparse. Depending on how we extrapolate the 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from 0.09 to 0.44 pressure scale heights, which is broadly compatible with observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

X-RAY OBSERVATIONS OF BOW SHOCKS AROUND RUNAWAY O STARS. THE CASE OF ζ OPH AND BD+43°3654

Non-thermal radiation has been predicted within bow shocks around runaway stars by recent theoretical works. We present X-ray observations toward the runaway stars ζ Oph by Chandra and Suzaku and of BD+43°3654 by XMM-Newton to search for the presence of non-thermal X-ray emission. We found no evidence of non-thermal emission spatially coincident with the bow shocks; nonetheless, diffuse emission was detected in the vicinity of ζ Oph. After a careful analysis of its spectral characteristics, we conclude that this emission has a thermal nature with a plasma temperature of T ≈ 2 × 10{sup 6} K. The cometary shape of this emission seems to be in line with recent predictions of radiation-hydrodynamic models of runaway stars. The case of BD+43°3654 is puzzling, as non-thermal emission has been reported in a previous work for this source.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2022

A Census of Young Stellar Objects in Two Line-of-Sight Star-Forming Regions Toward IRAS 22147+5948 in the Outer Galaxy

Context. Star formation in the outer Galaxy, namely, outside of the Solar circle, has not been extensively studied in part due to the low CO brightness of the molecular clouds linked with the negative metallicity gradient. Recent infrared surveys provide an overview of dust emission in large sections of the Galaxy, but they suffer from cloud confusion and poor spatial resolution at far-infrared wavelengths. Aims. We aim to develop a methodology to identify and classify young stellar objects (YSOs) in star-forming regions in the outer Galaxy and use it to resolve a long-standing disparity in terms of the distance and evolutionary status of IRAS 22147+5948. Methods. We used a support vector machine learning algorithm to complement standard color–color and color–magnitude diagrams in our search for YSOs in the IRAS 22147 region, based on publicly available data from the Spitzer Mapping of the Outer Galaxy survey. The agglomerative hierarchical clustering algorithm was used to identify clusters. Then the physical properties of individual YSOs were calculated. The distances were determined using CO 1–0 from the Five College Radio Astronomy Observatory survey. Results. We identified 13 Class I and 13 Class II YSO candidates using the color–color diagrams, along with an additional 2 and 21 sources, respectively, using the applied machine learning techniques. The spectral energy distributions of 23 sources were modeled with a star and a passive disk, corresponding to Class II objects. The models of three sources include envelopes that are typical for Class I objects. The objects were grouped into two clusters located at a distance of 2:2 kpc and 5 clusters at 5:6 kpc. The spatial extent of CO, radio continuum, and dust emission confirms the origin of YSOs in two distinct star-forming regions along a similar line of sight. Conclusions. The outer Galaxy may serve as a unique laboratory for exploring star formation across environments, on the condition that complementary methods and ancillary data are used to properly account for cloud confusion and distance uncertainties.

Agata Karska↗
DOE OSTI2012

STELLAR DIAMETERS AND TEMPERATURES. II. MAIN-SEQUENCE K- AND M-STARS

We present interferometric angular diameter measurements of 21 low-mass, K- and M-dwarfs made with the CHARA Array. This sample is enhanced by adding a collection of radius measurements published in the literature to form a total data set of 33 K-M-dwarfs with diameters measured to better than 5%. We use these data in combination with the Hipparcos parallax and new measurements of the star's bolometric flux to compute absolute luminosities, linear radii, and effective temperatures for the stars. We develop empirical relations for {approx}K0 to M4 main-sequence stars that link the stellar temperature, radius, and luminosity to the observed (B - V), (V - R), (V - I), (V - J), (V - H), and (V - K) broadband color index and stellar metallicity [Fe/H]. These relations are valid for metallicities ranging from [Fe/H] = -0.5 to +0.1 dex and are accurate to {approx}2%, {approx}5%, and {approx}4% for temperature, radius, and luminosity, respectively. Our results show that it is necessary to use metallicity-dependent transformations in order to properly convert colors into stellar temperatures, radii, and luminosities. Alternatively, we find no sensitivity to metallicity on relations we construct to the global properties of a star omitting color information, e.g., temperature-radius and temperature-luminosity. Thus, we are able to empirically quantify to what order the star's observed color index is impacted by the stellar iron abundance. In addition to the empirical relations, we also provide a representative look-up table via stellar spectral classifications using this collection of data. Robust examinations of single star temperatures and radii compared to evolutionary model predictions on the luminosity-temperature and luminosity-radius planes reveal that models overestimate the temperatures of stars with surface temperatures <5000 K by {approx}3%, and underestimate the radii of stars with radii <0.7 R{sub Sun} by {approx}5%. These conclusions additionally suggest that the models over account for the effects that the stellar metallicity may have on the astrophysical properties of an object. By comparing the interferometrically measured radii for the single star population to those of eclipsing binaries, we find that for a given mass, single and binary star radii are indistinguishable. However, we also find that for a given radius, the literature temperatures for binary stars are systematically lower compared to our interferometrically derived temperatures of single stars by {approx}200 to 300 K. The nature of this offset is dependent on the validation of binary star temperatures, where bringing all measurements to a uniform and correctly calibrated temperature scale is needed to identify any influence stellar activity may have on the physical properties of a star. Lastly, we present an empirically determined H-R diagram using fundamental properties presented here in combination with those in Boyajian et al. for a total of 74 nearby, main-sequence, A- to M-type stars, and define regions of habitability for the potential existence of sub-stellar mass companions in each system.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Models for various aspects of dwarf novae and nova-like stars

The first attempts to explain the nature of dwarf novae were based on the assumption of single-star phenomena, in which emission lines were assumed to be caused by circumstellar gas shells. The outburst behavior was tentatively ascribed to the kind of (also not understood) mechanism leading to nova outbursts. The realization that some, and possibly all, dwarf novae and nova-like stars (and novae) are binaries eventually led to models which bore more and more similarities to the modern interpretation on the basis of the Roche model. Not all cataclysmic variables are known binaries. In fact, with respect to the entire number of known objects, the proven binaries are still the minority, but all the brightest variables are in fact known to binaries. Not a single system is known which exhibits the usual characteristics of a cataclysmic variable and at the same time can be declared with certainty to be a single star. Two systems are known, the dwarf nova EY Cyg and the recurrent nova V1017 Sgr, in which, in spite of intensive search, no radial velocity variations have been found; but they still exhibit composite spectra consisting of a bright continuum, an emission spectrum, and a cool absorption spectrum. If the Roche model is correct, it is to be expected that a small percentage of objects is viewed pole-on, so orbital motions do not make themselves felt as Doppler shifts of spectral lines. So even these two systems support the hypothesis that all cataclysmic variables (with the possible exception of symbiotic stars) are binaries. In cataclysmic variables, it seems that the brightness changes observed in dwarf novae and nova-like stars in the optical and the UV are due directly to changes in the accretion disks. The study and understanding of accretion disks in these systems can bear potentially valuable consequences for many other fields in astronomy. The observed spectra of dwarf novae and nova-like stars comprise a fairly large range: pure emission spectra, pure absorption spectra, a mixture of both, asymmetric line profiles, very different slopes of the continuous flux distribution -- and one single system may exhibit all of these features at different times. Agreement and disagreement between computed and observed spectra should show whether or not the Roche model is applicable and where it probably will have to be modified and improved. Except for their outburst behavior and its immediate consequences, novae, dwarf novae, and nova-like stars cannot be physically distinguished from each other.

Ladous, Constanze↗
NASA NTRS2019

Star Tracker Accuracy Improvement and Optimization for Attitude Measurement in Three-Axis

High precision attitude measurement systems obviate the need for the beacon from the receiver making it possible for the spacecraft to beam a laser communications signal to a ground station without the ground station advertising its location. The research presented targets new detection and estimation methods to improve the accuracy in locating stars on a focal plane detector, and an understanding of the effects of changes in the optics design parameters and aberration, including defocus, on the navigation solution itself. This understanding can lead to an optimization of the attitude solution with respect to those optics realm parameter changes. The methodology discussed includes the development of a model of a current star tracker system. Using this model, multiple algorithms are implemented, including a multi-hypothesis method (MHT), to detect and estimate the position of the stars on the focal plane detector. It will be shown that using the MHT for detection and estimation, a greater accuracy can be found for each star estimation from more traditional detection and estimation algorithms. The approach then uses the model to develop statistics of the star tracker and the attitude estimation outputs to understand the accuracy, or variance, of the system's attitude solution. This solution is repeated for a range of defocus aberration, and a lower limit to the variance of the attitude solution is shown. A Cramer Rao lower bound solution is derived for the star tracker system and the results are compared to the Monte Carlo analysis from the model and shown to correlate very well. The approach uses a star image not as a Gaussian spot on the focal plane as done in previous work, and use of an image that includes the effects of aberrations of the optic system, and the effects of under-sampling and noise from the focal plane detector as well. Analysis includes exploring a star tracker's accuracy improvement through the combination of focus error and under-sampling effects alone, possibly contradicting conventional wisdom and approaches.

Lichter, Michael J.↗
DOE OSTI2014

Short-term variability of X-rays from accreting neutron star Vela X-1. II. Monte Carlo modeling

We develop a Monte Carlo Comptonization model for the X-ray spectrum of accretion-powered pulsars. Simple, spherical, thermal Comptonization models give harder spectra for higher optical depth, while the observational data from Vela X-1 show that the spectra are harder at higher luminosity. This suggests a physical interpretation where the optical depth of the accreting plasma increases with the mass accretion rate. We develop a detailed Monte Carlo model of the accretion flow, including the effects of the strong magnetic field (∼10{sup 12} G), both in geometrically constraining the flow into an accretion column and in reducing the cross section. We treat bulk-motion Comptonization of the infalling material as well as thermal Comptonization. These model spectra can match the observed broadband Suzaku data from Vela X-1 over a wide range of mass accretion rates. The model can also explain the so-called 'low state' in which the luminosity decreases by an order of magnitude. Here, thermal Comptonization should be negligible, so the spectrum is instead dominated by bulk-motion Comptonization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2010

Mimicking diffuse supernova antineutrinos with the sun as a source

Measuring the {nu}-bar{sub e} component of the cosmic diffuse supernova neutrino background (DSNB) is the next ambitious goal for low-energy neutrino astronomy. The largest flux is expected in the lowest accessible energy bin. However, for E {<=} 15 MeV a possible signal can be mimicked by a solar {nu}-bar{sub e} flux that originates from the usual {sup 8}B neutrinos by spin-flavor oscillations. We show that such an interpretation is possible within the allowed range of neutrino electromagnetic transition moments and solar turbulent field strengths and distributions. Therefore, an unambiguous detection of the DSNB requires a significant number of events at E {>=} 15 MeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2023

Modeling the Chromosphere and Transition Region of Planet-hosting Star GJ 436

Ahead of upcoming space missions intending to conduct observations of low-mass stars in the ultraviolet (UV) spectral region it becomes imperative to simultaneously conduct atmospheric modeling from the UV to the visible (VIS) and near-infrared (NIR). Investigations on extended spectral regions will help to improve the overall understanding of the diversity of spectral lines arising from very different atmospheric temperature regions. Here we investigate atmosphere models with a chromosphere and transition region for the M2.5V star GJ 436, which hosts a close-in Hot Neptune. The atmosphere models are guided by observed spectral features from the UV to the VIS/NIR originating in the chromosphere and transition region of GJ 436. High-resolution observations from the Hubble Space Telescope and Calar Alto high-Resolution search for M dwarfs with Exo-earths with Near-infrared and optical Echelle Spectrographs (CARMENES) are used to obtain an appropriate model spectrum for the investigated M dwarf. We use a large set of atomic species considered in nonlocal thermodynamic equilibrium conditions within our PHOENIX model computations to approximate the physics within the low-density atmospheric regions. In order to obtain an overall match for the nonsimultaneous observations, it is necessary to apply a linear combination of two model spectra, where one of them better reproduces the UV lines while the other better represents the lines from the VIS/NIR range. This is needed to adequately handle different activity states across the observations.

stars↗
DOE OSTI2011

THE STELLAR-ACTIVITY-ROTATION RELATIONSHIP AND THE EVOLUTION OF STELLAR DYNAMOS

We present a sample of 824 solar and late-type stars with X-ray luminosities and rotation periods. This is used to study the relationship between rotation and stellar activity and derive a new estimate of the convective turnover time. From an unbiased subset of this sample the power-law slope of the unsaturated regime, L{sub X} /L{sub bol}{proportional_to}Ro{sup {beta}}, is fit as {beta} = -2.70 {+-} 0.13. This is inconsistent with the canonical {beta} = -2 slope to a confidence of 5{sigma}, and argues for an additional term in the dynamo number equation. From a simple scaling analysis this implies {Delta}{Omega}/{Omega}{proportional_to}{Omega}{sup 0.7}, i.e., the differential rotation of solar-type stars gradually declines as they spin down. Supersaturation is observed for the fastest rotators in our sample and its parametric dependencies are explored. Significant correlations are found with both the corotation radius and the excess polar updraft, the latter theory providing a stronger dependence and being supported by other observations. We estimate mass-dependent empirical thresholds for saturation and supersaturation and map out three regimes of coronal emission. Late F-type stars are shown never to pass through the saturated regime, passing straight from supersaturated to unsaturated X-ray emission. The theoretical threshold for coronal stripping is shown to be significantly different from the empirical saturation threshold (Ro < 0.13), suggesting it is not responsible. Instead we suggest that a different dynamo configuration is at work in stars with saturated coronal emission. This is supported by a correlation between the empirical saturation threshold and the time when stars transition between convective and interface sequences in rotational spin-down models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

DISCOVERY OF A COMPANION CANDIDATE IN THE HD 169142 TRANSITION DISK AND THE POSSIBILITY OF MULTIPLE PLANET FORMATION

We present L'- and J-band high-contrast observations of HD 169142, obtained with the Very Large Telescope/NACO AGPM vector vortex coronagraph and the Gemini Planet Imager, respectively. A source located at 0.''156 ± 0.''032 north of the host star (P.A. = 7.°4 ± 11.°3) appears in the final reduced L' image. At the distance of the star (∼145 pc), this angular separation corresponds to a physical separation of 22.7 ± 4.7 AU, locating the source within the recently resolved inner cavity of the transition disk. The source has a brightness of L' = 12.2 ± 0.5 mag, whereas it is not detected in the J band (J >13.8 mag). If its L' brightness arose solely from the photosphere of a companion and given the J – L' color constraints, it would correspond to a 28-32 M {sub Jupiter} object at the age of the star, according to the COND models. Ongoing accretion activity of the star suggests, however, that gas is left in the inner disk cavity from which the companion could also be accreting. In this case, the object could be lower in mass and its luminosity enhanced by the accretion process and by a circumplanetary disk. A lower-mass object is more consistent with the observed cavity width. Finally, the observations enable us to place an upper limit on the L'-band flux of a second companion candidate orbiting in the disk annular gap at ∼50 AU, as suggested by millimeter observations. If the second companion is also confirmed, HD 169142 might be forming a planetary system, with at least two companions opening gaps and possibly interacting with each other.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A constraint on the formation timescale of the young open cluster NGC 2264: lithium abundance of pre-main sequence stars

The timescale of cluster formation is an essential parameter in order to understand the formation process of star clusters. Pre-main sequence (PMS) stars in nearby young open clusters reveal a large spread in brightness. If the spread were considered to be a result of a real spread in age, the corresponding cluster formation timescale would be about 5–20 Myr. Hence it could be interpreted that star formation in an open cluster is prolonged for up to a few tens of Myr. However, difficulties in reddening correction, observational errors, and systematic uncertainties introduced by imperfect evolutionary models for PMS stars can result in an artificial age spread. Alternatively, we can utilize Li abundance as a relative age indicator of PMS star to determine the cluster formation timescale. The optical spectra of 134 PMS stars in NGC 2264 have been obtained with MMT/Hectochelle. The equivalent widths have been measured for 86 PMS stars with a detectable Li line (3500<T{sub eff}[K]⩽6500). Li abundance under the condition of local thermodynamic equilibrium (LTE) was derived using the conventional curve of growth method. After correction for non-LTE effects, we find that the initial Li abundance of NGC 2264 is A(Li)=3.2±0.2. From the distribution of the Li abundances, the underlying age spread of the visible PMS stars is estimated to be about 3–4 Myr and this, together with the presence of embedded populations in NGC 2264, suggests that the cluster formed on a timescale shorter than 5 Myr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

New Improved Indirect Measurement of the {sup 19}F( p , α ){sup 16}O Reaction at Energies of Astrophysical Relevance

Fluorine abundance determination is of great importance in stellar physics to understand s-elements production and mixing processes in asymptotic giant branch (AGB) stars. Up to now, theoretical models overproduce F abundances in AGB stars with respect to the observed values, thus calling for further investigation of the reactions involving fluorine. In particular, the {sup 19}F( p , α ){sup 16}O reaction is the main destruction channel of fluorine at the bottom of the convective envelope in AGB stars, an H-rich environment where it can experience temperatures high enough to determine its destruction, owing to additional mixing processes. In this paper the Trojan horse method (THM) was used to extract the {sup 19}F( p , α {sub 0}){sup 16}O S-factor in the energy range of astrophysical interest ( E {sub cm} ≈ 0–1 MeV). This is the most relevant channel at the low temperatures (few 10{sup 7} K) characterizing the bottom of the convective envelope, according to current knowledge. A previous indirect experiment using the THM has observed three resonances in the energy regions below E {sub cm} ≈ 450 keV. These energies correspond to typical AGB temperatures, thus implying a significant increase in the reaction rate. Statistics are scarce for performing an accurate separation between resonances, preventing one from drawing a quantitative conclusion about their total widths and spin parities. Before THM measurement, only extrapolations were available below about 500 keV, showing a non-resonant behavior that sharply contradicts the trend of the astrophysical factor at higher energies. A new experiment has been performed to verify the measured TH astrophysical factor and to perform more accurate spectroscopy of the involved resonances.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

Adiabatic properties of pulsating DA white dwarfs. I - The treatment of the Brunt-Vaisala frequency and the region of period formation

The fundamental issue of the region of period formaton in a degenerate star is examined, with special attention given to the treatment of the Brunt-Vaisala frequency. It is shown that, in order to obtain reliable numerical results in degenerate stellar models, the Brunt-Vaisala frequency must be appropriately transformed, because it is defined in terms of a difference between two numbers which become nearly equal in highly degenerate matter, causing serious numerical problems and systematic errors. An alternative expression is derived, which is valid for multicomponent nonideal partially degenerate and partially ionized plasmas such as those encountered in white dwarf envelopes. This expression is used to compute the period structure of the same white dwarf considered by Pesnell (1987). It is shown that the implicit numerical differencing used in the Lagrangian pulsation code of Pesnell leads to very serious difficulties when used with models of degenerate stars.

Brassard, P.↗
NASA NTRS1990

Gamma-ray bursts from extinct neutron stars

The paper concentrates on disks around old extinct pulsars, that can produce gamma-ray bursts owing to viscous evolution of the disk bringing it into the near vicinity of the neutron star, with runaway ionization of the disk and simultaneous precipitation of this plasma onto the neutron star. An old extinct pulsar is modeled as a magnetized neutron star circled by a ring of cold dense matter with an orbital period approximately equal to the rotational period of the neutron star. The numerical estimates produced are found to be consistent with the observed properties of gamma-ray bursters.

Michel, F. C.↗
DOE OSTI2017

Two Coexisting Families of Compact Stars: Observational Implications for Millisecond Pulsars

It is usually thought that a single equation of state (EoS) model “correctly” represents cores of all compact stars. Here we emphasize that two families of compact stars, viz., neutron stars and strange stars, can coexist in nature, and that neutron stars can get converted to strange stars through the nucleation process of quark matter in the stellar center. From our fully general relativistic numerical computations of the structures of fast-spinning compact stars, known as millisecond pulsars, we find that such a stellar conversion causes a simultaneous spin-up and decrease in gravitational mass of these stars. This is a new type of millisecond pulsar evolution through a new mechanism, which gives rise to relatively lower mass compact stars with higher spin rates. This could have an implication for the observed mass and spin distributions of millisecond pulsars. Such a stellar conversion can also rescue some massive, spin-supported millisecond pulsars from collapsing into black holes. Besides, we extend the concept of critical mass M{sub cr} for the neutron star sequence to the case of fast-spinning neutron stars, and point out that neutron star EoS models cannot be ruled out by the stellar mass measurement alone. Finally, we emphasize the additional complexity for constraining EoS models, for example, by stellar radius measurements using X-ray observations, if two families of compact stars coexist.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

THE METALLICITY OF THE CM DRACONIS SYSTEM

The CM Draconis system comprises two eclipsing mid-M dwarfs of nearly equal mass in a 1.27 day orbit. This well-studied eclipsing binary has often been used for benchmark tests of stellar models, since its components are among the lowest mass stars with well-measured masses and radii ({approx}< 1% relative precision). However, as with many other low-mass stars, non-magnetic models have been unable to match the observed radii and effective temperatures for CM Dra at the 5%-10% level. To date, the uncertain metallicity of the system has complicated comparison of theoretical isochrones with observations. In this Letter, we use data from the SpeX instrument on the NASA Infrared Telescope Facility to measure the metallicity of the system during primary and secondary eclipses, as well as out of eclipse, based on an empirical metallicity calibration in the H and K near-infrared (NIR) bands. We derive an [Fe/H] = -0.30 {+-} 0.12 that is consistent across all orbital phases. The determination of [Fe/H] for this system constrains a key dimension of parameter space when attempting to reconcile model isochrone predictions and observations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1995

Surface temperature of a magnetized neutron star and interpretation of the ROSAT data. 1: Dipole fields

We model the temperature distribution at the surface of a magnetized neutron star and study the effects on the observed X-ray spectra and light curves. Generalrelativistic effects, i.e., redshift and lensing, are fully taken into account. Atmospheric effects on the emitted spectral flux are not included: we consider only blackbody emission at the local effective temperature. In this first paper we restrict ourselves to dipole fields. General features are studied and compared with the ROSAT data from the pulsars 0833 - 45 (Vela), 0656 + 14, 0630 + 178 (Geminga), and 1055 - 52, the four cases for which there is strong evidence that thermal radiation from the stellar surface is detected. The composite spectra we obtain are not very different from a blackbody spectrum at the star's effective temperature. We conclude that, as far as blackbody spectra are considered, temperature estimates using single-temperature models give results practically identical to our composite models. The change of the (composite blackbody) spectrum with the star's rotational phase is also not very large and may be unobservable inmost cases. Gravitational lensing strongly suppresses the light curve pulsations. If a dipole field is assumed, pulsed fractions comparable to the observed ones can be obtained only with stellar radii larger than those which are predicted by current models of neutron star struture, or with low stellar masses. Moreover, the shapes of the theoretical light curves with dipole fields do not correspond to the observations. The use of magnetic spectra may raise the pulsed fraction sufficiently but will certainly make the discrepancy with the light curve shapes worse: dipole fields are not sufficient to interpret the data. Many neutron star models with a meson condensate or hypersons predict very small radii, and hence very strong lensing, which will require highly nondipolar fields to be able to reproduce the observed pulsed fractions, if possible at all: this may be a new tool to constrain the size of neutron stars. The pulsed fractions obtained in all our models increase with photon energy: the strong decrease observed in Geminga at energies 0.3-0.5 keV is definitely a genuine effect of the magnetic field on the spectrum in contradistinction to the magnetic effects on the surface temperature considered her. Thus, a detailed analysis of thermal emission from the four pulsars we consider will require both complex surface field configurations and the inclusion of magnetic effects in the atmosphere (i.e., on the emitted spectrum).

Page, Dany↗