Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “circular data”

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,459 records · Page 81

Multiscale modeling of HBV infection integrating intra- and intercellular viral propagation to analyze extracellular viral markers

Chronic infection with hepatitis B virus (HBV) is caused by the persistence of closed circular DNA (cccDNA) in the nucleus of infected hepatocytes. Despite available therapeutic anti-HBV agents, eliminating the cccDNA remains challenging. Thus, quantifying and understanding the dynamics of cccDNA are essential for developing effective treatment strategies and new drugs. However, such study requires repeated liver biopsy to measure the intrahepatic cccDNA, which is basically not accepted because liver biopsy is potentially morbid and not common during hepatitis B treatment. We here aimed to develop a noninvasive method for quantifying cccDNA in the liver using surrogate markers in peripheral blood. We constructed a multiscale mathematical model that explicitly incorporates both intracellular and intercellular HBV infection processes. The model, based on age-structured partial differential equations, integrates experimental data from in vitro and in vivo investigations. By applying this model, we roughly predicted the amount and dynamics of intrahepatic cccDNA within a certain range using specific viral markers in serum samples, including HBV DNA, HBsAg, HBeAg, and HBcrAg. Our study represents a significant step towards advancing the understanding of chronic HBV infection. The noninvasive quantification of cccDNA using our proposed method holds promise for improving clinical analyses and treatment strategies. By comprehensively describing the interactions of all components involved in HBV infection, our multiscale mathematical model provides a valuable framework for further research and the development of targeted interventions.

59 BASIC BIOLOGICAL SCIENCES↗

Evidence for volcanism in NW Ishtar Terra, Venus

Venera 15/16 radar data for an area in NW Ishtar Terra, Venus, show an area with moderate radar return and a smooth textured surface which embays low lying areas of the surrounding mountainous terrain. Although this unit may be an extension of the lava plains of Lakshmi Planum to the southeast, detailed study suggests a separate volcanic center in NW Ishtar Terra. Lakshmi Planum, on the Ishtar Terra highland, exhibits major volcanic and tectonic features. On the Venera radar image radar brightness is influenced by slope and roughness; radar-facing slopes (east-facing) and rough surfaces (approx. 8 cm average relief) are bright, while west-facing slopes and smooth surfaces are dark. A series of semi-circular features, apparently topographic depressions, do not conform in orientation to major structural trends in this region of NW Ishtar Terra. The large depression in NW Ishtar Terra is similar to the calderas of Colette and Sacajawea Paterae, as all three structures are large irregular depressions. NW Ishtar Terra appears to be the site of a volcanic center with a complex caldera structure, possibly more than one eruptive vent, and associated lobed flows at lower elevations. The morphologic similarity between this volcanic center and those of Colette and Sacajawea suggests that centralized eruptions have been the dominant form of volcanism in Ishtar. The location of this volcanic center at the intersection of two major compressional mountain belts and the large size of the calders (with an inferred larg/deep magma source) support a crustal thickening/melting rather than a hot-spot origin for these magmas.

Gaddis, L.↗

Errors in reconstruction of dichroic X-ray orientation tomography due to polarization rotation of the incident beam

Dichroic X-ray tomography is a technique in which the crystal orientation or magnetization of a sample is resolved in three dimensions. The best-known uses of this technique are for observation of magnetic moments via circular dichroism, using left- and right-handed circularly polarized X-ray beams. Another variant uses linear dichroism to resolve the crystal orientation. In both these techniques, it is assumed that the absorption of X-rays along a path inside a material can be computed as a line integral of a local absorption coefficient along the ray path. For linear dichroism, this assumption is inaccurate because the polarization of the beam changes along the propagation direction when the optic axis of the material is not aligned along the polarization. In this work, a finite-element Maxwell solver is used to simulate tomography and reconstructions. The propagation effect can lead to significant errors in the reconstructed orientations. These errors may be mitigated by taking data at additional angles or by operating at energies at which the dichroism is weak. An iterative approach is proposed which may allow accurate reconstruction with fewer data than would otherwise be required.

X-ray linear dichroism↗

Flare stars and solar bursts - High resolution in time and frequency

