Context. Variability caused by stellar activity represents a challenge to the discovery and characterization of terrestrial exoplanets and complicates the interpretation of atmospheric planetary signals. Aims. We aim to use a detailed modeling tool to reproduce the effect of active regions on radial velocity measurements, which aids the identification of the key parameters that have an impact on the induced variability. Methods. We analyzed the effect of stellar activity on radial velocities as a function of wavelength by simulating the impact of the properties of spots, shifts induced by convective motions, and rotation. We focused our modeling effort on the active star YZ CMi (GJ 285), which was photometrically and spectroscopically monitored with CARMENES and the Telescopi Joan Or\'o. Results. We demonstrate that radial velocity curves at different wavelengths yield determinations of key properties of active regions, including spot filling factor, temperature contrast, and location, thus solving the degeneracy between them. Most notably, our model is also sensitive to convective motions. Results indicate a reduced convective shift for M dwarfs when compared to solar-type stars (in agreement with theoretical extrapolations) and points to a small global convective redshift instead of blueshift. Conclusions. Using a novel approach based on simultaneous chromatic radial velocities and light curves, we can set strong constraints on stellar activity, including an elusive parameter such as the net convective motion effect.

We report the confirmation of a transiting planet around the bright, inactive M0.5 V star TOI-1235 (TYC 4384-1735-1, V = 11.5 mag), whose transit signal was detected in the photometric time series of Sectors 14, 20, and 21 of the TESS space mission. We confirm the planetary nature of the transit signal, which has a period of 3.44 d, by using precise radial velocity measurements with CARMENES and HARPS-N spectrographs. A comparison of the properties derived for TOI-1235 b's with theoretical models reveals that the planet has a rocky composition, with a bulk density slightly higher than Earth's. In particular, we measure a mass of M_p = 5.9+/-0.6 M_Earth and a radius of R_p = 1.69+/-0.08 R_Earth, which together result in a density of rho_p = 6.7+1.3-1.1 g/cm3. When compared with other well-characterized exoplanetary systems, the particular combination of planetary radius and mass puts our discovery in the radius gap, a transition region between rocky planets and planets with significant atmospheric envelopes, with few known members. While the exact location of the radius gap for M dwarfs is still a matter of debate, our results constrain it to be located at around 1.7 R_Earth or larger at the insolation levels received by TOI-1235 b (~60 S_Earth), which makes it an extremely interesting object for further studies of planet formation and atmospheric evolution.

The high-resolution, dual channel, visible and near-infrared spectrograph CARMENES offers exciting opportunities for stellar and exoplanetary research on M dwarfs. In this work we address the challenge of reaching the highest radial velocity precision possible with a complex, actively cooled, cryogenic instrument, such as the near-infrared channel. We describe the performance of the instrument and the work flow used to derive precise Doppler measurements from the spectra. The capability of both CARMENES channels to detect small exoplanets is demonstrated with the example of the nearby M5.0 V star CD Cet (GJ 1057), around which we announce a super-Earth ($4.0\pm0.4\,M_\oplus$) companion on a 2.29 d orbit.

The He I infrared (IR) line at a vacuum wavelength of 10833 A is a diagnostic for the investigation of atmospheres of stars and planets orbiting them. For the first time, we study the behavior of the He I IR line in a set of chromospheric models for M-dwarf stars, whose much denser chromospheres may favor collisions for the level population over photoionization and recombination, which are believed to be dominant in solar-type stars. For this purpose, we use published PHOENIX models for stars of spectral types M2 V and M3 V and also compute new series of models with different levels of activity following an ansatz developed for the case of the Sun. We perform a detailed analysis of the behavior of the He I IR line within these models. We evaluate the line in relation to other chromospheric lines and also the influence of the extreme ultraviolet (EUV) radiation field. The analysis of the He I IR line strengths as a function of the respective EUV radiation field strengths suggests that the mechanism of photoionization and recombination is necessary to form the line for inactive models, while collisions start to play a role in our most active models. Moreover, the published model set, which is optimized in the ranges of the Na I D2, H$\alpha$, and the bluest Ca II IR triplet line, gives an adequate prediction of the He I IR line for most stars of the stellar sample. Because especially the most inactive stars with weak He I IR lines are fit worst by our models, it seems that our assumption of a 100% filling factor of a single inactive component no longer holds for these stars.

High resolution transit spectroscopy has proven to be a reliable technique for the characterization of the chemical composition of exoplanet atmospheres. Taking advantage of the broad spectral coverage of the CARMENES spectrograph, we initiated a survey aimed at characterizing a broad range of planetary systems. Here, we report our observations of three transits of \tplanet with CARMENES in search of \het\ absorption. On one of the nights, the He~{\sc i} region was heavily contaminated by OH$^-$ telluric emission and, thus, it was not useful for our purposes. The remaining two nights had a very different signal-to-noise ratio (S/N) due to weather. They both indicate the presence of \het\ absorption in the transmission spectrum of \tplanet, although a statistically valid detection can only be claimed for the night with higher S/N. For that night, we retrieved a 1.5$\pm$0.3\% absorption depth, translating into a $R_p(\lambda)/R_p = 1.15\pm 0.14$ at this wavelength. Spectro-photometric light curves for this same night also indicate the presence of extra absorption during the planetary transit with a consistent absorption depth. The \het\ absorption is modeled in detail using a radiative transfer code, and the results of our modeling efforts are compared to the observations. We find that the mass-loss rate, \mlr, is confined to a range of 3\,$\times\,10^{10}$\,\gs\ for $T$ = 6000\,K to 10\,$\times\,10^{10}$\,\gs\ for $T$ = 9000\,K. We discuss the physical mechanisms and implications of the He~{\sc i} detection in \tplanet and put it in context as compared to similar detections and non-detections in other Neptune-size planets. We also present improved stellar and planetary parameter determinations based on our visible and near-infrared observations.

