National Solar Observatory

Full list of authors: Price-Whelan, Adrian M.; Lim, Pey Lian; Earl, Nicholas; Starkman, Nathaniel; Bradley, Larry; Shupe, David L.; Patil, Aarya A.; Corrales, Lia; Brasseur, C. E.; Noethe, Maximilian; Donath, Axel; Tollerud, Erik; Morris, Brett M.; Ginsburg, Adam; Vaher, Eero; Weaver, Benjamin A.; Tocknell, James; Jamieson, William; van Kerkwijk, Marten H.; Robitaille, Thomas P.; Merry, Bruce; Bachetti, Matteo; Gunther, H. Moritz; Aldcroft, Thomas L.; Alvarado-Montes, Jaime A.; Archibald, Anne M.; Bodi, Attila; Bapat, Shreyas; Barentsen, Geert; Bazan, Juanjo; Biswas, Manish; Boquien, Mederic; Burke, D. J.; Cara, Daria; Cara, Mihai; Conroy, Kyle E.; Conseil, Simon; Craig, Matthew W.; Cross, Robert M.; Cruz, Kelle L.; D'Eugenio, Francesco; Dencheva, Nadia; Devillepoix, Hadrien A. R.; Dietrich, Jorg P.; Eigenbrot, Arthur Davis; Erben, Thomas; Ferreira, Leonardo; Foreman-Mackey, Daniel; Fox, Ryan; Freij, Nabil; Garg, Suyog; Geda, Robel; Glattly, Lauren; Gondhalekar, Yash; Gordon, Karl D.; Grant, David; Greenfield, Perry; Groener, Austen M.; Guest, Steve; Gurovich, Sebastian; Handberg, Rasmus; Hart, Akeem; Hatfield-Dodds, Zac; Homeier, Derek; Hosseinzadeh, Griffin; Jenness, Tim; Jones, Craig K.; Joseph, Prajwel; Kalmbach, J. Bryce; Karamehmetoglu, Emir; Kaluszynski, Mikolaj; Kelley, Michael S. P.; Kern, Nicholas; Kerzendorf, Wolfgang E.; Koch, Eric W.; Kulumani, Shankar; Lee, Antony; Ly, Chun; Ma, Zhiyuan; MacBride, Conor; Maljaars, Jakob M.; Muna, Demitri; Murphy, N. A.; Norman, Henrik; O'Steen, Richard; Oman, Kyle A.; Pacifici, Camilla; Pascual, Sergio; Pascual-Granado, J.; Patil, Rohit R.; Perren, Gabriel, I; Pickering, Timothy E.; Rastogi, Tanuj; Roulston, Benjamin R.; Ryan, Daniel F.; Rykoff, Eli S.; Sabater, Jose; Sakurikar, Parikshit; Salgado, Jesus; Sanghi, Aniket; Saunders, Nicholas; Savchenko, Volodymyr; Schwardt, Ludwig; Seifert-Eckert, Michael; Shih, Albert Y.; Jain, Anany Shrey; Shukla, Gyanendra; Sick, Jonathan; Simpson, Chris; Singanamalla, Sudheesh; Singer, Leo P.; Singhal, Jaladh; Sinha, Manodeep; Sipocz, Brigitta M.; Spitler, Lee R.; Stansby, David; Streicher, Ole; Sumak, Jani; Swinbank, John D.; Taranu, Dan S.; Tewary, Nikita; Tremblay, Grant R.; De Val-Borro, Miguel; Vasovic, Zlatan; Van Kooten, Samuel J.; Verma, Shresth; Cardoso, Jose Vinicius de Miranda; Williams, Peter K. G.; Wilson, Tom J.; Winkel, Benjamin; Wood-Vasey, W. M.; Xue, Rui; Yoachim, Peter; Zhang, Chen; Zonca, Andrea; Astropy Project Contributors; TARDIS Collaboration; Astropy Coordination Comm.--This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

