Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrostaticity”

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 775 records · Page 43

Inside the supernova: A powerful convective engine

We present an extensive study of the inception of supernova explosions by following the evolution of the cores of two massive stars (15 and 25 Solar mass) in multidimension. Our calculations begin at the onset of core collapse and stop several hundred milliseconds after the bounce, at which time successful explosions of the appropriate magnitude have been obtained. Similar to the classical delayed explosion mechanism of Wilson, the explosion is powered by the heating of the envelope due to neutrinos emitted by the protoneutron star as it radiates the gravitational energy liberated by the collapse. However, as was shown by Herant, Benz, & Colgate, this heating generates strong convection outside the neutrinosphere, which we demonstrate to be critical to the explosion. By breaking a purely stratified hydrostatic equilibrium, convection moves the nascent supernova away from a delicate radiative equilibrium between neutrino emission and absorption, Thus, unlike what has been observed in one-dimensional calculations, explosions are rendered quite insensitive to the details of the physical input parameters such as neutrino cross sections or nuclear equation of state parameters. As a confirmation, our comparative one-dimensional calculations with identical microphysics, but in which convection cannot occur, lead to dramatic failures. Guided by our numerical results, we have developed a paradigm for the supernova explosion mechanism. We view a supernova as an open cycle thermodynamic engine in which a reservoir of low-entropy matter (the envelope) is thermally coupled and physically connected to a hot bath (the protoneutron star) by a neutrino flux, and by hydrodynamic instabilities. This paradigm does not invoke new or modified physics over previous treatments, but relies on compellingly straightforward thermodynamic arguments. It provides a robust and self-regulated explosion mechanism to power supernovae that is effective under a wide range of physical parameters.

Herant, Marc↗

Loop models of low coronal structures observed by the Normal Incidence X-Ray Telescope (NIXT)

The X-ray pictures obtained with the Normal Incidence X-Ray Telescope (NIXT), apart from the ubiquitous coronal loops well known from previous X-ray observations, show a new and peculiar morphology: in many active regions there are wide and apparently low-lying areas of intense emission which resemble H alpha plages. By means of hydrostatic models of coronal arches, we analyze the distribution of temperature, density, emission measure, and plasma emissivity in the spectral band to which NIXT is sensitive, and we show that the above morphology can be explained by the characteristics of high pressure loops having a thin region of high surface brightness at the base. We therefore propose that this finding might help to identify high-pressure X-ray emitting coronal regions in NIXT images, and it is in principle applicable to any imaging instrument which has high sensitivity to 10(exp 4) - 10(exp 6) K plasma within a narrow coronal-temperature passband. As a more general result of this study, we propose that the comparison of NIXT observations with models of stationary loops might provide a new diagnostic: the determination of the loop plasma pressure from measurements of brightness distribution along the loop.

Peres, G.↗

Dynamics of emerging active region flux loops

The buoyant rise of a magnetic flux loop arising from a single perturbed segment of a toroidal flux ring lying slightly beneath the base of the convection zone is studied by way of numerical simulations. We have considered flux loop evolution assuming both solid-body rotation, and differential rotation consistent with recent results from helioseismology. Our major results are presented, and we offer some speculations on the decay of active regions, based on the results of our studies. We speculate that as plasma in the tube attempts to establish hydrostatic equilibrium along the field lines after the flux emergence has taken place, the tube field strength at some intermediate depths below the surface becomes sufficiently small at the surface portions of the tube (which have cooled and undergone convective collapse) become dynamically disconnected from those portions near the base of the convection zone. The surface proportions of the emerged flux tubes are then transported by motions near the photosphere, such as supergranular convection and meridional flow.

Fan, Y.↗

The r-process and neutrino-heated supernova ejecta