We present the discovery and characterisation of two transiting planets observed by the Transiting Exoplanet Survey Satellite (TESS) orbiting the nearby (d ~ 22 pc), bright (J ~ 9 mag) M3.5 dwarf LP 729-54 (TOI-732). We confirm both planets and their association with LP 729-54 via ground-based photometry and determine their masses using precise radial velocities measured with the CARMENES spectrograph. Precise stellar parameters determined from CARMENES high resolution spectra confirm that LP 729-54 is a mid-M dwarf with an effective temperature of T_eff = 3360 +\- 51 K, a surface gravity of log(g) = 4.81 +/- 0.04 (cgs), and an iron abundance of [Fe/H] = 0.09 +/- 0.16 dex, with an inferred mass of M_star = 0.379 +/- 0.016 M_sun and a radius of R_star = 0.382 +/- 0.012 R_sun. The ultra-short period planet LP 729-54 b (P_b = 0.77 d) with a radius of 1.35^{+0.06}_{-0.06} R_earth, a mass of 2.34^{+0.24}_{-0.23} M_earth, and a bulk density of 5.24^{+0.94}_{-0.81} g cm^{-3} joins the population of Earth-size planets with rocky, terrestrial composition. The outer planet, LP 729-54 c, with an orbital period of 12.25 d, radius of 2.42^{+0.10}_{-0.10} R_earth, mass of 6.29^{+0.63}_{-0.61} M_earth, and mean density of 2.45^{+0.44}_{-0.37} g cm^{-3} belongs to the population of dense sub-Neptunes. With the two planets located on opposite sides of the radius gap, this planetary system is an excellent target for testing planetary formation, evolution and atmospheric models. In particular, LP 729-54 c is an ideal object for atmospheric studies with the James Webb Space Telescope.

Statistical analyses from exoplanet surveys around low-mass stars indicate that super-Earth and Neptune-mass planets are more frequent than gas giants around such stars, in agreement with core accretion theory of planet formation. Using precise radial velocities derived from visual and near-infrared spectra, we report the discovery of a giant planet with a minimum mass of 0.46 Jupiter masses in an eccentric 204-day orbit around the very low-mass star GJ 3512. Dynamical models show that the high eccentricity of the orbit is most likely explained from planet-planet interactions. The reported planetary system challenges current formation theories and puts stringent constraints on the accretion and migration rates of planet formation and evolution models, indicating that disc instability may be more efficient in forming planets than previously thought.

Context. Teegarden's Star is the brightest and one of the nearest ultra-cool dwarfs in the solar neighbourhood. For its late spectral type (M7.0V), the star shows relatively little activity and is a prime target for near-infrared radial velocity surveys such as CARMENES. Aims. As part of the CARMENES search for exoplanets around M dwarfs, we obtained more than 200 radial-velocity measurements of Teegarden's Star and analysed them for planetary signals. Methods. We find periodic variability in the radial velocities of Teegarden's Star. We also studied photometric measurements to rule out stellar brightness variations mimicking planetary signals. Results. We find evidence for two planet candidates, each with $1.1M_\oplus$ minimum mass, orbiting at periods of 4.91 and 11.4 d, respectively. No evidence for planetary transits could be found in archival and follow-up photometry. Small photometric variability is suggestive of slow rotation and old age. Conclusions. The two planets are among the lowest-mass planets discovered so far, and they are the first Earth-mass planets around an ultra-cool dwarf for which the masses have been determined using radial velocities.

We present precise photospheric parameters of 282 M dwarfs determined from fitting the most recent version of PHOENIX models to high-resolution CARMENES spectra in the visible (0.52 - 0.96 $\mu$m) and near-infrared wavelength range (0.96 - 1.71 $\mu$m). With its aim to search for habitable planets around M dwarfs, several planets of different masses have been detected. The characterization of the target sample is important for the ability to derive and constrain the physical properties of any planetary systems that are detected. As a continuation of previous work in this context, we derived the fundamental stellar parameters effective temperature, surface gravity, and metallicity of the CARMENES M-dwarf targets from PHOENIX model fits using a $\chi^2$ method. We calculated updated PHOENIX stellar atmosphere models that include a new equation of state to especially account for spectral features of low-temperature stellar atmospheres as well as new atomic and molecular line lists. We show the importance of selecting magnetically insensitive lines for fitting to avoid effects of stellar activity in the line profiles. For the first time, we directly compare stellar parameters derived from multiwavelength range spectra, simultaneously observed for the same star. In comparison with literature values we show that fundamental parameters derived from visible spectra and visible and near-infrared spectra combined are in better agreement than those derived from the same spectra in the near-infrared alone.