The Interface Region Imaging Spectrograph (IRIS) small explorer spacecraft provides simultaneous spectra and images of the photosphere, chromosphere, transition region, and corona with 0.33 – 0.4 arcsec spatial resolution, two-second temporal resolution, and 1 km s−1 velocity resolution over a field-of-view of up to 175 arcsec × 175 arcsec. IRIS was launched into a Sun-synchronous orbit on 27 June 2013 using a Pegasus-XL rocket and consists of a 19-cm UV telescope that feeds a slit-based dual-bandpass imaging spectrograph. IRIS obtains spectra in passbands from 1332 – 1358 Å, 1389 – 1407 Å, and 2783 – 2834 Å, including bright spectral lines formed in the chromosphere (Mg ii h 2803 Å and Mg ii k 2796 Å) and transition region (C ii 1334/1335 Å and Si iv 1394/1403 Å). Slit-jaw images in four different passbands (C ii 1330, Si iv 1400, Mg ii k 2796, and Mg ii wing 2830 Å) can be taken simultaneously with spectral rasters that sample regions up to 130 arcsec × 175 arcsec at a variety of spatial samplings (from 0.33 arcsec and up). IRIS is sensitive to emission from plasma at temperatures between 5000 K and 10 MK and will advance our understanding of the flow of mass and energy through an interface region, formed by the chromosphere and transition region, between the photosphere and corona. This highly structured and dynamic region not only acts as the conduit of all mass and energy feeding into the corona and solar wind, it also requires an order of magnitude more energy to heat than the corona and solar wind combined. The IRIS investigation includes a strong numerical modeling component based on advanced radiative–MHD codes to facilitate interpretation of observations of this complex region. Approximately eight Gbytes of data (after compression) are acquired by IRIS each day and made available for unrestricted use within a few days of the observation.

Data from the Global Oscillation Network Group (GONG) project and other helioseismic experiments provide a test for models of stellar interiors and for the thermodynamic and radiative properties, on which the models depend, of matter under the extreme conditions found in the sun. Current models are in agreement with the helioseismic inferences, which suggests, for example, that the disagreement between the predicted and observed fluxes of neutrinos from the sun is not caused by errors in the models. However, the GONG data reveal subtle errors in the models, such as an excess in sound speed just beneath the convection zone. These discrepancies indicate effects that have so far not been correctly accounted for; for example, it is plausible that the sound-speed differences reflect weak mixing in stellar interiors, of potential importance to the overall evolution of stars and ultimately to estimates of the age of the galaxy based on stellar evolution calculations.

The splitting of the frequencies of the global resonant acoustic modes of the Sun by large-scale Ñows and rotation permits study of the variation of angular velocity) with both radius and latitude within the turbulent convection zone and the deeper radiative interior. The nearly uninterrupted Doppler imaging observations, provided by the Solar Oscillations Investigation (SOI) using the Michelson Doppler Imager (MDI) on the Solar and Heliospheric Observatory (SOHO) spacecraft positioned at the Lagrangian point in continuous sunlight, yield oscillation power spectra with very high signal-to-L 1noise ratios that allow frequency splittings to be determined with exceptional accuracy. This paper reports on joint helioseismic analyses of solar rotation in the convection zone and in the outer part of the radiative core. Inversions have been obtained for a medium-l mode set (involving modes of angular degree l extending to about 250) obtained from the Ðrst 144 day interval of SOI-MDI observations in 1996. Drawing inferences about the solar internal rotation from the splitting data is a subtle process. By applying more than one inversion technique to the data, we get some indication of what are the more robust and less robust features of our inversion solutions. Here we have used seven di†erent inversion methods. To test the reliability and sensitivity of these methods, we have performed a set of controlled experiments utilizing artiÐcial data. This gives us some conÐdence in the inferences we can draw from

