Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “normal stars”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Evolution of dynamo-generated magnetic fields in accretion disks around compact and young stars

Geometrically thin, optically thick, turbulent accretion disks are believed to surround many stars. Some of them are the compact components of close binaries, while the others are throught to be T Tauri stars. These accretion disks must be magnetized objects because the accreted matter, whether it comes from the companion star (binaries) or from a collapsing molecular cloud core (single young stars), carries an embedded magnetic field. In addition, most accretion disks are hot and turbulent, thus meeting the condition for the MHD turbulent dynamo to maintain and amplify any seed field magnetic field. In fact, for a disk's magnetic field to persist long enough in comparison with the disk viscous time it must be contemporaneously regenerated because the characteristic diffusion time of a magnetic field is typically much shorter than a disk's viscous time. This is true for most thin accretion disks. Consequently, studying magentic fields in thin disks is usually synonymous with studying magnetic dynamos, a fact that is not commonly recognized in the literature. Progress in studying the structure of many accretion disks was achieved mainly because most disks can be regarded as two-dimensional flows in which vertical and radial structures are largely decoupled. By analogy, in a thin disk, one may expect that vertical and radial structures of the magnetic field are decoupled because the magnetic field diffuses more rapidly to the vertical boundary of the disk than along the radius. Thus, an asymptotic method, called an adiabatic approximation, can be applied to accretion disk dynamo. We can represent the solution to the dynamo equation in the form B = Q(r)b(r,z), where Q(r) describes the field distribution along the radius, while the field distribution across the disk is included in the vector function b, which parametrically depends on r and is normalized by the condition max (b(z)) = 1. The field distribution across the disk is established rapidly, while the radial distribution Q(r) evolves on a considerably longer timescale. It is this evolution that is the subject of this paper.

Stepinski, Tomasz F.↗
NASA NTRS1996

Fault Diagnosis of Power Systems Using Intelligent Systems

The power system operator's need for a reliable power delivery system calls for a real-time or near-real-time Al-based fault diagnosis tool. Such a tool will allow NASA ground controllers to re-establish a normal or near-normal degraded operating state of the EPS (a DC power system) for Space Station Alpha by isolating the faulted branches and loads of the system. And after isolation, re-energizing those branches and loads that have been found not to have any faults in them. A proposed solution involves using the Fault Diagnosis Intelligent System (FDIS) to perform near-real time fault diagnosis of Alpha's EPS by downloading power transient telemetry at fault-time from onboard data loggers. The FDIS uses an ANN clustering algorithm augmented with a wavelet transform feature extractor. This combination enables this system to perform pattern recognition of the power transient signatures to diagnose the fault type and its location down to the orbital replaceable unit. FDIS has been tested using a simulation of the LeRC Testbed Space Station Freedom configuration including the topology from the DDCU's to the electrical loads attached to the TPDU's. FDIS will work in conjunction with the Power Management Load Scheduler to determine what the state of the system was at the time of the fault condition. This information is used to activate the appropriate diagnostic section, and to refine if necessary the solution obtained. In the latter case, if the FDIS reports back that it is equally likely that the faulty device as 'start tracker #1' and 'time generation unit,' then based on a priori knowledge of the system's state, the refined solution would be 'star tracker #1' located in cabinet ITAS2. It is concluded from the present studies that artificial intelligence diagnostic abilities are improved with the addition of the wavelet transform, and that when such a system such as FDIS is coupled to the Power Management Load Scheduler, a faulty device can be located and isolated from the rest of the system. The benefit of these studies provides NASA with the ability to quickly restore the operating status of a space station from a critical state to a safe degraded mode, thereby saving costs in experimentation rescheduling, fault diagnostics, and prevention of loss-of-life.

Momoh, James A.↗
DOE OSTI2026

$\Sigma ^{+}$ production in pp collisions at $\sqrt{\textit{s}}=13$ TeV

