Search NASA⌕ Search

SEARCH · Search NASA

Results for “astronomy computing”

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 361 records · Page 20

UVSTAR: An imaging spectrograph with telescope for the Shuttle Hitchhiker-M platform

UVSTAR is an EUV spectral imager intended as a facility instrument devoted to solar system and astronomy studies. It covers the wavelength range of 500 to 1250 A, with sufficient spectral resolution to separate emission lines and to form spectrally resolved images of extended plasma sources. Targets include the Io plasma torus at Jupiter, hot stars, planetary nebulae and bright galaxies. UVSTAR consists of a pair of telescopes and concave grating spectrographs that cover the overlapping spectral ranges of 500-900 and 850-1250 A. The telescopes use two 30 cm diameter off-axis paraboloids having focal length of 1.5 m. An image of the target is formed at the entrance slits of the two concave grating spectrographs. The gratings provide dispersion and re-image the slits at the detectors, intensified CCD's. The readout format of the detectors can be chosen by computer, and three slit widths are selectable to adapt the instrument to specific tasks. UVSTAR has internal gimbals which allow rotation of plus or minus 3 deg about each of two axes. Dedicated finding and tracking telescopes will acquire and track the target after rough pointing is achieved by orienting the Orbiter. Responsibilities for implementation and utilization of UVSTAR are shared by groups in Italy and the U.S. UVSTAR is scheduled for flight in early 1995, timed for an opportunity to observe the Jovian system.

Stalio, Roberto↗

Population-level Dark Energy Constraints from Strong Gravitational Lensing using Simulation-Based Inference

In this work, we present a scalable approach for inferring the dark energy equation-of-state parameter ($w$) from a population of strong gravitational lens images using Simulation-Based Inference (SBI). Strong gravitational lensing offers crucial insights into cosmology, but traditional Monte Carlo methods for cosmological inference are computationally prohibitive and inadequate for processing the thousands of lenses anticipated from future cosmic surveys. New tools for inference, such as SBI using Neural Ratio Estimation (NRE), address this challenge effectively. By training a machine learning model on simulated data of strong lenses, we can learn the likelihood-to-evidence ratio for robust inference. Our scalable approach enables more constrained population-level inference of $w$ compared to individual lens analysis, constraining $w$ to within $1\sigma$.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Problem-tailored Simulation of Energy Transport on Noisy Quantum Computers

The transport of conserved quantities like spin and charge is fundamental to characterizing the behavior of quantum many-body systems. Numerically simulating such dynamics is generically challenging, which motivates the consideration of quantum computing strategies. However, the relatively high gate errors and limited coherence times of today's quantum computers pose their own challenge, highlighting the need to be frugal with quantum resources. In this work we report simulations on quantum hardware of infinite-temperature energy transport in the mixed-field Ising chain, a paradigmatic many-body system that can exhibit a range of transport behaviors at intermediate times. We consider a chain with L = 12 sites and find results broadly consistent with those from ideal circuit simulators over 90 Trotter steps, containing up to 990 entangling gates. To obtain these results, we use two key problem-tailored insights. First, we identify a convenient basis – the Pauli Y basis – in which to sample the infinite-temperature trace and provide theoretical and numerical justifications for its efficiency relative to, e.g., the computational basis. Second, in addition to a variety of problem-agnostic error mitigation strategies, we employ a renormalization strategy that compensates for global nonconservation of energy due to device noise. We discuss the applicability of the proposed sampling approach beyond the mixed-field Ising chain and formulate a variational method to search for a sampling basis with small sample-to-sample fluctuations for an arbitrary Hamiltonian. This opens the door to applying these techniques in more general models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Measurement of the small-scale 3D Lyman- α forest power spectrum

Small-scale correlations measured in the Lyman-α (Lyα) forest encode information about the intergalactic medium and the primordial matter power spectrum. In this article, we present and implement a simple method to measure the 3-dimensional power spectrum, P 3D , of the Lyα forest at wavenumbers k corresponding to small, ~ Mpc scales. In order to estimate P 3D from sparsely and unevenly distributed data samples, we rely on averaging 1-dimensional Fourier Transforms, as previously carried out to estimate the 1-dimensional power spectrum of the Lyα forest, P 1D . Further, this methodology exhibits a very low computational cost. We confirm the validity of this approach through its application to Nyx cosmological hydrodynamical simulations. Subsequently, we apply our method to the eBOSS DR16 Lyα forest sample, providing as a proof of principle, a first P 3D measurement averaged over two redshift bins z = 2.2 and z = 2.4. This work highlights the potential for forthcoming P 3D measurements, from upcoming large spectroscopic surveys, to untangle degeneracies in the cosmological interpretation of P 1D .

