Here we use synthetic data to explore the performance of forward models and inverse methods for helioseismic holography. Specifically, this work presents the first comprehensive test of inverse modeling for flows using lateral-vantage (deep-focus) holography. We derive sensitivity functions in the Born approximation. We then use these sensitivity functions in a series of forward models and inversions of flows from a publicly available magnetohydrodynamic quiet-Sun simulation. The forward travel times computed using the kernels generally compare favorably with measurements obtained by applying holography, in a lateral-vantage configuration, on a 15-hour time series of artificial Dopplergrams extracted from the simulation. Inversions for the horizontal flow components are able to reproduce the flows in the upper 3Mm of the domain, but are compromised by noise at greater depths.

Wave propagation through sunspots involves conversion between waves of acoustic and magnetic character. In addition, the thermal structure of sunspots is very different than that of the quiet Sun. As a consequence, the interpretation of local helioseismic measurements of sunspots has long been a challenge. With the aim of understanding these measurements, we carry out numerical simulations of wave propagation through sunspots. Helioseismic holography measurements made from the resulting simulated wavefields show qualitative agreement with observations of real sunspots. We use additional numerical experiments to determine, separately, the influence of the thermal structure of the sunspot and the direct effect of the sunspot magnetic field. We use the ray approximation to show that the travel-time shifts in the thermal (non-magnetic) sunspot model are primarily produced by changes in the wave path due to the Wilson depression rather than variations in the wave speed. This shows that inversions for the subsurface structure of sunspots must account for local changes in the density. In some ranges of horizontal phase speed and frequency there is agreement (within the noise level in the simulations) between the travel times measured in the full magnetic sunspot model and the thermal model. If this conclusion proves to be robust for a wide range of models, it would suggest a path towards inversions for sunspot structure.

We continue the investigation of how to use the divergence-free condition to resolve the azimuthal ambiguity present in vector magnetogram data. In previous articles, by Crouch, Barnes, and Leka (Solar Physics, 260, 271, 2009) and Crouch (Solar Physics, 282, 107, 2013), all methods used an expression for the divergence of the magnetic field that involves differentiation of quantities that depend on the choice of azimuthal angle. As a result, all heights used to approximate line-of-sight derivatives should generally be disambiguated simultaneously. In this article, we investigate a set of methods that use an expression for the divergence that involves differentiation of quantities that do not depend on the choice of azimuthal angle. This results in an expression for the divergence that can be used to disambiguate each height independently. We test two methods using synthetic and find that the two-step, hybrid method, adapted to disambiguate each height independently, generally produces reasonable results. Moreover, the time required to compute solutions is substantially decreased in comparison to the corresponding method that disambiguates all relevant heights simultaneously.

The helioseismic properties of the wave scattering generated by monolithic and spaghetti sunspots are analyzed by means of numerical simulations. In these computations, an incident f or p1 mode travels through the sunspot model, which produces absorption and phase shift of the waves. The scattering is studied by inspecting the wavefield, computing travel-time shifts, and performing Fourier-Hankel analysis. The comparison between the results obtained for both sunspot models reveals that the differences in the absorption coefficient can be detected above noise level. The spaghetti model produces an steep increase of the phase shift with the degree of the mode at short wavelengths, while mode-mixing is more efficient for the monolithic model. These results provide a clue for what to look for in solar observations to discern the constitution of sunspots between the proposed monolithic and spaghetti models.

We use numerial simulations to study the absorption and phase shift of surface-gravity waves caused by groups of magnetic flux tubes. The dependence of the scattering coefficients with the distance between the tubes and their positions is analyzed for several cases with two or three flux tubes embedded in a quiet Sun atmosphere. The results are compared with those obtained neglecting completely or partially multiple scattering effects. We show that multiple scattering has a significant impact on the absorption measurements and tends to reduce the phase shift. We also consider more general cases of ensembles of randomly distributed flux tubes, and we have evaluated the effects on the scattering measurements of changing the number of tubes included in the bundle and the average distance between flux tubes. We find that for the longest wavelength incoming waves multiple scattering enhances the absorption, and its efficiency increases with the number of flux tubes and the reduction of the distance between them.

We investigate least-squares fitting methods for estimating the winding rate of field lines about the axis of twisted magnetic-flux tubes. These methods estimate the winding rate by finding the values for a set of parameters that correspond to the minimum of the discrepancy between magnetic-field measurements and predictions from a twisted flux-tube model. For the flux-tube model used in the fitting, we assume that the magnetic field is static, axisymmetric, and does not vary in the vertical direction. Using error-free, synthetic vector magnetic-field data constructed with models for twisted magnetic-flux tubes, we test the efficacy of fitting methods at recovering the true winding rate. Furthermore, we demonstrate how assumptions built into the flux-tube models used for the fitting influence the accuracy of the winding-rate estimates. We identify the radial variation of the winding rate within the flux tube as one assumption that can have a significant impact on the winding-rate estimates. We show that the errors caused by making a fixed, incorrect assumption about the radial variation of the winding rate can be largely avoided by fitting directly for the radial variation of the winding rate. Other assumptions that we investigate include the lack of variation of the field in the azimuthal and vertical directions in the flux tube model used for the fitting, and the inclination, curvature, and location of the flux-tube axis. When the observed field deviates substantially from the flux-tube model used for the fitting, we find that the winding-rate estimates can be unreliable. We conclude that the magnetic-flux tube models used in this investigation are probably too simple to yield reliable estimates for the winding rate of the field lines in solar magnetic structures in general, unless additional information is available to justify the choice of flux-tube model used for the fitting.