Coronal loops on the sun and nearby stars are studied using observations at 20-cm wavelength with high resolution in time and frequency. Observations of the dwarf M star AD Leonis with high time resolution led to the discovery of a quasi-periodic train of circularly polarized spikes with a mean periodicity of 32 + or - 5 ms and a total duration of 150 ms. Observations at closely spaced wavelengths using the VLA revealed a narrow-band structure in solar bursts and in the slowly varying radiation of the dwarf M star YZ Canis Minoris. This narrow-band emission may be attributed to electron-cyclotron maser radiation. Maser action at the second or first harmonic of the gyrofrequency implies magnetic field strengths of 250 and 500 G, respectively. Hence, observations with high resolution in time and frequency suggest coherent processs in the coronas of the sun and dwarf M stars. It is noted that the scientific potential suggested by the present data can be realized by the development of a solar-stellar synthesis radiotelescope. This instrument would be dedicated to solar and stellar observations with high angular, temporal, and frequency resolution.

Lang, K. R.↗

Progress in the prediction of unsteady heat transfer on turbines blades

Progress toward developing a general method for predicting unsteady heat transfer on turbine blades subject to blade-passing frequencies and Reynolds numbers relevant to the Space Shuttle Main Engine (SSME) is discussed. The method employs an invisid/viscous interactive procedure which has been tested extensively for steady subsonic and transonic external airfoil problems. One such example is shown. The agreement with experimental data and with Navier-Stokes calculations yields confidence in the method. The technique is extended to account for wake generated unsteadiness. The flow reversals around the stagnation point caused by the nonuniform onset velocity are accounted for by using the Characteristic Box scheme developed by Cebeci and Stewartson. The coupling between the inviscid and viscous methods is achieved by using a special procedure, which, with a novel inverse finite-difference boundary-layer method, allows the calculations to be performed for a wide range of flow conditions, including separation. Preliminary results are presented for the stagnation region of turbine blades for both laminar and turbulent flows. A laminar model problem corresponding to a flow on a circular cylinder which experiences the periodic passing of wakes from turbine blades is presented to demonstrate the ability of the method to calculate flow reversals around the stagnation region.

Cebeci, T.↗

Predicted performance of InP solar cells in Cassegrainian and slats space concentrator arrays at 20 to 100 AM0, 80 to 100 C

Researchers have calculated the expected performance dependence of near-optimally designed shallow homojunction n+pp+InP solar cells on incident intensities up 200 AM0 and temperatures up to 100 C (373K). Both circular and rectangular cells have been considered, the former for use in a Cassegrainian concentrator array at 100 AM0, 80 to 100 C and the latter for use in a Slats type concentrator array at 20 AM0 80 to 100 C. Calculation of the temperature dependence of the performance parameters I sub sc, V sub oc, FF and eta was done by first verifying that the use of the measured temperature variation of I sub sc, of the best published value of the temperature dependence of the bandgap of InP, and of the temperature dependences of the lifetimes and mobilities of electrons and holes the same as in equivalently doped GaAs, gave calculated results that closely matched measured data on the temperature variation of I sub sc, V sub oc, and FF of four existing InP cells at 1 AM0. It was then assumed that the same temperature dependences of I sub sc, the bandgap and lifetimes and mobilities would hold in the near-optimally designed cells at the higher concentrations.

Goradia, Chandra↗

Measurements of the Influence of Integral Length Scale on Stagnation Region Heat Transfer

The purpose was twofold: first, to determine if a length scale existed that would cause the greatest augmentation in stagnation region heat transfer for a given turbulence intensity and second, to develop a prediction tool for stagnation heat transfer in the presence of free stream turbulence. Toward this end, a model with a circular leading edge was fabricated with heat transfer gages in the stagnation region. The model was qualified in a low turbulence wind tunnel by comparing measurements with Frossling's solution for stagnation region heat transfer in a laminar free stream. Five turbulence generating grids were fabricated; four were square mesh, biplane grids made from square bars. Each had identical mesh to bar width ratio but different bar widths. The fifth grid was an array of fine parallel wires that were perpendicular to the axis of the cylindrical leading edge. Turbulence intensity and integral length scale were measured as a function of distance from the grids. Stagnation region heat transfer was measured at various distances downstream of each grid. Data were taken at cylinder Reynolds numbers ranging from 42,000 to 193,000. Turbulence intensities were in the range 1.1 to 15.9 percent while the ratio of integral length scale to cylinder diameter ranged from 0.05 to 0.30. Stagnation region heat transfer augmentation increased with decreasing length scale. An optimum scale was not found. A correlation was developed that fit heat transfer data for the square bar grids to within +4 percent. The data from the array of wires were not predicted by the correlation; augmentation was higher for this case indicating that the degree of isotropy in the turbulent flow field has a large effect on stagnation heat transfer. The data of other researchers are also compared with the correlation.

Vanfossen, G. James↗