79 ASTRONOMY AND ASTROPHYSICS↗

Fortran for the nineties

Fortran has largely enjoyed prominence for the past few decades as the computer programming language of choice for numerically intensive scientific, engineering, and process control applications. Fortran's well understood static language syntax has allowed resulting parsers and compiler optimizing technologies to often generate among the most efficient and fastest run-time executables, particularly on high-end scalar and vector supercomputers. Computing architectures and paradigms have changed considerably since the last ANSI/ISO Fortran release in 1978, and while FORTRAN 77 has more than survived, it's aged features provide only partial functionality for today's demanding computing environments. The simple block procedural languages have been necessarily evolving, or giving way, to specialized supercomputing, network resource, and object-oriented paradigms. To address these new computing demands, ANSI has worked for the last 12-years with three international public reviews to deliver Fortran 90. Fortran 90 has superseded and replaced ISO FORTRAN 77 internationally as the sole Fortran standard; while in the US, Fortran 90 is expected to be adopted as the ANSI standard this summer, coexisting with ANSI FORTRAN 77 until at least 1996. The development path and current state of Fortran will be briefly described highlighting the many new Fortran 90 syntactic and semantic additions which support (among others): free form source; array syntax; new control structures; modules and interfaces; pointers; derived data types; dynamic memory; enhanced I/O; operator overloading; data abstraction; user optional arguments; new intrinsics for array, bit manipulation, and system inquiry; and enhanced portability through better generic control of underlying system arithmetic models. Examples from dynamical astronomy, signal and image processing will attempt to illustrate Fortran 90's applicability to today's general scalar, vector, and parallel scientific and engineering requirements and object oriented programming paradigms. Time permitting, current work proceeding on the future development of Fortran 2000 and collateral standards will be introduced.

Himer, J. T.↗

Predictive Thermal Control (PTC) Technology to Enable Thermally Stable Telescopes: First Two Year Status

Predictive Thermal Control Technology (PTCT) development project is a multiyear effort initiated in Fiscal Year (FY) 2017, to mature the Technology Readiness Level (TRL) of critical technologies required to enable ultra-thermally-stable ultraviolet/optical/infrared (UVOIR) space telescope primary-mirror assemblies for ultra-high-contrast observations of exoplanets. Key accomplishments of 2017 to 2019 include: creating a high-fidelity STOP (Structural, Thermal, and Optical Performance) model of the AMTD-2 (Advance Mirror Technology Development) 1.5-meter Ultra-Low Expansion (ULE®) mirror (manufactured by Harris Corp) by merging 3D X-Ray computed tomography data of the ‘as-built’ mirror and coefficient of thermal expansion (CTE) data maps for each of the 18 core elements; partially validating this model by measuring the mirror’s response to bulk temperature changes and lateral thermal gradients; designed and built (with PTC partner Harris Corp) a 1.5-meter enclosure with 26 actively-control thermal zones; and defined specifications for a potential 4-meter primary mirror thermal enclosure for the Habitable Exoplanet (HabEx) Imager mission.

Astrophysics↗

Exact evaluation of large-charge correlation functions in nonrelativistic conformal field theory

The large-charge master field which generates all n -point correlation functions with an insertion of large charge Q in nonrelativistic conformal field theory is obtained. This field is used to compute Schrödinger-invariant n -point correlation functions of large-charge operators via a direct evaluation of the path integral. Conformal dimensions are found to agree with calculations based on the state-operator correspondence. The master field solution exhibits an emergent harmonic trap whose frequency is a function of the Euclidean time. The large-charge effective action with operator insertions describes a droplet of superfluid matter whose spatial size scales with the time separation of sources. The solution is used to compute Schrödinger symmetry breaking corrections in the large-charge effective field theory (EFT) due to a finite scattering length in the fundamental theory of fermions near unitarity. The scaling of these effects in the large-charge power counting scheme is established, and the size of the effects is quantified using input from quantum Monte Carlo simulations of the near-unitary gas, as well as from the large- N expansion at large charge. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗

