Showing 1–26 of 26 results for author: Frazin, R

High contrast imaging at the photon noise limit with self-calibrating WFS/C systems

Authors: Olivier Guyon, Barnaby Norris, Marc-Antoine Martinod, Kyohoon Ahn, Peter Tuthill, Jared Males, Alison Wong, Nour Skaf, Thayne Currie, Kelsey Miller, Steven P. Bos, Julien Lozi, Vincent Deo, Sebastien Vievard, Ruslan Belikov, Kyle van Gorkom, Benjamin Mazin, Michael Bottom, Richard Frazin, Alexander Rodack, Tyler Groff, Nemanja Jovanovic, Frantz Martinache

Abstract: High contrast imaging (HCI) systems rely on active wavefront control (WFC) to deliver deep raw contrast in the focal plane, and on calibration techniques to further enhance contrast by identifying planet light within the residual speckle halo. Both functions can be combined in an HCI system and we discuss a path toward designing HCI systems capable of calibrating residual starlight at the fundamen… ▽ More High contrast imaging (HCI) systems rely on active wavefront control (WFC) to deliver deep raw contrast in the focal plane, and on calibration techniques to further enhance contrast by identifying planet light within the residual speckle halo. Both functions can be combined in an HCI system and we discuss a path toward designing HCI systems capable of calibrating residual starlight at the fundamental contrast limit imposed by photon noise. We highlight the value of deploying multiple high-efficiency wavefront sensors (WFSs) covering a wide spectral range and spanning multiple optical locations. We show how their combined information can be leveraged to simultaneously improve WFS sensitivity and residual starlight calibration, ideally making it impossible for an image plane speckle to hide from WFS telemetry. We demonstrate residual starlight calibration in the laboratory and on-sky, using both a coronagraphic setup, and a nulling spectro-interferometer. In both case, we show that bright starlight can calibrate residual starlight. △ Less

Submitted 4 October, 2021; v1 submitted 28 September, 2021; originally announced September 2021.

Comments: 11 pages, 6 Figures, 1 Table; Proc. SPIE in press

SCExAO, a testbed for developing high-contrast imaging technologies for ELTs

Authors: Kyohoon Ahn, Olivier Guyon, Julien Lozi, Sébastien Vievard, Vincent Deo, Nour Skaf, Ruslan Belikov, Steven P. Bos, Michael Bottom, Thayne Currie, Richard Frazin, Kyle V. Gorkom, Tyler D. Groff, Sebastiaan Y. Haffert, Nemanja Jovanovic, Hajime Kawahara, Takayuki Kotani, Jared R. Males, Frantz Martinache, Benjamin A. Mazin, Kelsey Miller, Barnaby Norris, Alexander Rodack, Alison Wong

Abstract: To directly detect exoplanets and protoplanetary disks, the development of high accuracy wavefront sensing and control (WFS&C) technologies is essential, especially for ground-based Extremely Large Telescopes (ELTs). The Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) instrument is a high-contrast imaging platform to discover and characterize exoplanets and protoplanetary disks. It also serv… ▽ More To directly detect exoplanets and protoplanetary disks, the development of high accuracy wavefront sensing and control (WFS&C) technologies is essential, especially for ground-based Extremely Large Telescopes (ELTs). The Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) instrument is a high-contrast imaging platform to discover and characterize exoplanets and protoplanetary disks. It also serves as a testbed to validate and deploy new concepts or algorithms for high-contrast imaging approaches for ELTs, using the latest hardware and software technologies on an 8-meter class telescope. SCExAO is a multi-band instrument, using light from 600 to 2500 nm, and delivering a high Strehl ratio (>80% in median seeing in H-band) downstream of a low-order correction provided by the facility AO188. Science observations are performed with coronagraphs, an integral field spectrograph, or single aperture interferometers. The SCExAO project continuously reaches out to the community for development and upgrades. Existing operating testbeds such as the SCExAO are also unique opportunities to test and deploy the new technologies for future ELTs. We present and show a live demonstration of the SCExAO capabilities (Real-time predictive AO control, Focal plane WFS&C, etc) as a host testbed for the remote collaborators to test and deploy the new WFS&C concepts or algorithms. We also present several high-contrast imaging technologies that are under development or that have already been demonstrated on-sky. △ Less