We determine the radii and masses of 293 nearby, bright M dwarfs of the CARMENES survey. This is the first time that such a large and homogeneous high-resolution (R>80 000) spectroscopic survey has been used to derive these fundamental stellar parameters. We derived the radii using Stefan-Boltzmann's law. We obtained the required effective temperatures $T_{\rm eff}$ from a spectral analysis and we obtained the required luminosities L from integrated broadband photometry together with the Gaia DR2 parallaxes. The mass was then determined using a mass-radius relation that we derived from eclipsing binaries known in the literature. We compared this method with three other methods: (1) We calculated the mass from the radius and the surface gravity log g, which was obtained from the same spectral analysis as $T_{\rm eff}$. (2) We used a widely used infrared mass-magnitude relation. (3) We used a Bayesian approach to infer stellar parameters from the comparison of the absolute magnitudes and colors of our targets with evolutionary models. Between spectral types M0V and M7V our radii cover the range $0.1\,R_{\normalsize\odot}<R<0.6\,R_{\normalsize\odot}$ with an error of 2-3% and our masses cover $0.09\,{\mathcal M}_{\normalsize\odot}<{\mathcal M}<0.6\,{\mathcal M}_{\normalsize\odot}$ with an error of 3-5%. We find good agreement between the masses determined with these different methods for most of our targets. Only the masses of very young objects show discrepancies. This can be well explained with the assumptions that we used for our methods.

Chromospheric modeling of observed differences in stellar activity lines is imperative to fully understand the upper atmospheres of late-type stars. We present one-dimensional parametrized chromosphere models computed with the atmosphere code PHOENIX using an underlying photosphere of 3500 K. The aim of this work is to model chromospheric lines of a sample of 50 M2-3 dwarfs observed in the framework of the CARMENES, the Calar Alto high-Resolution search for M dwarfs with Exo-earths with Near-infrared and optical Echelle Spectrographs, exoplanet survey. The spectral comparison between observed data and models is performed in the chromospheric lines of Na I D2, H$\alpha$, and the bluest Ca II infrared triplet line to obtain best-fit models for each star in the sample. We find that for inactive stars a single model with a VAL C-like temperature structure is sufficient to describe simultaneously all three lines adequately. Active stars are rather modeled by a combination of an inactive and an active model, also giving the filling factors of inactive and active regions. Moreover, the fitting of linear combinations on variable stars yields relationships between filling factors and activity states, indicating that more active phases are coupled to a larger portion of active regions on the surface of the star.

We report the detection of a Neptune-mass exoplanet around the M4.0 dwarf GJ 4276 (G 232-070) based on radial velocity (RV) observations obtained with the CARMENES spectrograph. The RV variations of GJ 4276 are best explained by the presence of a planetary companion that has a minimum mass of $m_{\rm b}\sin i \approx 16\, M_\oplus$ on a $P_{\rm b}=13.35$ day orbit. The analysis of the activity indicators and spectral diagnostics exclude stellar induced RV perturbations and prove the planetary interpretation of the RV signal. We show that a circular single-planet solution can be excluded by means of a likelihood ratio test. Instead, we find that the RV variations can be explained either by an eccentric orbit or interpreted as a pair of planets on circular orbits near a period ratio of 2:1. Although the eccentric single-planet solution is slightly preferred, our statistical analysis indicates that none of these two scenarios can be rejected with high confidence using the RV time series obtained so far. Based on the eccentric interpretation, we find that GJ 4276 b is the most eccentric ($e_{\rm b} = 0.37$) exoplanet around an M dwarf with such a short orbital period known today.

At a distance of 1.8 parsecs, Barnard's star (Gl 699) is a red dwarf with the largest apparent motion of any known stellar object. It is the closest single star to the Sun, second only to the alpha Centauri triple stellar system. Barnard's star is also among the least magnetically active red dwarfs known and has an estimated age older than our Solar System. Its properties have made it a prime target for planet searches employing techniques such as radial velocity, astrometry, and direct imaging, all with different sensitivity limits but ultimately leading to disproved or null results. Here we report that the combination of numerous measurements from high-precision radial velocity instruments reveals the presence of a low-amplitude but significant periodic signal at 233 days. Independent photometric and spectroscopic monitoring, as well as the analysis of instrumental systematic effects, show that this signal is best explained as arising from a planetary companion. The candidate planet around Barnard's star is a cold super-Earth with a minimum mass of 3.2 Earth masses orbiting near its snow-line. The combination of all radial velocity datasets spanning 20 years additionally reveals a long-term modulation that could arise from a magnetic activity cycle or from a more distant planetary object. Because of its proximity to the Sun, the proposed planet has a maximum angular separation of 220 milli-arcseconds from Barnard's star, making it an excellent target for complementary direct imaging and astrometric observations.

We announce the discovery of two planetary companions orbiting around the low mass stars Ross 1020 (GJ 3779, M4.0V) and LP 819-052 (GJ 1265, M4.5V). The discovery is based on the analysis of CARMENES radial velocity observations in the visual channel as part of its survey for exoplanets around M dwarfs. In the case of GJ 1265, CARMENES observations were complemented with publicly available Doppler measurements from HARPS. The datasets reveal one planetary companion for each star that share very similar properties: minimum masses of $8.0\pm0.5$ M$_{\oplus}$ and $7.4\pm0.5$ M$_{\oplus}$ in low-eccentricity orbits with periods of $3.023\pm0.001$ d and $3.651\pm0.001$ d for GJ 3779 b and GJ 1265 b, respectively. The periodic signals around three days found in the radial velocity data have no counterpart in any spectral activity indicator. Besides, we collected available photometric data for the two host stars, which confirm that the additional Doppler variations found at periods around 95 d can be attributed to the rotation of the stars. The addition of these planets in a mass-period diagram of known planets around M dwarfs suggests a bimodal distribution with a lack of short-period low-mass planets in the range of 2-5 M$_{\oplus}$. It also indicates that super-Earths (> 5 M$_{\oplus}$) currently detected by radial velocity and transit techniques around M stars are usually found in systems dominated by a single planet.

