Source author record

W. Schmidt

W. Schmidt appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

77works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

77 published item(s)

preprint2026arXiv

Water Production of Interstellar Comet 3I/ATLAS from SOHO/SWAN Observations after Perihelion

The Solar Wind ANisotropies (SWAN) all-sky hydrogen Lyman-alpha camera on the Solar and Heliosphere Observatory (SOHO) observed the hydrogen coma of interstellar comet 3I/ATLAS, also called C/2025 N1 (ATLAS), beginning on November 6, 2025, 9 days after perihelion. Water production rates were calculated from each image of 3I/ATLAS using the methodology of Makinen & Combi (2005, Icarus 177, 217) and fluorescence rates, g-factors, calculated using the daily solar Lyman-alpha fluxes from the LASP database (https://lasp.colorado.edu/lisird/data) corrected for solar rotation and for the comet heliocentric velocity. The method has been used for over 90 comet apparitions (Combi 2022; Combi et al. 2019). A water production rate of 3.17 x 10^29 s^-1 was found on November 6 when the comet was at a heliocentric distance of 1.40 au and at a sufficient solar elongation angle. It decreased over time after that, down to 1-2 x 10^28 s^-1 around 40 days post-perihelion (December 9).

preprint2022arXiv

The magnetic drivers of campfires seen by the Polarimetric and Helioseismic Imager (PHI) on Solar Orbiter

The Extreme Ultraviolet Imager (EUI) on board the Solar Orbiter (SO) spacecraft observed small extreme ultraviolet (EUV) bursts, termed campfires, that have been proposed to be brightenings near the apexes of low-lying loops in the quiet-Sun atmosphere. The underlying magnetic processes driving these campfires are not understood. During the cruise phase of SO and at a distance of 0.523\,AU from the Sun, the Polarimetric and Helioseismic Imager on Solar Orbiter (SO/PHI) observed a quiet-Sun region jointly with SO/EUI, offering the possibility to investigate the surface magnetic field dynamics underlying campfires at a spatial resolution of about 380~km. In 71\% of the 38 isolated events, campfires are confined between bipolar magnetic features, which seem to exhibit signatures of magnetic flux cancellation. The flux cancellation occurs either between the two main footpoints, or between one of the footpoints of the loop housing the campfire and a nearby opposite polarity patch. In one particularly clear-cut case, we detected the emergence of a small-scale magnetic loop in the internetwork followed soon afterwards by a campfire brightening adjacent to the location of the linear polarisation signal in the photosphere, that is to say near where the apex of the emerging loop lays. The rest of the events were observed over small scattered magnetic features, which could not be identified as magnetic footpoints of the campfire hosting loops. The majority of campfires could be driven by magnetic reconnection triggered at the footpoints, similar to the physical processes occurring in the burst-like EUV events discussed in the literature. About a quarter of all analysed campfires, however, are not associated to such magnetic activity in the photosphere, which implies that other heating mechanisms are energising these small-scale EUV brightenings.

preprint2022arXiv

The on-ground data reduction and calibration pipeline for SO/PHI-HRT

The ESA/NASA Solar Orbiter space mission has been successfully launched in February 2020. Onboard is the Polarimetric and Helioseismic Imager (SO/PHI), which has two telescopes, a High Resolution Telescope (HRT) and the Full Disc Telescope (FDT). The instrument is designed to infer the photospheric magnetic field and line-of-sight velocity through differential imaging of the polarised light emitted by the Sun. It calculates the full Stokes vector at 6 wavelength positions at the Fe I 617.3 nm absorption line. Due to telemetry constraints, the instrument nominally processes these Stokes profiles onboard, however when telemetry is available, the raw images are downlinked and reduced on ground. Here the architecture of the on-ground pipeline for HRT is presented, which also offers additional corrections not currently available on board the instrument. The pipeline can reduce raw images to the full Stokes vector with a polarimetric sensitivity of $10^{-3}\cdot I_{c}$ or better.

preprint2020arXiv

Neutron inelastic scattering study of rare-earth orthoferrite HoFeO$_3$

By the single crystal inelastic neutron scattering the orthoferrite HoFeO3 was studied. We show that the spin dynamics of the Fe subsystem does not change through the spin-reorientation transitions. The observed spectrum of magnetic excitations was analyzed in the frames of linear spin-wave theory. Within this approach the antiferromagnetic exchange interactions of nearest neighbors and next nearest neighbors were obtained for Fe subsystem. Parameters of Dzyaloshinskii-Moriya interactions at Fe subsystem were refined. The temperature dependence of the gap in Fe spin-wave spectrum indicates the temperature evolution of the anisotropy parameters. The estimations for the values of Fe-Ho and Ho-Ho exchange interaction were made as well.

preprint2020arXiv

Power spectrum of turbulent convection in the solar photosphere

The solar photosphere provides us with a laboratory for understanding turbulence in a layer where the fundamental processes of transport vary rapidly and a strongly superadiabatic region lies very closely to a subadiabatic layer. Our tools for probing the turbulence are high-resolution spectropolarimetric observations such as have recently been obtained with the two sunrise missions, and numerical simulations. Our aim is to study photospheric turbulence with the help of Fourier power spectra that we compute from observations and simulations. We also attempt to explain some properties of the photospheric overshooting flow with the help of its governing equations and simulations. We find that quiet-Sun observations and smeared simulations exhibit a power-law behavior in the subgranular range of their Doppler velocity power spectra with an index of$~\approx -2$. The unsmeared simulations exhibit a power-law index of$~\approx -2.25$. The smearing considerably reduces the extent of the power-law-like portion of the spectra. Therefore, the limited spatial resolution in some observations might eventually result in larger uncertainties in the estimation of the power-law indices. The simulated vertical velocity power spectra as a function of height show a rapid change in the power-law index from the solar surface to $300$~km above it. A scale-dependent transport of the vertical momentum occurs. At smaller scales, the vertical momentum is more efficiently transported sideways than at larger scales. This results in less vertical velocity power transported upward at small scales than at larger scales and produces a progressively steeper vertical velocity power law below $180$ km. Above this height, the gravity work progressively gains importance at all scales, making the atmosphere progressively more hydrostatic and resulting in a gradually less steep power law.

preprint2016arXiv

A nonlinear structural subgrid-scale closure for compressible MHD Part II: a priori comparison on turbulence simulation data

Even though compressible plasma turbulence is encountered in many astrophysical phenomena, its effect is often not well understood. Furthermore, direct numerical simulations are typically not able to reach the extreme parameters of these processes. For this reason, large-eddy simulations (LES), which only simulate large and intermediate scales directly, are employed. The smallest, unresolved scales and the interactions between small and large scales are introduced by means of a subgrid-scale (SGS) model. We propose and verify a new set of nonlinear SGS closures for future application as an SGS model in LES of compressible magnetohydrodynamics (MHD). We use 15 simulations (without explicit SGS model) of forced, isotropic, homogeneous turbulence with varying sonic Mach number $\mathrm{M_s} = 0.2$ to $20$ as reference data for the most extensive \textit{a priori} tests performed so far in literature. In these tests we explicitly filter the reference data and compare the performance of the new closures against the most widely tested closures. These include eddy-viscosity and scale-similarity type closures with different normalizations. Performance indicators are correlations with the turbulent energy and cross-helicty flux, the average SGS dissipation, the topological structure and the ability to reproduce the correct magnitude and direction of the SGS vectors. We find that only the new nonlinear closures exhibit consistently high correlations (median value \textgreater$0.8$) with the data over the entire parameter space and outperform the other closures in all tests. Moreover, we show that these results are independent of resolution and chosen filter scale. Additionally, the new closures are effectively coefficient-free with a deviation of less than $20\%$.

preprint2016arXiv

Active region fine structure observed at 0.08 arcsec resolution

The various mechanisms of magneto-convective energy transport determines the structure of sunspots and active regions. We characterise the appearance of light bridges and other fine-structure details and elaborate on their magneto-convective nature. We present speckle-reconstructed images taken with the broad-band imager at the 1.5 m GREGOR telescope in the 486nm and 589nm bands. We estimate the spatial resolution from the noise characteristics of the image bursts and obtain 0.08" at 589nm. We describe structure details in individual best images as well as the temporal evolution of selected features. We find branched dark lanes extending along thin (~1") light bridges in sunspots at various heliocentric angles. In thick (~2") light bridges the branches are disconnected from the central lane and have a Y shape with a bright grain toward the umbra. The images reveal that light bridges exist on varying intensity levels and that their small-scale features evolve on timescales of minutes. Faint light bridges show dark lanes outlined by the surrounding bright features. Dark lanes are very common and are also found in the boundary of pores. They have a characteristic width of 0.1" or smaller. Intergranular dark lanes of that width are seen in active region granulation. While central dark lanes in thin light bridges are elevated and associated with an density increase above upflows, the dark lane branches correspond to locations of downflows and are depressed relative to the adjacent bright plasma. Thick light bridges with central dark lanes show no projection effect. They have a flat elevated plateau that falls off steeply at the umbral boundary. There, Y-shaped filaments form as they do in the inner penumbra. This indicates the presence of inclined magnetic fields, meaning that the umbral magnetic field is wrapped around the convective light bridge.

preprint2016arXiv

Deep probing of the photospheric sunspot penumbra: no evidence for magnetic field-free gaps

Some models for the topology of the magnetic field in sunspot penumbrae predict the existence of field-free or dynamically weak-field regions in the deep Photosphere. To confirm or rule out the existence of weak-field regions in the deepest photospheric layers of the penumbra. The magnetic field at $\logτ_5=0$ is investigated by means of inversions of spectropolarimetric data of two different sunspots located very close to disk center with a spatial resolution of approximately 0.4-0.45 arcsec. The data have been recorded using the GRIS instrument attached to the 1.5-meters GREGOR solar telescope at El Teide observatory. It includes three Fe I lines around 1565 nm, whose sensitivity to the magnetic field peaks at half a pressure-scale-height deeper than the sensitivity of the widely used Fe I spectral line pair at 630 nm. Prior to the inversion, the data is corrected for the effects of scattered light using a deconvolution method with several point spread functions. At $\logτ_5=0$ we find no evidence for the existence of regions with dynamically weak ($B<500$ Gauss) magnetic fields in sunspot penumbrae. This result is much more reliable than previous investigations done with Fe I lines at 630 nm. Moreover, the result is independent of the number of nodes employed in the inversion, and also independent of the point spread function used to deconvolve the data, and does not depend on the amount of straylight (i.e. wide-angle scattered light) considered.

preprint2016arXiv

Electronic and magnetic nano phase separation in cobaltates La$_{2-x}$Sr$_{x}$CoO$_4$

The single-layer perovskite cobaltates have attracted enormous attention due to the recent observation of hour-glass shaped magnetic excitation spectra which resemble the ones of the famous high-temperature superconducting cuprates. Here, we present an overview of our most recent studies of the spin and charge correlations in floating-zone grown cobaltate single crystals. We find that frustration and a novel kind of electronic and magnetic nano phase separation are intimately connected to the appearance of the hour-glass shaped spin excitation spectra. We also point out the difference between nano phase separation and conventional phase separation.

preprint2016arXiv

Fitting peculiar spectral profiles in He I 10830 Å absorption features

The new generation of solar instruments provides better spectral, spatial, and temporal resolution for a better understanding of the physical processes that take place on the Sun. Multiple-component profiles are more commonly observed with these instruments. Particularly, the He I 10830 Å triplet presents such peculiar spectral profiles, which give information on the velocity and magnetic fine structure of the upper chromosphere. The purpose of this investigation is to describe a technique to efficiently fit the two blended components of the He I 10830 Å triplet, which are commonly observed when two atmospheric components are located within the same resolution element. The observations used in this study were taken on 2015 April 17 with the very fast spectroscopic mode of the GREGOR Infrared Spectrograph (GRIS) attached to the 1.5-meter GREGOR solar telescope, located at the Observatorio del Teide, Tenerife, Spain. We apply a double-Lorentzian fitting technique using Levenberg-Marquardt least-squares minimization. This technique is very simple and much faster than inversion codes. Line-of-sight Doppler velocities can be inferred for a whole map of pixels within just a few minutes. Our results show sub- and supersonic downflow velocities of up to 32 km/s for the fast component in the vicinity of footpoints of filamentary structures. The slow component presents velocities close to rest.

preprint2016arXiv

Flows in and around active region NOAA12118 observed with the GREGOR solar telescope and SDO/HMI

Accurate measurements of magnetic and velocity fields in and around solar active regions are key to unlocking the mysteries of the formation and the decay of sunspots. High spatial resolution image and spectral sequences with a high cadence obtained with the GREGOR solar telescope give us an opportunity to scrutinize 3-D flow fields with local correlation tracking and imaging spectroscopy. We present GREGOR early science data acquired in 2014 July - August with the GREGOR Fabry-Pérot Interferometer and the Blue Imaging Channel. Time-series of blue continuum (? 450.6 nm) images of the small active region NOAA 12118 were restored with the speckle masking technique to derive horizontal proper motions and to track the evolution of morphological changes. In addition, high-resolution observations are discussed in the context of synoptic data from the Solar Dynamics Observatory.

preprint2016arXiv

Horizontal flow fields in and around a small active region-- The transition period between flux emergence and decay

Aims. Combining high-resolution and synoptic observations aims to provide a comprehensive description of flux emergence at photospheric level and of the growth process that eventually leads to a mature active region. Methods. Small active region NOAA 12118 was observed on 2014 July 18 with the 1.5-meter GREGOR solar telescope on 2014 July 18. High-resolution time-series of blue continuum and G-band images acquired in the blue imaging channel (BIC) of the GREGOR Fabry-Pérot Interferometer (GFPI) were complemented by LOS magnetograms and continuum images obtained with the HMI onboard the SDO. Horizontal proper motions and horizontal plasma velocities were computed with local correlation tracking (LCT) and the differential affine velocity estimator, respectively. Morphological image processing was employed to measure the photometric/magnetic area, magnetic flux, and the separation profile of the EFR during its evolution. Results. The computed growth rates for photometric area, magnetic area, and magnetic flux are about twice as high as the respective decay rates. The space-time diagram using HMI magnetograms of five days traces a leaf-like structure, which is determined by the initial separation of the two polarities, a rapid expansion phase, a time when the spread stalls, and a period when the region slowly shrinks again. The separation rate of 0.26 km\s is highest in the initial stage, and it decreases when the separation comes to a halt. Horizontal plasma velocities computed at four evolutionary stages indicate a changing pattern of inflows. In LCT maps we find persistent flow patterns such as outward motions in the outer part of the two major pores, a diverging feature near the trailing pore marking the site of upwelling plasma and flux emergence, and low velocities in the interior of pores. We detected many elongated rapidly expanding granules between the two major polarities.

preprint2016arXiv

Incommensurate spin correlations in highly oxidized cobaltates La$_{2-x}$Sr$_{x}$CoO$_{4}$

We observe quasi-static incommensurate magnetic peaks in neutron scattering experiments on layered cobalt oxides La2-xSrxCoO4 with high Co oxidation states that have been reported to be paramagnetic. This enables us to measure the magnetic excitations in this highly hole-doped incommensurate regime and compare our results with those found in the low-doped incommensurate regime that exhibit hourglass magnetic spectra. The hourglass shape of magnetic excitations completely disappears given a high Sr doping. Moreover, broad low-energy excitations are found, which are not centered at the incommensurate magnetic peak positions but around the quarter-integer values that are typically exhibited by excitations in the checkerboard charge ordered phase. Our findings suggest that the strong inter-site exchange interactions in the undoped islands are critical for the emergence of hourglass spectra in the incommensurate magnetic phases of La2-xSrxCoO4.

preprint2016arXiv

Magnetic field reconstruction based on sunspot oscillations

The magnetic field of a sunspot guides magnetohydrodynamic waves toward higher atmospheric layers. In the upper photosphere and lower chromosphere, wave modes with periods longer than the acoustic cut-off period become evanescent. The cut-off period essentially changes due to the atmospheric properties, e.g., increases for larger zenith inclinations of the magnetic field. In this work, we aim at introducing a novel technique of reconstructing the magnetic field inclination on the basis of the dominating wave periods in the sunspot chromosphere and upper photosphere. On 2013 August 21st, we observed an isolated, circular sunspot (NOAA11823) for 58 min in a purely spectroscopic multi-wavelength mode with the Interferometric Bidimensional Spectro-polarimeter (IBIS) at the Dunn Solar Telescope. By means of a wavelet power analysis, we retrieved the dominating wave periods and reconstructed the zenith inclinations in the chromosphere and upper photosphere. The results are in good agreement with the lower photospheric HMI magnetograms. The sunspot's magnetic field in the chromosphere inclines from almost vertical (0$^{\circ}$) in the umbra to around 60$^{\circ}$ in the outer penumbra. With increasing altitude in the sunspot atmosphere, the magnetic field of the penumbra becomes less inclined. We conclude that the reconstruction of the magnetic field topology on the basis of sunspot oscillations yields consistent and conclusive results. The technique opens up a new possibility to infer the magnetic field inclination in the solar chromosphere.

preprint2016arXiv

Magnetic fields of opposite polarity in sunspot penumbrae

Context. A significant part of the penumbral magnetic field returns below the surface in the very deep photosphere. For lines in the visible, a large portion of this return field can only be detected indirectly by studying its imprints on strongly asymmetric and three-lobed Stokes V profiles. Infrared lines probe a narrow layer in the very deep photosphere, providing the possibility of directly measuring the orientation of magnetic fields close to the solar surface. Aims. We study the topology of the penumbral magnetic field in the lower photosphere, focusing on regions where it returns below the surface. Methods. We analyzed 71 spectropolarimetric datasets from Hinode and from the GREGOR infrared spectrograph. We inferred the quality and polarimetric accuracy of the infrared data after applying several reduction steps. Techniques of spectral inversion and forward synthesis were used to test the detection algorithm. We compared the morphology and the fractional penumbral area covered by reversed-polarity and three-lobed Stokes V profiles for sunspots at disk center. We determined the amount of reversed-polarity and three-lobed Stokes V profiles in visible and infrared data of sunspots at various heliocentric angles. From the results, we computed center-to-limb variation curves, which were interpreted in the context of existing penumbral models. Results. Observations in visible and near-infrared spectral lines yield a significant difference in the penumbral area covered by magnetic fields of opposite polarity. In the infrared, the number of reversed-polarity Stokes V profiles is smaller by a factor of two than in the visible. For three-lobed Stokes V profiles the numbers differ by up to an order of magnitude.

preprint2016arXiv

Probing deep photospheric layers of the quiet Sun with high magnetic sensitivity

Context. Investigations of the magnetism of the quiet Sun are hindered by extremely weak polarization signals in Fraunhofer spectral lines. Photon noise, straylight, and the systematically different sensitivity of the Zeeman effect to longitudinal and transversal magnetic fields result in controversial results in terms of the strength and angular distribution of the magnetic field vector. Aims. The information content of Stokes measurements close to the diffraction limit of the 1.5 m GREGOR telescope is analyzed. We took the effects of spatial straylight and photon noise into account. Methods. Highly sensitive full Stokes measurements of a quiet-Sun region at disk center in the deep photospheric Fe I lines in the 1.56 μm region were obtained with the infrared spectropolarimeter GRIS at the GREGOR telescope. Noise statistics and Stokes V asymmetries were analyzed and compared to a similar data set of the Hinode spectropolarimeter (SOT/SP). Simple diagnostics based directly on the shape and strength of the profiles were applied to the GRIS data. We made use of the magnetic line ratio technique, which was tested against MHD simulations. Results. About 80% of the GRIS spectra of a very quiet solar region show polarimetric signals above a 3σ level. Area and amplitude asymmetries agree well with small-scale surface dynamo MHD simulations. The magnetic line ratio analysis reveals ubiquitous magnetic regions in the ten to hundred Gauss range with some concentrations of kilo-Gauss fields. Conclusions. The GRIS spectropolarimetric data at a spatial resolution of 0.40" are so far unique in the combination of high spatial resolution scans and high magnetic field sensitivity. Nevertheless, the unavoidable effect of spatial straylight and the resulting dilution of the weak Stokes profiles means that inversion techniques still bear a high risk of misinterpretating the data.

preprint2016arXiv

Slipping reconnection in a solar flare observed in high resolution with the GREGOR solar telescope

A small flare ribbon above a sunspot umbra in active region 12205 was observed on November 7, 2014, at 12:00 UT in the blue imaging channel of the 1.5 m GREGOR telescope, using a 1 A Ca II H interference filter. Context observations from the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO), the Solar Optical Telescope (SOT) onboard Hinode, and the Interface Region Imaging Spectrograph (IRIS) show that this ribbon is part of a larger one that extends through the neighboring positive polarities and also participates in several other flares within the active region. We reconstructed a time series of 140 seconds of Ca II H images by means of the multiframe blind deconvolution method, which resulted in spatial and temporal resolutions of 0.1 arcsec and 1 s. Light curves and horizontal velocities of small-scale bright knots in the observed flare ribbon were measured. Some knots are stationary, but three move along the ribbon with speeds of 7-11 km/s. Two of them move in the opposite direction and exhibit highly correlated intensity changes, which provides evidence of a slipping reconnection at small spatial scales.