The measurement of $\Sigma ^{+}$ production in pp collisions at $\sqrt{s}=13$ TeV is presented. The measurement is performed at midrapidity in both minimum-bias and high-multiplicity pp collisions at $\sqrt{s}=13$ TeV. The $\Sigma ^{+}$ is reconstructed via its weak-decay topology in the decay channel $\Sigma ^{+} \rightarrow \mathrm{{p}} + \pi ^{0}$ with $\pi ^{0} \rightarrow \gamma + \gamma$. In a novel approach, the neutral pion is reconstructed by combining photons that convert in the detector material with photons measured in the calorimeters. The transverse-momentum $(p_{\textrm{T}})$ distributions of the $\Sigma ^{+}$ and its rapidity densities dN/dy in both event classes are reported. The $p_{\textrm{T}}$ spectrum in minimum-bias collisions is compared to QCD-inspired event generators. The ratio of $\Sigma ^{+}$ to previously measured Λ baryons is in good agreement with calculations from the Statistical Hadronization Model. The high efficiency and purity of the novel reconstruction method for $\Sigma ^{+}$ presented here will enable future studies of the interaction of $\Sigma ^{+}$ with protons in the context of femtoscopic measurements, which could be crucial for understanding the equation of state of neutron stars.

Abualrob, I. J. [University of Houston] (ORCID:000↗
NASA NTRS2017

Looking at A 0535+26 at Low Luminosities with NuSTAR

We report on two NuSTAR observations of the high-mass X-ray binary A 0535+26 taken toward the end of its normal 2015 outburst at very low 3-50 keV luminosities of approximately 1.4 times x 10 (sup 36) ergs per second and approximately 5 times x 10 (sup 35) ergs per second, which are complemented by nine Swift observations. The data clearly confirm indications seen in earlier data that the source's spectral shape softens as it becomes fainter. The smooth exponential rollover at high energies seen in the first observation evolves to a much more abrupt steepening of the spectrum at 20-30 keV. The continuum evolution can be nicely described with emission from a magnetized accretion column, modeled using the compmag model modified by an additional Gaussian emission component for the fainter observation. Between the two observations, the optical depth changes from 0.75 plus or minus 0.04 to 0.56 plus 0.01 (sup) minus 0.04 (sub), the electron temperature remains constant, and there is an indication that the column decreases in radius. Since the energy-resolved pulse profiles remain virtually unchanged in shape between the two observations, the emission properties of the accretion column reflect the same accretion regime. This conclusion is also confirmed by our result that the energy of the cyclotron resonant scattering feature (CRSF) at approximately 45 keV is independent of the luminosity, implying that the magnetic field in the region in which the observed radiation is produced is the same in both observations. Finally, we also constrain the evolution of the continuum parameters with the rotational phase of the neutron star. The width of the CRSF could only be constrained for the brighter observation. Based on Monte Carlo simulations of CRSF formation in single accretion columns, its pulse phase dependence supports a simplified fan beam emission pattern. The evolution of the CRSF width is very similar to that of the CRSF depth, which is, however, in disagreement with expectations.

Ballhausen, Ralf↗
DOE OSTI2022

The 3D Kinematics of the Orion Nebula Cluster: NIRSPEC-AO Radial Velocities of the Core Population

The kinematics and dynamics of stellar and substellar populations within young, still-forming clusters provide valuable information for constraining theories of formation mechanisms. Using Keck II NIRSPEC+AO data, we have measured radial velocities for 56 low-mass sources within 4' of the core of the Orion Nebula Cluster (ONC). We also remeasure radial velocities for 172 sources observed with SDSS/APOGEE. These data are combined with proper motions measured using HST ACS/WFPC2/WFC3IR and Keck II NIRC2, creating a sample of 135 sources with all three velocity components. The velocities measured are consistent with a normal distribution in all three components. We measure intrinsic velocity dispersions of (${\sigma }_{{v}_{\alpha }}$, ${\sigma }_{{v}_{\delta }}$, ${\sigma }_{{v}_{r}}$) = (1.64 ± 0.12, 2.03 ± 0.13, ${2.56}_{-0.17}^{+0.16}$) km s –1 . Our computed intrinsic velocity dispersion profiles are consistent with the dynamical equilibrium models from Da Rio et al. (2014) in the tangential direction but not in the line-of-sight direction, possibly indicating that the core of the ONC is not yet virialized, and may require a nonspherical potential to explain the observed velocity dispersion profiles. We also observe a slight elongation along the north–south direction following the filament, which has been well studied in previous literature, and an elongation in the line-of-sight to tangential velocity direction. These 3D kinematics will help in the development of realistic models of the formation and early evolution of massive clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2020

The Fast Radio Burst Luminosity Function and Death Line in the Low-Twist Magnetar Model

We explore the burst energy distribution of fast radio bursts(FRBs)in the low-twist magnetar model ofWadiasingh & Timokhin(WT19). Motivated by the power-lawfluence distributions of FRB121102, we proposean elementary model for the FRB luminosity function of individual repeaters with an inversion protocol thatdirectly relates the power-law distribution index of magnetar short burstfluences to that for FRBs. The protocolindicates that the FRB energy scales virtually linearly with crust/field dislocation amplitude, if magnetar shortbursts prevail in the magnetoelastic regime. Charge starvation in the magnetosphere during bursts(required inWT19)for individual repeaters implies the predicted burstfluence distribution is narrow,3 decades for yieldingstrains and oscillation frequencies feasible in magnetar crusts. Requiring magnetic confinement and chargestarvation, we obtain a death line for FRBs, which segregates magnetars from the normal pulsar population,suggesting only the former will host recurrent FRBs. We convolve the burst energy distribution for individualmagnetars to define the distribution of luminosities in evolved magnetar populations. The broken power-lawluminosity function’s low-energy character depends on the population model, while the high-energy index tracesthat of individual repeaters. Independent of the evolved population, the broken power-law isotropic-equivalentenergy/luminosity function peaks at∼1037–1040erg with a low-energy cutoff at∼1037erg. Lastly, we considerthe localfluence distribution of FRBs andfind that it can constrain the subset of FRB-producing magnetarprogenitors. Our model suggests that improvements in sensitivity may reveal aflattening of the global FRBfluencedistribution and saturation in FRB rates.

Cosmology↗
DOE OSTI2026

Probing gluon saturation with forward di-hadron correlations in proton-nucleus collisions

We present a detailed numerical investigation of semi-inclusive forward di-hadron production in proton–nucleus collisions employing the Color Glass Condensate effective theory. We focus on the regime where di-hadrons are produced nearly back-to-back in the transverse plane, thereby justifying a transverse-momentum-dependent factorization approach in terms of small-x gluon distributions. Our computation integrates several key elements: i) non-linear rapidity evolution via the Balitsky–Kovchegov equation with running coupling, ii) both perturbative and non-perturbative Sudakov resummation, and iii) a phenomenologically constrained model for the initial conditions for small-x gluon distributions. We compare this phenomenological framework to experimental data from the STAR Collaboration on azimuthal correlations in forward di-pion production in both proton–proton and proton–gold collisions. We analyze the systematic theoretical uncertainties associated with the saturation scales of nuclei at the initial scale for rapidity evolution and with those associated with the hadronization process. Finally, we make predictions for the kinematics anticipated to be covered by the ALICE Forward Calorimeter (FoCal) upgrade at the Large Hadron Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1994