A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be-+12.390.981.30km and-+2.0720.0660.067Me respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is=-+TlogK5.99100.060.05([])for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters(within the prior support)and can thus be computed from the posterior samples.

Millisecond pulsars↗

New Voyager radio spectrograms of Uranus

New, high-resolution spectrograms of the Voyager-2 radio observations at Uranus were produced from the original, six-second Planetary Radio Astronomy (PRA) data and these show a number of new features which were not obvious in previous versions. Among these new features are the detailed structure of the so-called broadband-bursty (b-bursty) emissions, unexpected sloping striations in the smooth high-frequency (SHF) component, and the overlap of these two components during the first rotation after closest approach. In addition, a slightly different planetary rotation rate from the b-bursty emissions, was found, and at the initial onset of the SHF component, what appears to be the shadow of a Uranian plasmasphere. These new spectrograms were prepared using a special dithering algorithm to show signal strengths as gray shadings, and the data were also manually cleaned to suppress noise and interference. This produced spectrograms of exceptional quality and certain details of their production on a stand-alone personal computer are also discussed.

Calvert, W.↗

Ray tracing of Jovian kilometric radiation

Results of computer ray tracing of Jovian kilometric radiation from 56.2 kHz to 1 MHz in a model Jovian magnetosphere with an Io torus are presented. Ray tracing calculations indicate that the Io torus presents a propagation barrier to the radiation and that the Jovian kilometric radiation must be generated in the L-O mode from a source near Jupiter on field lines passing through the Io torus. One effect of the Io torus is to refract the rays away from the magnetic equator forming a shadow zone at radial distances beyond the torus. In general, at radial distances greater than 10 Jovian radii, as the wave frequency increases (greater than 200 kHz) so does the magnetic latitude of the shadow zone. These and other features of the ray tracing calculations are in good qualitative agreement with the observations from the plasma wave receiver and planetary radio astronomy experiment on board both Voyagers 1 and 2.

Green, J. L.↗

Emission from small dust particles in diffuse and molecular cloud medium

Infrared Astronomy Satellite (IRAS) observations of the whole galaxy has shown that long wavelength emission (100 and 60 micron bands) can be explained by thermal emission from big grains (approx 0.1 micron) radiating at their equilibrium temperature when heated by the InterStellar Radiation Field (ISRF). This conclusion has been confirmed by continuum sub-millimeter observations of the galactic plane made by the EMILIE experiment at 870 microns (Pajot et al. 1986). Nevertheless, shorter wavelength observations like 12 and 25 micron IRAS bands, show an emission from the galactic plane in excess with the long wavelength measurements which can only be explained by a much hotter particles population. Because dust at equilibrium cannot easily reach high temperatures required to explain this excess, this component is thought to be composed of very small dust grains or big molecules encompassing thermal fluctuations. Researchers present here a numerical model that computes emission, from Near Infrared Radiation (NIR) to Sub-mm wavelengths, from a non-homogeneous spherical cloud heated by the ISRF. This model fully takes into account the heating of dust by multi-photon processes and back-heating of dust in the Visual/Infrared Radiation (VIS-IR) so that it is likely to describe correctly emission from molecular clouds up to large A sub v and emission from dust experiencing temperature fluctuations. The dust is a three component mixture of polycyclic aromatic hydrocarbons, very small grains, and classical big grains with independent size distributions (cut-off and power law index) and abundances.

Bernard, J. P.↗

Heavy dark matter in white dwarfs: multiple-scattering capture and thermalization

We present an improved treatment for the scattering of heavy dark matter from the ion constituents of a white dwarf. In the heavy dark matter regime, multiple collisions are required for the dark matter to become gravitationally captured. Our treatment incorporates all relevant physical effects including the dark matter trajectories, nuclear form factors, and radial profiles for the white dwarf escape velocity and target number densities. Our capture rates differ by orders of magnitude from previous estimates, which have typically used approximations developed for dark matter scattering in the Earth. We also compute the time for the dark matter to thermalize in the center of the white dwarf, including in-medium effects such as phonon emission and absorption from the ionic lattice in the case where the star has a crystallized core. We find much shorter thermalization timescales than previously estimated, especially if the white dwarf core has crystallized. We illustrate the importance of our improved approach by determining the cross section required for accumulated asymmetric dark matter to self-gravitate.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Toward a Ginsparg-Wilson lattice Hamiltonian