CFD Simulations of the IHF Arc-Jet Flow: Compression-Pad/Separation Bolt Wedge Tests

This paper reports computational analyses in support of two wedge tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using two different wedge models, each placed in a free jet downstream of a corresponding different conical nozzle in the Ames 60-MW Interaction Heating Facility. Each panel test article included a metallic separation bolt imbedded in Orion compression-pad and heatshield materials, resulting in a circular protuberance over a flat plate. The protuberances produce complex model flowfields, containing shock-shock and shock-boundary layer interactions, and multiple augmented heating regions on the test plate. As part of the test calibration runs, surface pressure and heat flux measurements on water-cooled calibration plates integrated with the wedge models were also obtained. Surface heating distributions on the test articles as well as arc-jet test environment parameters for each test configuration are obtained through computational fluid dynamics simulations, consistent with the facility and calibration measurements. The present analysis comprises simulations of the nonequilibrium flowfield in the facility nozzle, test box, and flowfield over test articles, and comparisons with the measured calibration data.

CFD↗

American Made Challenges Battery Voucher Program Cooperative Research and Development Agreement (Cooperative Research and Development Final Report, CRADA Number CRD-21-17533)

Renewance is a Phase II winner of the U.S. Department of Energy Lithium-ion Battery Recycling Prize. The Prize is designed to incentivize American entrepreneurs to develop and demonstrate processes that, when scaled, have the potential to profitably capture 90% of all discarded or spent lithium-based batteries (LIB) in the Unites States for eventual recovery of key materials for re-introductions into the U.S. supply chain. The objective of this work is to enable a more efficient evaluation of battery sources for second life applications prior to ultimately being recycled, through evaluation of chemistry characteristics, projected battery lifetime, and application history. This work will develop the capability to identify groups of batteries that may be useful for second life and reduce the cost of end-of-life (EOL) LIB evaluation and repurposing. To meet the objective, NREL will use existing and new data to create a refined algorithm that could be used to evaluate batches of batteries for potential reuse based on manufacturing date and historical use characteristics. Based on current battery market prices and compiled literature data, a starting-point estimate of the market value of the batteries for reuse based on expected lifetime will be included in the algorithm. With the projected surge in LIB demand, battery second life is a new area ripe for development and investment from companies like Renewance. With so few large format batteries reaching EOL to date, this is a new market with a variety of areas for optimization and adding value. This work with Renewance is an example of how existing expertise in battery degradation at NREL can be used to reduce the cost of shifting a battery into a second life application. With these cost reductions, this work is also facilitating the development of a battery circular economy in the United States. A robust circular economy can maximize the utilization of critical metals demanded by battery technology such as nickel and cobalt while also reducing the costs of batteries in the marketplace for the many end-uses needed for the green energy transition. The supply of these metals is limited, and we face a supply chain shortage both domestically and globally unless we can ensure they are being used to their maximum potential. This research can improve the economics of a battery circular economy to make it a more likely path for EOL batteries with critical metals. CRADA benefit to DOE, Participant, and US Taxpayer: assists laboratory in achieving programmatic scope competencies, uses the laboratory's core competencies.

25 ENERGY STORAGE↗

Structure of the Diphtheria Toxin at Acidic pH: Implications for the Conformational Switching of the Translocation Domain

Diphtheria toxin, an exotoxin secreted by Corynebacterium that causes disease in humans by inhibiting protein synthesis, enters the cell via receptor-mediated endocytosis. The subsequent endosomal acidification triggers a series of conformational changes, resulting in the refolding and membrane insertion of the translocation (T-)domain and ultimately leading to the translocation of the catalytic domain into the cytoplasm. Here, we use X-ray crystallography along with circular dichroism and fluorescence spectroscopy to gain insight into the mechanism of the early stages of pH-dependent conformational transition. For the first time, we present the high-resolution structure of the diphtheria toxin at a mildly acidic pH (5–6) and compare it to the structure at neutral pH (7). We demonstrate that neither catalytic nor receptor-binding domains change their structure upon this acidification, while the T-domain undergoes a conformational change that results in the unfolding of the TH2–3 helices. Surprisingly, the TH1 helix maintains its conformation in the crystal of the full-length toxin even at pH 5. This contrasts with the evidence from the new and previously published data, obtained by spectroscopic measurements and molecular dynamics computer simulations, which indicate the refolding of TH1 upon the acidification of the isolated T-domain. The overall results imply that the membrane interactions of the T-domain are critical in ensuring the proper conformational changes required for the preparation of the diphtheria toxin for the cellular entry.