Context. The CARMENES spectrograph is surveying ~300 M dwarf stars in search for exoplanets. Among the target stars, spectroscopic binary systems have been discovered, which can be used to measure fundamental properties of stars. Aims. Using spectroscopic observations, we determine the orbital and physical properties of nine new double-line spectroscopic binary systems by analysing their radial velocity curves. Methods. We use two-dimensional cross-correlation techniques to derive the radial velocities of the targets, which are then employed to determine the orbital properties. Photometric data from the literature are also analysed to search for possible eclipses and to measure stellar variability, which can yield rotation periods. Results. Out of the 342 stars selected for the CARMENES survey, 9 have been found to be double-line spectroscopic binaries, with periods ranging from 1.13 to ~8000 days and orbits with eccentricities up to 0.54. We provide empirical orbital properties and minimum masses for the sample of spectroscopic binaries. Absolute masses are also estimated from mass-luminosity calibrations, ranging between ~0.1 and ~0.6 Msol . Conclusions. These new binary systems increase the number of double-line M dwarf binary systems with known orbital parameters by 15%, and they have lower mass ratios on average.

Despite their activity, low-mass stars are of particular importance for the search of exoplanets by the means of Doppler spectroscopy, as planets with lower masses become detectable. We report on the discovery of a planetary companion around HD 180617, a bright J = 5.58 mag, low-mass M = 0.45 M_{sun} star of spectral type M2.5 V. The star, located at a distance of 5.9 pc, is the primary of the high proper motion binary system containing vB 10, a star with one of the lowest masses known in most of the twentieth century. Our analysis is based on new radial velocity (RV) measurements made at red-optical wavelengths provided by the high-precision spectrograph CARMENES, which was designed to carry out a survey for Earth-like planets around M dwarfs. The available CARMENES data are augmented by archival Doppler measurements from HIRES and HARPS. Altogether, the RVs span more than 16 years. The modeling of the RV variations, with a semi-amplitude of K = 2.85-0.25/+0.16m/s yields a Neptune-like planet with a minimum mass of 12.2-1.4/+1.0 M_{Earth} on a 105.90-0.10/+0.09d circumprimary orbit, which is partly located in the host star's habitable zone. The analysis of time series of common activity indicators does not show any dependence on the detected RV signal. The discovery of HD 180617 b not only adds information to a currently hardly filled region of the mass-period diagram of exoplanets around M dwarfs, but the investigated system becomes the third known binary consisting of M dwarfs and hosting an exoplanet in an S-type configuration. Its proximity makes it an attractive candidate for future studies.

The CARMENES radial velocity (RV) survey is observing 324 M dwarfs to search for any orbiting planets. In this paper, we present the survey sample by publishing one CARMENES spectrum for each M dwarf. These spectra cover the wavelength range 520--1710nm at a resolution of at least $R > 80,000$, and we measure its RV, H$\alpha$ emission, and projected rotation velocity. We present an atlas of high-resolution M-dwarf spectra and compare the spectra to atmospheric models. To quantify the RV precision that can be achieved in low-mass stars over the CARMENES wavelength range, we analyze our empirical information on the RV precision from more than 6500 observations. We compare our high-resolution M-dwarf spectra to atmospheric models where we determine the spectroscopic RV information content, $Q$, and signal-to-noise ratio. We find that for all M-type dwarfs, the highest RV precision can be reached in the wavelength range 700--900nm. Observations at longer wavelengths are equally precise only at the very latest spectral types (M8 and M9). We demonstrate that in this spectroscopic range, the large amount of absorption features compensates for the intrinsic faintness of an M7 star. To reach an RV precision of 1ms$^{-1}$ in very low mass M dwarfs at longer wavelengths likely requires the use of a 10m class telescope. For spectral types M6 and earlier, the combination of a red visual and a near-infrared spectrograph is ideal to search for low-mass planets and to distinguish between planets and stellar variability. At a 4m class telescope, an instrument like CARMENES has the potential to push the RV precision well below the typical jitter level of 3-4ms$^{-1}$.

