Researcher profile

Richard A. Frazin

Richard A. Frazin contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

Empirical Green's Function Approach for Utilizing Millisecond Focal and Pupil Plane Telemetry in Exoplanet Imaging

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 $λ/ 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).

preprint2016arXiv

Finding the Dark Hole with the Lights On: A New Approach to Focal Plane Wavefront Sensing

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.

preprint2016arXiv

Time-Dependent Tomographic Reconstruction of the Solar Corona

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.

preprint2015arXiv

Direct Imaging of Exoplanets Without Background Subtraction: Implications for ELTs

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.

preprint2014arXiv

Simultaneous ultra-high contrast imaging and determination of time-dependent, non-common path aberrations in the presence of detector noise

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 $α\cos(k \cdot r - ωt + \vartheta)$ is considered. It is shown that, when sampled at millisecond time-scales, the image plane data are sensitive to $\mbox{arg}(α)$, $\vartheta$ and $ω$, 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.

preprint2013arXiv

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

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 $λ$\ of $1.1 \, μ$, 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λ/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.

preprint2012arXiv

Newly Discovered Global Temperature Structures in the Quiet Sun at Solar Minimum

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.

preprint2012arXiv

The First Stray Light Corrected EUV Images of Solar Coronal Holes

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 Å). 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 $β$\ and ion-electron collision rates.

preprint2010arXiv

The WHI Corona from Differential Emission Measure Tomography

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)..