preprint2016arXiv

Spectropolarimetric observations of an arch filament system with the GREGOR solar telescope

Arch filament systems occur in active sunspot groups, where a fibril structure connects areas of opposite magnetic polarity, in contrast to active region filaments that follow the polarity inversion line. We used the GREGOR Infrared Spectrograph (GRIS) to obtain the full Stokes vector in the spectral lines Si I 1082.7 nm, He I 1083.0 nm, and Ca I 1083.9 nm. We focus on the near-infrared calcium line to investigate the photospheric magnetic field and velocities, and use the line core intensities and velocities of the helium line to study the chromospheric plasma. The individual fibrils of the arch filament system connect the sunspot with patches of magnetic polarity opposite to that of the spot. These patches do not necessarily coincide with pores, where the magnetic field is strongest. Instead, areas are preferred not far from the polarity inversion line. These areas exhibit photospheric downflows of moderate velocity, but significantly higher downflows of up to 30 km/s in the chromospheric helium line. Our findings can be explained with new emerging flux where the matter flows downward along the fieldlines of rising flux tubes, in agreement with earlier results.

preprint2016arXiv

Three-dimensional structure of a sunspot light bridge

Active regions are the most prominent manifestations of solar magnetic fields; their generation and dissipation are fundamental problems in solar physics. Light bridges are commonly present during sunspot decay, but a comprehensive picture of their role in the removal of photospheric magnetic field is still missing. We study the three dimensional configuration of a sunspot and in particular its light bridge during one of the last stages of its decay. We present the magnetic and thermodynamical stratification inferred from full Stokes inversions of the photospheric Si I 10827 Å and Ca I 10839 Å lines obtained with the GREGOR Infrared Spectrograph of the GREGOR telescope at Observatorio del Teide, Tenerife, Spain. The analysis is complemented by a study of continuum images covering the disk passage of the active region, which are provided by the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory. The sunspot shows a light bridge with penumbral continuum intensity that separates the central umbra from a smaller umbra. We find that in this region the magnetic field lines form a canopy with lower magnetic field strength in the inner part. The photospheric light bridge is dominated by gas pressure (high-$β$), as opposed to the surrounding umbra where the magnetic pressure is higher. A convective flow is observed in the light bridge. This flow is able to bend the magnetic field lines and to produce field reversals. The field lines close above the light bridge and become as vertical and strong as in the surrounding umbra. We conclude that it develops because of two highly magnetized regions which come closer during the sunspot evolution.