The new CARMENES instrument comprises two high-resolution and high-stability spectrographs that are used to search for habitable planets around M dwarfs in the visible and near-infrared regime via the Doppler technique. Characterising our target sample is important for constraining the physical properties of any planetary systems that are detected. The aim of this paper is to determine the fundamental stellar parameters of the CARMENES M-dwarf target sample from high-resolution spectra observed with CARMENES. We also include several M-dwarf spectra observed with other high-resolution spectrographs, that is CAFE, FEROS, and HRS, for completeness. We used a {chi}^2 method to derive the stellar parameters effective temperature T_eff, surface gravity log g, and metallicity [Fe/H] of the target stars by fitting the most recent version of the PHOENIX-ACES models to high-resolution spectroscopic data. These stellar atmosphere models incorporate a new equation of state to describe spectral features of low-temperature stellar atmospheres. Since T_eff, log g, and [Fe/H] show degeneracies, the surface gravity is determined independently using stellar evolutionary models. We derive the stellar parameters for a total of 300 stars. The fits achieve very good agreement between the PHOENIX models and observed spectra. We estimate that our method provides parameters with uncertainties of {sigma} T_eff = 51 K, {sigma} log g = 0.07, and {sigma} [Fe/H] = 0.16, and show that atmosphere models for low-mass stars have significantly improved in the last years. Our work also provides an independent test of the new PHOENIX-ACES models, and a comparison for other methods using low-resolution spectra. In particular, our effective temperatures agree well with literature values, while metallicities determined with our method exhibit a larger spread when compared to literature results.

Context: The main goal of the CARMENES survey is to find Earth-mass planets around nearby M-dwarf stars. Seven M-dwarfs included in the CARMENES sample had been observed before with HIRES and HARPS and either were reported to have one short period planetary companion (GJ15A, GJ176, GJ436, GJ536 and GJ1148) or are multiple planetary systems (GJ581 and GJ876). Aims: We aim to report new precise optical radial velocity measurements for these planet hosts and test the overall capabilities of CARMENES. Methods: We combined our CARMENES precise Doppler measurements with those available from HIRES and HARPS and derived new orbital parameters for the systems. Bona-fide single planet systems are fitted with a Keplerian model. The multiple planet systems were analyzed using a self-consistent dynamical model and their best fit orbits were tested for long-term stability. Results: We confirm or provide supportive arguments for planets around all the investigated stars except for GJ15A, for which we find that the post-discovery HIRES data and our CARMENES data do not show a signal at 11.4 days. Although we cannot confirm the super-Earth planet GJ15Ab, we show evidence for a possible long-period ($P_{\rm c}$ = 7025$_{-629}^{+972}$ d) Saturn-mass ($m_{\rm c} \sin i$ = 51.8$_{-5.8}^{+5.5}M_\oplus$) planet around GJ15A. In addition, based on our CARMENES and HIRES data we discover a second planet around GJ1148, for which we estimate a period $P_{\rm c}$ = 532.6$_{-2.5}^{+4.1}$ d, eccentricity $e_{\rm c}$ = 0.34$_{-0.06}^{+0.05}$ and minimum mass $m_{\rm c} \sin i$ = 68.1$_{-2.2}^{+4.9}M_\oplus$. Conclusions: The CARMENES optical radial velocities have similar precision and overall scatter when compared to the Doppler measurements conducted with HARPS and HIRES. We conclude that CARMENES is an instrument that is up to the challenge of discovering rocky planets around low-mass stars.

We report on the first star discovered to host a planet detected by radial velocity (RV) observations obtained within the CARMENES survey for exoplanets around M dwarfs. HD 147379 ($V = 8.9$ mag, $M = 0.58 \pm 0.08$ M$_{\odot}$), a bright M0.0V star at a distance of 10.7 pc, is found to undergo periodic RV variations with a semi-amplitude of $K = 5.1\pm0.4$ m s$^{-1}$ and a period of $P = 86.54\pm0.06$ d. The RV signal is found in our CARMENES data, which were taken between 2016 and 2017, and is supported by HIRES/Keck observations that were obtained since 2000. The RV variations are interpreted as resulting from a planet of minimum mass $m_{\rm p}\sin{i} = 25 \pm 2$ M$_{\oplus}$, 1.5 times the mass of Neptune, with an orbital semi-major axis $a = 0.32$ au and low eccentricity ($e < 0.13$). HD 147379b is orbiting inside the temperate zone around the star, where water could exist in liquid form. The RV time-series and various spectroscopic indicators show additional hints of variations at an approximate period of 21.1d (and its first harmonic), which we attribute to the rotation period of the star.

We present a brief overview of a splinter session on determining the metallicity of low-mass dwarfs that was organized as part of the Cool Stars 16 conference. We review contemporary spectroscopic and photometric techniques for estimating metallicity in low-mass dwarfs and discuss the importance of measuring accurate metallicities for studies of Galactic and chemical evolution using subdwarfs, creating metallicity benchmarks for brown dwarfs, and searching for extrasolar planets that are orbiting around low--mass dwarfs. In addition, we present the current understanding of the effects of metallicity on stellar evolution and atmosphere models and discuss some of the limitations that are important to consider when comparing theoretical models to data.

The core of scientific research is turning new ideas into reality. From the school science fair to the search for the secrets of dark energy, high-quality research consists of scientific investigation constrained within the scope of a well-defined project. Large or small, generously funded or just scraping by,scientific projects use time, money, and information to turn ideas into plans, plans into action, and action into results. While we, as a community, do much to educate students in the techniques of research, we do not systematically train students in the nature and organization of scientific projects or in the techniques of project management. We propose a two-pronged attack to address this issue in the next decade. First, to generate a broad base of future scientists who have a basic familiarity with the ideas of projects, we propose that the community develop standards for the content of a project design and management course in astronomy and astrophysics. Second, to train future scientists to assume leadership roles in new investigations in astronomy and astrophysics, we propose that the community develop standards for graduate programs in the area of research project leadership.