As a neutron star is formed by the collapse of the iron core of a massive star, its Kelvin-Helmholtz evolution is characterized by the release of gravitational binding energy as neutrinos. The interaction of these neutrinos with heated material above the neutron star generates a hot bubble in an atmosphere that is nearly in hydrostatic equilibrium and heated, after approximately 10 s, to an entropy of S/N(sub AS)k greater than or approximately = 400. The neutron-to-proton ratio for material moving outward through this bubble is set by the balance between neutrino and antineutrino capture on nucleons. Because the electron antineutrino spectrum at this time is hotter than the electron neutrino spectrum, the bubble is neutron-rich (0.38 less than or approximately = Y(sub e) less than or approximately = 0.47). Previous work using a schematic model has shown that these conditions are well suited to the production of heavy elements by the r-process. In this paper we have advanced the numerical modeling of a 20 solar mass 'delayed' supernova explosion to the point that we can follow the detailed evolution of material moving through the bubble at the late times appropiate to r-process nucleosynthesis. The supernova model predicts a final kinetic energy for the ejecta of 1.5 x 10(exp 51) ergs and leaves behind a remnant with a baryon mass of 1.50 solar mass (and a gravitational mass of 1.445 solar mass). We follow the thermodynamic and compositional evolution of 40 trajectories in rho(t), T(t), Y(sub e)(t) for a logarithmic grid of mass elements for the last approximately = 0.03 solar mass to be ejected by the proto-neutron star down to the last less than 10(exp -6) solar mass of material expelled at up to approximately = 18 s after core collapse. We find that an excellent fit to the solar r-process abundance distribution is obtained with no adjustable parameters in the nucleosynthesis calculations. Moreover, the abundances are produced in the quantities required to account for the present Galactic abundances. However, at earlier times, this one-dimensional model ejects too much material with entropies S/N(sub A)k approximately 50 and Y(sub e) approximately 0.46. This leads to an acceptable over production of N = 50 nuclei, particularly Sr-88, Y-89, and Zr-90, relative to their solar abundances. We speculate on various means to avoid the early overproduction and/or ejection of N = 50 isotonic nuclei while still producing and ejecting the correct amount of r-process material.

Woosley, S. E.↗

Super-Eddington radiation transfer in soft gamma repeaters

Bursts from soft gamma repeaters (SGRs) have been shown to be super-Eddington by a factor of 1000 and have been persuasively associated with compact objects. Super-Eddington radiation transfer on the surface of a strongly magnetic (greater than or equal to 10(exp 13) G) neutron star is studied and related to the observational constraints on SGRs. In strong magnetic fields, Thompson scattering is suppressed in one polarization state, so super-Eddington fluxes can be radiated while the plasma remains in hydrostatic equilibrium. We discuss a model which offers a somewhat natural explanation for the observation that the energy spectra of bursts with varying intensity are similar. The radiation produced is found to be linearly polarized to one part in 1000 in a direction determined by the local magnetic field, and intensity variations between bursts are understood as a change in the radiating area on the source. The net polarization is inversely correlated with burst intensity. Further, it is shown that for radiation transfer calculations in limit of superstrong magnetic fields, it is sufficient to solve the radiation transfer for the low opacity state rather than the coupled equations for both. With this approximation, standard stellar atmosphere techniques are utilized to calculate the model energy spectrum.

Ulmer, Andrew↗

Disruption of a helmet streamer by photospheric shear

Helmet streamers on the Sun have been observed to be the site of coronal mass ejections, dynamic events that eject coronal plasma and magnetic fields into the solar wind. We develop a two-dimensional (azimuthally symmetric) helmet streamer configuration by computing solutions of the time-dependent magnetohydrodynamic (MHD) equations, and we investigate the evolution of the configuration when photospheric shearing motions are imposed. We find that the configuration disrupts when a critical shear is exceeded, ejecting a plasmoid into the solar wind. The results are similar to the case of a sheared dipole magnetic field in a hydrostatic atmosphere (Mikic & Linker 1994). However, the presence of the outflowing solar wind makes the disruption significantly more energetic when a helmet streamer is sheared. Our resutls suggest that shearing of helmet streamers may initiate coronal mass ejections.

Linker, Jon A.↗

HST eclipse mapping of dwarf nova OY Carinae in quiescence: An 'Fe II curtain' with Mach approx. = 6 velocity dispersion veils the white dwarf