preprint2016arXiv

Upper Chromospheric Magnetic Field of a Sunspot Penumbra: Observations of Fine Structure

The fine-structure of magnetic field of a sunspot penumbra in the upper chromosphere is to be explored and compared to that in the photosphere. High spatial resolution spectropolarimetric observations were recorded with the 1.5-meter GREGOR telescope using the GREGOR Infrared Spectrograph (GRIS). The observed spectral domain includes the upper chromospheric He I triplet at 1083.0 nm and the photospheric Si I 1082.7 nm and Ca I 1083.3 nm spectral lines. The upper chromospheric magnetic field is obtained by inverting the He I triplet assuming a Milne-Eddington type model atmosphere. A height dependent inversion was applied to the Si I 1082.7 nm and Ca I 1083.3 nm lines to obtain the photospheric magnetic field. We find that the inclination of the magnetic field shows variations in the azimuthal direction both in the photosphere, but also in the upper chromosphere. The chromospheric variations remarkably well coincide with the variations in the inclination of the photospheric field and resemble the well-known spine and inter-spine structure in the photospheric layers of penumbrae. The typical peak-to-peak variations in the inclination of the magnetic field in the upper chromosphere is found to be 10-15 degree, i.e., roughly half the variation in the photosphere. In contrast, the magnetic field strength of the observed penumbra does not show variations on small spatial scales in the upper chromosphere. Thanks to the high spatial resolution observations possible with the GREGOR telescope at 1.08 microns, we find that the prominent small-scale fluctuations in the magnetic field inclination, which are a salient part of the property of sunspot penumbral photospheres, also persist in the chromosphere, although at somewhat reduced amplitudes. Such a complex magnetic configuration may facilitate penumbral chromospheric dynamic phenomena, such as penumbral micro-jets or transient bright dots.

preprint2015arXiv

Hour-glass magnetic excitations induced by nanoscopic phase separation in cobalt oxides La$_{2-x}$Sr$_x$CoO$_4$

The magnetic excitations in the cuprate superconductors might be essential for an understanding of high-temperature superconductivity. In these cuprate superconductors the magnetic excitation spectrum resembles an hour-glass and certain resonant magnetic excitations within are believed to be connected to the pairing mechanism which is corroborated by the observation of a universal linear scaling of superconducting gap and magnetic resonance energy. So far, charge stripes are widely believed to be involved in the physics of hour-glass spectra. Here we study an isostructural cobaltate that also exhibits an hour-glass magnetic spectrum. Instead of the expected charge stripe order we observe nano phase separation and unravel a microscopically split origin of hour-glass spectra on the nano scale pointing to a connection between the magnetic resonance peak and the spin gap originating in islands of the antiferromagnetic parent insulator. Our findings open new ways to theories of magnetic excitations and superconductivity in cuprate superconductors.

preprint2015arXiv

The photospheric solar oxygen project: III. Investigation of the centre-to-limb variation of the 630nm [OI]-NiI blend

The solar photospheric abundance of oxygen is still a matter of debate. For about ten years some determinations have favoured a low oxygen abundance which is at variance with the value inferred by helioseismology. Among the oxygen abundance indicators, the forbidden line at 630nm has often been considered the most reliable even though it is blended with a NiI line. In Papers I and Paper II of this series we reported a discrepancy in the oxygen abundance derived from the 630nm and the subordinate [OI] line at 636nm in dwarf stars, including the Sun. Here we analyse several, in part new, solar observations of the the centre-to-limb variation of the spectral region including the blend at 630nm in order to separate the individual contributions of oxygen and nickel. We analyse intensity spectra observed at different limb angles in comparison with line formation computations performed on a CO5BOLD 3D hydrodynamical simulation of the solar atmosphere. The oxygen abundances obtained from the forbidden line at different limb angles are inconsistent if the commonly adopted nickel abundance of 6.25 is assumed in our local thermodynamic equilibrium computations. With a slightly lower nickel abundance, A(Ni)~6.1, we obtain consistent fits indicating an oxygen abundance of A(O)=8.73+/-0.05. At this value the discrepancy with the subordinate oxygen line remains. The derived value of the oxygen abundance supports the notion of a rather low oxygen abundance in the solar hotosphere. However, it is disconcerting that the forbidden oxygen lines at 630 and 636nm give noticeably different results, and that the nickel abundance derived here from the 630nm blend is lower than expected from other nickel lines.

preprint2014arXiv

Cosmological Fluid Mechanics with Adaptively Refined Large Eddy Simulations

We investigate turbulence generated by cosmological structure formation by means of large eddy simulations using adaptive mesh refinement. In contrast to the widely used implicit large eddy simulations, which resolve a limited range of length scales and treat the effect of turbulent velocity fluctuations below the grid scale solely by numerical dissipation, we apply a subgrid-scale model for the numerically unresolved fraction of the turbulence energy. For simulations with adaptive mesh refinement, we utilize a new methodology that allows us to adjust the scale-dependent energy variables in such a way that the sum of resolved and unresolved energies is globally conserved. We test our approach in simulations of randomly forced turbulence, a gravitationally bound cloud in a wind, and the Santa Barbara cluster. To treat inhomogeneous turbulence, we introduce an adaptive Kalman filtering technique that separates turbulent velocity fluctuations on resolved length scales from the non-turbulent bulk flow. From the magnitude of the fluctuating component and the subgrid-scale turbulence energy, a total turbulent velocity dispersion of several 100 km/s is obtained for the Santa Barbara cluster, while the low-density gas outside the accretion shocks is nearly devoid of turbulence. The energy flux through the turbulent cascade and the dissipation rate predicted by the subgrid-scale model correspond to dynamical time scales around 5 Gyr, independent of numerical resolution.

preprint2014arXiv

Influence of adaptive mesh refinement and the hydro solver on shear-induced mass stripping in a minor-merger scenario

We compare two different codes for simulations of cosmological structure formation to investigate the sensitivity of hydrodynamical instabilities to numerics, in particular, the hydro solver and the application of adaptive mesh refinement (AMR). As a simple test problem, we consider an initially spherical gas cloud in a wind, which is an idealized model for the merger of a subcluster or galaxy with a big cluster. Based on an entropy criterion, we calculate the mass stripping from the subcluster as a function of time. Moreover, the turbulent velocity field is analyzed with a multi-scale filtering technique. We find remarkable differences between the commonly used PPM solver with directional splitting in the Enzo code and an unsplit variant of PPM in the Nyx code, which demonstrates that different codes can converge to systematically different solutions even when using uniform grids. For the test case of an unbound cloud, AMR simulations reproduce uniform-grid results for the mass stripping quite well, although the flow realizations can differ substantially. If the cloud is bound by a static gravitational potential, however, we find strong sensitivity to spurious fluctuations which are induced at the cutoff radius of the potential and amplified by the bow shock. This gives rise to substantial deviations between uniform-grid and AMR runs performed with Enzo, while the mass stripping in Nyx simulations of the subcluster is nearly independent of numerical resolution and AMR. Although many factors related to numerics are involved, our study indicates that unsplit solvers with advanced flux limiters help to reduce grid effects and to keep numerical noise under control, which is important for hydrodynamical instabilities and turbulent flows.

preprint2014arXiv

Magnetic fields during the formation of supermassive black holes

Observations of quasars at $\rm z> 6$ report the existence of a billion solar mass black holes. Comprehending their formation in such a short time scale is a matter of ongoing research. One of the most promising scenarios to assemble supermassive black holes is a monolithic collapse of protogalactic gas clouds in atomic cooling halos with $\rm T_{vir} \geq 10^{4} K$. In this article, we study the amplification and impact of magnetic fields during the formation of seed black holes in massive primordial halos. We perform high resolution cosmological magnetohydrodynamics simulations for four distinct halos and follow their collapse for a few free-fall times until the simulations reach a peak density of $\rm 7 \times 10^{-10} g/cm^{3}$. Our findings show that irrespective of the initial seed field, the magnetic field strength reaches a saturated state in the presence of strong accretion shocks. Under such conditions, the growth time becomes very short and amplification occurs rapidly within a small fraction of the free-fall time. We find that the presence of such strong magnetic fields provides additional support against gravity and helps in suppressing fragmentation. Massive clumps of a few hundred solar masses are formed at the end of our simulations and high accretion rates of $\rm 1 M_{\odot}/yr$ are observed. We expect that in the presence of such accretion rates, the clumps will grow to form supermassive stars of $\rm \sim 10^{5} M_{\odot}$. Overall, the role of the magnetic fields seems supportive for the formation of massive black holes.

preprint2014arXiv

Magnetic structure of the Eu2+ moments in superconducting EuFe2(As1-xPx)2 with x = 0.19