Shock propagation and the generation of magnetohydrodynamic wave fields in inhomogeneous molecular clouds

We develop a simple one-dimensional model for the interaction of a steady, thin, planar shock wave with a nonrigid cloud which may be in motion relative to the surrounding medium, and we apply the model to shocks impinging on, and propagating through, molecular clouds. Both 'adiabatic' (gamma = 5/3) and radiative (gamma = 1) shocks are considered and we allow for the presence of a uniform magnetic field directed either parallel or perpendicular to the shock normal. The former field orientation is equivalent to the hydrodynamic case, and the latter involves only fast MHD shocks. We focus on the manner in which such shocks can generate internal kinetic motions in the cloud on a range of size and density scales through the direct acceleration of cores and clumps by shocks transmitted into them and through the generation of an MHD wavefield via the reflection of the incident shock at clump boundaries. We find that stronger incident Mach numbers and smaller density contrasts lead to more efficient cloud acceleration, as do isothermal intercloud shocks and small intercloud magnetic field strengths. The acceleration efficiency is insensitive to the adiabatic index and the magnetic field strength in the cloud itself. For typical parameter choices, the direct acceleration of clouds and clumps by strong shocks is found to be substantial and could at least in part account for their observed velocity dispersions. If the shocks are moderately weak, the final velocity of the cloud is linearly related to its initial velocity, with higher acceleration giving shallower slopes (i.e., final velocity distributions which are less sensitive to the initial distribution). Compared to the kinetic energy of the postshock cloud, the energy given to the wavefield at each encounter is small, and the heating of the interclump medium by the dissipation of this wavefield is found to be insufficient to balance the cooling rate in the cloud as a whole (although it may be important in particular regions), even if this medium is warm, unless it is also extremely tenuous (n approximately less than 0.1/cu cm). Nevertheless, the correction for the velocity imparted to the cloud leads to a substantial increase in the critical incident Mach numner for wave emission over that reported by Spitzer for the rigid case. The implications of our model for shock-induced star formation are discussed briefly.

Miesch, Mark S.↗
DOE OSTI2023