We investigate how the azimuthal ambiguity in solar vector magnetogram data can be resolved by using the divergence-free property of magnetic fields. In a previous article, by Crouch, Barnes, and Leka (Solar Phys. 260, 271, 2009), error-free synthetic data were used to test several methods that each make a different assumption about how the divergence-free property can be used to resolve the ambiguity. In this paper this testing is continued with an examination of the effects of Poisson photon noise and limited instrumental spatial resolution. We find that all currently available methods based on the divergence-free property can produce undesirable results when photon noise or unresolved structure are present in the data. We perform a series of experiments aimed at improving the performance of the global minimisation method, which is the most promising of the methods. We present a two-step approach that produces reasonable results in tests using synthetic data. The first step of this approach involves the global minimisation of a combination of the absolute value of the approximation for the divergence and a smoothness constraint, which is designed to minimise the difference between the magnetic field in neighbouring pixels. In the second step, pixels with measurements known to be strongly affected by photon noise are revisited with a smoothing algorithm that also seeks to minimise the difference between the magnetic field in neighbouring pixels.

The scattering of f-modes by magnetic tubes is analyzed using three-dimensional numerical simulations. An f-mode wave packet is propagated through a solar atmosphere embedded with three different flux tube models which differ in radius and total magnetic flux. A quiet Sun simulation without a tube present is also performed as a reference. Waves are excited inside the flux tube and propagate along the field lines, and jacket modes are generated in the surroundings of the flux tube, carrying 40% as much energy as the tube modes. The resulting scattered wave is mainly an f-mode composed of a mixture of m=0 and m=+/-1 modes. The amplitude of the scattered wave approximately scales with the magnetic flux. A small amount of power is scattered into the p_1-mode. We have evaluated the absorption and phase shift from a Fourier-Hankel decomposition of the photospheric vertical velocities. They are compared with the results obtained from the emsemble average of 3400 small magnetic elements observed in high-resolution MDI Doppler datacubes. The comparison shows that the observed dependence of the phase shift with wavenumber can be matched reasonably well with the simulated flux tube model. The observed variation of the phase-shifts with the azimuthal order $m$ appears to depend on details of the ensemble averaging, including possible motions of the magnetic elements and asymmetrically shaped elements.

We address points recently discussed in Georgoulis (2011) in reference to Leka et al. (2009b). Most importantly, we find that the results of Georgoulis (2011) support a conclusion of Leka et al. (2009b): that limited spatial resolution and the presence of unresolved magnetic structures can challenge ambiguity- resolution algorithms. Moreover, the findings of both Metcalf et al. (2006) and Leka et al. (2009b) are confirmed in Georgoulis (2011): a method's performance can be diminished when the observed field fails to conform to that method's assumptions. The implication of boundaries in models of solar magnetic structures is discussed; we confirm that the distribution of the field components in the model used in Leka et al. (2009b) is closer to what is observed on the Sun than what is proposed in Georgoulis (2011). It is also shown that method does matter with regards to simulating limited spatial resolution and avoiding an inadvertent introduction of bias. Finally, the assignment of categories to data- analysis algorithms is revisited; we argue that assignments are only useful and elucidating when used appropriately.

We investigate how the divergence-free property of magnetic fields can be exploited to resolve the azimuthal ambiguity present in solar vector magnetogram data, by using line-of-sight and horizontal heliographic derivative information as approximated from discrete measurements. Using synthetic data we test several methods that each make different assumptions about how the divergence-free property can be used to resolve the ambiguity. We find that the most robust algorithm involves the minimisation of the absolute value of the divergence summed over the entire field of view. Away from disk centre this method requires the sign and magnitude of the line-of-sight derivatives of all three components of the magnetic field vector.

We present observations of microlensing event MACHO-98-BLG-35 which reached a peak magnification factor of almost 80. These observations by the Microlensing Planet Search (MPS) and the MOA Collaborations place strong constraints on the possible planetary system of the lens star and show intriguing evidence for a low mass planet with a mass fraction $4\times 10^{-5} \leq \epsilon \leq 2\times 10^{-4}$. A giant planet with $\epsilon = 10^{-3}$ is excluded from 95% of the region between 0.4 and 2.5 $R_E$ from the lens star, where $R_E$ is the Einstein ring radius of the lens. This exclusion region is more extensive than the generic "lensing zone" which is $0.6 - 1.6 R_E$. For smaller mass planets, we can exclude 57% of the "lensing zone" for $\epsilon = 10^{-4}$ and 14% of the lensing zone for $\epsilon = 10^{-5}$. The mass fraction $\epsilon = 10^{-5}$ corresponds to an Earth mass planet for a lensing star of mass $\sim 0.3 \msun$. A number of similar events will provide statistically significant constraints on the prevalence of Earth mass planets. In order to put our limits in more familiar terms, we have compared our results to those expected for a Solar System clone averaging over possible lens system distances and orientations. We find that such a system is ruled out at the 90% confidence level. A copy of the Solar System with Jupiter replaced by a second Saturn mass planet can be ruled out at 70% confidence. Our low mass planetary signal (few Earth masses to Neptune mass) is significant at the $4.5\sigma$ confidence level. If this planetary interpretation is correct, the MACHO-98-BLG-35 lens system constitutes the first detection of a low mass planet orbiting an ordinary star without gas giant planets.