Abstract We present cosmological results from the measurement of baryon acoustic oscillations (BAO) in galaxy, quasar and Lyman- α forest tracers from the first year of observations from the Dark Energy Spectroscopic Instrument (DESI), to be released in the DESI Data Release 1. DESI BAO provide robust measurements of the transverse comoving distance and Hubble rate, or their combination, relative to the sound horizon, in seven redshift bins from over 6 million extragalactic objects in the redshift range 0.1 < z < 4.2. To mitigate confirmation bias, a blind analysis was implemented to measure the BAO scales. DESI BAO data alone are consistent with the standard flat ΛCDM cosmological model with a matter density Ω m =0.295±0.015. Paired with a baryon density prior from Big Bang Nucleosynthesis and the robustly measured acoustic angular scale from the cosmic microwave background (CMB), DESI requires H 0 =(68.52±0.62) km s -1 Mpc -1 . In conjunction with CMB anisotropies from Planck and CMB lensing data from Planck and ACT, we find Ω m =0.307± 0.005 and H 0 =(67.97±0.38) km s -1 Mpc -1 . Extending the baseline model with a constant dark energy equation of state parameter w , DESI BAO alone require w =-0.99 +0.15 -0.13 . In models with a time-varying dark energy equation of state parametrised by w 0 and w a , combinations of DESI with CMB or with type Ia supernovae (SN Ia) individually prefer w 0 > -1 and w a < 0. This preference is 2.6 σ for the DESI+CMB combination, and persists or grows when SN Ia are added in, giving results discrepant with the ΛCDM model at the 2.5 σ , 3.5 σ or 3.9 σ levels for the addition of the Pantheon+, Union3, or DES-SN5YR supernova datasets respectively. For the flat ΛCDM model with the sum of neutrino mass ∑ m ν free, combining the DESI and CMB data yields an upper limit ∑ m ν < 0.072 (0.113) eV at 95% confidence for a ∑ m ν > 0 (∑ m ν > 0.059) eV prior. These neutrino-mass constraints are substantially relaxed if the background dynamics are allowed to deviate from flat ΛCDM.

Alfvén waves, transverse incompressible magnetic oscillations, have been proposed as a possible mechanism to heat the Sun's corona to millions of degrees by transporting convective energy from the photosphere into the diffuse corona. We report the detection of Alfvén waves in intensity, line-of-sight velocity, and linear polarization images of the solar corona taken using the FeXIII 1074.7-nanometer coronal emission line with the Coronal Multi-Channel Polarimeter (CoMP) instrument at the National Solar Observatory, New Mexico. Ubiquitous upward propagating waves were seen, with phase speeds of 1 to 4 megameters per second and trajectories consistent with the direction of the magnetic field inferred from the linear polarization measurements. An estimate of the energy carried by the waves that we spatially resolved indicates that they are too weak to heat the solar corona; however, unresolved Alfvén waves may carry sufficient energy.

Helioseismology requires nearly continuous observations of the oscillations of the solar surface for long periods of time in order to obtain precise measurements of the sun's normal modes of oscillation. The GONG project acquires velocity images from a network of six identical instruments distributed around the world. The GONG network began full operation in October 1995. It has achieved a duty cycle of 89 percent and reduced the magnitude of spectral artifacts by a factor of 280 in power, compared with single-site observations. The instrumental noise is less than the observed solar background.

▪ Abstract Radio emission from solar flares offers a number of unique diagnostic tools to address long-standing questions about energy release, plasma heating, particle acceleration, and particle transport in magnetized plasmas. At millimeter and centimeter wavelengths, incoherent gyrosynchrotron emission from electrons with energies of tens of kilo electron volts to several mega electron volts plays a dominant role. These electrons carry a significant fraction of the energy released during the impulsive phase of flares. At decimeter and meter wavelengths, coherent plasma radiation can play a dominant role. Particularly important are type III and type III–like radio bursts, which are due to upward- and downward-directed beams of nonthermal electrons, presumed to originate in the energy release site. With the launch of Yohkoh and the Compton Gamma-Ray Observatory, the relationship between radio emission and energetic photon emissions has been clarified. In this review, recent progress on our Dunderstanding of radio emission from impulsive flares and its relation to X-ray emission is discussed, as well as energy release in flare-like phenomena (microflares, nanoflares) and their bearing on coronal heating.