High resolution synthetic stellar libraries are of fundamental importance for the preparation of the Gaia Mission. We present new sets of spectral stellar libraries covering two spectral ranges: 300 --1100 nm at 0.1 nm resolution, and 840 -- 890 nm at 0.001 nm resolution. These libraries span a large range in atmospheric parameters, from super-metal-rich to very metal-poor (-5.0 $<$[Fe/H]$<$+1.0), from cool to hot (\teff=3000--50000 K) stars, including peculiar abundance variations. The spectral resolution, spectral type coverage and number of models represent a substantial improvement over previous libraries used in population synthesis models and in atmospheric analysis.

Two peculiar stars showing an apparent extremely broadened and strong NaI D absorption have been discovered in surveys for cool white dwarfs by Oppenheimer et al. (2001) and Harris et al. (SDSS, 2003). We discuss the nature of these objects using PHOENIX atmosphere models for metal-poor brown dwarfs/very low mass stars, and new white dwarf LTE and NLTE models for hydrogen- and helium-dominated atmospheres with metals. These include complete molecular formation in chemical equilibrium and a model for the alkali resonance line broadening based on the damping profiles of Allard et al. (2003), as well as new molecular line opacities for metal hydrides. First results of our calculations indicate good agreement with a hydrogen-dominated WD atmosphere with a Na abundance roughly consistent with a state of high accretion. We analyse deviations of the abundances of Na, K, Mg and Ca from the cosmic pattern and comment on implications of these results for standard accretion scenarios.

We present a complete line list for the X 1Sigma+ electronic ground state of LiCl computed using fully quantum-mechanical techniques. This list includes transition energies and oscillator strengths in the spectral region 0.3-39,640.7 cm-1 for all allowed rovibrational transitions in absorption within the electronic ground state. The calculations were performed using an accurate hybrid potential constructed from a spectral inversion fit of experimental data and from recent multi-reference single- and double-excitation configuration interaction calculations. The line list was incorporated into the stellar atmosphere code PHOENIX to compute spectra for a range of young to old T dwarf models. The possibility of observing a signature of LiCl in absorption near 15.8 microns is addressed and the proposal to use this feature to estimate the total lithium elemental abundance for these cool objects is discussed.

Brown dwarfs bridge the gap between the stellar and planetary mass regimes. Evolving from conditions very similar to the lowest-mass stars, the atmospheres of older brown dwarfs closely resemble those expected in close-in extrasolar giant planets, and with cooler BDs still being discovered, more and more approach the properties of gas giants at wider separation. Interpreting the spectra of BDs is therefore a crucial step towards understanding and predicting the spectral and thermal properties of EGPs. Essential properties of substellar atmospheres are massive molecular line-blanketing and the condensation of species with decreasing Teff, changing the chemical equilibrium composition and causing absorption from dust grains. More complex details involve the distribution of dust clouds over the surface giving rise to temporal variability, and possible deviations from chemical equilibrium conditions. In the case of close-in EGPs and some BDs in binary systems, the effect of irradiation from the primary significantly affects the spectral properties and thermal evolution.

We briefly describe the current version of the PHOENIX code. We then present some illustrative results from the modeling of Type Ia and Type II supernovae, hot stars, and irradiated giant planets. Good fits to observations can be obtained, when account is taken for spherically symmetric, line-blanketed, static or expanding atmospheres.

The opacity due to photodissociation of 24MgH is investigated in the atmospheres of cool stars. The lowest two electronic transitions A 2Pi -- X 2Sigma+ and B' 2Sigma+ -- X 2 Sigma+ are considered where the cross sections for the latter were published previously (Weck, Stancil, & Kirby 2002) while the former are presented in this work. Model atmospheres calculated with the PHOENIX code are used to investigate the effect of the photodissociation opacity on spectra of cool stars. The A -- X photodissociation cross sections are obtained using a combination of ab initio and experimentally derived potential curves and dipole transition moments. Partial cross sections have been evaluated over the accessible wavelength range 1770-4560 Angstrom for all rotational transitions from the vibrational levels v''=0-11. Assuming a Boltzmann distribution of the rovibrational levels of the X 2Sigma+ state, LTE photodissociation cross sections are presented for temperatures between 1000 and 5000 K. Shape resonances, arising from rotational predissociation of quasi-bound levels of the A 2Pi state near threshold, characterize the LTE photodissociation cross sections. A sum rule is proposed as a check on the accuracy of the photodissociation calculations.

A new, complete, theoretical rotational and vibrational line list for the A-X electronic transition in MgH is presented. The list includes transition energies and oscillator strengths for all possible allowed transitions and was computed using the best available theoretical potential energies and dipole transition moment function with the former adjusted to account for experimental data. The A-X line list, as well as new line lists for the B'-X and the X-X (pure rovibrational) transitions, were included in comprehensive stellar atmosphere models for M, L, and T dwarfs and solar-type stars. The resulting spectra, when compared to models lacking MgH, show that MgH provides significant opacity in the visible between 4400 and 5600 Angstrom. Further, comparison of the spectra obtained with the current line list to spectra obtained using the line list constructed by Kurucz (1993) show that the Kurucz list significantly overestimates the opacity due to MgH particularly for the bands near 5150 and 4800 Angstrom with the discrepancy increasing with decreasing effective temperature.

