Helioseismic Holography of an Artificial Submerged Sound Speed Perturbation and Implications for the Detection of Pre-Emergence Signatures of Active Regions

We use a publicly available numerical wave-propagation simulation of Hartlep et al. 2011 to test the ability of helioseismic holography to detect signatures of a compact, fully submerged, 5% sound-speed perturbation placed at a depth of 50 Mm within a solar model. We find that helioseismic holography as employed in a nominal "lateral-vantage" or "deep-focus" geometry employing quadrants of an annular pupil is capable of detecting and characterizing the perturbation. A number of tests of the methodology, including the use of a plane-parallel approximation, the definition of travel-time shifts, the use of different phase-speed filters, and changes to the pupils, are also performed. It is found that travel-time shifts made using Gabor-wavelet fitting are essentially identical to those derived from the phase of the Fourier transform of the cross-covariance functions. The errors in travel-time shifts caused by the plane-parallel approximation can be minimized to less than a second for the depths and fields of view considered here. Based on the measured strength of the mean travel-time signal of the perturbation, no substantial improvement in sensitivity is produced by varying the analysis procedure from the nominal methodology in conformance with expectations. The measured travel-time shifts are essentially unchanged by varying the profile of the phase-speed filter or omitting the filter entirely. The method remains maximally sensitive when applied with pupils that are wide quadrants, as opposed to narrower quadrants or with pupils composed of smaller arcs. We discuss the significance of these results for the recent controversy regarding suspected pre-emergence signatures of active regions

Helioseismology of Sunspots: A Case Study of NOAA Region 9787

Various methods of helioseismology are used to study the subsurface properties of the sunspot in NOAA Active Region 9787. This sunspot was chosen because it is axisymmetric, shows little evolution during 20-28 January 2002, and was observed continuously by the MDI/SOHO instrument. (...) Wave travel times and mode frequencies are affected by the sunspot. In most cases, wave packets that propagate through the sunspot have reduced travel times. At short travel distances, however, the sign of the travel-time shifts appears to depend sensitively on how the data are processed and, in particular, on filtering in frequency-wavenumber space. We carry out two linear inversions for wave speed: one using travel-times and phase-speed filters and the other one using mode frequencies from ring analysis. These two inversions give subsurface wave-speed profiles with opposite signs and different amplitudes. (...) From this study of AR9787, we conclude that we are currently unable to provide a unified description of the subsurface structure and dynamics of the sunspot.Comment: 28 pages, 18 figure

Predicting the Amplitude of a Solar Cycle Using the North-South Asymmetry in the Previous Cycle: II. An Improved Prediction for Solar Cycle~24

Recently, using Greenwich and Solar Optical Observing Network sunspot group data during the period 1874-2006, (Javaraiah, MNRAS, 377, L34, 2007: Paper I), has found that: (1) the sum of the areas of the sunspot groups in 0-10 deg latitude interval of the Sun's northern hemisphere and in the time-interval of -1.35 year to +2.15 year from the time of the preceding minimum of a solar cycle n correlates well (corr. coeff. r=0.947) with the amplitude (maximum of the smoothed monthly sunspot number) of the next cycle n+1. (2) The sum of the areas of the spot groups in 0-10 deg latitude interval of the southern hemisphere and in the time-interval of 1.0 year to 1.75 year just after the time of the maximum of the cycle n correlates very well (r=0.966) with the amplitude of cycle n+1. Using these relations, (1) and (2), the values 112 + or - 13 and 74 + or -10, respectively, were predicted in Paper I for the amplitude of the upcoming cycle 24. Here we found that in case of (1), the north-south asymmetry in the area sum of a cycle n also has a relationship, say (3), with the amplitude of cycle n+1, which is similar to (1) but more statistically significant (r=0.968) like (2). By using (3) it is possible to predict the amplitude of a cycle with a better accuracy by about 13 years in advance, and we get 103 + or -10 for the amplitude of the upcoming cycle 24. However, we found a similar but a more statistically significant (r=0.983) relationship, say (4), by using the sum of the area sum used in (2) and the north-south difference used in (3). By using (4) it is possible to predict the amplitude of a cycle by about 9 years in advance with a high accuracy and we get 87 + or - 7 for the amplitude of cycle 24.Comment: 21 pages, 7 figures, Published in Solar Physics 252, 419-439 (2008

Stokes Diagnostis of 2D MHD-simulated Solar Magnetogranulation