To address quantum computation of quantities in quantum chromodynamics (QCD) for which chiral symmetry is important, it would be useful to have the Hamiltonian for a fermion satisfying the Ginsparg-Wilson (GW) equation. I work with a solution to the GW equation which is fractional linear in time derivatives. The resulting Hamiltonian is nonlocal and has ghosts, but is free of doublers and has the correct continuum limit. This construction works in general odd spatial dimensions, and I provide an explicit expression for the Hamiltonian in one spatial dimension. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

A Model for Removing Fringes from STIS Slitless Spectra

We have developed a model that allows us to defringe slitless 2-dimensional spectra taken with the Space Telescope Imaging Spectrograph (STIS). Like all thin chip CCDs used in astronomy, the STIS CCD detector acts like an interferometer in the presence of monochromatic light, producing bright and dark fringes where there is constructive and destructive interference. These fringes are especially troublesome for spectra in the near-IR (> 7000A), reaching a peak amplitude of about 20% near 9500 A. Removing fringes from STIS spectra taken with a slit is simply a matter of dividing by an appropriate "fringe flat". A suitable flat is obtained by observing a continuum calibration source through the same slit. This however is not possible for slitless spectra, whose wavelength mapping onto the CCD detector varies with position of the object in the field. Beginning with an approximate knowledge of the STIS CCD structure, we have used 50 continuum source spectral flats taken at various central wavelengths to constrain the thickness of the detection layer at each pixel. Rom this result, we can compute a fringe flat for any astronomical source, no matter how it is positioned in the field. Flats made by this method reduce the fringing amplitude in astronomical spectra by a factor of about 5. The model fringe flat is a function of the index of refraction of each of four layers in the CCD, the thickness of each layer, and the wavelength of the light hitting each pixel. We show how the CCD structural parameters are incorporated into this "fitting function", and describe the procedure used to solve for the free parameters (the thickness of each layer at each pixel). The effectiveness of this technique will be demonstrated using slitless spectra taken as part of the STIS parallel program.

Malumuth, E.M.↗

Thermodynamics of the near-extremal Kerr spacetime

We examine the thermodynamics of a near-extremal Kerr black hole, and demonstrate that the geometry behaves as an ordinary quantum system with a vanishingly small degeneracy at low temperatures. This is in contrast with the classical analysis, which instead predicts a macroscopic entropy for the extremal Kerr black hole. Our results follow from a careful analysis of the gravitational path integral. Specifically, the low temperature canonical partition function behaves as Z ~ T$\frac{3}{2}$ e S 0 +c log S 0 , with S 0 the classical degeneracy and c a numerical coefficient we compute. This is in line with the general expectations for non-supersymmetric near-extremal black hole thermodynamics, as has been clarified in the recent past, although cases without spherical symmetry have not yet been fully analyzed until now. We also point out some curious features relating to the rotational zero modes of the near-extremal Kerr black hole background that affects the coefficient c. This raises a puzzle when considering similar black holes in string theory. Our results generalize to other rotating black holes, as we briefly exemplify.

79 ASTRONOMY AND ASTROPHYSICS↗

Cosmic Ray Mediated Thermal Fronts in the Warm-Hot Circumgalactic Medium

We investigate the 1D plane-parallel front connecting the warm ($10^4$ K) and hot ($10^6$ K) phases of the circumgalactic medium (CGM), focusing on the influence of cosmic rays (CRs) in shaping these transition layers. We find that cosmic rays dictate the thermal balance while other fluxes (thermal conduction, radiative cooling, and gas flow) adjust to compensate. We compute column density ratios for selected transition temperature ions and compare them with observational data. While most of our models fail to reproduce the observations, a few are successful, although we make no claims for their uniqueness. Some of the discrepancies may indicate challenges in capturing the profiles in cooler, photoionized regions, as has been suggested for by previous efforts to model thermal transition layers.

79 ASTRONOMY AND ASTROPHYSICS↗