The Origin and Evolution of the Normal Type Ia SN 2018aoz with Infant-phase Reddening and Excess Emission

SN 2018aoz is a Type Ia SN with a B-band plateau and excess emission in infant-phase light curves ≲1 day after the first light, evidencing an over-density of surface iron-peak elements as shown in our previous study. Here, we advance the constraints on the nature and origin of SN 2018aoz based on its evolution until the nebular phase. Near-peak spectroscopic features show that the SN is intermediate between two subtypes of normal Type Ia: core normal and broad line. The excess emission may be attributable to the radioactive decay of surface iron-peak elements as well as the interaction of ejecta with either the binary companion or a small torus of circumstellar material. Nebular-phase limits on Hα and He ι favor a white dwarf companion, consistent with the small companion size constrained by the low early SN luminosity, while the absence of [O ι] and He ι disfavors a violent merger of the progenitor. Of the two main explosion mechanisms proposed to explain the distribution of surface iron-peak elements in SN 2018aoz, the asymmetric Chandrasekhar-mass explosion is less consistent with the progenitor constraints and the observed blueshifts of nebular-phase [Fe ιι] and [Ni ιι]. The helium-shell double-detonation explosion is compatible with the observed lack of C spectral features, but current 1D models are incompatible with the infant-phase excess emission, $B_{max}$--$V_{max}$ color, and weak strength of nebular-phase [Ca ιι]. Although the explosion processes of SN 2018aoz still need to be more precisely understood, the same processes could produce a significant fraction of Type Ia SNe that appear to be normal after ~1 day.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1995

Hubble Space Telescope observations of the dwarf Nova Z Chamaeleontis through two eruption cycles

We have obtained the first high-speed photometry of the eclipsing dwarf nova Z Cha at ultraviolet wavelengths with the Hubble Space Telescope (HST). We observed the eclipse roughly every 4 days over two cycles of the normal eruptions of Z Cha, giving a uniquely complete coverage of its outburst cycle. The accretion disk dominated the ultraviolet light curve of Z Cha at the peak of an eruption; the white dwarf, the bright spot on the edge of the disk, and the boundary layer were all invisible. We were able to obtain an axisymmetric map of the accretion disk at this time only by adopting a flared disk with an opening angle of approximately 8 deg. The run of brightness temperature with radius in the disk at the peak of the eruption was too flat to be consistent with a steady state, optically thick accretion disk. The local rate of mass flow through the disk was approximately 5 x 10(exp -10) solar masses/yr near the center of the disk and approximately 5 x 10(exp -9) solar masses/yr near the outer edge. The white dwarf, the accretion disk, and the boundary layer were all significant contributors to the ultraviolet flux on the descending branches of the eruptions. The temperature of the white dwarf during decline was 18,300 K less than T(sub wd) less than 21,800 K, which is significantly greater than at minimum light. Six days after the maximum of an eruption Z Cha has faded to near minimum light at ultraviolet wavelenghts, but was still approximately 70% brighter at minimum light in the B band. About one-quarter of the excess flux in the B band came from the accretion disk. Thus, the accretion disk faded and became invisible at ultraviolet wavelengths before it faded at optical wavelenghts. The disk did, however, remain optically thick and obscured the lower half of the white dwarf at ultraviolet and possibly at optical wavelenghts for 2 weeks after the eruption ended. By the third week after eruptiuons the eclipse looked like a simple occultation of an unobscured, spherical white dwarf by a dark secondary star. The center of the accretion disk was, therfore, optically thin at ultraviolet wavelenghts and the boundary layer was too faint to be visible.

Robinson, E. L.↗
NASA NTRS2011

Broadband Achromatic Phase Shifter for a Nulling Interferometer

