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] oai:arXiv.org:1905.04275 [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] oai:arXiv.org:1811.05096 [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] oai:arXiv.org:1808.09632 [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] oai:arXiv.org:1804.01011 [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] oai:arXiv.org:1708.09044 [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] oai:arXiv.org:1706.00365 [pdf] - 1584143

Energy Input Flux in the Global Quiet-Sun Corona

Comments:

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] oai:arXiv.org:1608.04616 [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] oai:arXiv.org:1607.08586 [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] oai:arXiv.org:1607.06308 [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] oai:arXiv.org:1601.08239 [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] oai:arXiv.org:1602.05486 [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] oai:arXiv.org:1502.01702 [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] oai:arXiv.org:1503.02473 [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] oai:arXiv.org:1412.8288 [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] oai:arXiv.org:1404.7453 [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] oai:arXiv.org:1211.4914 [pdf] - 1157949

Utilization of the Wavefront Sensor and Short-Exposure Images for
Simultaneous Estimation of Quasi-static Aberration and Exoplanet Intensity

Comments:

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] oai:arXiv.org:1207.6661 [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] oai:arXiv.org:1110.3052 [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] oai:arXiv.org:1012.5953 [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] oai:arXiv.org:astro-ph/0509472 [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.