Hubble Space Telescope (HST) observations of the eclipsing dwarf nova OY Car in its quiescent state are used to isolate the ultraviolet spectrum (1150-2500 A at 9.2 A Full Width at Half Maximum (FWHM) resolution) of the white dwarf, the accretion disk, and the bright spot. The white dwarf spectrum has a Stark-broadened photospheric L(alpha) absorption, but is veiled by a forest of blended Fe II features that we attribute to absorption by intervening disk material. A fit gives T(sub w) approx. = 16.5 x 10(exp 3) K for the white dwarf with a solar-abundance, log g = 8 model atmosphere, and T approx. = 10(exp 4) K, n(sub e) approx. = 10(exp 13)/cu cm, N(sub H) approx. = 10(exp 22) sq cm, and velocity dispersion delta V approx. = 60 km/s for the veil of homogeneous solar-abundance local thermodynamic equilibrium (LTE) gas. The veil parameters probably measure characteristic physical conditions in the quiescent accretion disk or its chromosphere. The large velocity dispersion is essential for a good fit; it lowers (chi square)/778 from 22 to 4. Keplerian shear can produce the velocity dispersion if the veiling gas is located at R approx. = 5 R(sub W) with (delta R)/R approx. = 0.3, but this model leaves an unobscured view to the upper hemisphere of the white dwarf, incompatible with absorptions that are up to 80% deep. The veiling gas may be in the upper atmosphere of the disk near its outer rim, but we then require supersonic (Mach approx. = 6) but sub-Keplerian (delta V/V(sub Kep) approx. = 0.07) velocity disturbances in this region to produce both the observed radial velocity dispersion and vertical motions sufficient to elevate the gas to z/R = cos i = 0.12. Such motions might be driven by the gas stream, since it may take several Kepler periods to reestablish the disk's vertical hydrostatic equilibrium. The temperature and column density of the gas we see as Fe II absorption in the ultraviolet are similar to what is required to produce the strong Balmer jump and line emissions seen in optical spectra of OY Car and similar quiescent dwarf novae. The outer accretion disk is detected at mid-eclipse with a spectrum that rises from 0.05 to 0.3 mJy between 2000 and 2500 A, consistent with combinations of cool blackbodies, blended Fe II emission lines, and Balmer continuum emission. The total disk flux density is 0.5 mJy at 2500 A, and this shallow disk eclipse implies a roughly flat surface brightness distribution. The bright spot, somewhat bluer than the disk, has a flux density rising from 0.05 to 0.15 mJy between 1600 and 2500 A. The C IV emission line has a broad shallow eclipse, but the radial velocity variations observed during the eclipse do not clearly distinguish between a disk or wind origin. The only possible indications of boundary layer emission are fast UV flares that appear to arise from near the central object -- not from the bright spot.

Horne, Keith↗

X-ray-emitting gas surrounding the spiral galaxy NGC 891

We observed the edge-on spiral galaxy NGC 891 with the Position Sensitive Proportional Counter (PSPC) on Roentgen Satellite (ROSAT) to search for how extraplanar gas expected in the galactic fountain model. Diffuse X-ray emission surrounds the disk with a Half Width at Half Maximum (HWHM) for the surface brightness perpendicular to the disk of 50 sec (2.4 kpc) and a radial extent of approximately 6.5 kpc, both of which are similar in extent to the extended H(alpha) and radio halo component; the implied density scale height for the hot gas is 7 kpc. The spectrum is best fitted with a hard stellar component and a soft diffuse gas component of temperature 3.6 x 10(exp 6) K. The density of this gas is 2 x 10(exp -3)/cu cm, the luminosity is 4.4 x 10(exp 39) ergs/s, the mass is 1 x 10(exp 8) solar mass, and the pressure (P/k) is 1.4 10(exp 4) K/cu cm. These data are consistent with this gas participating in a galactic fountain, where the material approaches hydrostatic equilibrium before cooling at a rate of 0.12 solar mass/yr. The cooled material may be responsible for some of the H(alpha) emission.

Bregman, Joel N.↗

Diabatic forcing and intialization with assimilation of cloud water and rainwater in a forecast model

In this study, diabatic forcing, and liquid water assimilation techniques are tested in a semi-implicit hydrostatic regional forecast model containing explicit representations of grid-scale cloud water and rainwater. Diabatic forcing, in conjunction with diabatic contributions in the initialization, is found to help the forecast retain the diabatic signal found in the liquid water or heating rate data, consequently reducing the spinup time associated with grid-scale precipitation processes. Both observational Special Sensor Microwave/Imager (SSM/I) and model-generated data are used. A physical retrieval method incorporating SSM/I radiance data is utilized to estimate the 3D distribution of precipitating storms. In the retrieval method the relationship between precipitation distributions and upwelling microwave radiances is parameterized, based upon cloud ensemble-radiative model simulations. Regression formulae relating vertically integrated liquid and ice-phase precipitation amounts to latent heating rates are also derived from the cloud ensemble simulations. Thus, retrieved SSM/I precipitation structures can be used in conjunction with the regression-formulas to infer the 3D distribution of latent heating rates. These heating rates are used directly in the forecast model to help initiate Tropical Storm Emily (21 September 1987). The 14-h forecast of Emily's development yields atmospheric precipitation water contents that compare favorably with coincident SSM/I estimates.