The magnetic structure of the Eu2+ moments in the superconducting EuFe2(As1-xPx)2 sample with x = 0.19 has been determined using neutron scattering. We conclude that the Eu2+ moments are aligned along the c direction below T_C = 19.0(1) K with an ordered moment of 6.6(2) mu_B in the superconducting state. An impurity phase similar to the underdoped phase exists within the bulk sample which orders antiferromagnetically below T_N = 17.0(2) K. We found no indication of iron magnetic order, nor any incommensurate magnetic order of the Eu2+ moments in the sample.

preprint2014arXiv

The formation and disintegration of magnetic bright points observed by Sunrise/IMaX

The evolution of the physical parameters of magnetic bright points (MBPs) located in the quiet Sun (mainly in the interwork) during their lifetime is studied. First we concentrate on the detailed description of the magnetic field evolution of three MBPs. This reveals that individual features follow different, generally complex, and rather dynamic scenarios of evolution. Next we apply statistical methods on roughly 200 observed MBP evolutionary tracks. MBPs are found to be formed by the strengthening of an equipartition field patch, which initially exhibits a moderate downflow. During the evolution, strong downdrafts with an average velocity of 2.4 km/s set in. These flows, taken together with the concurrent strengthening of the field, suggest that we are witnessing the occurrence of convective collapses in these features, although only 30% of them reach kG field strengths. This fraction might turn out to be larger when the new 4 m class solar telescopes are operational as observations of MBPs with current state of the art instrumentation could still be suffering from resolution limitations. Finally, when the bright point disappears (although the magnetic field often continues to exist) the magnetic field strength has dropped to the equipartition level and is generally somewhat weaker than at the beginning of the MBP's evolution. Noteworthy is that in about 10% of the cases we observe in the vicinity of the downflows small-scale strong (exceeding 2 km/s) intergranular upflows related spatially and temporally to these downflows.

preprint2013arXiv

Black hole formation in the early universe

Supermassive black holes with up to a $\rm 10^{9}~M_{\odot}$ dwell in the centers of present-day galaxies, and their presence has been confirmed at z $\geq$ 6. Their formation at such early epochs is still an enigma. Different pathways have been suggested to assemble supermassive black holes in the first billion years after the Big Bang. Direct collapse has emerged as a highly plausible scenario to form black holes as it provides seed masses of $\rm 10^{5}-10^{6}~M_{\odot}$. Gravitational collapse in atomic cooling haloes with virial temperatures T$_{vir} \geq 10^{4}$~K may lead to the formation of massive seed black holes in the presence of an intense background UV flux. Turbulence plays a central role in regulating accretion and transporting angular momentum. We present here the highest resolution cosmological large-eddy simulations to date which track the evolution of high-density regions on scales of $0.25$~AU beyond the formation of the first peak, and study the impact of subgrid-scale turbulence. The peak density reached in these simulations is $\rm 1.2 \times 10^{-8}~g~cm^{-3}$. Our findings show that while fragmentation occasionally occurs, it does not prevent the growth of a central massive object resulting from turbulent accretion and occasional mergers. The central object reaches $\rm \sim 1000~M_{\odot}$ within $4$ free-fall times, and we expect further growth up to $\rm 10^{6}~M_{\odot}$ through accretion in about 1 million years. The direct collapse model thus provides a viable pathway of forming high-mass black holes at early cosmic times.

preprint2013arXiv

First high-resolution images of the Sun in the 2796 Å Mg II k line

We present the first high-resolution solar images in the Mg II k 2796 Å line. The images, taken through a 4.8 Å broad interference filter, were obtained during the second science flight of SUNRISE in June 2013 by the SuFI instrument. The Mg II k images display structures that look qualitatively very similar to images taken in the core of Ca II H. The Mg II images exhibit reversed granulation (or shock waves) in the internetwork regions of the quiet Sun, at intensity contrasts that are similar to those found in Ca II H. Very prominent in Mg II are bright points, both in the quiet Sun and in plage regions, particularly near disk center. These are much brighter than at other wavelengths sampled at similar resolution. Furthermore, Mg II k images also show fibril structures associated with plage regions. Again, the fibrils are similar to those seen in Ca II H images, but tend to be more pronounced, particularly in weak plage.

preprint2013arXiv

Is magnetic reconnection the cause of supersonic upflows in granular cells ?

In a previous work, we reported on the discovery of supersonic magnetic upflows on granular cells in data from the {\sc Sunrise}/IMaX instrument. In the present work we investigate the physical origin of these events employing data of the same instrument but with higher spectral sampling. By means of the inversion of Stokes profiles we are able to recover the physical parameters (temperature, magnetic field, line-of-sight velocity, etc) present in the solar photosphere at the time of these events. The inversion is performed in a Monte-Carlo-like fashion, that is, repeating it many times with different initializations and retaining only the best result. We find that many of the events are characterized by a reversal in the polarity of the magnetic field along the vertical direction in the photosphere, accompanied by an enhancement in the temperature and by supersonic line-of-sight velocities. In about half of the studied events, large blue-shifted and red-shifted line-of-sight velocities coexist above/below each other. These features can be explained in terms of magnetic reconnection, where the energy stored in the magnetic field is released in the form of kinetic and thermal energy when magnetic field lines of opposite polarities coalesce. However, the agreement with magnetic reconnection is not perfect and therefore, other possible physical mechanisms might also play a role.

preprint2013arXiv

Local support against gravity in magneto-turbulent fluids

Comparisons of the integrated thermal pressure support of gas against its gravitational potential energy lead to critical mass scales for gravitational instability such as the Jeans and the Bonnor-Ebert masses, which play an important role in analysis of many physical systems, including the heuristics of numerical simulations. In a strict theoretical sense, however, neither the Jeans nor the Bonnor-Ebert mass are meaningful when applied locally to substructure in a self-gravitating turbulent medium. For this reason, we investigate the local support by thermal pressure, turbulence, and magnetic fields against gravitational compression through an approach that is independent of these concepts. At the centre of our approach is the dynamical equation for the divergence of the velocity field. We carry out a statistical analysis of the source terms of the local compression rate (the negative time derivative of the divergence) for simulations of forced self-gravitating turbulence in periodic boxes with zero, weak, and moderately strong mean magnetic fields (measured by the averages of the magnetic and thermal pressures). We also consider the amplification of the magnetic field energy by shear and by compression. Thereby, we are able to demonstrate that the support against gravity is dominated by thermal pressure fluctuations, although magnetic pressure also yields a significant contribution. The net effect of turbulence in the highly supersonic regime, however, is to enhance compression rather than supporting overdense gas even if the vorticity is very high. This is incommensurate with the support of the highly dynamical substructures in magneto-turbulent fluids being determined by local virial equilibria of volume energies without surface stresses.

preprint2013arXiv

Magnetic anisotropic energy gap and low energy spin wave excitation in antiferromagnetic block phase of K2Fe4Se5

Neutron scattering experiments were performed to investigate magnetic order and magnetic excitations in ternary iron chalcogenide K2Fe4Se5. The formation of a superlattice structure below 580 K together with the decoupling between the Fe-vacancy order-disorder transition and the antiferromagnetic order transition appears to be a common feature in the A2Fe4Se5 family. The study of spin dynamics of K2Fe4Se5 reveals two distinct energy gaps at the magnetic Brillouin zone center, which indicates the presence of magnetic anisotropy and the decrease of local symmetry due to electronic and orbital anisotropy. The low-energy spin wave excitations of K2Fe4Se5 can be properly described by linear spin wave theory within a Heisenberg model. Compared to iron pnictides, K2Fe4Se5 exhibits a more two-dimensional magnetism as characterized by large differences not only between out-of-plane and in-plane spin wave velocities, but also between out-of-plane and in-plane exchange interactions.

preprint2013arXiv

The characteristic black hole mass resulting from direct collapse in the early universe

Black holes of a billion solar masses are observed in the infant universe a few hundred million years after the Big Bang. The direct collapse of protogalactic gas clouds in primordial halos with $\rm T_{vir} \geq 10^{4} K$ provides the most promising way to assemble massive black holes. In this study, we aim to determine the characteristic mass scale of seed black holes and the time evolution of the accretion rates resulting from the direct collapse model. We explore the formation of supermassive black holes via cosmological large eddy simulations (LES) by employing sink particles and following their evolution for twenty thousand years after the formation of the first sink. As the resulting protostars were shown to have cool atmospheres in the presence of strong accretion, we assume here that UV feedback is negligible during this calculation. We confirm this result in a comparison run without sinks. Our findings show that black hole seeds with characteristic mass of $\rm 10^{5} M_{\odot}$ are formed in the presence of strong Lyman Werner flux which leads to an isothermal collapse. The characteristic mass is a about two times higher in LES compared to the implicit large eddy simulations (ILES). The accretion rates increase with time and reach a maximum value of 10 $\rm M_{\odot}/yr$ after $\rm 10^{4}$ years. Our results show that the direct collapse model is clearly feasible as it provides the expected mass of the seed black holes.

preprint2013arXiv

The formation of massive Pop III stars in the presence of turbulence

Population III stars forming in the infant universe at z=30 heralded the end of the cosmic dark ages. They are presumed to be assembled in so-called minihaloes with virial temperatures of a few thousand K where collapse is triggered by molecular hydrogen cooling. A central question concerns their final masses, and whether fragmentation occurs during their formation. While studies employing Lagrangian codes suggest fragmentation via a self-gravitating disk, recent high resolution simulations indicated that disk formation is suppressed. Here we report the first high-resolution large-eddy simulations performed with the Eulerian grid-based code Enzo following the evolution beyond the formation of the first peak, to investigate the accretion of the central massive clump and potential fragmentation. For a total of 3 halos, we see that a disk forms around the first clump. The central clump reaches $\sim10$ solar masses after 40 years, while subsequent accretion is expected at a rate of $10^{-2}$ solar masses per year. In one of these halos, additional clumps form as a result of fragmentation which proceeds at larger scales. We note that subgrid-scale turbulence yields relevant contributions to the stability of the protostellar disks. We conclude that the first protostar may reach masses up to $\rm 40-100 M_{\odot}$, which are only limited by the effect of radiative feedback.

preprint2013arXiv

The structure of the thermally bistable and turbulent atomic gas in the local interstellar medium

This paper is a numerical study of the condensation of the warm neutral medium (WNM) into cold neutral medium (CNM) structures under the effect of turbulence and thermal instability. Using low resolution simulations we explored the impact of the WNM initial density and properties of the turbulence (stirring in Fourier with a varying mix of solenoidal and compressive modes) on the cold gas formation. Two sets of initial conditions which match the observations were selected to produce high resolution simulations (1024^3) allowing to study in details the properties of the produced dense structures. For typical values of the density, pressure and velocity dispersion of the WNM in the solar neighborhood, the turbulent motions of the HI can not provoque the phase transition from WNM to CNM, whatever their amplitude and their distribution in solenoidal and compressive modes. On the other hand we show that a quasi-isothermal increase in WNM density of a factor of 2 to 4 is enough to induce the phase transition, leading to the transition of about 40 percent of the gas to the cold phase within 1 Myr. Given the observed properties of the HI in the local ISM, the WNM and individual CNM structures in the local ISM are sub or transsonic and their dynamics are tightly interwoven. The velocity field bears the evidence of subsonic turbulence with a 2D power spectrum following the Kolmogorov law as P(k) \sim k^{-8/3} while the density is highly contrasted with a singificantly shallower power spectrum, reminiscent of what is observed in the cold ISM. Supra-thermal line width observed for CNM might be the result of relative velocity between cold structures. Finally, the cold structures denser than 5 cm^{-3} reproduce well the laws M \sim L^{2.25-2.28} and sigma(v) \sim 0.5-0.8L^{1/3} generally observed in molecular clouds.