Using a sample of 69,919 red giants from the SDSS-III/APOGEE Data Release 12, we measure the distribution of stars in the [/Fe] versus [Fe/H] plane and the metallicity distribution functions (MDFs) across an unprecedented volume of the Milky Way disk, with radius 3 < R < 15 kpc and height z 2 < | | kpc. Stars in the inner disk (R < 5 kpc) lie along a single track in [/Fe] versus [Fe/H], starting with -enhanced, metal-poor stars and ending at [/Fe] 0 and [Fe/H] +0.4. At larger radii we find two distinct sequences in [/Fe] versus

We measure the large-scale cross-correlation of quasars with the Lyα forest absorption, using over 164,000 quasars from Data Release 11 of the SDSS-III Baryon Oscillation Spectroscopic Survey. We extend the previous study of roughly 60,000 quasars from Data Release 9 to larger separations, allowing a measurement of the Baryonic Acoustic Oscillation (BAO) scale along the line of sight c /( H ( z = 2.36) r s ) = 9.0±0.3 and across the line of sight D A ( z = 2.36)/ r s = 10.8±0.4, consistent with CMB and other BAO data. Using the best fit value of the sound horizon from Planck data ( r s = 147.49 Mpc), we can translate these results to a measurement of the Hubble parameter of H ( z = 2.36) = 226±8 km s −1 Mpc −1 and of the angular diameter distance of D A ( z = 2.36) = 1590±60 Mpc. The measured cross-correlation function and an update of the code to fit the BAO scale (baofit) are made publicly available.

A multilevel accelerated lambda iteration (MALI) method for radiative transfer calculations with partial frequency redistribution (PRD) is presented. The method, which is based on Rybicki & Hummer's complete frequency redistribution (CRD) formalism with full preconditioning, consistently accounts for overlapping radiative transitions. Its extension to PRD is implemented in a very natural way through the use of the Ψ operator operating on the emissivity rather than the commonly used Λ operator, which operates on the source function. Apart from requiring an additional inner computational loop to evaluate the PRD emission-line profiles with fixed population numbers, implementation of the presented method requires only a trivial addition of computer code. Since the presented method employs a diagonal operator, it is easily extended to different geometries. Currently, it has been implemented for one-, two-, and three-dimensional Cartesian grids and spherical symmetry. In all cases, the speed of convergence with PRD is very similar to that in CRD, with the former sometimes even surpassing the latter. Sample calculations exhibiting the favorable convergence behavior of the PRD code are presented in the case of the Ca II H and K lines, the Mg II h and k lines, and the hydrogen Lyα and Lyβ lines in a one-dimensional solar model and the Ca II resonance lines in a two-dimensional flux-sheet model.

We have detected changes in the rotation of the sun near the base of its convective envelope, including a prominent variation with a period of 1.3 years at low latitudes. Such helioseismic probing of the deep solar interior has been enabled by nearly continuous observation of its oscillation modes with two complementary experiments. Inversion of the global-mode frequency splittings reveals that the largest temporal changes in the angular velocity Omega are of the order of 6 nanohertz and occur above and below the tachocline that separates the sun's differentially rotating convection zone (outer 30% by radius) from the nearly uniformly rotating deeper radiative interior beneath. Such changes are most pronounced near the equator and at high latitudes and are a substantial fraction of the average 30-nanohertz difference in Omega with radius across the tachocline at the equator. The results indicate variations of rotation close to the presumed site of the solar dynamo, which may generate the 22-year cycles of magnetic activity.

Abstract The goal of the SunPy project is to facilitate and promote the use and development of community-led, free, and open source data analysis software for solar physics based on the scientific Python environment. The project achieves this goal by developing and maintaining the sunpy core package and supporting an ecosystem of affiliated packages. This paper describes the first official stable release (version 1.0) of the core package, as well as the project organization and infrastructure. This paper concludes with a discussion of the future of the SunPy project.