59 BASIC BIOLOGICAL SCIENCES↗

The Escape Velocity Profile of the Milky Way from Gaia DR3

Abstract The escape velocity profile of the Milky Way offers a crucial and independent measurement of its underlying mass distribution and dark matter (DM) properties. Using a sample of stars from the third data release of Gaia with 6D kinematics and strict quality cuts, we obtain an escape velocity profile of the Milky Way from 4 to 11 kpc in Galactocentric radius. To infer the escape velocity in radial bins, we model the tail of the stellar speed distribution with both traditional power-law models and a new functional form that we introduce. While power-law models tend to rely on extrapolation to high speeds, we find our new functional form gives the most faithful representation of the observed distribution. Using this for the escape velocity profile, we constrain the properties of the Milky Way’s DM halo modeled as a Navarro–Frenck–White profile. Combined with constraints from the circular velocity at the solar position, we obtain a concentration and mass of c 200 c DM = 13.9 − 4.3 + 6.2 and M 200 c DM = 0.55 − 0.14 + 0.15 × 10 12 M ⊙ . This corresponds to a total Milky Way mass of M 200 c = 0.64 − 0.14 + 0.15 × 10 12 M ⊙ , which is on the low end of the historic range of the galaxy’s mass, but in line with other recent estimates.

Astronomy & Astrophysics↗

Spectral simulation of unsteady compressible flow past a circular cylinder

An unsteady compressible viscous wake flow past a circular cylinder was successfully simulated using spectral methods. A new approach in using the Chebyshev collocation method for periodic problems is introduced. It was further proved that the eigenvalues associated with the differentiation matrix are purely imaginary, reflecting the periodicity of the problem. It was been shown that the solution of a model problem has exponential growth in time if improper boundary conditions are used. A characteristic boundary condition, which is based on the characteristics of the Euler equations of gas dynamics, was derived for the spectral code. The primary vortex shedding frequency computed agrees well with the results in the literature for Mach = 0.4, Re = 80. No secondary frequency is observed in the power spectrum analysis of the pressure data.

Don, Wai-Sun↗

Cloud Computing Methods for Near Rectilinear Halo Orbit Trajectory Design

Complicated mission design problems require innovative computational solutions. As spacecraft depart from a proposed Gateway in a Near Rectilinear Halo Orbit (NRHO), recontact analysis is required to avoid risk of collision and ensure safe operations. Escape dynamics from NRHOs are governed by multiple gravitational bodies, yielding a trajectory design space that is exhaustively large. This paper summarizes the recontact analysis for departure from the NRHO and describes how the Deep Space Trajectory Explorer (DSTE) trajectory design software incorporates high performance cloud computing to compute and visualize the orbit design space. Recent focus on exploration missions to cislunar space has kindled accelerated interest in multibody orbit solutions. Trajectory analysis in the presence of multiple gravity fields is complex, and innovative computational tools are needed to simplify complicated design spaces, to generate large quantities of data quickly, and to visualize the output for user accessibility. The Gateway mission is a prime example. The Gateway1 is proposed as a human outpost in deep space. The current baseline orbit for the Gateway is a Near Rectilinear Halo Orbit (NRHO) near the Moon.2 The NRHO exists in a regime that experiences the gravitational effects of the Earth and the Moon simultaneously, complicating orbit analysis. The mission design process benefits greatly from updated computational tools for multibody missions like the Gateway. As an example, consider the problem of assessing the risk of collision in an NRHO. As a staging location to missions to the lunar surface and beyond the Earth-Moon system, the Gateway will experience spacecraft and other objects regularly arriving and departing. Departing objects potentially include spent logistics modules, visiting crew vehicles, debris objects, wastewater particles, and cubesats. Each departure is governed by the dynamics of the Gateway orbit and the surrounding dynamical environment. Over time, any unmaintained object in such an orbit eventually departs due to the small instabilities associated with the NRHOs. A separation maneuver speeds the departure from the NRHO, but the effects of the maneuver on the spacecraft behavior depend on the location, magnitude, and direction of the burn. Escape dynamics from the NRHO with regard to these maneuver options open up an enormous potential trajectory design space where subtle changes in input can produce dramatically large changes in the results. Any departing object must avoid recontacting the Gateway as it leaves the lunar vicinity, and a recontact analysis thus involves a significant number of computations and extensive output data. To explore the dynamics of this extensive design space, the Deep Space Trajectory Explorer3 (DSTE) trajectory design software incorporates new High Performance Computing (HPC) services and novel interactive visualizations. This paper details the HPC and cloud infrastructure techniques that are implemented in the DSTE, applying the new capabilities to analysis of recontact risk with the Gateway in NRHO. NEAR RECTILINEAR HALO ORBITS The Gateway is planned to fly in a lunar NRHO as its baseline orbit. The NRHO families of orbits are subsets of the larger halo families, which originate from planar orbits near the L1 and L2 libration points; the Earth-Moon L2 halo family appears in Figure 1. Each halo orbit is perfectly periodic in the Circular Restricted 3-Body Problem (CR3BP) and becomes a quasi-periodic orbit in a higher fidelity ephemeris force model. The NRHOs are defined as those members of the halo family with bounded stability properties;2 they pass near the Moon at perilune and are nearly polar. Families exist with apolunes located both above the lunar north pole and above the lunar south pole; the Gateway is planned to reside in a southern L2 NRHO in a 9:2 resonance with the lunar synodic period. The 9:2 NRHO is characterized by a period of about 6.5 days, a perilune radius of about 3,500 km, and an apolune radius of about 71,000 km; it is strongly affected by the gravity of both the Earth and the Moon simultaneously. This NRHO offers extended communications with assets on the south pole of the Moon,4 as well as low-cost orbit maintenance and attitude control,5 favorable eclipse avoidance properties,6 and inexpensive transfers from Earth and to other destinations.5,7 The NRHO portion of the southern L2 halo family is highlighted in black in Figure 1, and the 9:2 NRHO appears in blue.

