sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Frazin, R. A.

Normalized to: Frazin, R.

20 article(s) in total. 87 co-authors, from 1 to 5 common article(s). Median position in authors list is 1,0.

[1]  [pdf] - 1881118
Directly Imaging Rocky Planets from the Ground
Comments: 8 pages, 1 figure, Astro2020 Science White Paper
Submitted: 2019-05-10
Over the past three decades instruments on the ground and in space have discovered thousands of planets outside the solar system. These observations have given rise to an astonishingly detailed picture of the demographics of short-period planets, but are incomplete at longer periods where both the sensitivity of transit surveys and radial velocity signals plummet. Even more glaring is that the spectra of planets discovered with these indirect methods are either inaccessible (radial velocity detections) or only available for a small subclass of transiting planets with thick, clear atmospheres. Direct detection can be used to discover and characterize the atmospheres of planets at intermediate and wide separations, including non-transiting exoplanets. Today, a small number of exoplanets have been directly imaged, but they represent only a rare class of young, self-luminous super-Jovian-mass objects orbiting tens to hundreds of AU from their host stars. Atmospheric characterization of planets in the <5 AU regime, where radial velocity (RV) surveys have revealed an abundance of other worlds, is technically feasible with 30-m class apertures in combination with an advanced AO system, coronagraph, and suite of spectrometers and imagers. There is a vast range of unexplored science accessible through astrometry, photometry, and spectroscopy of rocky planets, ice giants, and gas giants. In this whitepaper we will focus on one of the most ambitious science goals --- detecting for the first time habitable-zone rocky (<1.6 R_Earth) exoplanets in reflected light around nearby M-dwarfs
[2]  [pdf] - 1784067
Treating Wavefront Measurement Error in Estimation of Non-Common Path Aberration for Direct Imaging of Exoplanets
Comments: 12 pages, no figures
Submitted: 2018-11-12
One of the major difficulties limiting ground-based direct imaging of exoplanets with adaptive optics is quasi-static speckles in the science camera (SC) that obscure the planetary image. These speckles are caused by aberrations, called non-common path aberrations (NCPA), that are not corrected in the adaptive optics loop, and all attempts to subtract them in post-processing have been problematic. The method of Frazin (2013) (F13) uses simultaneous millisecond telemetry from wavefront sensor (WFS) and the SC to estimate the both the NCPA and the exoplanet image in a self-consistent manner. Rodack et al. (2018) proposed correcting for the NCPA in real-time while on-sky using the F13 estimation method, and called this procedure the "Real-Time Frazin Algorithm." The original regression model underlying the F13 method did not account for uncertainty in the WFS measurements, and this cannot be done with standard statistical methodology since these uncertainties manifest themselves in the independent variables (i.e., they cannot be treated as another source of noise in the SC data). Further, simulations show that simply using the noisy wavefront measurements without accounting for their uncertainties leads to estimates of the NCPA with unacceptably large bias. Here, the source of this bias is explained in terms of an "errors in variables" statistical model. Then, the method of F13 is generalized to account for WFS measurement error using a new sequential estimation technique that treats the nonlinear coupling between NCPA, WFS measurements and the error covariance of the WFS measurements. This new technique keeps a running estimate of the NCPA, the exoplanet image and their joint covariance matrix. The sequential implementation of the method should make it computationally efficient enough to be suitable for on-sky correction of the NCPA as well as off-line analysis.
[3]  [pdf] - 1741340
Direct Imaging in Reflected Light: Characterization of Older, Temperate Exoplanets With 30-m Telescopes
Comments: (March 2018) Submitted to the Exoplanet Science Strategy committee of the NAS
Submitted: 2018-08-29
Direct detection, also known as direct imaging, is a method for discovering and characterizing the atmospheres of planets at intermediate and wide separations. It is the only means of obtaining spectra of non-transiting exoplanets. Characterizing the atmospheres of planets in the <5 AU regime, where RV surveys have revealed an abundance of other worlds, requires a 30-m-class aperture in combination with an advanced adaptive optics system, coronagraph, and suite of spectrometers and imagers - this concept underlies planned instruments for both TMT (the Planetary Systems Imager, or PSI) and the GMT (GMagAO-X). These instruments could provide astrometry, photometry, and spectroscopy of an unprecedented sample of rocky planets, ice giants, and gas giants. For the first time habitable zone exoplanets will become accessible to direct imaging, and these instruments have the potential to detect and characterize the innermost regions of nearby M-dwarf planetary systems in reflected light. High-resolution spectroscopy will not only illuminate the physics and chemistry of exo-atmospheres, but may also probe rocky, temperate worlds for signs of life in the form of atmospheric biomarkers (combinations of water, oxygen and other molecular species). By completing the census of non-transiting worlds at a range of separations from their host stars, these instruments will provide the final pieces to the puzzle of planetary demographics. This whitepaper explores the science goals of direct imaging on 30-m telescopes and the technology development needed to achieve them.
[4]  [pdf] - 1678569
Efficient, Nonlinear Phase Estimation with the Non-Modulated Pyramid Wavefront Sensor
Comments: 14 pages, 6 figures. arXiv admin note: text overlap with arXiv:1708.09044
Submitted: 2018-04-02
The sensitivity of the the pyramid wavefront sensor (PyWFS) has made it a popular choice for astronomical adaptive optics (AAO) systems, and it is at its most sensitive when it is used without modulation of the input beam. In non-modulated mode, the device is highly nonlinear. Hence, all PyWFS implementations on current AAO systems employ modulation to make the device more linear. The upcoming era of 30-m class telescopes and the demand for ultra-precise wavefront control stemming from science objectives that include direct imaging of exoplanets make using the PyWFS without modulation desirable. This article argues that nonlinear estimation based on Newton's method for nonlinear optimization can be useful for mitigating the effects of nonlinearity in the non-modulated PyWFS. The proposed approach requires all optical modeling to be pre-computed, which has the advantage of avoiding real-time simulations of beam propagation. Further, the required real-time calculations are amenable to massively parallel computation. Numerical experiments simulate a currently operational PyWFS. A singular value analysis shows that the common practice of calculating two "slope" images from the four PyWFS pupil images discards critical information and is unsuitable for the non-modulated PyWFS simulated here. Instead, this article advocates estimators that use the raw pixel values not only from the four geometrical images of the pupil, but from surrounding pixels as well. The simulations indicate that nonlinear estimation can be effective when the Strehl ratio of the input beam is greater than 0.3, and the improvement relative to linear estimation tends to increase at larger Strehl ratios. At Strehl ratios less than about 0.5, the performances of both the nonlinear and linear estimators are relatively insensitive to noise, since they are dominated by nonlinearity error.
[5]  [pdf] - 1587665
Fast, Nonlinear Phase Estimation with the Non-Modulated Pyramid Wavefront Sensor at Low Strehl Ratio
Comments: 9 pages, 4 figures
Submitted: 2017-08-29
Most adaptive optics (AO) systems using pyramid wavefront sensors (PyWFS) to estimate the phase of the pupil field use mechanical modulation of the beam in order to increase the dynamic range in low-order modes so the PyWFS can usefully operate at low Strehl ratio. The tradeoff for this approach is reduced sensitivity, which, in turn, makes it difficult to attain a high Strehl ratio once the loop has been closed. We propose an algorithm that increases the dynamic range of the PyWFS without modulation. The proposed algorithm achieves this in two ways: 1) it allows the PyWFS to be treated with any desired optical modeling algorithms, and 2) it employs Newton's method for nonlinear optimization to create an estimator that is more accurate than the corresponding linear estimator. Numerical simulations show that nonlinear optimization can make more accurate estimates of the phase of the wavefront than the corresponding linear estimator for Strehl ratios of the input beam that are greater than about 0.2. As the input Strehl ratio increases, so does the advantage of the nonlinear estimator over the linear one. For example, when the input beam had a Strehl ratio of 0.4 (corresponding to a standard deviation of the phase of about 0.96 radians), the linear estimator error had a standard deviation of about 0.65 radians, while the nonlinear estimator error had a standard deviation of about 0.27 radians (this depends only weakly on the noise level, assuming there is enough signal for the PyWFS to work in the linear regime). The new algorithm can be implemented in massively parallel modes, since the required calculations have almost no inter-dependencies. It is suggested that the required computations can be performed quickly enough for the purposes of AO on modern computer systems.
[6]  [pdf] - 1584143
Energy Input Flux in the Global Quiet-Sun Corona
Submitted: 2017-06-01
We present first results of a novel technique that provides, for the first time, constraints on the energy input flux at the coronal base ($r\sim\,1.025\,{\rm R}_\odot$) of the quiet-Sun at a global scale. By combining differential emission measure tomography (DEMT) of EUV images, with global models of the coronal magnetic field, we estimate the energy input flux at the coronal base that is required to maintain thermodynamically stable structures. The technique is described in detail and first applied to data provided by the Extreme Ultraviolet Imager (EUVI) instrument, on board the Solar TErrestrial RElations Observatory (STEREO) mission, and the Atmospheric Imaging Assembly (AIA) instrument, on board the Solar Dynamics Observatory (SDO) mission, for two solar rotations with different levels of activity. Our analysis indicates that the typical energy input flux at the coronal base of magnetic loops in the quiet-Sun is in the range $\sim\,0.5-2.0\times 10^5\,{\rm (erg\,sec^{-1}\,{\rm cm}^{-2})}$, depending on the structure size and level of activity. A large fraction of this energy input, or even its totality, could be accounted for by Alfv\'en waves, as shown by recent independent observational estimates derived from determinations of the non-thermal broadening of spectral lines in the coronal base of quiet-Sun regions. This new tomography product will be useful for validation of coronal heating models in magnetohydrodinamic simulations of the global corona.
[7]  [pdf] - 1458619
Finding the Dark Hole with the Lights On: A New Approach to Focal Plane Wavefront Sensing
Comments: 14 pages, 2 figures
Submitted: 2016-08-16, last modified: 2016-08-17
In direct imaging of exoplanets from space, achieving the required dynamic range (i.e., planet-to-star contrast in brightness) currently relies on coronagraphic technology combined with active control of one or more deformable mirrors (DMs) to create a dark region in the image plane, sometimes called a "dark hole." While many algorithms have been proposed for this purpose, all of them employ focal plane wavefront sensing (FPWS) in order to calculate the optimal DM configuration to create the desired dark hole. All current algorithms are limited by their own success in that, as the dark hole is achieved, the FPWS procedure becomes shot-noise limited due to he low intensity in the dark hole. This article proposes a FPWS procedure that allows determination of the optimal DM configuration without relying on information obtained when the DM is near the optimal configuration. This article gives regression procedures for FPWS that do not assume the DM step size is small, which should allow two important improvements to the control loop: 1) performing informative FPWS observations with DM configurations that are sufficiently distant from the optimal dark hole configuration to mitigate shot-noise limitations, and 2) more accurately predicting the DM configuration that will achieve the desired objective in the dark hole control loop. In order to treat this more challenging FPWS problem, two different representations are presented. The first of these, is called the empirical Green's function (EGF), is easy to implement, and has a block-diagonal matrix structure that is well-suited to parallel processing. The other representation, based on an explicit aberration expansion (EAE) requires the regression to estimate a smaller number of parameters than the EGF, but leads to a dense matrix structure. The EGF and EAE methods both simultaneously estimate the planetary image.
[8]  [pdf] - 1521271
Empirical Green's Function Approach for Utilizing Millisecond Focal and Pupil Plane Telemetry in Exoplanet Imaging
Comments: 9 pages, 1 figure
Submitted: 2016-07-28
Millisecond focal plane telemetry is now becoming practical due to a new generation of near-IR detector arrays with sub-electron noise that are capable of kHz readout rates. Combining these data with those simultaneously available from the wavefront sensing system allows the possibility of self-consistently determining the optical aberrations (the cause of quasi-static speckles) and the planetary image. This approach may be especially advantageous for finding planets within about 3 $\lambda / D$ of the star where differential imaging is ineffective. As shown in a recent article by the author (J. Opt. Soc. Am. A., 33, 712, 2016), one must account for unknown aberrations in several non-conjugate planes of the optical system, which, in turn, requires ability to computational propagate the field between these planes. These computations are likely to be difficult to implement and expensive. Here, a far more convenient alternative based on empirical Green's functions is provided. It is shown that the empirical Green's function (EGF), which accounts for all multi-planar, non-common path aberrations, and results in a much more tractable and highly parallel computational problem. It is also shown that the EGF can be generalized to treat polarization, resulting in the empirical Green's tensor (EGT).
[9]  [pdf] - 1489885
Time-Dependent Tomographic Reconstruction of the Solar Corona
Comments: 24 pages, 18 figures
Submitted: 2016-07-21
Solar rotational tomography (SRT) applied to white-light coronal images observed at multiple aspect angles has been the preferred approach for determining the three-dimensional (3D) electron density structure of the solar corona. However, it is seriously hampered by the restrictive assumption that the corona is time-invariant which introduces significant errors in the reconstruction. We first explore several methods to mitigate the temporal variation of the corona by decoupling the "fast-varying" inner corona from the "slow-moving" outer corona using multiple masking (either by juxtaposition or recursive combination) and radial weighting. Weighting with a radial exponential profile provides some improvement over a classical reconstruction but only beyond 3 Rsun. We next consider a full time-dependent tomographic reconstruction involving spatio-temporal regularization and further introduce a co-rotating regularization aimed at preventing concentration of reconstructed density in the plane of the sky. Crucial to testing our procedure and properly tuning the regularization parameters is the introduction of a time-dependent MHD model of the corona based on observed magnetograms to build a time-series of synthetic images of the corona. Our procedure, which successfully reproduces the time-varying model corona, is finally applied to a set of of 53 LASCO-C2 pB images roughly evenly spaced in time from 15 to 29 March 2009. Our procedure paves the way to a time-dependent tomographic reconstruction of the coronal electron density to the whole set of LASCO-C2 images presently spanning 20 years.
[10]  [pdf] - 1392823
A Statistical Framework for the Utilization of Simultaneous Pupil Plane and Focal Plane Telemetry for Exoplanet Imaging, Part I: Accounting for Aberrations in Multiple Planes
Comments: 21 pages, 3 figures
Submitted: 2016-01-29, last modified: 2016-03-03
A new generation of telescopes with mirror diameters of 20 m or more, called extremely large telescopes (ELTs) has the potential to provide unprecedented imaging and spectroscopy of exo-planetary systems, if the difficulties in achieving the extremely high dynamic range required to differentiate the planetary signal from the star can be overcome to a sufficient degree. Fully utilizing the potential of ELTs for exoplanet imaging will likely require simultaneous and self-consistent determination of both the planetary image and the unknown aberrations in multiple planes of the optical system, using statistical inference based on the wavefront sensor and science camera data streams. This approach promises to overcome the most important systematic errors inherent in the various schemes based on differential imaging, such as ADI and SDI. This paper is the first in a series on this subject, in which a formalism is established for the exoplanet imaging problem, setting the stage for the statistical inference methods to follow in the future. Every effort has been made to be rigorous and complete, so that validity of approximations to be made later can be assessed. Here, the polarimetric image is expressed in terms of aberrations in the various planes of a polarizing telescope with an adaptive optics system. Further, it is shown that current methods that utilize focal plane sensing to correct the speckle field, e.g., electric field conjugation, rely on the tacit assumption that aberrations on multiple optical surfaces can be represented as aberration on a single optical surface, ultimately limiting their potential effectiveness for ground-based astronomy.
[11]  [pdf] - 1360454
A Statistical Framework for Utilization of Simultaneous Pupil Plane and Focal Plane Telemetry for Exoplanet Imaging, Part II: The Science Camera Image as a Function of the Wavefront Sensor Field
Comments: 18 pages
Submitted: 2016-02-17
In an effort to transcend the limitations of differential imaging of exoplanets in the era of extremely large telescopes (ELTs), the first paper in this series established a rigorous, fully polarimetric framework for determining the science camera (SC) image given a turbulent wavefront and unknown aberrations in multiple planes the optical system. This article builds on the structure developed in Paper I in order to rigorously express the polarimetric SC image in terms of the field impinging on the wavefront sensor (WFS), thereby providing a direct connection between the measurements made in both subsystems. This formulation allows the SC image to be written as a function of the WFS measurements, including the following unknown quantities which can, in principle, be estimated via statistical inference: the non-common path aberration (NCPA), WFS gain errors, aberrations downstream of the beamsplitter, and the planetary image. It is demonstrated that WFS bias error is mathematically equivalent to NCPA. Thus, with the ability to treat WFS bias and gains, the method should not be overly sensitive to WFS calibration problems. Importantly, this formulation includes stochastic processes that represent noisy measurement of the SC image, noisy WFS measurements, and high-frequency components of the wavefront to which the WFS is insensitive. It is shown that wavefront error due to noise in the WFS measurements has a convenient semi-analytical representation in terms of the WFS measurement operator's singular functions. Further, the first and second order statistics of these processes are specified, thereby setting the stage for the application of statistical inference methods to be describe in later papers in this series.
[12]  [pdf] - 1251902
Direct Imaging of Exoplanets Without Background Subtraction: Implications for ELTs
Comments: This article contains little new work. 12 pages, 5 figures. arXiv admin note: substantial text overlap with arXiv:1404.7453
Submitted: 2015-02-05, last modified: 2015-07-28
The ultra-high contrast capability required to form images of other solar systems is arguably the highest-profile challenge in astronomy today. The current high-contrast imaging efforts all require background subtraction to separate the planetary image from the image of the host star. Background estimation is difficult due to the presence of non-common path aberrations (NCPAs) that change with time. The only major source of information that is not being utilized by current efforts is the random encoding of the planetary image and the NCPAs by the atmosphere on millisecond time-scales. Here, a method that utilizes this information in order to avoid background subtraction altogether is proposed. This new paradigm will allow simultaneous estimation of the time-dependent NCPAs and the planetary image via rigorous statistical inference procedures. These procedures are fully compatible with other information sources, such as diurnal field rotation and spectral diversity. Given the open-ended nature of the background subtraction issues, the ideas explained herein may well the key to imaging habitable planets with Extremely Large Telescopes (ELTs). Fully exploiting the information content of millisecond exposures will require significant design modifications of the ELT wavefront sensors and science camera systems, if ultra-high contrast imaging is to be priority.
[13]  [pdf] - 1288333
Multimodal Differential Emission Measure in the Solar Corona
Comments: 52 pages, 18 figures
Submitted: 2015-03-09
The Atmospheric Imaging Assembly (AIA) telescope on board the Solar Dynamics Observatory (SDO) provides coronal EUV imaging over a broader temperature sensitivity range than the previous generations of instruments (EUVI, EIT, and TRACE). Differential emission measure tomography (DEMT) of the solar corona based on AIA data is presented here for the first time. The main product of DEMT is the three-dimensional (3D) distribution of the local differential emission measure (LDEM). While in previous studies, based on EIT or EUVI data, there were 3 available EUV bands, with a sensitivity range $\sim 0.60 - 2.70$ MK, the present study is based on the 4 cooler AIA bands (aimed at studying the quiet sun), sensitive to the range $\sim 0.55 - 3.75$ MK. The AIA filters allow exploration of new parametric LDEM models. Since DEMT is better suited for lower activity periods, we use data from Carrington Rotation 2099, when the Sun was in its most quiescent state during the AIA mission. Also, we validate the parametric LDEM inversion technique by applying it to standard bi-dimensional (2D) differential emission measure (DEM) analysis on sets of simultaneous AIA images, and comparing the results with DEM curves obtained using other methods. Our study reveals a ubiquitous bimodal LDEM distribution in the quiet diffuse corona, which is stronger for denser regions. We argue that the nanoflare heating scenario is less likely to explain these results, and that alternative mechanisms, such as wave dissipation appear better supported by our results.
[14]  [pdf] - 1128971
A Steady-State Picture of Solar Wind Acceleration and Charge State Composition Derived from a Global Wave-Driven MHD Model
Comments: Submitted to the Astrophysical Journal
Submitted: 2014-12-29
The higher charge states found in slow ($<$400km s$^{-1}$) solar wind streams compared to fast streams have supported the hypothesis that the slow wind originates in closed coronal loops, and released intermittently through reconnection. Here we examine whether a highly ionized slow wind can also form along steady and open magnetic field lines. We model the steady-state solar atmosphere using AWSoM, a global magnetohydrodynamic model driven by Alfv{\'e}n waves, and apply an ionization code to calculate the charge state evolution along modeled open field lines. This constitutes the first charge states calculation covering all latitudes in a realistic magnetic field. The ratios $O^{+7}/O^{+6}$ and $C^{+6}/C^{+5}$ are compared to in-situ Ulysses observations, and are found to be higher in the slow wind, as observed; however, they are under-predicted in both wind types. The modeled ion fractions of S, Si, and Fe are used to calculate line-of-sight intensities, which are compared to EIS observations above a coronal hole. The agreement is partial, and suggests that all ionization rates are under-predicted. Assuming the presence of suprathermal electrons improved the agreement with both EIS and Ulysses observations; importantly, the trend of higher ionization in the slow wind was maintained. The results suggest there can be a sub-class of slow wind that is steady and highly ionized. Further analysis shows it originates from coronal hole boundaries (CHB), where the modeled electron density and temperature are higher than inside the hole, leading to faster ionization. This property of CHBs is global, and observationally supported by EUV tomography.
[15]  [pdf] - 1209240
Simultaneous ultra-high contrast imaging and determination of time-dependent, non-common path aberrations in the presence of detector noise
Comments: 7 pages, 3 figures, to be published in Proc. SPIE
Submitted: 2014-04-29, last modified: 2014-05-05
Ground-based ultra-high contrast imaging, as required for direct imaging of exoplanets and other solar systems, is limited by difficulty of separating the planetary emission from the effects of optical aberrations that are not compensated by the adaptive optics (AO) system, so-called "non-common path aberrations" (NCPAs). Simultaneous ($\sim$ millisecond) exposures by the science camera and the AO system enable the use of "phase diversity" to estimate both the NCPAs and the scene via a processing procedure first described by the author (R. Frazin 2013, ApJ, 767, article id. 21). This method is fully compatible with more standard concepts used in long-exposure high-contrast imaging, such as angular differential imaging and spectral deconvolution. Long-exposure methods find time-dependent NCPAs, such as those caused by vibrations, particularly challenging. Here, an NCPA of the form of $\alpha \cos(k \cdot r - \omega t + \vartheta)$ is considered. It is shown that, when sampled at millisecond time-scales, the image plane data are sensitive to $\mbox{arg}(\alpha)$, $\vartheta$ and $\omega$, and, therefore such NCPAs can be simultaneously estimated with the scene. Simulations of observations with ms exposure times are reported. These simulations include substantial detector noise and a sinusoidal NCPA that places a speckle exactly at the location of a planet. Simulations show that the effects of detector noise can be mitigated by mixing exposures of various lengths, allowing estimation of the planet's brightness.
[16]  [pdf] - 1157949
Utilization of the Wavefront Sensor and Short-Exposure Images for Simultaneous Estimation of Quasi-static Aberration and Exoplanet Intensity
Submitted: 2012-11-20, last modified: 2013-02-05
This paper provides a framework for the incorporation of the wavefront sensor measurements in the context of observing modes in which the science camera takes millisecond exposures. In this formulation, the wavefront sensor measurements provide a means to jointly estimate the static speckle and the planetary signal. The ability to estimate planetary intensities in as little as few seconds has the potential to greatly improve the efficiency of exoplanet search surveys. Unlike currently used methods, in which increasing the observation time beyond a certain threshold is useless, this method produces estimates whose error covariances decrease more quickly than inversely proportional to the observation time. This is due to the fact that the estimates of the quasi-static aberrations are informed by a new random (but approximately known) wavefront every millisecond. The method can be extended to include angular (due to diurnal field rotation) and spectral diversity. Numerical experiments are performed with wavefront data from the AEOS Adaptive Optics System sensing at 850 nm. These experiments assume a science camera wavelength $\lambda$\ of $1.1 \, \mu$, that the measured wavefronts are exact, and a Gaussian approximation of shot-noise. The effects of detector read-out noise and other issues are left to future investigations. A number of static aberrations are introduced, including one with a spatial frequency exactly corresponding the planet location, which was at a distance of $\approx 3\lambda/D$\ from the star. Using only 4 seconds of simulated observation time, a planetary intensity, of $\approx 1$ photon/ms, a stellar intensity of $\approx 10^5$\ photons/ms (contrast ratio $10^5$), the short-exposure estimation method recovers the amplitudes static aberrations with a 1% accuracy, and the planet brightness with a with 20% accuracy.
[17]  [pdf] - 1125155
Newly Discovered Global Temperature Structures in the Quiet Sun at Solar Minimum
Comments: Accepted for publication in The Astrophysical Journal, waiting for the full biblio info
Submitted: 2012-07-27
Magnetic loops are building blocks of the closed-field corona. While active region loops are readily seen in images taken at EUV and X-ray wavelengths, quiet Sun loops are seldom identifiable and therefore difficult to study on an individual basis. The first analysis of solar minimum (Carrington Rotation 2077) quiet Sun (QS) coronal loops utilizing a novel technique called the Michigan Loop Diagnostic Technique (MLDT) is presented. This technique combines Differential Emission Measure Tomography (DEMT) and a potential field source surface (PFSS) model, and consists of tracing PFSS field lines through the tomographic grid on which the Local Differential Emission Measure (LDEM) is determined. As a result, the electron temperature Te and density Ne at each point along each individual field line can be obtained. Using data from STEREO/EUVI and SOHO/MDI, the MLDT identifies two types of QS loops in the corona: so-called "up" loops in which the temperature increases with height, and so-called "down" loops in which the temperature decreases with height. Up loops are expected, however, down loops are a surprise, and furthermore, they are ubiquitous in the low-latitude corona. Up loops dominate the QS at higher latitudes. The MLDT allows independent determination of the empirical pressure and density scale heights, and the differences between the two remain to be explained. The down loops appear to be a newly discovered property of the solar minimum corona that may shed light on the physics of coronal heating. The results are shown to be robust to the calibration uncertainties of the EUVI instrument.
[18]  [pdf] - 622363
The First Stray Light Corrected EUV Images of Solar Coronal Holes
Comments: Accepted to Astrophysical Journal Letters
Submitted: 2011-10-13, last modified: 2012-03-07
Coronal holes are the source regions of the fast solar wind, which fills most of the solar system volume near the cycle minimum. Removing stray light from extreme ultraviolet (EUV) images of the Sun's corona is of high astrophysical importance, as it is required to make meaningful determinations of temperatures and densities of coronal holes. EUV images tend to be dominated by the component of the stray light due to the long-range scatter caused by microroughness of telescope mirror surfaces, and this component has proven very difficult to measure in pre-flight characterization. In-flight characterization heretofore has proven elusive due to the fact that the detected image is simultaneously nonlinear in two unknown functions: the stray light pattern and the true image which would be seen by an ideal telescope. Using a constrained blind deconvolution technique that takes advantage of known zeros in the true image provided by a fortuitous lunar transit, we have removed the stray light from solar images seen by the EUVI instrument on STEREO-B in all four filter bands (171, 195, 284, and 304 \AA). Uncertainty measures of the stray light corrected images, which include the systematic error due to misestimation of the scatter, are provided. It is shown that in EUVI, stray light contributes up to 70% of the emission in coronal holes seen on the solar disk, which has dramatic consequences for diagnostics of temperature and density and therefore estimates of key plasma parameters such as the plasma $\beta$\ and ion-electron collision rates.
[19]  [pdf] - 1042804
The WHI Corona from Differential Emission Measure Tomography
Comments: 13 Figures
Submitted: 2010-12-29
A three dimensional (3D) tomographic reconstruction of the local differential emission measure (LDEM) of the global solar corona during the whole heliosphere interval (WHI, Carrington rotation CR-2068) is presented, based on STEREO/EUVI images. We determine the 3D distribution of the electron density, mean temperature, and temperature spread, in the range of heliocentric heights 1.03 to 1.23 Rsun. The reconstruction is complemented with a potential field source surface (PFSS) magnetic-field model. The streamer core, streamer legs, and subpolar regions are analyzed and compared to a similar analysis previously performed for CR-2077, very near the absolute minimum of the Solar Cycle 23. In each region, the typical values of density and temperature are similar in both periods. The WHI corona exhibits a streamer structure of relatively smaller volume and latitudinal extension than during CR-2077, with a global closed-to-open density contrast about 6% lower, and a somewhat more complex morphology. The average basal electron density is found to be about 2.23 and 1.08 x 10^8 cm^-3, in the streamer core and subpolar regions, respectively. The electron temperature is quite uniform over the analyzed height range, with average values of about 1.13 and 0.93 MK, in the streamer core and subpolar regions, respectively. Within the streamer closed region, both periods show higher temperatures at mid-latitudes and lower temperatures near the equator. Both periods show beta>1 in the streamer core and beta<1 in the surrounding open regions, with CR-2077 exhibiting a stronger contrast. Hydrostatic fits to the electron density are performed, and the scale height is compared to the LDEM mean electron temperature. Within the streamer core, the results are consistent with an isothermal hydrostatic plasma regime, with the temperatures of ions and electrons differing by up to about 10% .. (continues)..
[20]  [pdf] - 272935
Time-Dependent Reconstruction of Non-Stationary Objects with Tomographic or Interferometric Measurements
Comments: 12 pages, 3 figures
Submitted: 2005-09-15
In a number of astrophysical applications one tries to determine the two-dimensional or three-dimensional structure of an object from a time series of measurements. While most methods used for reconstruction assume that object is static, the data are often acquired over a time interval during which the object may change significantly. This problem may be addressed with time-dependent reconstruction methods such as Kalman filtering which models the temporal evolution of the unknown object as a random walk that may or may not have a deterministic component. Time-dependent reconstructions of a hydrodynamical simulation from its line-integral projections are presented. In these examples standard reconstructions based on the static assumption are poor while the Kalman based reconstructions are of good quality. Implications for various astrophysical applications, including tomography of the solar corona and radio aperture synthesis, are discussed.