The recent announcement that sodium absorption has been observed in the atmosphere of HD209458b, the only EGP observed to transit its parent star, is the first direct detection of an EGP atmosphere. We explore the possibility that neutral sodium is {\em not} in local thermodynamic equilibrium (LTE) in the outer atmosphere of irradiated EGPs and that the sodium concentration may be underestimated by models that make the LTE assumption. Our results indicate that it may not be necessary to invoke excessive photoionization, low metallicity, or even high altitude clouds to explain the observations.

We present a spectral analysis of a sample of late L dwarfs. We use our latest model atmospheres and synthetic spectra and optical and K band spectra to determine effective temperatures. We derive effective temperatures of 1400 K to 1700 K for L8 to L6 dwarfs. The analysis demonstrates that our recent models that rain out the formed dust completely are applicable to optical spectra of late L dwarfs and that more consistent models are needed for intermediate L dwarfs and for infrared spectra. We compare the results for the effective temperatures with the temperatures of the onset of Methane formation. Our models predict Methane absorption at 3.3micron to occur at about 400 K higher temperatures than Methane absorption at 2.2micron. This is consistent with our data and previous observations which show Methane absorption at 3.3micron but not at 2.2micron in late L dwarfs.

We have studied the stellar population and internal structure of the core of the Ursa Minor dwarf spheroidal galaxy, using a sample of stars selected to be members based on their proper motions. In agreement with previous studies, we find Ursa Minor to be dominated by an old, metal-poor stellar population. A small number of stars with high membership probabilities lie redward of the red giant branch. The brightest (V <= 18) such stars are known to be Carbon stars, rather than metal-rich first-ascent giants. A number of stars with high membership probabilities lie blueward of the red giant branch, and are more luminous than the horizontal branch. We speculate that these are post-horizontal branch stars. There may also be one or two stars in the post-AGB phase. Spectroscopy of the candidate post-HB and post-AGB stars is required to determine their nature. We recover the internal substructure in Ursa Minor that has been noted by several authors in the last 15 years. Using a variety of two- and three-dimensional statistical tests, we conclude that this substructure is statistically significant at the 0.005 level. There is no evidence that the regions of density excess have stellar populations that differ from the main body of Ursa Minor. The crossing time for a typical density excess is only ~5 million years. They are therefore clearly not due to intermediate age star-forming bursts. We conclude that they are instead due to tidal interactions between the Galaxy and Ursa Minor.

We present opacity sampling model atmospheres, synthetic spectra and colors for brown dwarfs and very low mass stars in two limiting case of dust grain formation: 1) inefficient gravitational settling i.e. the dust is distributed according to the chemical equilibrium predictions, 2) efficient gravitational settling i.e. the dust forms and depletes refractory elements from the gas, but their opacity does not affect the thermal structure. The models include the formation of over 600 gas phase species, and 1000 liquids and crystals, and the opacities of 30 different types of grains including corundum (Al$_2$O$_3$), the magnesium aluminum spinel MgAl$_2$O$_4$, iron, enstatite (MgSiO$_3$), forsterite (Mg$_2$SiO$_4$), amorphous carbon, SiC, and a number of calcium silicates. The models extend from the beginning of the grain formation regime well into the condensation regime of water ice ($\teff= 3000 - 100$ K) and encompasses the range of $\log g= 2.5 - 6.0$ at solar metallicity. We find that silicate dust grains can form abundantly in the outer atmospheric layers of red and brown dwarfs with spectral type later than M8. The greenhouse effects of dust opacities provide a natural explanation for the peculiarly red spectroscopic distribution of the latest M dwarfs and young brown dwarfs. The grainless (Cond) models on the other hand, correspond closely to methane brown dwarfs such as Gliese 229B. We also recover that the $\lambda$5891,5897\AA Na I D and $\lambda$7687,7701\AA K I resonance doublets plays a critical role in T dwarfs where their red wing define the pseudo-continuum from the $I$ to the $Z$ bandpass.

We present analyses of high resolution and medium resolution spectra of early L dwarfs. We used our latest set of model atmospheres to reproduce and analyze the observed features. We can model the optical flux and the atomic line profiles with the best accuracy to date. The models used to reproduce the observations include dust condensation and dust opacities. Compared to previous studies using older models we find that our dust treatment is much improved. The derived parameters for the objects are well in the expected range for old very low mass objects. This is also supported by the absence of Li in most of the objects. For the objects showing Li we can be almost certain that those are brown dwarfs. However, a spectral analysis in general, and this one in particular can only very roughly determine mass and age.

We have obtained 1.0-2.5um spectra at R~600 of 14 disk dwarfs with spectral types M6 to L7. For four of the dwarfs we have also obtained infrared spectra at R~3000 in narrow intervals. In addition, we present new L' photometry for four of the dwarfs in the sample, which allows improved determinations of their bolometric luminosities. We resolve the L-dwarf Denis-P J 0205-1159 into an identical pair of objects separated by 0.35". The spectra, with the published energy distribution for one other dwarf, are compared to synthetic spectra generated by upgraded model atmospheres. Good matches are found for 2200> Teff K>1900 (spectral types around M9 to L3), but discrepancies exist at Teff> 2300 K (M8) and for Teff<1800 K (L4-L7). At the higher temperatures the mismatches are due to incompleteness in the water vapor linelist. At the lower temperatures the disagreement is probably due to our treatment of dust: we assume a photospheric distribution in equilibrium with the gas phase. We derive effective temperatures for the sample from the comparison with synthetic spectra and also by comparing our observed total intrinsic luminosities to structural model calculations (which are mostly independent of the atmosphere but are dependent on the unknown masses and ages of the targets). The two derivations agree to ~200 K except for the faintest object in the sample where the discrepancy is larger. Agreement with other temperature determinations is also ~200 K, except for the L7 dwarf.