Submitted 28 September, 2021; v1 submitted 27 September, 2021; originally announced September 2021.

Comments: 13 pages, 8 figures, 2021 SPIE Optics+Photonics

Millisecond Exoplanet Imaging, II: Regression Equations and Technical Discussion

Authors: Richard A Frazin, Alexander T Rodack

Abstract: The leading difficulty in achieving the contrast necessary to directly image exoplanets and associated structures (eg. protoplanetary disks) at wavelengths ranging from the visible to the infrared are quasi-static speckles, and they are hard to distinguish from planets at the necessary level of precision. The source of the quasi-static speckles is hardware aberrations that are not compensated by t… ▽ More The leading difficulty in achieving the contrast necessary to directly image exoplanets and associated structures (eg. protoplanetary disks) at wavelengths ranging from the visible to the infrared are quasi-static speckles, and they are hard to distinguish from planets at the necessary level of precision. The source of the quasi-static speckles is hardware aberrations that are not compensated by the adaptive optics system. These aberrations are called non-common path aberrations (NCPA). In 2013, Frazin showed how, in principle, simultaneous millisecond (ms) telemetry from the wavefront sensor (WFS) and the science camera behind a stellar coronagraph can be used as input into a regression scheme that simultaneously and self-consistently estimates the NCPA and the sought-after image of the planetary system (the exoplanet image). The physical principle underlying the regression method is rather simple: the wavefronts, which are measured by the WFS, modulate the speckles caused by the NCPA and therefore can be used as probes of the optical system. The most important departure from realism in the author's 2013 article was the assumption that the WFS made error-free measurements. The simulations in Part I provide results on the joint regression on the NCPA and the exoplanet image from three different methods, called the ideal, the naive, and the bias-corrected estimators. The ideal estimator is not physically realizable but is a useful as a benchmark for simulation studies, but the other two are, at least in principle. This article provides the regression equations for all three of these estimators as well as a supporting technical discussion. Briefly, the naive estimator simply uses the noisy WFS measurements without any attempt to account for the errors, and the bias-corrected estimator uses statistical knowledge of the wavefronts to treat errors in the WFS measurements. △ Less

Submitted 30 June, 2021; v1 submitted 13 May, 2021; originally announced May 2021.

Comments: 13 pages, 2 figures, submitted to JOSAA

Millisecond Exoplanet Imaging, I: Method and Simulation Results

Authors: Alexander T Rodack, Richard A Frazin, Jared R Males, Olivier Guyon

Abstract: One of the top remaining science challenges in astronomical optics is the direct imaging and characterization of extrasolar planets and planetary systems. Directly imaging exoplanets from ground-based observatories requires combining high-order adaptive optics with a stellar coronagraph observing at wavelengths ranging from the visible to the mid-IR. A limiting factor in achieving the required con… ▽ More One of the top remaining science challenges in astronomical optics is the direct imaging and characterization of extrasolar planets and planetary systems. Directly imaging exoplanets from ground-based observatories requires combining high-order adaptive optics with a stellar coronagraph observing at wavelengths ranging from the visible to the mid-IR. A limiting factor in achieving the required contrast (planet-to-star intensity ratio) is quasi-static speckles, caused largely by non-common path aberrations (NCPA) in the coronagraph. Starting with a realistic simulator of a telescope with an AO system and a coronagraph, this article provides simulations of several closely related millisecond regression models requiring inputs of the measured wavefronts and science camera images. The simplest regression model, called the naive estimator, does not treat the noise and other sources of information loss in the WFS. The naive estimator provided a useful estimate of the NCPA of $\sim$ 0.5 radian RMS, with an accuracy of $\sim$ 0.06 radian RMS in one minute of simulated sky time on a magnitude 8 star. The bias-corrected estimator generalizes the regression model to account for the noise and information loss in the WFS. A simulation of the bias-corrected estimator with four minutes of sky time included an NCPA of $\sim 0.05 \,$ radian RMS and an extended exoplanet scene. The joint regression of the bias-corrected estimator simultaneously achieved an NCPA estimate with an accuracy of $\sim 5\times10^{-3} \,$radian and contrast of $\sim 10^{-5}$ on the exoplanet scene. In addition, the estimate of the exoplanet image was completely free of the subtraction artifacts that always plague differential imaging. The estimate of the exoplanet image obtained by the joint regression was nearly identical to the image obtained by subtraction of a perfectly known point-spread function. △ Less