Nulling interferometry is a technique for imaging exoplanets in which light from the parent star is suppressed using destructive interference. Light from the star is divided into two beams and a phase shift of radians is introduced into one of the beams. When the beams are recombined, they destructively interfere to produce a deep null. For monochromatic light, this is implemented by introducing an optical path difference (OPD) between the two beams equal to lambda/2, where lambda is the wavelength of the light. For broadband light, however, a different phase shift will be introduced at each wavelength and the two beams will not effectively null when recombined. Various techniques have been devised to introduce an achromatic phase shift a phase shift that is uniform across a particular bandwidth. One popular technique is to use a series of dispersive elements to introduce a wavelength-dependent optical path in one or both of the arms of the interferometer. By intelligently choosing the number, material and thickness of a series of glass plates, a nearly uniform, arbitrary phase shift can be introduced between two arms of an interferometer. There are several constraints that make choosing the number, type, and thickness of materials a difficult problem, such as the size of the bandwidth to be nulled. Several solutions have been found for bandwidths on the order of 20 to 30 percent (Delta(lambda)/lambda(sub c)) in the mid-infrared region. However, uniform phase shifts over a larger bandwidth in the visible regime between 480 to 960 nm (67 percent) remain difficult to obtain at the tolerances necessary for exoplanet detection. A configuration of 10 dispersive glass plates was developed to be used as an achromatic phase shifter in nulling interferometry. Five glass plates were placed in each arm of the interferometer and an additional vacuum distance was also included in the second arm of the interferometer. This configuration creates a phase shift of pi radians with an average error of 5.97 x 10(exp -8) radians and standard deviation of 3.07 x 10(exp -4) radians. To reduce ghost reflections and interference effects from neighboring elements, the glass plates are tilted such that the beam does not strike each plate at normal incidence. Reflections will therefore walk out of the system and not contribute to the intensity when the beams are recombined. Tilting the glass plates, however, introduces several other problems that must be mitigated: (1) the polarization of a beam changes when refracted at an interface at non-normal incidence; (2) the beam experiences lateral chromatic spread as it traverses multiple glass plates; (3) at each surface, wavelength- dependent intensity losses will occur due to reflection. For a fixed angle of incidence, each of these effects must be balanced between each arm of the interferometer in order to ensure a deep null. The solution was found using a nonlinear optimization routine that minimized an objective function relating phase shift, intensity difference, chromatic beam spread, and polarization difference to the desired parameters: glass plate material and thickness. In addition to providing a uniform, broadband phase shift, the configuration achieves an average difference in intensity transmission between the two arms of the interferometer of 0.016 percent with a standard deviation of 3.64 x 10(exp -4) percent, an average difference in polarization between the two arms of the interferometer of 5.47 x 10(exp -5) percent with a standard deviation of 1.57 x 10(exp -6) percent, and an average chromatic beam shift between the two arms of the interferometer of -47.53 microns with a wavelength-by-wavelength spread of 0.389 microns.

Bolcar, Matthew R.↗
DOE OSTI2025

Correlations between maximum mass of neutron stars and the nuclear matter properties and the constraints from PSR J0740+6620 and GW190814

The correlations between the maximum mass of neutron stars (NSs) with the properties of nuclear matter at saturation density ρ 0 , i.e., the incompressibility K 0 , the symmetry energy E sym (ρ 0 ), and its slope L(ρ 0 ), the isoscalar effective mass m$^{*}_{s0}$, and the parameter f I , have been investigated using several sets of Skyrme interactions based on the momentum-dependent (MD) Skyrme interactions SAMi-J27 (soft), SAMi-J31 (stiff), and SAMi-J35 (super-stiff), as well as the momentum-independent (MI) Skyrme interaction SkT5. Only one physical quantity is varied at a time for a given set of Skyrme interactions. It is found that the maximum mass of NSs has (1) a strong positive linear correlation with K 0 for stiff and super-stiff equation of state (EOS) and no correlation with K 0 for soft EOS; (2) a negative logarithm correlation with the symmetry energy E sym (ρ 0 ) for all the EOSs; (3) a positive logarithm correlation with the slope L(ρ 0 ) for all the EOSs; (4) a positive logarithm correlation with the parameter f I for MD EOS; (5) a negative power-law correlation with m$^{*}_{s0}$ for MD EOS. Based on these established correlations, the constraints to the nuclear matter properties from the mass of PSR J0740+6620 and GW190814 are investigated. Furthermore, we find that while rotation enhances the maximum mass of NSs, it does not effect the underlying correlations between the maximum mass and the aforementioned nuclear matter parameters.