We present a grid of spherically symmetric model atmospheres for young pre-MS stars. This grid spans the parameter range $2000\K \leq \Teff \leq 6800\K$ and $2.0 \leq \logg \leq 3.5$ for $M=0.1\Msun$, appropriate for low mass stars and brown dwarfs. A major improvement is the replacement of TiO and \water line lists with the newer line list calculated by the NASA-AMES group, for TiO (about 175 million lines of 5 isotopes) and for \water (about 350 million lines in 2 isotopes). We provide the model structures, spectra and broad-band colors in standard filters in electronic form.

We present a technique to treat systems with very many levels, like molecules, in non-LTE. This method is based on a superlevel formalism coupled with rate operator splitting. Superlevels consist of many individual levels that are assumed to be in LTE relative to each other. The usage of superlevels reduces the dimensionality of the rate equations dramatically and, thereby, makes the problem computationally more easily treatable. Our superlevel formalism retains maximum accuracy by using direct opacity sampling (dOS) when calculating the radiative transitions and the opacities. We developed this method in order to treat molecules in cool dwarf model calculations in non-LTE. Cool dwarfs have low electron densities and a radiation field that is far from a black body radiation field, both properties may invalidate the conditions for the common LTE approximation. Therefore, the most important opacity sources, the molecules, need to be treated in non-LTE. As a case study we applied our method to carbon monoxide. We find that our method gives accurate results since the conditions for the superlevel method are very well met for molecules. Due to very high collisional cross sections with hydrogen, and the high densities of H_2 the population of CO itself shows no significant deviation from LTE.

We present high-spatial resolution HST and adaptive optics observations, and high-sensitivity ISO (ISOCAM & ISOPHOT) observations of a sample of X-ray selected weak-line (WTTS) and post (PTTS) T Tauri stars located in the nearby Chamaeleon T and Scorpius-Centaurus OB associations. HST/NICMOS and adaptive optics observations aimed at identifying substellar companions at separations >=30 A.U. from the primary stars. No such objects were found within 300 A.U. of any of the target stars, and a number of faint objects at larger separations can very likely be attributed to a population of background stars. ISOCAM observations of 5 to 15 Myr old WTTS and PTTS in ScoCen reveal infrared excesses which are clearly above photospheric levels, and which have a spectral index intermediate between that of younger (1 to 5 Myr) T Tauri stars in Cha and that of pure stellar photospheres. The difference in the spectral index of the older PTTS in ScoCen compared to the younger classical and WTTS in Cha can be attributed to a deficiency of smaller size (0.1 to 1mu) dust grains relative to larger size (~5mu) dust grains in the disks of the PTTS. The lack of small dust grains is either due to the environment (effect of nearby O stars and supernova explosions) or due to disk evolution. If the latter is the case, it would hint that circumstellar disks start to get dust depleted at an age between 5 to 15 Myr. Dust depletion is very likely related to the build-up of larger particles (ultimately rocks and planetesimals) and thus an indicator for the onset of the period of planet formation.

We use a preliminary version of our ``NextGen'' grid of cool star model atmospheres to compute synthetic line profiles which fit high resolution Keck spectra of the cool M~dwarf VB10 satisfactorily well. We show that the parameters derived from the Keck data are consistent with the parameters derived from lower resolution spectra with larger wavelength coverage. We discuss the treatment of van der Waals broadening in cool, molecular (mostly ${\rm H_2}$) dominated stellar atmospheres. The line profiles are dominated by van der Waals pressure broadening and are a sensitive indicator for the gravity and metallicity. Therefore, the high-resolution Keck spectra are useful for determining the parameters of M dwarfs. There is some ambiguity between the metallicity and gravity. For VB10, we find from the high-resolution spectra that $ 5.0 < \log(g) < 5.5$ and $0 < \left[\case{{\rm M}}{{\rm H}}\right] < +0.5$ for an adopted fixed effective temperature of 2700~K (Schweitzer, 1995), which is consistent with recent interior calculations (Baraffe etal, 1995).

The atmospheres of M stars are dominated by a small number of very strong molecular compounds (H$_2$O, TiO, H$_2$, CO, VO). Most of the hydrogen is locked in molecular H$_2$, most of the carbon in CO; and H$_2$O, TiO and VO opacities define a pseudo-continuum covering the entire flux distribution of these stars. The optical ``continuum'' is due to TiO vibrational bands which are often used as temperature indicators for these stars. These may be the depth of the bands relative to the troughs in between them; or the depth of the VO bands; or of the atomic lines relative to the local ``continuum''; or even the strength of the infrared water bands; all of these depend on the strength of the TiO bands and the amount of flux-redistribution to longer wavelengths exerted by them. Departures from LTE of the Ti I atom, and thus the concentration of the important TiO molecule, could, therefore, have severe and measurable consequences on the atmospheric structure and spectra of these stars. In this paper we discuss NLTE effects of Ti I in fully self-consistent models for a few representative M/Brown dwarf and M giant model atmospheres and spectra.