We present the first results from a speckle imaging survey of stars classified as candidate exoplanet host stars discovered by the Kepler mission. We use speckle imaging to search for faint companions or closely aligned background stars that could contribute flux to the Kepler light curves of their brighter neighbors. Background stars are expected to contribute significantly to the pool of false positive candidate transiting exoplanets discovered by the Kepler mission, especially in the case that the faint neighbors are eclipsing binary stars. Here, we describe our Kepler follow-up observing program, the speckle imaging camera used, our data reduction, and astrometric and photometric performance. Kepler stars range from R = 8 to 16 and our observations attempt to provide background non-detection limits 5-6 mag fainter and binary separations of ~0.05-2.0 arcsec. We present data describing the relative brightness, separation, and position angles for secondary sources, as well as relative plate limits for non-detection of faint nearby stars around each of 156 target stars. Faint neighbors were found near 10 of the stars.

We measure the power spectrum, P_F(k,z), of the transmitted flux in the Ly-alpha forest using 3035 high redshift quasar spectra from the Sloan Digital Sky Survey. This sample is almost two orders of magnitude larger than any previously available data set, yielding statistical errors of ~0.6% and ~0.005 on, respectively, the overall amplitude and logarithmic slope of P_F(k,z). This unprecedented statistical power requires a correspondingly careful analysis of the data and of possible systematic contaminations in it. For this purpose we reanalyze the raw spectra to make use of information not preserved by the standard pipeline. We investigate the details of the noise in the data, resolution of the spectrograph, sky subtraction, quasar continuum, and metal absorption. We find that background sources such as metals contribute significantly to the total power and have to be subtracted properly. We also find clear evidence for SiIII correlations with the Ly-alpha forest and suggest a simple model to account for this contribution to the power. While it is likely that our newly developed analysis technique does not eliminate all systematic errors in the P_F(k,z) measurement below the level of the statistical errors, our tests indicate that any residual systematics in the analysis are unlikely to affect the inference of cosmological parameters from P_F(k,z). These results should provide an essential ingredient for all future attempts to constrain modeling of structure formation, cosmological parameters, and theories for the origin of primordial fluctuations.

The SDSS-III/Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey operated from 2011-2014 using the APOGEE spectrograph, which collects high-resolution (R 22,500), near-IR (1.51-1.70 m) spectra with a multiplexing (300 fiber-fed objects) capability. We describe the survey data products that are publicly available, which include catalogs with radial velocity, stellar parameters, and 15 elemental abundances for over 150,000 stars, as well as the more than 500,000 spectra from which these quantities are derived. Calibration relations for the stellar parameters (T eff , g log , [M/H], [/M]) and abundances (C, N, O, Na, Mg, Al, Si, S, K, Ca, Ti, V, Mn, Fe, Ni) are presented and discussed. The internal scatter of the abundances within clusters indicates that abundance precision is generally between 0.05 and 0.09 dex across a broad temperature range; it is smaller for some elemental abundances within more limited ranges and at high signal-to-noise ratio. We assess the accuracy of the abundances using comparison of mean cluster metallicities with literature values, APOGEE observations of the solar spectrum and of Arcturus, comparison of individual star abundances with other measurements, and consideration of the locus of derived parameters and abundances of the entire sample, and find that it is challenging to determine the absolute abundance scale; external accuracy may be good to 0.1-0.2 dex. Uncertainties may be larger at cooler temperatures (T eff 4000 K < ). Access to the public data release and data products is described, and some guidance for using the data products is provided.

Splitting of the sun's global oscillation frequencies by large-scale flows can be used to investigate how rotation varies with radius and latitude within the solar interior. The nearly uninterrupted observations by the Global Oscillation Network Group (GONG) yield oscillation power spectra with high duty cycles and high signal-to-noise ratios. Frequency splittings derived from GONG observations confirm that the variation of rotation rate with latitude seen at the surface carries through much of the convection zone, at the base of which is an adjustment layer leading to latitudinally independent rotation at greater depths. A distinctive shear layer just below the surface is discernible at low to mid-latitudes.