Zhou, M. [Shaanxi Normal University, Xi’an (China)↗
DOE OSTI2021

Auxiliary function approach for determining symmetry energy at suprasaturation densities

Nuclear symmetry energy E sym (ρ) at density ρ is normally expanded or simply parameterized as a function of χ = (ρ – ρ 0 )/3ρ 0 in the form of Esym(ρ) ≈ S + Lχ + 2 –1 K sym χ 2 + 6 –1 J sym χ 3 + … using its magnitude S, slope L, curvature K sym and skewness J sym at the saturation density ρ0 of nuclear matter. Much progress has been made in recent years in constraining especially the S and L parameters using various terrestrial experiments and astrophysical observations. Furthermore, such kind of expansions/parameterizations do not converge at supra-saturation densities where χ is not small enough, hindering an accurate determination of high-density E sym (ρ) even if its characteristic parameters at ρ0 are all well determined by experiments/observations. By expanding the E sym (ρ) in terms of a properly chosen auxiliary function Π sym (χ,Θ sym ) with a parameter Θsym fixed accurately by an experimental E sym (ρ r ) value at a reference density ρr, we show that the shortcomings of the χ-expansion can be completely removed or significantly reduced in determining the high-density behavior of E sym (ρ). In particular, using two significantly different auxiliary functions, we show that the new approach effectively incorporates higher χ-order contributions and converges to the same E sym (ρ) much faster than the conventional χ-expansion at densities ≲3ρ 0 . Several quantitative demonstrations using Monte Carlo simulations are given.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2016

The First Neptune Analog or Super-Earth with a Neptune-Like Orbit: MOA-2013-BLG-605Lb

We present the discovery of the first Neptune analog exoplanet or super-Earth with a Neptune-like orbit, MOA- 2013-BLG-605Lb. This planet has a mass similar to that of Neptune or a super-Earth and it orbits at 9 approximately 14 times the expected position of the snow line, a(sub snow), which is similar to Neptuneʼs separation of 11 a(sub snow) from the Sun. The planet/host-star mass ratio is q = (3.6 +/- 0.7) × 10(exp −4) and the projected separation normalized by the Einstein radius is s = 2.39 +/- 0.05. There are three degenerate physical solutions and two of these are due to a new type of degeneracy in the microlensing parallax parameters, which we designate "the wide degeneracy." The three models have (i) a Neptune-mass planet with a mass of M(sub p) = 21(+6/-7)(M) orbiting a low-mass M-dwarf with a mass of M(sub h) = 0.19(+0.05/-0.06 (solar mass)), (ii) a mini-Neptune with M(sub p) = 7.9(+1.8/-1.5)(M)) orbiting a brown dwarf host with M(sub h) = 0.068(+0.019/-0.011(solar mass)), and (iii) a super-Earth with M(sub p) = 3.2(+0.5/-0.3(M)) orbiting a low-mass brown dwarf host with M(sub h) = 0.025(+0.005/-0.004)(solar mass)), which is slightly favored. The 3D planet-host separations are 4.6(+4.7/-1.2)au, 2.1(+1.0/-0.2)au, and 0.94(+0.67/-0.02)au, which are 8.9(+10.5/-1.4)m 12(+7/-1), or 14(+11/-1) times larger than a(sub snow) for these models, respectively. Keck adaptive optics observations confirm that the lens is faint. This discovery suggests that low-mass planets with Neptune-like orbits are common. Therefore processes similar to the one that formed Neptune in our own solar system or cold super-Earths may be common in other solar systems.

Einstein radius↗
DOE OSTI2020

An Asymmetric Eclipse Seen toward the Pre-main-sequence Binary System V928 Tau

K2 observations of the weak-lined T Tauri binary V928 Tau A and B show the detection of a single, asymmetric eclipse, which may be due to a previously unknown substellar companion eclipsing one component of the binary with an orbital period >66 days. Over an interval of about 9 hr, one component of the binary dims by around 60%, returning to its normal brightness about 5 hr later. From modeling of the eclipse shape, we find evidence that the eclipsing companion may be surrounded by a disk or a vast ring system. The modeled disk has a radius of 0.9923 ± 0.0005 R∗, with an inclination of 56.°78 ± 0.°03, a tilt of 41.°22 ± 0.°05, an impact parameter of -0.2506 ± 0.0002 R∗, and an opacity of 1.00. The occulting disk must also move at a transverse velocity of 6.637 ± 0.002 R∗ day -1 , which, depending on whether it orbits V928 Tau A or B, corresponds to approximately 73.53 or 69.26 km s -1 . A search in ground-based archival data reveals additional dimming events, some of which suggest periodicity, but no unambiguous period associated with the eclipse observed by K2. We present a new epoch of astrometry that is used to further refine the orbit of the binary, presenting a new lower bound of 67 yr, and constraints on the possible orbital periods of the eclipsing companion. The binary is also separated by 18' (~2250 au) from the lower-mass CFHT-BD-Tau 7, which is likely associated with V928 Tau A and B. We also present new high-dispersion optical spectroscopy that we use to characterize the unresolved stellar binary.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2005

Influence of Coronal Abundance Variations