Submitted 29 June, 2021; v1 submitted 13 May, 2021; originally announced May 2021.

Comments: 16 pages, 18 Figures, 4 Tables, submitted to JOSAA

Laboratory Demonstration of Spatial Linear Dark Field Control For Imaging Extrasolar Planets in Reflected Light

Authors: Thayne Currie, Eugene Pluzhnik, Olivier Guyon, Ruslan Belikov, Kelsey Miller, Steven Bos, Jared Males, Dan Sirbu, Charlotte Bond, Richard Frazin, Tyler Groff, Brian Kern, Julien Lozi, Benjamin Mazin, Bijan Nemati, Barnaby Norris, Hari Subedi, Scott Will

Abstract: Imaging planets in reflected light, a key focus of future NASA missions and ELTs, requires advanced wavefront control to maintain a deep, temporally correlated null of stellar halo -- i.e. a dark hole -- at just several diffraction beam widths. Using the Ames Coronagraph Experiment testbed, we present the first laboratory tests of Spatial Linear Dark Field Control (LDFC) approaching raw contrasts… ▽ More Imaging planets in reflected light, a key focus of future NASA missions and ELTs, requires advanced wavefront control to maintain a deep, temporally correlated null of stellar halo -- i.e. a dark hole -- at just several diffraction beam widths. Using the Ames Coronagraph Experiment testbed, we present the first laboratory tests of Spatial Linear Dark Field Control (LDFC) approaching raw contrasts ($\sim$ 5$\times$10$^{-7}$) and separations (1.5--5.2 $λ$/D) needed to image jovian planets around Sun-like stars with space-borne coronagraphs like WFIRST-CGI and image exo-Earths around low-mass stars with future ground-based 30m class telescopes. In four separate experiments and for a range of different perturbations, LDFC largely restores (to within a factor of 1.2--1.7) and maintains a dark hole whose contrast is degraded by phase errors by an order of magnitude. Our implementation of classical speckle nulling requires a factor of 2--5 more iterations and 20--50 DM commands to reach contrasts obtained by spatial LDFC. Our results provide a promising path forward to maintaining dark holes without relying on DM probing and in the low-flux regime, which may improve the duty cycle of high-contrast imaging instruments, increase the temporal correlation of speckles, and thus enhance our ability to image true solar system analogues in the next two decades. △ Less

Submitted 28 July, 2020; originally announced July 2020.

Comments: 13 pages, 7 figures, accepted for publication in Publications of the Astronomical Society of the Pacific

Thermodynamic Structure of the Solar Corona: Tomographic Reconstructions and MHD Modeling

Authors: D. G. Lloveras, A. M. Vásquez, F. A. Nuevo, C. Mac Cormack, N. Sachdeva, W. Manchester IV, B. Van der Holst, R. A. Frazin

Abstract: We carry out a study of the global three-dimensional (3D) structure of the electron density and temperature of the quiescent inner solar corona ($r<1.25 R_\odot$) by means of tomographic reconstructions and magnetohydrodynamic simulations. We use differential emission measure tomography (DEMT) and the Alfvén Wave Solar Model (AWSoM), in their latest versions. Two target rotations were selected fro… ▽ More We carry out a study of the global three-dimensional (3D) structure of the electron density and temperature of the quiescent inner solar corona ($r<1.25 R_\odot$) by means of tomographic reconstructions and magnetohydrodynamic simulations. We use differential emission measure tomography (DEMT) and the Alfvén Wave Solar Model (AWSoM), in their latest versions. Two target rotations were selected from the solar minimum between solar cycles (SCs) 23 and 24 and the declining phase of SC 24. We report in quantitative detail on the 3D thermodynamic structure of the core and outer layers of the streamer belt, and of the high latitude coronal holes (CH), as revealed by the DEMT analysis. We report on the presence of two types of structures within the streamer belt, loops with temperature decreasing/increasing with height (dubbed down/up loops), as reported first in previous DEMT studies. We also estimate the heating energy flux required at the coronal base to keep these structures stable, found to be or order $10^5 erg\, cm^{-2} s^{-1}$, consistently with previous DEMT and spectroscopic studies. We discuss how these findings are consistent with coronal dissipation of Alfvén waves. We compare the 3D results of DEMT and AWSoM in distinct magnetic structures. We show that the agreement between the products of both techniques is the best so far, with an overall agreement $\lesssim 20\%$, depending on the target rotation and the specific coronal region. In its current implementation the ASWsoM model can not reproduce down loops though. Also, in the source region of the fast and slow components of the solar wind, the electron density of the AWSoM model increases with latitude, opposite to the trend observed in DEMT reconstructions. △ Less