Phillips, Sean M.↗

Comparative Habitability of Transiting Exoplanets

Exoplanet habitability is traditionally assessed by comparing a planet’s semimajor axis to the location of its host star’s “habitable zone,” the shell around a star for which Earth-like planets can possess liquid surface water. The Kepler space telescope has discovered numerous planet candidates near the habitable zone, and many more are expected from missions such as K2, TESS, and PLATO. These candidates often require significant follow-up observations for validation, so prioritizing planets for habitability from transit data has become an important aspect of the search for life in the universe. We propose a method to compare transiting planets for their potential to support life based on transit data, stellar properties and previously reported limits on planetary emitted flux. For a planet in radiative equilibrium, the emitted flux increases with eccentricity, but decreases with albedo. As these parameters are often unconstrained, there is an “eccentricity-albedo degeneracy” for the habitability of transiting exoplanets. Our method mitigates this degeneracy, includes a penalty for large-radius planets, uses terrestrial mass–radius relationships, and, when available, constraints on eccentricity to compute a number we call the “habitability index for transiting exoplanets” that represents the relative probability that an exoplanet could support liquid surface water. We calculate it for Kepler objects of interest and find that planets that receive between 60% and 90% of the Earth’s incident radiation, assuming circular orbits, are most likely to be habitable. Finally, we make predictions for the upcoming TESS and James Webb Space Telescope missions.

Planet and satellite surfaces↗

6 centimeter observations of solar bursts with 6-inch resolution

During May 1974 data were collected on nine 6-cm solar bursts originating in two active regions. One-dimensional fan-beam scans were obtained every 30 seconds in order to study burst structure and evolution. Maximum burst flux was in the 0.34-12.0 sfu range. At maximum intensity the estimated brightness temperatures were in the 10 to the 6th to 2 x 10 to the 7th K range. Angular sizes ranged from 7-23 arcsec. In most cases an expansion of the burst core after maximum intensity was observed along with a drift in the position of the burst core. Four bursts were associated important H-alpha flares, and in all cases burst maximum occurred before flare maximum. The bursts were circularly polarized in the extraordinary mode sense during maximum intensity. Since the polarization was the same over the extent of the burst source, it is suggested that if the burst is associated with loop sources, the emission is associated with one leg of the loop.

Alissandrakis, C. E.↗

Multi-waveband observations of colliding galaxies