preprint2012arXiv

3D simulations of globules and pillars formation around HII regions: turbulence and shock curvature

We investigate the interplay between the ionization radiation from massive stars and the turbulence inside the surrounding molecular gas thanks to 3D numerical simulations. We used the 3D hydrodynamical code HERACLES to model an initial turbulent medium that is ionized and heated by an ionizing source. Three different simulations are performed with different mean Mach numbers (1, 2 and 4). A non-equilibrium model for the ionization and the associated thermal processes was used. This revealed to be crucial when turbulent ram pressure is of the same order as the ionized-gas pressure. The density structures initiated by the turbulence cause local curvatures of the dense shell formed by the ionization compression. When the curvature of the shell is sufficient, the shell collapse on itself to form a pillar while a smaller curvature leads to the formation of dense clumps that are accelerated with the shell and therefore remain in the shell during the simulation. When the turbulent ram pressure of the cold gas is sufficient to balance the ionized-gas pressure, some dense-gas bubbles have enough kinetic energy to penetrate inside the ionized medium, forming cometary globules. This suggests a direct relation in the observations between the presence of globules and the relative importance of the turbulence compared to the ionized-gas pressure. The probability density functions present a double peak structure when the turbulence is low relative to the ionized-gas pressure. This could be used in observations as an indication of the turbulence inside molecular clouds.

preprint2012arXiv

A retrospective of the GREGOR solar telescope in scientific literature

In this review, we look back upon the literature, which had the GREGOR solar telescope project as its subject including science cases, telescope subsystems, and post-focus instruments. The articles date back to the year 2000, when the initial concepts for a new solar telescope on Tenerife were first presented at scientific meetings. This comprehensive bibliography contains literature until the year 2012, i.e., the final stages of commissioning and science verification. Taking stock of the various publications in peer-reviewed journals and conference proceedings also provides the "historical" context for the reference articles in this special issue of Astronomische Nachrichten/Astronomical Notes.

preprint2012arXiv

High resolution studies of massive primordial haloes

Atomic cooling haloes with virial temperatures $\rm T_{vir} \geq 10^{4}$ K are the most plausible sites for the formation of the first galaxies and the first intermediate mass black holes. It is therefore important to assess whether one can obtain robust results concerning their main properties from numerical simulations. A major uncertainty is the presence of turbulence, which is barely resolved in cosmological simulations. We explore the latter both by pursuing high-resolution simulations with up to 64 cells per Jeans length and by incorporating a subgrid-scale turbulence model to account for turbulent pressure and viscosity on unresolved scales. We find that the main physical quantities in the halo, in particular the density, temperature and energy density profile, are approximately converged. However, the morphologies in the central 500 AU change significantly with increasing resolution and appear considerably more turbulent. In a systematic comparison of three different haloes, we further found that the turbulence subgrid-scale model gives rise to more compact central structures, and decreases the amount of vorticity. Such compact morphologies may in particular favor the accretion onto the central object.

preprint2012arXiv

Magnetically driven superconductivity in CeCu2Si2

The origin of unconventional superconductivity, including high-temperature and heavy-fermion superconductivity, is still a matter of controversy. Spin excitations instead of phonons are thought to be responsible for the formation of Cooper pairs. Using inelastic neutron scattering, we present the first in-depth study of the magnetic excitation spectrum in momentum and energy space in the superconducting and the normal states of CeCu2Si2. A clear spin excitation gap is observed in the superconducting state. We determine a lowering of the magnetic exchange energy in the superconducting state, in an amount considerably larger than the superconducting condensation energy. Our findings identify the antiferromagnetic excitations as the major driving force for superconducting pairing in this prototypical heavy-fermion compound located near an antiferromagnetic quantum critical point.

preprint2012arXiv

Resolving the internal magnetic structure of the solar network

We analyze the spectral asymmetry of Stokes V (circularly polarized) profiles of an individual network patch in the quiet Sun observed by Sunrise/IMaX. At a spatial resolution of 0.15"-0.18", the network elements contain substructure which is revealed by the spatial distribution of Stokes V asymmetries. The area asymmetry between the red and blue lobes of Stokes V increases from nearly zero at the core of the structure to values close to unity at its edges (one-lobed profiles). Such a distribution of the area asymmetry is consistent with magnetic fields expanding with height, i.e., an expanding magnetic canopy (which is required to fulfill pressure balance and flux conservation in the solar atmosphere). Inversion of the Stokes I and V profiles of the patch confirms this picture, revealing a decreasing field strength and increasing height of the canopy base from the core to the periphery of the network patch. However, the non-roundish shape of the structure and the presence of negative area and amplitude asymmetries reveal that the scenario is more complex than a canonical flux tube expanding with height surrounded by downflows.

preprint2012arXiv

Supersonic Magnetic Flows in the Quiet Sun

In this contribution we describe some recent observations of high-speed magnetized flows in the quiet Sun granulation. These observations were carried out with the Imaging Magnetograph eXperiment (IMaX) onboard the stratospheric balloon {\sc Sunrise}, and possess an unprecedented spatial resolution and temporal cadence. These flows were identified as highly shifted circular polarization (Stokes $V$) signals. We estimate the LOS velocity responsible for these shifts to be larger than 6 km s$^{-1}$, and therefore we refer to them as {\it supersonic magnetic flows}. The average lifetime of the detected events is 81.3 s and they occupy an average area of about 23\,000 km$^2$. Most of the events occur within granular cells and correspond therefore to upflows. However some others occur in intergranular lanes or bear no clear relation to the convective velocity pattern. We analyze a number of representative examples and discuss them in terms of magnetic loops, reconnection events, and convective collapse.

preprint2012arXiv

Variation in sunspot properties between 1999 and 2011 as observed with the Tenerife Infrared Polarimeter

We study the variation in the magnetic field strength and the umbral intensity of sunspots during the declining phase of the solar cycle no.23 and in the beginning of cycle no.24. We analyze a sample of 183 sunspots observed from 1999 until 2011 with the Tenerife Infrared Polarimeter at the German Vacuum Tower Telescope. The magnetic field strength is derived from the Zeeman splitting of the Stokes-V signal in one near-infrared spectral line, either Fe I 1564.8 nm, Fe I 1089.6 nm, or Si I 1082.7 nm. This avoids the effects of the unpolarized stray light from the field-free quiet Sun surroundings. The minimum umbral continuum intensity and umbral area are also measured. We find that there is a systematic trend for sunspots in the late stage of the solar cycle no.23 to be weaker, i.e., to have a smaller maximum magnetic field strength than those at the start of the cycle. The decrease in the field strength with time of about 94 G/yr is well beyond the statistical fluctuations that would be expected because of the larger number of sunspots close to cycle maximum (14 G/yr). In the same time interval, the continuum intensity of the umbra increases with a rate of 1.3 (+- 0.4)% of Ic/yr, while the umbral area does not show any trend above the statistical variance. Sunspots in the new cycle no.24 show higher field strengths and lower continuum intensities than those at the end of cycle no.23, interrupting the trend. Sunspots have an intrinsically weaker field strength and brighter umbrae at the late stages of solar cycles compared to their initial stages, without any significant change in their area. The abrupt increase in field strength in sunspots of the new cycle suggests that the cyclic variations are dominating over any long-term trend that continues across cycles. We find a slight decrease in field strength and an increase in intensity as a long-term trend across the cycles.

preprint2011arXiv

Magnetic field emergence in mesogranular-sized exploding granules observed with SUNRISE/IMaX data

We report on magnetic field emergences covering significant areas of exploding granules. The balloon-borne mission SUNRISE provided high spatial and temporal resolution images of the solar photosphere. Continuum images, longitudinal and transverse magnetic field maps and Dopplergrams obtained by IMaX onboard SUNRISE are analyzed by Local Correlation Traking (LCT), divergence calculation and time slices, Stokes inversions and numerical simulations are also employed. We characterize two mesogranular-scale exploding granules where $\sim$ 10$^{18}$ Mx of magnetic flux emerges. The emergence of weak unipolar longitudinal fields ($\sim$100 G) start with a single visible magnetic polarity, occupying their respective granules' top and following the granular splitting. After a while, mixed polarities start appearing, concentrated in downflow lanes. The events last around 20 min. LCT analyses confirm mesogranular scale expansion, displaying a similar pattern for all the physical properties, and divergence centers match between all of them. We found a similar behaviour with the emergence events in a numerical MHD simulation. Granule expansion velocities are around 1 \kms while magnetic patches expand at 0.65 \kms. One of the analyzed events evidences the emergence of a loop-like structure. Advection of the emerging magnetic flux features is dominated by convective motion resulting from the exploding granule due to the magnetic field frozen in the granular plasma. Intensification of the magnetic field occurs in the intergranular lanes, probably because of being directed by the downflowing plasma.

preprint2011arXiv

The frontier between Small-scale bipoles and Ephemeral Regions in the solar photosphere: Emergence and Decay of an Intermediate-scale bipole observed with IMaX/SUNRISE

We report on the photospheric evolution of an intermediate-scale (~4 Mm footpoint separation) magnetic bipole, from emergence to decay, observed in the quiet Sun at high spatial 0".3 and temporal (33 s) resolution. The observations were acquired by the IMaX imaging magnetograph during the first science flight of the Sunrise balloon-borne solar observatory. The bipole flux content is 6 x 10^17 Mx, representing a structure bridging the gap between granular scale bipoles and the smaller ephemeral regions. Footpoints separate at a speed of 3.5 km s-1 and reach a maximum distance of 4.5 Mm before the field dissolves. The evolution of the bipole is revealed to be very dynamic: we found a proper motion of the bipole axis and detected a change of the azimuth angle of 90° in 300 seconds. The overall morphology and behaviour are in agreement with previous analyses of bipolar structures emerging at granular scale, but we also found several similarities with emerging flux structures at larger scale. The flux growth rate is 2.6 x 15 Mx s-1, while the mean decay rate is one order of magnitude smaller. We describe in some detail the decay phase of the bipole footpoints which includes break up into smaller structures, interaction with pre-existing fields leading to cancellation but appears to be dominated by an as-yet unidentified diffusive process that removes most of the flux with an exponential flux decay curve. The diffusion constant (8 x 10^2 km^2 s-1) associated with this decay is similar to the values used to describe the large scale diffusion in flux transport models.