Raymond, William H.↗

A simulation of the intracluster medium with feedback from cluster galaxies

We detail method and report first results from a three-dimensional hydrodynamical and N-body simulation of the formation and evolution of a Coma-sized cluster of galaxies, with the intent of studying the history of the hot, X-ray emitting intracluster medium. Cluster gas, galaxies, and dark matter are included in the model. The galaxies and dark matter fell gravitational forces; the cluster gas also undergoes hydrodynamical effects such as shock heating and PdV work. For the first time in three dimensions, we include modeling of ejection of processed gas from the simulated galaxies by winds, including heating and heavy element enrichment. For comparison, we employ a `pure infall' simulation using the same initial conditions but with no galaxies or winds. We employ an extreme ejection history for galactic feedback in order to define the boundary of likely models. As expected, feedback raises the entropy of the intracluster gas, preventing it from collapsing to densities as high as those attained in the infall model. The effect is more pronounced in subclusters formed at high redshift. The cluster with feedback is always less X-ray luminous, but experiences more rapid luminosity evolution, than the pure infall cluster. Even employing an extreme ejection model, the final gas temperature is only approximately 15% larger than in the infall model. The radial temperature profile is very nearly isothermal within 1.5 Mpc. The cluster galaxies in the feedback model have a velocity dispersion approximately 15% lower than the dark matter. This results in the true ratio of specific energies in galaxies to gas being less than one, beta(sub spec) approximately 0.7. The infall model predicts beta(sub spec) approximately 1.2. Large excursions in these values occur over time, following the complex dynamical history of the cluster. The morphology of the X-ray emission is little affected by feedback. The emission profiles of both clusters are well described by the standard beta-model with beta(sub fit) approximately equal to 0.7 - 0.9. X-ray mass estimates based on the assumptions of hydrostatic equilibrium and the applicability of the beta-model are quite accurate in both cases. A strong, radial iron abundance gradient is present, which develops as a consequence of the steepening of the galaxy density profile over time. Spectroscopic observations using nonimaging detectors with wide (approximately 45 min) fields of view dramatically smear the gradient. Observations with arcminute resolution, made available with the ASCA satellite, would readily resolve the gradient.

Metzler, Christopher A.↗

Interactions between gravity waves and cold air outflows in a stably stratified uniform flow

Interactions between gravity waves and cold air outflows in a stably stratified uniform flow forced by various combinations of prescribed heat sinks and sources are studied using a hydrostatic two-dimensional nonlinear numerical model. The formation time for the development of a stagnation point or reversed flow at the surface is not always directly proportional to the Froude number when wave reflections exist from upper levels. A density current is able to form by the wave-otuflow interaction, even though the Froude number is greater than a critical value. This is the result of the wave-outflow interaction shifting the flow response to a different location in the characteristic parameter space. A density current is able to form or be destroyed due to the wave-outflow interaction between a traveling gravity wave and cold air outflow. This is proved by performing experiments with a steady-state heat sink and an additional transient heat source. In a quiescent fluid, a region of cold air, convergence, and upward motion is formed after the collision between two outflows produced by two prescribed heat sinks. After the collision, the individual cold air outflows lose their own identity and merge into a single, stationary, cold air outflow region. Gravity waves tend to suppress this new stationary cold air outflow after the collision. The region of upward motion associated with the collision is confined to a very shallow layer. In a moving airstream, a density current produced by a heat sink may be suppressed or enhanced nonlinearly by an adjacent heat sink due to the wave-outflow interaction.

Lin, Yuh-Lang↗

ROSAT PSPC observations of the early-type galaxies NGC 507 and NGC 499: Central cooling and mass determination

We present the results of a deep observation of NGC 507 and NGC 499 with the ROSAT Position Sensitive Proportional Counter (PSPC). The X-ray emission of NGC 507 is extended at least out to 1000 sec (458 kpc at a distance of 94.5 Mpc). The radial profile of X-ray surface brightness goes as Sigma(sub x) is approximately r(exp -1.8) outside the core region. The radial profile is a function of energy such that the softer X-rays have a smaller core radius and a flatter slope. Spectral analysis reveals that the emission temperature, with an average of 1 keV, peaks at an intermediate radius of 2-3 min and falls toward the center (possibly decreases outward as well). The absorption column density is consistent with the Galactic line-of-sight value. The X-ray emission of NGC 499 is extended to 300 sec and suggests a similarly cooler core. The cooler cores of NGC 507 and NGC 499 are strong evidence of the presence of cooling flows in these galaxies. Assuming hydrostatic equilibrium outside the cooling radius, the estimated mass-to-light ratio of NGC 507 is 97 +/- 16 within 458 kpc, indicative of the presence of a heavy halo. Similarly, the mass-to-light ratio of NGC 499 is 89 +/- 14 within 137 kpc. Near the edge of the X-ray-emitting region of NGC 507 we detect 19 soft, unresolved sources. These sources do not have optical counterparts and are significantly in excess of the expected number of background serendipitous sources. We speculate that they may represent cooling clumps in the halo of NGC 507. If there are many undetected cooling clumps distributed at large radii, then the radial profile of the X-ray surface brightness does not directly reflect the potential, adding uncertainty to the measurement of the binding mass; the gas mass could also be overestimated.

Kim, Dong-Woo↗

An X-ray temperature map of Abell 754: A major merger

We present the first two-dimensional X-ray temperature map of the rich cluster of galaxies Abell 754. We also present an X-ray surface brightness map with improved spatial resolution and sensitivity compared with previous maps. Both the temperature map and the surface brightness map show that A754 is in the throes of a violent merger; it is probably far from hydrostatic equilibrium.

Henry, J. Patrick↗

The formation of protostellar disks. 2: Disks around intermediate-mass stars

Hydrodynamical calculations of the evolution of a collapsing, rotating axisymmetric 10 solar masses molecular clump, including the effects of radiative acceleration but without magnetic fields, are represented. The initial cloud is assumed to be uniformly rotating, centrally condensed sphere with rho is proportional to r(exp -2). Several cases are considered, in which both the overall clump size and the total amount of angular momentum are varied. The calculations show how a warm, quasi-hydrostatic disk surrounding a central unresolved core of only a few solar masses forms and grows in size and mass. The disk is encased in two distinct accretion shock fronts, both of which are several scale heights above the equatorial plane. At the end of the calculation of our standard case, the central unresolved region is found to have a mass of 2.7 solar masses and a ratio of rotational to gravitational energy of approximately 0.45, sufficiently large to be unstable to nonaxisymmetric perturbations. In addition, the inner portions of the disk containing most of the mass are unstable according to the local Toomre criterion, implying that also in this region nonaxisymmetric perturbations will lead to rapid evolution. Under the assumption that gravitational torques would transport angular momentum out of this region, a central core of less than or approximately 8 solar masses with a stable disk of greater than or approximately = 2 solar masses should result. Frequency-dependent radiative transfer calculations of the standard case at selected ages show how the continuum spectrum of the structure depends on the disk's orientation and age and how the observed isophotal contours vary with wavelength. Because of the strong dependence on viewing angle, continuum spectra alone should not be used to estimate the evolutionary stage of development of these objects. Comparable results were obtained for the other cases considered.

Yorke, Harold W.↗

The transition regions of Capella

We have used the Goddard High Resolution Spectrometer (GHRS) to observe the spectoscopic binary system Capella (G8 III + G1 III). Exposures with the G140L, G140M, G160M, G200M, and echelle gratings provide emission line profiles with unprecedented signal-to-noise and spectral resolving power (lambda/Delta-lambda) up to 92,000. Multi-Gaussin fits to the line profiles show that the hotter star contributes 60%-70% of the total flux in the chromospheric O I and Mg II resonance lines, but about 90% of the flux in the Si III, Si IV, and C IV lines formed in the transition region at T less than or = 10(exp 5) K. We find clear evidence that the emission lines from the hotter star are systemtically redshifted relative to the photosphere with Doppler shifts of 5 +/- 1 km/s for the +9 +/- 3 km/s in the chromospheric Mg II and O I lines, respectively, increasing to +24 +/- 5 km/s for the transition region Si IV 1393.8A line. The multi-Gaussian fits to permitted transition region lines of SI III, Si IV, C IV, and N V indicate the presence of three components: moderately broad lines formed in the transition region of the hotter star (component H), narrow lines formed in the transition region of the cooler star (component C), and very broad lines that we think are formed in microflares on the hotter star (component B). The He II 1640.4 A feature has an broad profile, which indicates that it is formed by collisional excitation primarily from the hotter star, and a weak narrow component that we interpret as due to radiative recombination on the cooler star. We observed spin-forbidden emission lines of C III), O III), Si III), O IV), O V), and S IV) that are sensitive to electron density. Fainter members of the O IV) multiplet and all of the S IV) lines have never before been seen in any star than the Sun. We determine electron densities in the transition regions of the Capella stars using lines ratios of O IV) lines and emission measure analysis. The emission measures are self-consistent only when the fluxes from each emitting component are considered separately. In particular, the transition region abundance distributions appear to be different on the two stars, and the spin-forbidden lines were not detected in the B component. The emission measures for component H (the G1 star) are consistent with a constant electron pressure (or hydrostatic equilibrium) transition region with P(sub e) = 10(exp 15 + 0.1)/cu cm K and possibly solar 'coronal' abundances.

Linsky, Jeffrey L.↗

Submillimeter heterodyne spectroscopy and remote sensing of the upper atmosphere

Submillimeter-wavelength heterodyne spectroscopy by microwave limb sounding (MLS) from orbiting satellite is a measurement technique for studying Earth's stratosphere, mesosphere and lowerr thermosphere on a global scale. Development and deployment of such tools are timely, as stratospheric ozone shileds life from solar ultraviolet radiation but is depleted by pollution from industrial activities. MLS experiments on NASA's Upper Atmosphere Research Satellite (UARS) and Earth Observing System (EOS) are now being developed and implementend for global monitoring. Atmospheric thermal-emission spectra at millimeter and submillimeter wave-lengths are measured as the instrument field of view (FOV) is scanned through the limb from above. Atmospheric profiles of molecular abundances, temperature, pressure, wind and magnetic field can be determined from the measured emission spectra. Intensity of the emission can provide abundance and temperature. Measured linewidths, and emission from temperature-insensitive O2 lines, can provide pressure. Differentitaion of measured pressure with repsect ot measured height differential (obtaiend from the instrument FOV scan encoder) can also provide temperature through atmospheric hydrostatic equilibrium (which relates temperature to pressure and height differential). Doppler shifts of spectral lines can provide wind, and Zeeman splitting of the magnetic-dipole lines of O2 can provide magnetic field.

Waters, J. W.↗

Gravitational coefficients and internal structures of the icy Galilean satellites: An assessment of the Galileo Orbiter mission

Simulated radio tracking data from one flyby of Europa, two flybys of Ganymede, and two flybys of Callisto by the Galileo Orbiter yield estimates of the standard errors in the gravitational coefficient J(sub 2) of 44, 8.5, and 12 (in units of 10(exp -6)) for Europa, Ganymede, and Callisto, respectively; errors in C(sub 22) (in units of 10(exp -6)) are 12, 1.7, and 2.1. These errors are sufficiently small that the values of J(sub 2) and C(sub 22) to be measured by Galileo should suffice to determine if the ice and rock in the satellite interiors are uniformly mixed or separated, so long as the bodies are in hydrostatic equilibrium.

Schubert, G.↗

Model atmospheres for M (sub)dwarf stars. 1: The base model grid

We have calculated a grid of more than 700 model atmospheres valid for a wide range of parameters encompassing the coolest known M dwarfs, M subdwarfs, and brown dwarf candidates: 1500 less than or equal to T(sub eff) less than or equal to 4000 K, 3.5 less than or equal to log g less than or equal to 5.5, and -4.0 less than or equal to (M/H) less than or equal to +0.5. Our equation of state includes 105 molecules and up to 27 ionization stages of 39 elements. In the calculations of the base grid of model atmospheres presented here, we include over 300 molecular bands of four molecules (TiO, VO, CaH, FeH) in the JOLA approximation, the water opacity of Ludwig (1971), collision-induced opacities, b-f and f-f atomic processes, as well as about 2 million spectral lines selected from a list with more than 42 million atomic and 24 million molecular (H2, CH, NH, OH, MgH, SiH, C2, CN, CO, SiO) lines. High-resolution synthetic spectra are obtained using an opacity sampling method. The model atmospheres and spectra are calculated with the generalized stellar atmosphere code PHOENIX, assuming LTE, plane-parallel geometry, energy (radiative plus convective) conservation, and hydrostatic equilibrium. The model spectra give close agreement with observations of M dwarfs across a wide spectral range from the blue to the near-IR, with one notable exception: the fit to the water bands. We discuss several practical applications of our model grid, e.g., broadband colors derived from the synthetic spectra. In light of current efforts to identify genuine brown dwarfs, we also show how low-resolution spectra of cool dwarfs vary with surface gravity, and how the high-regulation line profile of the Li I resonance doublet depends on the Li abundance.

Allard, France↗