Submitted 14 April, 2020; originally announced April 2020.

Directly Imaging Rocky Planets from the Ground

Authors: B. Mazin, É. Artigau, V. Bailey, C. Baranec, C. Beichman, B. Benneke, J. Birkby, T. Brandt, J. Chilcote, M. Chun, L. Close, T. Currie, I. Crossfield, R. Dekany, J. R. Delorme, C. Dong, R. Dong, R. Doyon, C. Dressing, M. Fitzgerald, J. Fortney, R. Frazin, E. Gaidos, O. Guyon, J. Hashimoto , et al. (38 additional authors not shown)

Abstract: 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 spe… ▽ More 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 △ Less

Submitted 10 May, 2019; originally announced May 2019.

Comments: 8 pages, 1 figure, Astro2020 Science White Paper

Treating Wavefront Measurement Error in Estimation of Non-Common Path Aberration for Direct Imaging of Exoplanets

Authors: Richard A. Frazin

Abstract: 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.… ▽ More 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. △ Less

Submitted 12 November, 2018; originally announced November 2018.

Comments: 12 pages, no figures

Direct Imaging in Reflected Light: Characterization of Older, Temperate Exoplanets With 30-m Telescopes

Authors: Etienne Artigau, Rebecca A. Bernstein, Timothy Brandt, Jeffrey Chilcote, Laird Close, Ian Crossfield, Jacques-Robert Delorme, Courtney Dressing, Michael P. Fitzgerald, Jonathan Fortney, Andrew Howard, Richard Frazin, Nemanja Jovanovic, Quinn Konopacky, Julien Lozi, Jared R. Males, Christian Marois, Benjamin A. Mazin, Max A. Millar-Blanchaer, Katie M. Morzinski, Lewis Roberts, Eugene Serabyn, Gautam Vasisht, J. Kent Wallace, Ji Wang

Abstract: 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 combinatio… ▽ More 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. △ Less

Submitted 29 August, 2018; originally announced August 2018.

Comments: (March 2018) Submitted to the Exoplanet Science Strategy committee of the NAS

Efficient, Nonlinear Phase Estimation with the Non-Modulated Pyramid Wavefront Sensor

Authors: Richard A. Frazin

Abstract: 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… ▽ More 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. △ Less

Submitted 2 April, 2018; originally announced April 2018.

Comments: 14 pages, 6 figures. arXiv admin note: text overlap with arXiv:1708.09044

Journal ref: Journal of the Optical Society of America A, Vol. 35, No. 4. pp.594-607, 2018

Fast, Nonlinear Phase Estimation with the Non-Modulated Pyramid Wavefront Sensor at Low Strehl Ratio

Authors: Richard Frazin

Abstract: 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… ▽ More 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. △ Less

Submitted 29 August, 2017; originally announced August 2017.

Comments: 9 pages, 4 figures

Energy Input Flux in the Global Quiet-Sun Corona

Authors: C. Mac Cormack, A. M. Vásquez, M. López Fuentes, F. A. Nuevo, E. Landi, R. A. Frazin

Abstract: 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 requir… ▽ More 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én 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. △ Less

Submitted 1 June, 2017; originally announced June 2017.

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

Authors: Richard A. Frazin

Abstract: 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… ▽ More 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. △ Less

Submitted 17 August, 2016; v1 submitted 16 August, 2016; originally announced August 2016.

Comments: 14 pages, 2 figures

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

Authors: Richard A. Frazin

Abstract: 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… ▽ More 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). △ Less

Submitted 28 July, 2016; originally announced July 2016.