preprint2011arXiv

Turbulence Modelling and Stirring Mechanisms in the Cosmological Large-scale Structure

FEARLESS (Fluid mEchanics with Adaptively Refined Large Eddy SimulationS) is a numerical scheme for modelling subgrid-scale turbulence in cosmological adaptive mesh refinement simulations. In this contribution, the main features of this tool will be outlined. We discuss the application of this method to cosmological simulations of the large-scale structure. The simulations show that the production of turbulence has a different redshift dependence in the intra-cluster medium and the warm-hot intergalactic medium, caused by the distinct stirring mechanisms (mergers and shock interactions) acting in them. Some properties of the non-thermal pressure support in the two baryon phases are also described.

preprint2011arXiv

Turbulence production and turbulent pressure support in the intergalactic medium

The injection and evolution of turbulence in the intergalactic medium is studied by means of mesh-based hydrodynamical simulations, including a subgrid scale (SGS) model for small-scale unresolved turbulence. The simulations show that the production of turbulence has a different redshift dependence in the intracluster medium (ICM) and the warm-hot intergalactic medium (WHIM). We show that turbulence in the ICM is produced chiefly by merger-induced shear flows, whereas the production in the WHIM is dominated by shock interactions. Secondly, the effect of dynamical pressure support on the gravitational contraction has been studied. This turbulent support is stronger in the WHIM gas at baryon overdensities 1 < delta < 100, and less relevant for the ICM. Although the relative mass fraction of the gas with large vorticity is considerable (52% in the ICM), we find that for only about 10% in mass this is dynamically relevant, namely not associated to an equally large thermal pressure support. According to this result, a significant non-thermal pressure support counteracting the gravitational contraction is a localised characteristic in the cosmic flow, rather than a widespread feature.

preprint2011arXiv

Unnoticed magnetic field oscillations in the very quiet Sun revealed by Sunrise/IMaX

We present observational evidence for oscillations of magnetic flux density in the quiet areas of the Sun. The majority of magnetic fields on the solar surface have strengths of the order of or lower than the equipartition field (300-500 G). This results in a myriad of magnetic fields whose evolution is largely determined by the turbulent plasma motions. When granules evolve they squash the magnetic field lines together or pull them apart. Here we report on the periodic deformation of the shapes of features in circular polarization observed at high resolution with Sunrise. In particular, we note that the area of patches with constant magnetic flux oscillates with time, which implies that the apparent magnetic field intensity oscillates in antiphase. The periods associated to this oscillatory pattern is compatible with the granular life-time and change abruptly, which suggests that these oscillations might not correspond to characteristic oscillatory modes of magnetic structures, but to the forcing by granular motions. In one particular case, we find three patches around the same granule oscillating in phase, which means that the spatial coherence of these oscillations can reach 1600 km. Interestingly, the same kind of oscillatory phenomenon is found also in the upper photosphere.

preprint2011arXiv

Venus transit 2004: Illustrating the capability of exoplanet transmission spectroscopy

The transit of Venus in 2004 offered the rare possibility to remotely sense a well-known planetary atmosphere using ground-based observations for absorption spectroscopy. Transmission spectra of Venus' atmosphere were obtained in the near infrared using the Vacuum Tower Telescope (VTT) in Tenerife. Since the instrument was designed to measure the very bright photosphere of the Sun, extracting Venus' atmosphere was challenging. CO_2 absorption lines could be identified in the upper Venus atmosphere. Moreover, the relative abundance of the three most abundant CO_2 isotopologues could be determined. The observations resolved Venus' limb, showing Doppler-shifted absorption lines that are probably caused by high-altitude winds. This paper illustrates the ability of ground-based measurements to examine atmospheric constituents of a terrestrial planet atmosphere which might be applied in future to terrestrial extrasolar planets.

preprint2010arXiv

A Fluid-Dynamical Subgrid Scale Model for Highly Compressible Astrophysical Turbulence

We formulate and implement the Euler equations with SGS dynamics and provide numerical tests of an SGS turbulence energy model that predicts the turbulent pressure of unresolved velocity fluctuations and the rate of dissipation for highly compressible turbulence. We test closures for the turbulence energy cascade by filtering data from high-resolution simulations of forced isothermal and adiabatic turbulence. Optimal properties and an excellent correlation are found for a linear combination of the eddy-viscosity closure that is employed in LES of weakly compressible turbulence and a term that is non-linear in the Jacobian matrix of the velocity. Using this mixed closure, the SGS turbulence energy model is validated in LES of turbulence with stochastic forcing. It is found that the SGS model satisfies several important requirements: 1. The mean SGS turbulence energy follows a power law for varying grid scale. 2. The root mean square (RMS) Mach number of the unresolved velocity fluctuations is proportional to the RMS Mach number of the resolved turbulence, independent of the forcing. 3. The rate of dissipation and the turbulence energy flux are constant. Moreover, we discuss difficulties with direct estimates of the turbulent pressure and the dissipation rate on the basis of resolved flow quantities that have recently been proposed. In combination with the energy injection by stellar feedback and other unresolved processes, the proposed SGS model is applicable to a variety of problems in computational astrophysics. Computing the SGS turbulence energy, the treatment of star formation and stellar feedback in galaxy simulations can be improved. Further, we expect that the turbulent pressure on the grid scale affects the stability of gas against gravitational collapse.

preprint2010arXiv

Bright points in the quiet Sun as observed in the visible and near-UV by the balloon-borne observatory Sunrise

Bright points (BPs) are manifestations of small magnetic elements in the solar photosphere. Their brightness contrast not only gives insight into the thermal state of the photosphere (and chromosphere) in magnetic elements, but also plays an important role in modulating the solar total and spectral irradiance. Here we report on simultaneous high-resolution imaging and spectropolarimetric observations of BPs using Sunrise balloon-borne observatory data of the quiet Sun at disk center. BP contrasts have been measured between 214 nm and 525 nm, including the first measurements at wavelengths below 388 nm. The histograms of the BP peak brightness show a clear trend toward broader contrast distributions and higher mean contrasts at shorter wavelengths. At 214 nm we observe a peak brightness of up to five times the mean quiet-Sun value, the highest BP contrast so far observed. All BPs are associated with a magnetic signal, although in a number of cases it is surprisingly weak. Most of the BPs show only weak downflows, the mean value being 240 m/s, but some display strong down- or upflows reaching a few km/s.

preprint2010arXiv

Comparing the statistics of interstellar turbulence in simulations and observations: Solenoidal versus compressive turbulence forcing

We study two limiting cases of turbulence forcing in numerical experiments: solenoidal (divergence-free) forcing, and compressive (curl-free) forcing, and compare our results to observations reported in the literature. We solve the equations of hydrodynamics on grids with up to 1024^3 cells for purely solenoidal and purely compressive forcing. Eleven lower-resolution models with mixtures of both forcings are also analysed. We find velocity dispersion--size relations consistent with observations and independent numerical simulations, irrespective of the type of forcing. However, compressive forcing yields stronger turbulent compression at the same RMS Mach number than solenoidal forcing, resulting in a three times larger standard deviation of volumetric and column density probability distributions (PDFs). We conclude that the strong dependence of the density PDF on the type of forcing must be taken into account in any theory using the PDF to predict properties of star formation. We supply a quantitative description of this dependence. We find that different observed regions show evidence of different mixtures of compressive and solenoidal forcing, with more compressive forcing occurring primarily in swept-up shells.

preprint2010arXiv

Detection of large acoustic energy flux in the solar atmosphere

We study the energy flux carried by acoustic waves excited by convective motions at sub-photospheric levels. The analysis of high-resolution spectropolarimetric data taken with IMaX/Sunrise provides a total energy flux of ~ 6400--7700 Wm$^{-2}$ at a height of ~ 250 km in the 5.2-10 mHz range, i.e. at least twice the largest energy flux found in previous works. Our estimate lies within a factor of 2 of the energy flux needed to balance radiative losses from the chromosphere according to Anderson & Athay (1989) and revives interest in acoustic waves for transporting energy to the chromosphere. The acoustic flux is mainly found in the intergranular lanes but also in small rapidly-evolving granules and at the bright borders, forming dark dots and lanes of splitting granules.

preprint2010arXiv

Detection of vortex tubes in solar granulation from observations with Sunrise