The PI of this project was Jeff Scargle of NASA/Ames. Co-I's were Alma Connors of Eureka Scientific/Wellesley, and myself. Part of the work was subcontracted to Eureka Scientific via SAO, with Vinay Kashyap as PI. This project was originally assigned grant number NCC2-1206, and was later changed to NCC2-1350 for administrative reasons. The goal of the project was to obtain, derive, and develop statistical and data analysis tools that would be of use in the analyses of high-resolution, high-sensitivity data that are becoming available with new instruments. This is envisioned as a cross-disciplinary effort with a number of "collaborators" including some at SA0 (Aneta Siemiginowska, Peter Freeman) and at the Harvard Statistics department (David van Dyk, Rostislav Protassov, Xiao-li Meng, Epaminondas Sourlas, et al). We have developed a new tool to reliably measure the metallicities of thermal plasma. It is unfeasible to obtain high-resolution grating spectra for most stars, and one must make the best possible determination based on lower-resolution, CCD-type spectra. It has been noticed that most analyses of such spectra have resulted in measured metallicities that were significantly lower than when compared with analyses of high- resolution grating data where available (see, e.g., Brickhouse et al., 2000, ApJ 530,387). Such results have led to the proposal of the existence of so-called Metal Abundance Deficient, or "MAD" stars (e.g., Drake, J.J., 1996, Cool Stars 9, ASP Conf.Ser. 109, 203). We however find that much of these analyses may be systematically underestimating the metallicities, and using a newly developed method to correctly treat the low-counts regime at the high-energy tail of the stellar spectra (van Dyk et al. 2001, ApJ 548,224), have found that the metallicities of these stars are generally comparable to their photospheric values. The results were reported at the AAS (Sourlas, Yu, van Dyk, Kashyap, and Drake, 2000, BAAS 196, v32, #54.02), and at the conference on Statistical Challenges in Modem Astronomy (Sourlas, van Dyk, Kashyap, Drake, and Pease, 2003, SCMA 111, Eds. E.D.Feigelson, G.J.Babu, New York:Springer, p489-490). We also described the limitations of one of the most egregiously misused and misapplied statistical tests in astrophysical literature, the F-test for verifying model components (Protassov, van Dyk, Connors, Kashyap, and Siemiginowska, 2002, ApJ, 571,545). Indeed, a search through the ApJ archives turned up 170 papers in the 5 previous years that used the F-test explicitly in some form or the other, and with the vast majority of them not using it correctly! Indeed, looking at just 4 issues of the ApJ in 2001, we found 13 instances of its use, of which nine were demonstrably incorrect. Clearly, it is difficult to understate the importance of this issue. We also worked on speeding up Bayes Blocks and Sparse Bayes Blocks algorithms to make them more tractable for large searches. We also supported staistics students and postdocs in both explicit physics- model-based (spectra with tens of thousands of atomic lines) and "model-free" -- i.e. non-parametric or semi-parametric -- algorithms. Work on using more of the latter is just beginning; while using multi-scale methods for Poisson imaging has come to hition. In fact, "An Image Restoration Technique with Error Estimates", by D. Esch, A. Connors, M. Karovska, and D. van Dyk, was published by ApJ (Esch et a1.2004, ApJ, 610, 1213). The code has been delivered to M. Karovska for CXC; and is available for beta-testing upon request. The other large project we worked on was on the self-consistent modeling of logN-logs curves in the Poisson limit. logN-logs curves are a fundamental tool in the study of source populations, luminosity functions, and cosmological parameters. However, their determination is hampered by statistical effects such as the Eddington bias, incompleteness due to detection efficiency, faint source flux fluctuations, etc. We have develed a new and powerful method using the full Poisson machinery that allows us to model the logN-logs distribution of X-ray sources in a self-consistent manner. Because we properly account for all the above statistical effects, our modeling is valid over the full range of the data, and not just for strong sources, as is normally done. Using a Bayesian approach and modeling the fluxes with known functional forms such as simple or broken power-laws, and conditioning the expected photon counts on the fluxes, the background contamination, effective area, detector vignetting, and detection probability, we can delve deeply into the low counts regime and extend the usefulness of medium sensitivity surveys such as ChAMP by orders of magnitude. The built-in flexibility of the algorithm also allows a simultaneous analysis of multiple datasets. We have applied this analysis to a set a Chandra observations (Sourlas, Kashyap, Zezas, van Dyk, 2004, HEAD #8, #16.32)

Scargle, Jeffrey D.↗
DOE OSTI2025