Colliding galaxies represent a major challenge to both theorists and observers because of the large variety of phenomena which are expected to come into play during the interaction. Strong gravitational fluctuations may drive non-linear waves and instabilities throughout the stars and gas leading to enhanced star formation, nuclear activity and ultimately a mixing of the morphological components of the original galaxies. One relatively uncomplicated class of colliding galaxy where stellar waves play an important role in star formation are ring galaxies. Ring galaxies are probably formed when a companion galaxy passes through the center of a disk system driving circular waves through the disk (Lynds and Toomre 1976, Toomre 1978, Struck-Marcell 1990). Off-center collisions can generate non-circular waves and can be loosely described as banana-shaped although they may exhibit more complex forms as the waves expand into the disk. The propagation of such stellar and gaseous waves through the disk leads to enhanced star formation (e.g., Appleton and Struck-Marcell 1987a; Jeske 1986) and provides a unique probe of the response of the interstellar medium (ISM) to a propagating wave (see Appleton and Struck-Marcell 1987b). Here, the authors report results for 3 systems; the irregular ring Arp 143 (=VV 117); Wakamatsu's Seyfert ring (A0959-755; see Wakamatsu and Nishida 1987) and the brighter member of the pair of ring galaxies comprising of AM 1358-221. The most complete multi-wavelength data is for Arp 143. Optical charge coupled device (CCD) observations made with the 60 inch Palomar telescope at BV and r band, near-IR images at J (1.25 microns), H (1.65 microns) and k (2.2 microns) bands from the infrared camera (IRCAM) InSb array camera on the 3.8m United Kingdon Infrared Telescope (UKIRT) telescope and very large array (VLA) observations at 20cm in both the neutral hydrogen line and radio continuum are described. The observations of Wakamatsu's ring and AM 1358 were made only in the near-IR, and a comparison is made with available optical plate material.

Appleton, P. N.↗

Experimental Wind-Tunnel Investigation of the Transonic Damping-in-Pitch Characteristics of Two Wing-Body Combinations

The results of an experimental wind-tunnel investigation of the damping in pitch of two wing-body combinations are presented. The tests were conducted in the Ames 14-foot transonic wind tunnel over a Mach number range from 0.60 to 1.18. Reynolds numbers varied from 2.3 million to 5.5 million. One model with a triangular wing of aspect ratio 2 having NACA 0003-63 sections was oscillated at an amplitude of 1.5 and a frequency of 17 cycles per second. The second model with a straight, tapered wing of aspect ratio 3 having 3-percent biconvex circular-arc sections was oscillated at an amplitude of 1.0 deg and a frequency of 21 cycles per second. The tests were made with the models at a mean angle of attack of 0 deg. The models were oscillated with a dynamic balance that was actuated by an electrohydraulic servo valve. The results of this investigation indicate the usefulness of this new apparatus. The experimental results of a previous damping-in-pitch investigation conducted in the Ames 6- by 6-foot supersonic wind tunnel at Mach numbers from 1.2 to 1.7 are included along with the theoretical results for this Mach number range. In the region of Mach numbers available for comparison, good agreement is shown to exist between the data obtained in the two facilities, except for some inconsistency in the slopes of the curves at M = 1.2 for the triangular wing. The results of this investigation clearly show that for the models tested the maximum values of the damping in pitch occur at Mach numbers very close to 1.0, and that abrupt changes in the pitch damping are encountered near sonic velocity.

Emerson, Horace F.↗

Spitzer Secondary Eclipses of Qatar-1b

Aims: Previous secondary eclipse observations of the hot Jupiter Qatar-1B in the Ks band suggest that it may have an unusually high day-side temperature, indicative of minimal heat redistribution. There have also been indications that the orbit may be slightly eccentric, possibly forced by another planet in the system. We investigate the day-side temperature and orbital eccentricity using secondary eclipse observations with Spitzer. Methods: We observed the secondary eclipse with Spitzer/IRAC (Infrared Array Camera) in subarray mode, in both 3.6 and 4.5 micron wavelengths. We used pixel-level de-correlation to correct for Spitzer's intra-pixel sensitivity variations and thereby obtain accurate eclipse depths and central phases. Results: Our 3.6 micron eclipse depth is 0.149 plus or minus 0.051 percent and the 4.5 micron depth is 0.273 plus or minus 0.049 percent. Fitting a blackbody planet to our data and two recent Ks band eclipse depths indicates a brightness temperature of 1506 plus or minus 71 degrees Kelvin. Comparison to model atmospheres for the planet indicates that its degree of longitudinal heat redistribution is intermediate between fully uniform and day-side only. The day-side temperature of the planet is unlikely to be as high (1885 degrees Kelvin) as indicated by the ground-based eclipses in the Ks band, unless the planet's emergent spectrum deviates strongly from model atmosphere predictions. The average central phase for our Spitzer eclipses is 0.4984 plus or minus 0.0017, yielding [a model parameter] e cos omega equal to minus 0.0028 plus or minus 0.0027. Our results are consistent with a circular orbit, and we constrain e cos omega much more strongly than has been possible with previous observations.

Garhart, Emily↗