Comments: 9 pages, 1 figure

Time-Dependent Tomographic Reconstruction of the Solar Corona

Authors: Didier Vibert, Christelle Peillon, Philippe Lamy, Richard A. Frazin, Julien Wojak

Abstract: 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 expl… ▽ More 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. △ Less

Submitted 21 July, 2016; originally announced July 2016.

Comments: 24 pages, 18 figures

Journal ref: Astronomy and Computing, Volume 17, October 2016, Pages 144-162

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

Authors: Richard A. Frazin

Abstract: 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… ▽ More 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. △ Less

Submitted 17 February, 2016; originally announced February 2016.

Comments: 18 pages

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

Authors: Richard A. Frazin

Abstract: 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 potentia… ▽ More 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. △ Less

Submitted 3 March, 2016; v1 submitted 29 January, 2016; originally announced January 2016.

Comments: 21 pages, 3 figures

Multimodal Differential Emission Measure in the Solar Corona

Authors: Federico A. Nuevo, Alberto M. Vásquez, Enrico Landi, Richard Frazin

Abstract: 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-… ▽ More 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. △ Less

Submitted 9 March, 2015; originally announced March 2015.

Comments: 52 pages, 18 figures

Direct Imaging of Exoplanets Without Background Subtraction: Implications for ELTs

Authors: Richard A. Frazin

Abstract: 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 tim… ▽ More 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. △ Less

Submitted 28 July, 2015; v1 submitted 5 February, 2015; originally announced February 2015.

Comments: This article contains little new work. 12 pages, 5 figures. arXiv admin note: substantial text overlap with arXiv:1404.7453

A Steady-State Picture of Solar Wind Acceleration and Charge State Composition Derived from a Global Wave-Driven MHD Model

Authors: Rona Oran, Enrico Landi, Bart van der Holst, Susan T. Lepri, Alberto M. Vásquez, Federico. A. Nuevo, Richard Frazin, Ward B. Manchester IV, Igor V. Sokolov, Tamas I. Gombosi

Abstract: 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 AWSo… ▽ More 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{é}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. △ Less

Submitted 29 December, 2014; originally announced December 2014.

Comments: Submitted to the Astrophysical Journal

Journal ref: ApJ, 806, 55, 2015

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

Authors: Richard A. Frazin

Abstract: 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… ▽ More 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. △ Less

Submitted 5 May, 2014; v1 submitted 29 April, 2014; originally announced April 2014.

Comments: 7 pages, 3 figures, to be published in Proc. SPIE

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

Authors: Richard A. Frazin

Abstract: 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… ▽ More 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. △ Less

Submitted 5 February, 2013; v1 submitted 20 November, 2012; originally announced November 2012.

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

Authors: Zhenguang Huang, Richard A. Frazin, Enrico Landi, Ward B. Manchester IV, Alberto M. Vásquez, Tamas I. Gombosi

Abstract: 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 Diagno… ▽ More 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. △ Less

Submitted 27 July, 2012; originally announced July 2012.

Comments: Accepted for publication in The Astrophysical Journal, waiting for the full biblio info

The First Stray Light Corrected EUV Images of Solar Coronal Holes

Authors: Paul Shearer, Richard A. Frazin, Alfred O. Hero III, Anna C. Gilbert

Abstract: 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… ▽ More 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. △ Less

Submitted 7 March, 2012; v1 submitted 13 October, 2011; originally announced October 2011.

Comments: Accepted to Astrophysical Journal Letters

The WHI Corona from Differential Emission Measure Tomography

Authors: Alberto M. Vásquez, Zhenguang Huang, Ward B. Manchester IV, Richard A. Frazin

Abstract: 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.… ▽ More 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).. △ Less

Submitted 29 December, 2010; v1 submitted 29 December, 2010; originally announced December 2010.

Comments: 13 Figures

Time-Dependent Reconstruction of Non-Stationary Objects with Tomographic or Interferometric Measurements

Authors: R. A. Frazin, M. D. Butala, A. Kemball, F. Kamalabadi

Abstract: 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 reconstruct… ▽ More 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. △ Less

Submitted 15 September, 2005; originally announced September 2005.

Comments: 12 pages, 3 figures

Journal ref: Astrophys.J.635:L197-L200,2005