We study the properties of solar magnetic fields on scales less than the spatial resolution of solar telescopes. A synthetic infrared spectropolarimetric diagnostics based on a 2D MHD simulation of magnetoconvection is used for this. We analyze two time sequences of snapshots that likely represent two regions of the network fields with their immediate surrounding on the solar surface with the unsigned magnetic flux density of 300 and 140 G. In the first region we find from probability density functions of the magnetic field strength that the most probable field strength at logtau_5=0 is equal to 250 G. Weak fields (B < 500 G) occupy about 70% of the surface, while stronger fields (B 1000 G) occupy only 9.7% of the surface. The magnetic flux is -28 G and its imbalance is -0.04. In the second region, these parameters are correspondingly equal to 150 G, 93.3 %, 0.3 %, -40 G, and -0.10. We estimate the distribution of line-of-sight velocities on the surface of log tau_5=-1. The mean velocity is equal to 0.4 km/s in the first simulated region. The averaged velocity in the granules is -1.2 km/s and in the intergranules is 2.5 km/s. In the second region, the corresponding values of the mean velocities are equal to 0, -1.8, 1.5 km/s. In addition we analyze the asymmetry of synthetic Stokes-V profiles of the Fe I 1564.8 nm line. The mean values of the amplitude and area asymmetry do not exceed 1%. The spatially smoothed amplitude asymmetry is increased to 10% while the area asymmetry is only slightly varied.Comment: 24 pages, 12 figure

Local Helioseismology of Sunspots: Current Status and Perspectives (Invited Review)

Mechanisms of the formation and stability of sunspots are among the longest-standing and intriguing puzzles of solar physics and astrophysics. Sunspots are controlled by subsurface dynamics hidden from direct observations. Recently, substantial progress in our understanding of the physics of the turbulent magnetized plasma in strong-field regions has been made by using numerical simulations and local helioseismology. Both the simulations and helioseismic measurements are extremely challenging, but it becomes clear that the key to understanding the enigma of sunspots is a synergy between models and observations. Recent observations and radiative MHD numerical models have provided a convincing explanation to the Evershed flows in sunspot penumbrae. Also, they lead to the understanding of sunspots as self-organized magnetic structures in the turbulent plasma of the upper convection zone, which are maintained by a large-scale dynamics. Local helioseismic diagnostics of sunspots still have many uncertainties, some of which are discussed in this review. However, there have been significant achievements in resolving these uncertainties, verifying the basic results by new high-resolution observations, testing the helioseismic techniques by numerical simulations, and comparing results obtained by different methods. For instance, a recent analysis of helioseismology data from the Hinode space mission has successfully resolved several uncertainties and concerns (such as the inclined-field and phase-speed filtering effects) that might affect the inferences of the subsurface wave-speed structure of sunspots and the flow pattern. It becomes clear that for the understanding of the phenomenon of sunspots it is important to further improve the helioseismology methods and investigate the whole life cycle of active regions, from magnetic-flux emergence to dissipation.Comment: 34 pages, 18 figures, submitted to Solar Physic

Modeling the Subsurface Structure of Sunspots

While sunspots are easily observed at the solar surface, determining their subsurface structure is not trivial. There are two main hypotheses for the subsurface structure of sunspots: the monolithic model and the cluster model. Local helioseismology is the only means by which we can investigate subphotospheric structure. However, as current linear inversion techniques do not yet allow helioseismology to probe the internal structure with sufficient confidence to distinguish between the monolith and cluster models, the development of physically realistic sunspot models are a priority for helioseismologists. This is because they are not only important indicators of the variety of physical effects that may influence helioseismic inferences in active regions, but they also enable detailed assessments of the validity of helioseismic interpretations through numerical forward modeling. In this paper, we provide a critical review of the existing sunspot models and an overview of numerical methods employed to model wave propagation through model sunspots. We then carry out an helioseismic analysis of the sunspot in Active Region 9787 and address the serious inconsistencies uncovered by \citeauthor{gizonetal2009}~(\citeyear{gizonetal2009,gizonetal2009a}). We find that this sunspot is most probably associated with a shallow, positive wave-speed perturbation (unlike the traditional two-layer model) and that travel-time measurements are consistent with a horizontal outflow in the surrounding moat.Comment: 73 pages, 19 figures, accepted by Solar Physic

Stellar structure and compact objects before 1940: Towards relativistic astrophysics

Since the mid-1920s, different strands of research used stars as "physics laboratories" for investigating the nature of matter under extreme densities and pressures, impossible to realize on Earth. To trace this process this paper is following the evolution of the concept of a dense core in stars, which was important both for an understanding of stellar evolution and as a testing ground for the fast-evolving field of nuclear physics. In spite of the divide between physicists and astrophysicists, some key actors working in the cross-fertilized soil of overlapping but different scientific cultures formulated models and tentative theories that gradually evolved into more realistic and structured astrophysical objects. These investigations culminated in the first contact with general relativity in 1939, when J. Robert Oppenheimer and his students George Volkoff and Hartland Snyder systematically applied the theory to the dense core of a collapsing neutron star. This pioneering application of Einstein's theory to an astrophysical compact object can be regarded as a milestone in the path eventually leading to the emergence of relativistic astrophysics in the early 1960s.Comment: 83 pages, 4 figures, submitted to the European Physical Journal