We have investigated a time series of continuum intensity maps and corresponding Dopplergrams of granulation in a very quiet solar region at the disk center, recorded with the Imaging Magnetograph eXperiment (IMaX) on board the balloon-borne solar observatory Sunrise. We find that granules frequently show substructure in the form of lanes composed of a leading bright rim and a trailing dark edge, which move together from the boundary of a granule into the granule itself. We find strikingly similar events in synthesized intensity maps from an ab initio numerical simulation of solar surface convection. From cross sections through the computational domain of the simulation, we conclude that these `granular lanes' are the visible signature of (horizontally oriented) vortex tubes. The characteristic optical appearance of vortex tubes at the solar surface is explained. We propose that the observed vortex tubes may represent only the large-scale end of a hierarchy of vortex tubes existing near the solar surface.

preprint2010arXiv

Field induced spin reorientation and giant spin-lattice coupling in EuFe2As2

We have studied a EuFe2As2 single crystal by neutron diffraction under magnetic fields up to 3.5 T and temperatures down to 2 K. A field induced spin reorientation is observed in the presence of a magnetic field along both the a and c axes, respectively. Above critical field, the ground state antiferromagnetic configuration of Eu$^{2+}$ moments transforms into a ferromagnetic structure with moments along the applied field direction. The magnetic phase diagram for Eu magnetic sublattice in EuFe2As2 is presented. A considerable strain ($\sim$0.9%) is induced by the magnetic field, caused by the realignment of the twinning structure. Furthermore, the realignment of the twinning structure is found to be reversible with the rebound of magnetic field, which suggested the existence of magnetic shape-memory effect. The Eu moment ordering exhibits close relationship with the twinning structure. We argue that the Zeeman energy in combined with magnetic anisotropy energy is responsible for the observed spin-lattice coupling.

preprint2010arXiv

Forced turbulence in thermally bistable gas: A parameter study

Context: The thermal instability is one of the dynamical agents for turbulence in the diffuse interstellar medium, where both, turbulence and thermal instability interact in a highly non-linear manner. Aims: We study basic properties of turbulence in thermally bistable gas for variable simulation parameters. The resulting cold gas fractions can be applied as parameterisation in simulations on galactic scales. Methods: Turbulent flow is induced on large scales by means of compressive stochastic forcing in a periodic box. The compressible Euler equations with constant UV heating and a parameterised cooling function are solved on uniform grids. We investigate several values of the mean density of the gas and different magnitudes of the forcing. For comparison with other numerical studies, solenoidal forcing is applied as well. Results: After a transient phase, we observe that a state of statistically stationary turbulence is approached. Compressive forcing generally produces a two-phase medium, with a decreasing fraction of cold gas for increasing forcing strength. This behaviour can be explained on the basis of turbulent mixing. We also find indications for power-law tails of probability density functions of the gas density. Solenoidal forcing, on the other hand, appears to prevent the evolution into a two-phase-medium for certain parameter regions. Conclusions: The dynamics of thermally bistable turbulence shows a substantial sensitivity on the initial state and the forcing properties.

preprint2010arXiv

Fully resolved quiet-Sun magnetic flux tube observed with the Sunrise IMaX instrument

Until today, the small size of magnetic elements in quiet Sun areas has required the application of indirect methods, such as the line-ratio technique or multi-component inversions, to infer their physical properties. A consistent match to the observed Stokes profiles could only be obtained by introducing a magnetic filling factor that specifies the fraction of the observed pixel filled with magnetic field. Here, we investigate the properties of a small magnetic patch in the quiet Sun observed with the IMaX magnetograph on board the balloon-borne telescope Sunrise with unprecedented spatial resolution and low instrumental stray light. We apply an inversion technique based on the numerical solution of the radiative transfer equation to retrieve the temperature stratification and the field strength in the magnetic patch. The observations can be well reproduced with a one-component, fully magnetized atmosphere with a field strength exceeding 1 kG and a significantly enhanced temperature in the mid- to upper photosphere with respect to its surroundings, consistent with semi-empirical flux tube models for plage regions. We therefore conclude that, within the framework of a simple atmospheric model, the IMaX measurements resolve the observed quiet-Sun flux tube.

preprint2010arXiv

Mesogranulation and the solar surface magnetic field distribution

The relation of the solar surface magnetic field with mesogranular cells is studied using high spatial (~ 100 km) and temporal (~ 30 sec) resolution data obtained with the IMaX instrument aboard SUNRISE. First, mesogranular cells are identified using Lagrange tracers (corks) based on horizontal velocity fields obtained through Local Correlation Tracking. After ~ 20 min of integration, the tracers delineate a sharp mesogranular network with lanes of width below about 280 km. The preferential location of magnetic elements in mesogranular cells is tested quantitatively. Roughly 85% of pixels with magnetic field higher than 100 G are located in the near neighborhood of mesogranular lanes. Magnetic flux is therefore concentrated in mesogranular lanes rather than intergranular ones. Secondly, magnetic field extrapolations are performed to obtain field lines anchored in the observed flux elements. This analysis, therefore, is independent of the horizontal flows determined in the first part. A probability density function (PDF) is calculated for the distribution of distances between the footpoints of individual magnetic field lines. The PDF has an exponential shape at scales between 1 and 10 Mm, with a constant characteristic decay distance, indicating the absence of preferred convection scales in the mesogranular range. Our results support the view that mesogranulation is not an intrinsic convective scale (in the sense that it is not a primary energy-injection scale of solar convection), but also give quantitative confirmation that, nevertheless, the magnetic elements are preferentially found along mesogranular lanes.

preprint2010arXiv

Sunrise: instrument, mission, data and first results

The Sunrise balloon-borne solar observatory consists of a 1m aperture Gregory telescope, a UV filter imager, an imaging vector polarimeter, an image stabilization system and further infrastructure. The first science flight of Sunrise yielded high-quality data that reveal the structure, dynamics and evolution of solar convection, oscillations and magnetic fields at a resolution of around 100 km in the quiet Sun. After a brief description of instruments and data, first qualitative results are presented. In contrast to earlier observations, we clearly see granulation at 214 nm. Images in Ca II H display narrow, short-lived dark intergranular lanes between the bright edges of granules. The very small-scale, mixed-polarity internetwork fields are found to be highly dynamic. A significant increase in detectable magnetic flux is found after phase-diversity-related reconstruction of polarization maps, indicating that the polarities are mixed right down to the spatial resolution limit, and probably beyond.

preprint2010arXiv

SUNRISE/IMaX observations of convectively driven vortex flows in the Sun

We characterize the observational properties of the convectively driven vortex flows recently discovered on the quiet Sun, using magnetograms, Dopplergrams and images obtained with the 1-m balloon-borne Sunrise telescope. By visual inspection of time series, we find some 3.1e-3 vortices/(Mm^2 min), which is a factor of 1.7 larger than previous estimates. The mean duration of the individual events turns out to be 7.9 min, with a standard deviation of 3.2 min. In addition, we find several events appearing at the same locations along the duration of the time series (31.6 min). Such recurrent vortices show up in the proper motion flow field map averaged over the time series. The typical vertical vorticities are <= 6e-3 1/sec, which corresponds to a period of rotation of some 35 min. The vortices show a preferred counterclockwise sense of rotation, which we conjecture may have to do with the preferred vorticity impinged by the solar differential rotation.

preprint2010arXiv

Supersonic Magnetic Upflows in Granular Cells Observed with Sunrise/IMaX

Using the IMaX instrument on-board the Sunrise stratospheric balloon-telescope we have detected extremely shifted polarization signals around the Fe I 5250.217 Å spectral line within granules in the solar photosphere. We interpret the velocities associated with these events as corresponding to supersonic and magnetic upflows. In addition, they are also related to the appearance of opposite polarities and highly inclined magnetic fields. This suggests that they are produced by the reconnection of emerging magnetic loops through granular upflows. The events occupy an average area of 0.046 arcsec$^2$ and last for about 80 seconds, with larger events having longer lifetimes. These supersonic events occur at a rate of $1.3\times10^{-5}$ occurrences per second per arcsec$^{2}$.

preprint2010arXiv

Surface waves in solar granulation observed with {\sc Sunrise}

Solar oscillations are expected to be excited by turbulent flows in the intergranular lanes near the solar surface. Time series recorded by the IMaX instrument aboard the {\sc Sunrise} observatory reveal solar oscillations at high resolution, which allow studying the properties of oscillations with short wavelengths. We analyze two times series with synchronous recordings of Doppler velocity and continuum intensity images with durations of 32\thinspace min and 23\thinspace min, resp., recorded close to the disk center of the Sun to study the propagation and excitation of solar acoustic oscillations. In the Doppler velocity data, both the standing acoustic waves and the short-lived, high-degree running waves are visible. The standing waves are visible as temporary enhancements of the amplitudes of the large-scale velocity field due to the stochastic superposition of the acoustic waves. We focus on the high-degree small-scale waves by suitable filtering in the Fourier domain. Investigating the propagation and excitation of $f$- and $p_1$-modes with wave numbers $k > 1.4$\thinspace 1/Mm we find that also exploding granules contribute to the excitation of solar $p$-modes in addition to the contribution of intergranular lanes.

preprint2010arXiv

The Filter Imager SuFI and the Image Stabilization and Light Distribution System ISLiD of the Sunrise Balloon-Borne Observatory: Instrument Description

We describe the design of the Sunrise Filter Imager (SuFI) and the Image Stabilization and Light Distribution (ISLiD) unit onboard the Sunrise balloon borne solar observatory. This contribution provides the necessary information which is relevant to understand the instruments working principles, the relevant technical data, and the necessary information about calibration issues directly related to the science data.

preprint2010arXiv

The Imaging Magnetograph eXperiment (IMaX) for the Sunrise balloon-borne solar observatory

The Imaging Magnetograph eXperiment (IMaX) is a spectropolarimeter built by four institutions in Spain that flew on board the Sunrise balloon-borne telesocope in June 2009 for almost six days over the Arctic Circle. As a polarimeter IMaX uses fast polarization modulation (based on the use of two liquid crystal retarders), real-time image accumulation, and dual beam polarimetry to reach polarization sensitivities of 0.1%. As a spectrograph, the instrument uses a LiNbO3 etalon in double pass and a narrow band pre-filter to achieve a spectral resolution of 85 mAA. IMaX uses the high Zeeman sensitive line of Fe I at 5250.2 AA and observes all four Stokes parameters at various points inside the spectral line. This allows vector magnetograms, Dopplergrams, and intensity frames to be produced that, after reconstruction, reach spatial resolutions in the 0.15-0.18 arcsec range over a 50x50 arcsec FOV. Time cadences vary between ten and 33 seconds, although the shortest one only includes longitudinal polarimetry. The spectral line is sampled in various ways depending on the applied observing mode, from just two points inside the line to 11 of them. All observing modes include one extra wavelength point in the nearby continuum. Gauss equivalent sensitivities are four Gauss for longitudinal fields and 80 Gauss for transverse fields per wavelength sample. The LOS velocities are estimated with statistical errors of the order of 5-40 m/s. The design, calibration and integration phases of the instrument, together with the implemented data reduction scheme are described in some detail.

preprint2010arXiv

The Wave-Front Correction System for the Sunrise Balloon-Borne Solar Observatory

This paper describes the wave-front correction system developed for the Sunrise balloon telescope, and provides information about its in-flight performance. For the correction of low-order aberrations, a Correlating Wave-Front Sensor (CWS) was used. It consisted of a six-element Shack-Hartmann wave-front sensor (WFS), a fast tip-tilt mirror for the compensation of image motion, and an active telescope secondary mirror for focus correction. The CWS delivered a stabilized image with a precision of 0.04 arcsec (rms), whenever the coarse pointing was better than 90 arcsec peak-to-peak. The automatic focus adjustment maintained a focus stability of 0.01 waves in the focal plane of the CWS. During the 5.5 day flight, good image quality and stability was achieved during 33 hours, containing 45 sequences that lasted between 10 and 45 minutes.

preprint2010arXiv

Transverse component of the magnetic field in the solar photosphere observed by Sunrise

We present the first observations of the transverse component of photospheric magnetic field acquired by the imaging magnetograph Sunrise/IMaX. Using an automated detection method, we obtain statistical properties of 4536 features with significant linear polarization signal. Their rate of occurrence is 1-2 orders of magnitude larger than values reported by previous studies. We show that these features have no characteristic size or lifetime. They appear preferentially at granule boundaries with most of them being caught in downflow lanes at some point in their evolution. Only a small percentage are entirely and constantly embedded in upflows (16%) or downflows (8%).

preprint2010arXiv

Where the granular flows bend

Based on IMaX/Sunrise data, we report on a previously undetected phenomenon in solar granulation. We show that in a very narrow region separating granules and intergranular lanes the spectral line width of the Fe I 5250.2 A line becomes extremely small. We offer an explanation of this observation with the help of magneto-convection simulations. These regions with extremely small line widths correspond to the places where the granular flows bend from mainly upflow in granules to downflow in intergranular lanes. We show that the resolution and image stability achieved by IMaX/Sunrise are important requisites to detect this interesting phenomenon.

preprint2009arXiv

Adaptively refined large eddy simulations of clusters

We present a numerical scheme for modelling unresolved turbulence in cosmological adaptive mesh refinement codes. As a first application, we study the evolution of turbulence in the intra-cluster medium and in the core of a galaxy cluster. Simulations with and without subgrid scale model are compared in detail. Since the flow in the ICM is subsonic, the global turbulent energy contribution at the unresolved length scales is smaller than 1% of the internal energy. We find that the production of turbulence is closely correlated with merger events occurring in the cluster environment, and its dissipation locally affects the cluster energy budget. Because of this additional source of dissipation, the core temperature is larger and the density is smaller in the presence of subgrid scale turbulence than in the standard adiabatic run, resulting in a higher entropy core value.

preprint2009arXiv

Algorithmic comparisons of decaying, isothermal, supersonic turbulence

Contradicting results have been reported in the literature with respect to the performance of the numerical techniques employed for the study of supersonic turbulence. We aim at characterising the performance of different particle-based and grid-based techniques on the modelling of decaying supersonic turbulence. Four different grid codes (ENZO, FLASH, TVD, ZEUS) and three different SPH codes (GADGET, PHANTOM, VINE) are compared. We additionally analysed two calculations denoted as PHANTOM A and PHANTOM B using two different implementations of artificial viscosity. Our analysis indicates that grid codes tend to be less dissipative than SPH codes, though details of the techniques used can make large differences in both cases. For example, the Morris & Monaghan viscosity implementation for SPH results in less dissipation (PHANTOM B and VINE versus GADGET and PHANTOM A). For grid codes, using a smaller diffusion parameter leads to less dissipation, but results in a larger bottleneck effect (our ENZO versus FLASH runs). As a general result, we find that by using a similar number of resolution elements N for each spatial direction means that all codes (both grid-based and particle-based) show encouraging similarity of all statistical quantities for isotropic supersonic turbulence on spatial scales k<N/32 (all scales resolved by more than 32 grid cells), while scales smaller than that are significantly affected by the specific implementation of the algorithm for solving the equations of hydrodynamics. At comparable numerical resolution, the SPH runs were on average about ten times more computationally intensive than the grid runs, although with variations of up to a factor of ten between the different SPH runs and between the different grid runs. (abridged)

preprint2009arXiv

Highly Dispersive Spin Excitations in the Chain Cuprate Li2CuO2

We present an inelastic neutron scattering investigation of Li2CuO2 detecting the long sought quasi-1D magnetic excitations with a large dispersion along the CuO2-chains studied up to 25 meV. The total dispersion is governed by a surprisingly large ferromagnetic (FM) nearest-neighbor exchange integral J1=-228 K. An anomalous quartic dispersion near the zone center and a pronounced minimum near (0,0.11,0.5) r.l.u. (corresponding to a spiral excitation with a pitch angle about 41 degree point to the vicinity of a 3D FM-spiral critical point. The leading exchange couplings are obtained applying standard linear spin-wave theory. The 2nd neighbor inter-chain interaction suppresses a spiral state and drives the FM in-chain ordering below the Ne'el temperature. The obtained exchange parameters are in agreement with the results for a realistic five-band extended Hubbard Cu 3d O 2p model and L(S)DA+U predictions.

preprint2009arXiv

Solenoidal versus compressive turbulence forcing

We analyze the statistics and star formation rate obtained in high-resolution numerical experiments of forced supersonic turbulence, and compare with observations. We concentrate on a systematic comparison of solenoidal (divergence-free) and compressive (curl-free) forcing, which are two limiting cases of turbulence driving. Our results show that for the same RMS Mach number, compressive forcing produces a three times larger standard deviation of the density probability distribution. When self-gravity is included in the models, the star formation rate is more than one order of magnitude higher for compressive forcing than for solenoidal forcing.

preprint2009arXiv

Thermonuclear explosions of rapidly rotating white dwarfs - I. Deflagrations

Context: Turbulent deflagrations of Chandrasekhar mass White Dwarfs are commonly used to model Type Ia Supernova explosions. In this context, rapid rotation of the progenitor star is plausible but has so far been neglected. Aims: The aim of this work is to explore the influence of rapid rotation on the deflagration scenario. Methods: We use three dimensional hydrodynamical simulations to model turbulent deflagrations ignited within a variety of rapidly rotating CO WDs obeying rotation laws suggested by accretion studies. Results: We find that rotation has a significant impact on the explosion. The flame develops a strong anisotropy with a preferred direction towards the stellar poles, leaving great amounts of unburnt matter along the equatorial plane. Conclusions: The large amount of unburnt matter is contrary to observed spectral features of SNe Ia. Thus, rapid rotation of the progenitor star and the deflagration scenario are incompatible in order to explain SNe Ia.

preprint2009arXiv

Thermonuclear explosions of rapidly rotating white dwarfs - II. Detonations

Context: Superluminous type Ia supernovae (SNe Ia) may be explained by super-Chandrasekhar-mass explosions of rapidly rotating white dwarfs (WDs). In a preceding paper, we showed that the deflagration scenario applied to rapidly rotating WDs generates explosions that cannot explain the majority of SNe Ia. Aims: Rotation of the progenitor star allows super-Chandrasekhar-mass WDs to form that have a shallower density stratification. We use simple estimates of the production of intermediate and iron group elements in pure detonations of rapidly rotating WDs to assess their viability in explaining rare SNe Ia. Methods: We numerically construct WDs in hydrostatic equilibrium that rotate according to a variety of rotation laws. The explosion products are estimated by considering the density stratification and by evaluating the result of hydrodynamics simulations. Results: We show that a significant amount of intermediate mass elements is produced for theoretically motivated rotation laws, even for prompt detonations of WDs. Conclusions: Rapidly rotating WDs that detonate may provide an explanation of rare superluminous SNe Ia in terms of both burning species and explosion kinematics.

preprint2009arXiv

Turbulence in a 3D deflagration model for type Ia SNe: II. Intermittency and the deflagration-to-detonation transition probability

The delayed detonation model describes the observational properties of the majority of type Ia supernovae very well. Using numerical data from a three-dimensional deflagration model for type Ia supernovae, the intermittency of the turbulent velocity field and its implications on the probability of a deflagration-to-detonation (DDT) transition are investigated. From structure functions of the turbulent velocity fluctuations, we determine intermittency parameters based on the log-normal and the log-Poisson models. On the other hand, the analysis of the turbulent velocity fluctuations in the vicinity of the flame front by Roepke suggests a much higher probability of large velocity fluctuations on the grid scale in comparison to the log-normal intermittency model. Following Pan et al., we computed probability density functions for a DDT for the different distributions. Assuming that a DDT can occur in the stirred flame regime, as proposed by Woosley et al., the log-normal model would imply a delayed detonation between 0.7 and 0.8 seconds after the beginning of the deflagration phase for the multi-spot ignition scenario used in the simulation. However, the probability drops to virtually zero if a DDT is further constrained by the requirement that the turbulent velocity fluctuations reach about 500 km/s. Under this condition, delayed detonations are only possible if the distribution of the velocity fluctuations is not log-normal. From our calculations follows that the distribution obtained by Roepke allow for multiple DDTs around 0.8 seconds after ignition at a transition density close to 1x10^7 g/cm^3.

preprint2009arXiv

Turbulence in a three-dimensional deflagration model for Type Ia supernovae: I. Scaling properties

We analyze the statistical properties of the turbulent velocity field in the deflagration model for Type Ia supernovae. In particular, we consider the question of whether turbulence is isotropic and consistent with the Kolmogorov theory at small length scales. Using numerical data from a high-resolution simulation of a thermonuclear supernova explosion, spectra of the turbulence energy and velocity structure functions are computed. We show that the turbulent velocity field is isotropic at small length scales and follows a scaling law that is consistent with the Kolmogorov theory until most of the nuclear fuel is burned. At length scales greater than a certain characteristic scale, turbulence becomes anisotropic. Here, the radial velocity fluctuations follow the scaling law of the Rayleigh-Taylor instability, whereas the angular component still obeys Kolmogorov scaling. In the late phase of the explosion, this characteristic scale drops below the numerical resolution of the simulation. The analysis confirms that a subgrid-scale model for the unresolved turbulence energy is required for the consistent calculation of the flame speed in deflagration models of Type Ia supernovae, and that the assumption of isotropy on these scales is appropriate.

preprint2009arXiv

Turbulence modeling and the physics of the intra-cluster medium

FEARLESS (Fluid mEchanics with Adaptively Refined Large Eddy SimulationS) is a new numerical scheme arising from the combined use of subgrid scale (SGS) model for turbulence at the unresolved length scales and adaptive mesh refinement (AMR) for resolving the large scales. This tool is especially suitable for the study of turbulent flows in strongly clumped media. In this contribution, the main features of FEARLESS are briefly outlined. We then summarize the main results of FEARLESS cosmological simulations of galaxy cluster evolution. In clusters, the production of turbulence is closely correlated with merger events; for minor mergers, we find that turbulent dissipation affects the cluster energy budget only locally. The level of entropy in the cluster core is enhanced in FEARLESS simulations, in accord with a better modeling of the unresolved flow, and with its feedback on the resolved mixing in the ICM.

preprint2005arXiv

Level set simulations of turbulent thermonuclear deflagration in degenerate carbon and oxygen

We study the dynamics of thermonuclear flames propagating in fuel stirred by stochastic forcing. The fuel consists of carbon and oxygen in a state which is encountered in white dwarfs close to the Chandrasekhar limit. The level set method is applied to represent the flame fronts numerically. The computational domain for the numerical simulations is cubic, and periodic boundary conditions are imposed. The goal is the development of a suitable flame speed model for the small-scale dynamics of turbulent deflagration in thermonuclear supernovae. Because the burning process in a supernova explosion is transient and spatially inhomogeneous, the localised determination of subgrid scale closure parameters is essential. We formulate a semi-localised model based on the dynamical equation for the subgrid scale turbulence energy $k_{\mathrm{sgs}}$. The turbulent flame speed $s_{\mathrm{t}}$ is of the order $\sqrt{2k_{\mathrm{sgs}}}$. In particular, the subgrid scale model features a dynamic procedure for the calculation of the turbulent energy transfer from resolved toward subgrid scales, which has been successfully applied to combustion problems in engineering. The options of either including or suppressing inverse energy transfer in the turbulence production term are compared. In combination with the piece-wise parabolic method for the hydrodynamics, our results favour the latter option. Moreover, different choices for the constant of proportionality in the asymptotic flame speed relation, $s_{\mathrm{t}}\propto\sqrt{2k_{\mathrm{sgs}}}$, are investigated.

preprint2000arXiv

Vibrational States of Glassy and Crystalline Orthotherphenyl

Low-frequency vibrations of glassy and crystalline orthoterphenyl are studied by means of neutron scattering. Phonon dispersions are measured along the main axes of a single crystal, and the corresponding longitudinal and transversal sound velocities are obtained. For glassy and polycrystalline samples, a density of vibrational states is determined and cross-checked against other dynamic observables. In the crystal, low-lying zone-boundary modes lead to an excess over the Debye density of states. In the glass, the boson peak is located at even lower frequencies. With increasing temperature, both glass and crystal show anharmonicity.