Onset of Constituent Quark Number Scaling in Heavy-Ion Collisions at RHIC

Partonic collectivity is one of the necessary signatures for the formation of quark-gluon plasma in high-energy nuclear collisions. Number of constituent quarks (NCQ) scaling has been observed for hadron elliptic flow 𝑣2 in top energy nuclear collisions at the Relativistic Heavy Ion Collider and the LHC, and this has been theoretically suggested as strong evidence for partonic collectivity. In this Letter, a systematic analysis of 𝑣 2 of 𝜋 ± , 𝐾 ± , 𝐾$^{0}_{𝑆}$, 𝑝, and Λ in Au+Au collisions at $\sqrt{s_{NN}}$ = 3.2, 3.5, 3.9, and 4.5 GeV, with the STAR experiment at the Relativistic Heavy Ion Collider, is presented. NCQ scaling is markedly violated at 3.2 GeV, consistent with a hadronic-interaction dominated equation of state. However, as the collision energy increases, a gradual evolution to NCQ scaling is observed. In conclusion, this beam-energy dependence of 𝑣 2 for all hadrons studied provides evidence for the onset of dominant partonic interactions by $\sqrt{s_{NN}}$ = 4.5 GeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2007

Radiation Test Results for a MEMS Microshutter Operating at 60 K

The James Webb Space Telescope (JWST), the successor to the Hubble Space Telescope, is due to be launched in 2013 with the goal of searching the very distant Universe for stars that formed shortly after the Big Bang. Because this occurred so far back in time, the available light is strongly red-shifted, requiring the use of detectors sensitive to the infrared portion of the electromagnetic spectrum. HgCdTe infrared focal plane arrays, cooled to below 30 K to minimize noise, will be used to detect the faint signals. One of the instruments on JWST is the Near Infrared Spectrometer (NIRSPEC) designed to measure the infrared spectra of up to 100 separate galaxies simultaneously. A key component in NIRSPEC is a Micro-Electromechanical System (MEMS), a two-dimensional micro-shutter array (MSA) developed by NASA/GSFC. The MSA is inserted in front of the detector to allow only the light from the galaxies of interest to reach the detector and to block the light from all other sources. The MSA will have to operate at 30 K to minimize the amount of thermal radiation emitted by the optical components from reaching the detector array. It will also have to operate in the space radiation environment that is dominated by the MSA will be exposed to a large total ionizing dose of approximately 200 krad(Si). Following exposure to ionizing radiation, a variety of MEMS have exhibited performance degradation. MEMS contain moving parts that are either controlled or sensed by changes in electric fields. Radiation degradation can be expected for those devices where there is an electric field applied across an insulating layer that is part of the sensing or controlling structure. Ionizing radiation will liberate charge (electrons and holes) in the insulating layers, some of which may be trapped within the insulating layer. Trapped charge will partially cancel the externally applied electric field and lead to changes in the operation of the MEMS. This appears to be a general principle for MEMS. Knowledge of the above principle has raised the concern at NASA that the MSA might also exhibit degraded performance because, i) each shutter flap is a multilayer structure consisting of metallic and insulating layers and ii) the movement of the shutter flaps is partially controlled by the application of an electric field between the shutter flap and the substrate (vertical support grid). The whole mission would be compromised if radiation exposure were to prevent the shutters from opening and closing properly. energetic ionizing particles. Because it is located A unique feature of the MSA is that, as outside the spacecraft and has very little shielding, previously mentioned, it will have to operate at temperatures near 30 K. To date, there are no published reports on how very low temperatures (- 30K) affect the response of MEMS devices to total ionizing dose. Experiments on SiO2 structures at low temperatures (80 K) indicate that the electrons generated by the ionizing radiation are mobile and will move rapidly under the application of an external electric field. Holes, on the other hand, that would normally move in the opposite direction through the SiO2 via a "thermal hopping" process, are effectively immobile at low electric fields as they are trapped close to their generation sites. However, for sufficiently large electric fields (greater than 3 MV/cm) holes are able to move through the SiO2. The larger the field, the more rapidly the holes move. The separation of the electrons and holes leads to a reduced electric field within the insulating layer. To overcome this reduction in electric field, a greater external voltage will have to be applied that alters the normal operation of the device. This report presents the results of radiation testing of the MSA at 60 K. The temperature was higher than the targeted temperature because of a faulty electrical interconnect on the test board. Specifically, our goal was to determine whether the MSA would function propey after a TID of 200 krad(Si).

Rapchun, David A.↗