Modelling spectral and timing properties of accreting black holes: the hybrid hot flow paradigm

The general picture that emerged by the end of 1990s from a large set of optical and X-ray, spectral and timing data was that the X-rays are produced in the innermost hot part of the accretion flow, while the optical/infrared (OIR) emission is mainly produced by the irradiated outer thin accretion disc. Recent multiwavelength observations of Galactic black hole transients show that the situation is not so simple. Fast variability in the OIR band, OIR excesses above the thermal emission and a complicated interplay between the X-ray and the OIR light curves imply that the OIR emitting region is much more compact. One of the popular hypotheses is that the jet contributes to the OIR emission and even is responsible for the bulk of the X-rays. However, this scenario is largely ad hoc and is in contradiction with many previously established facts. Alternatively, the hot accretion flow, known to be consistent with the X-ray spectral and timing data, is also a viable candidate to produce the OIR radiation. The hot-flow scenario naturally explains the power-law like OIR spectra, fast OIR variability and its complex relation to the X-rays if the hot flow contains non-thermal electrons (even in energetically negligible quantities), which are required by the presence of the MeV tail in Cyg X-1. The presence of non-thermal electrons also lowers the equilibrium electron temperature in the hot flow model to <100 keV, making it more consistent with observations. Here we argue that any viable model should simultaneously explain a large set of spectral and timing data and show that the hybrid (thermal/non-thermal) hot flow model satisfies most of the constraints.Comment: 26 pages, 13 figures. To be published in the Space Science Reviews and as hard cover in the Space Sciences Series of ISSI - The Physics of Accretion on to Black Holes (Springer Publisher

Recovering Joys Law as a Function of Solar Cycle, Hemisphere, and Longitude

Bipolar active regions in both hemispheres tend to be tilted with respect to the East West equator of the Sun in accordance with Joys law that describes the average tilt angle as a function of latitude. Mt. Wilson observatory data from 1917 to 1985 are used to analyze the active-region tilt angle as a function of solar cycle, hemisphere, and longitude, in addition to the more common dependence on latitude. Our main results are as follows: i) We recommend a revision of Joys law toward a weaker dependence on latitude (slope of 0.13 to 0.26) and without forcing the tilt to zero at the Equator. ii) We determine that the hemispheric mean tilt value of active regions varies with each solar cycle, although the noise from a stochastic process dominates and does not allow for a determination of the slope of Joys law on an 11-year time scale. iii) The hemispheric difference in mean tilt angles, 1.1 degrees + 0.27, over Cycles 16 to 21 was significant to a three-sigma level, with average tilt angles in the northern and southern hemispheres of 4.7 degrees + 0.26 and 3.6 degrees + 0.27 respectively. iv) Area-weighted mean tilt angles normalized by latitude for Cycles 15 to 21 anticorrelate with cycle strength for the southern hemisphere and whole-Sun data, confirming previous results by Dasi-Espuig, Solanki, Krivova, et al. (2010, Astron. Astrophys. 518, A7). The northern hemispheric mean tilt angles do not show a dependence on cycle strength. vi) Mean tilt angles do not show a dependence on longitude for any hemisphere or cycle. In addition, the standard deviation of the mean tilt is 29 to 31 degrees for all cycles and hemispheres indicating that the scatter is due to the same consistent process even if the mean tilt angles vary.Comment: 13 pages, 4 figures, 3 table

Search for displaced vertices arising from decays of new heavy particles in 7 TeV pp collisions at ATLAS

We present the results of a search for new, heavy particles that decay at a significant distance from their production point into a final state containing charged hadrons in association with a high-momentum muon. The search is conducted in a pp-collision data sample with a center-of-mass energy of 7 TeV and an integrated luminosity of 33 pb^-1 collected in 2010 by the ATLAS detector operating at the Large Hadron Collider. Production of such particles is expected in various scenarios of physics beyond the standard model. We observe no signal and place limits on the production cross-section of supersymmetric particles in an R-parity-violating scenario as a function of the neutralino lifetime. Limits are presented for different squark and neutralino masses, enabling extension of the limits to a variety of other models.Comment: 8 pages plus author list (20 pages total), 8 figures, 1 table, final version to appear in Physics Letters