We analyze the SDSS Ly-alpha forest P_F(k,z) measurement to determine the linear theory power spectrum. Our analysis is based on fully hydrodynamic simulations, extended using hydro-PM simulations. We account for the effect of absorbers with damping wings, which leads to an increase in the slope of the linear power spectrum. We break the degeneracy between the mean level of absorption and the linear power spectrum without significant use of external constraints. We infer linear theory power spectrum amplitude Delta^2_L(k_p=0.009s/km,z_p=3.0)=0.452_{-0.057-0.116}^{+0.069+0.141} and slope n_eff=-2.321_{-0.047-0.102}^{+0.055+0.131} (possible systematic errors are included through nuisance parameters in the fit - a factor >~5 smaller errors would be obtained on both parameters if we ignored modeling uncertainties). The errors are correlated and not perfectly Gaussian, so we provide a chi^2 table to accurately describe the results. The result corresponds to sigma_8=0.85, n=0.94, for a LCDM model with Omega_m=0.3, Omega_b=0.04, and h=0.7, but is most useful in a combined fit with the CMB. The inferred curvature of the linear power spectrum and the evolution of its amplitude and slope with redshift are consistent with expectations for LCDM models, with the evolution of the slope, in particular, being tightly constrained. We use this information to constrain systematic contamination, e.g., fluctuations in the UV background. This paper should serve as a starting point for more work to refine the analysis, including technical improvements such as increasing the size and number of the hydrodynamic simulations, and improvements in the treatment of the various forms of feedback from galaxies and quasars.

Our knowledge of the long-term evolution of solar activity and of its primary modulation, the 11-year cycle, largely depends on a single direct observational record: the visual sunspot counts that retrace the last 4 centuries, since the invention of the astronomical telescope. Currently, this activity index is available in two main forms: the International Sunspot Number initiated by R. Wolf in 1849 and the Group Number constructed more recently by Hoyt and Schatten (Sol. Phys. 179:189–219, 1998a, 181:491–512, 1998b). Unfortunately, those two series do not match by various aspects, inducing confusions and contradictions when used in crucial contemporary studies of the solar dynamo or of the solar forcing on the Earth climate. Recently, new efforts have been undertaken to diagnose and correct flaws and biases affecting both sunspot series, in the framework of a series of dedicated Sunspot Number Workshops. Here, we present a global overview of our current understanding of the sunspot number calibration.

Abstract Transmission spectra are differential measurements that utilize stellar illumination to probe transiting exoplanet atmospheres. Any spectral difference between the illuminating light source and the disk-integrated stellar spectrum due to starspots and faculae will be imprinted in the observed transmission spectrum. However, few constraints exist for the extent of photospheric heterogeneities in M dwarfs. Here we model spot and faculae covering fractions consistent with observed photometric variabilities for M dwarfs and the associated 0.3–5.5 μ m stellar contamination spectra. We find that large ranges of spot and faculae covering fractions are consistent with observations and corrections assuming a linear relation between variability amplitude, and covering fractions generally underestimate the stellar contamination. Using realistic estimates for spot and faculae covering fractions, we find that stellar contamination can be more than 10× larger than the transit depth changes expected for atmospheric features in rocky exoplanets. We also find that stellar spectral contamination can lead to systematic errors in radius and therefore the derived density of small planets. In the case of the TRAPPIST-1 system, we show that TRAPPIST-1's rotational variability is consistent with spot covering fractions and faculae covering fractions . The associated stellar contamination signals alter the transit depths of the TRAPPIST-1 planets at wavelengths of interest for planetary atmospheric species by roughly 1–15× the strength of planetary features, significantly complicating JWST follow-up observations of this system. Similarly, we find that stellar contamination can lead to underestimates of the bulk densities of the TRAPPIST-1 planets of , thus leading to overestimates of their volatile contents.