BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Tiago M. D. Pereira

Tiago M. D. Pereira contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.SRastro-ph.IMphysics.ed-ph

Trust snapshot

Quick read

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

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

15 published item(s)

preprint2026arXiv

Spectral signatures of bright grains determine chromospheric heating

Chromospheric heating is an important ingredient in the energy budget of the solar atmosphere, which is challenging to quantify from observations. By using 3D radiative magnetohydrodynamic simulations of the solar atmosphere combined with non-LTE spectral synthesis, we estimated chromospheric heating from synthetic spectra and studied the spectral and temporal signatures of heating events. We performed k-means clustering on the Mg II h, Ca II H, and Ca II 8542 Å lines to identify representative profiles associated with elevated chromospheric heating and studied their atmospheric stratification. We find that locations with the strongest chromospheric heating show spectral signatures with strong emission. Profiles with strong emission in the blue wing of the lines (blue grains) are created by upward-propagating shock waves and have an order of magnitude higher heating in the chromosphere than the ambient heating. Profiles with strong emission in the red wing (red grains) also display heating that is an order of magnitude stronger than the baseline, but these spectra do not show a characteristic atmospheric stratification. Spectra classified as blue grains have a consistent temporal evolution, which is an oscillating sawtooth pattern in the line core and emission in the blue wing. However, spectra classified as red grains did not show a consistent temporal signature: Red wing emission from the simulations can appear spontaneously or be associated with an oscillation. While red and blue grain profiles account for around 3% of our synthetic spectra, they account for more than 12% of the total chromospheric heating in these simulations. By comparing two quiet Sun simulations, we find that the prevalence of bright grains is influenced by the magnetic field configuration, with a unipolar configuration showing fewer bright grains and consequently a lower share of heating from such events.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Effects of spatial resolution on inferences of atmospheric quantities from simulations

Small scale processes are thought to be important for the dynamics of the solar atmosphere. While numerical resolution fundamentally limits their inclusion in MHD simulations, real observations at the same nominal resolution should still contain imprints of sub-resolution effects. This means that the synthetic observables from a simulation of given resolution might not be directly comparable to real observables at the same resolution. It is thus of interest to investigate how inferences based on synthetic spectra from simulations with different numerical resolutions compare, and whether these differences persist after the spectra have been spatially degraded to a common resolution. We aim to compare synthetic spectra obtained from realistic 3D radiative magnetohydrodynamic (rMHD) simulations run at three different numerical resolutions from the same initial atmosphere, using very simple methods for inferring line-of-sight velocities and magnetic fields. Additionally we examine how the differing spatial resolution impacts the results retrieved from the STiC inversion code. We find that while the simple inferences for all three simulations reveal the same large-scale tendencies, the higher resolutions yield more fine-grained structures and more extreme line-of-sight velocities/magnetic fields in concentrated spots even after spatial smearing. We also see indications that the imprints of sub-resolution effects on the degraded spectra result in systematic errors in the inversions, and that these errors increase with the amount of sub-resolution effects included. Fortunately, however, we find that including successively more sub-resolution yields smaller additional effects; i.e. there is a clear trend of diminishing importance for progressively finer sub-resolution effects.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

SunnyNet: A neural network approach to 3D non-LTE radiative transfer

Context. Computing spectra from 3D simulations of stellar atmospheres when allowing for departures from local thermodynamic equilibrium (non-LTE) is computationally very intensive. Aims. We develop a machine learning based method to speed up 3D non-LTE radiative transfer calculations in optically thick stellar atmospheres. Methods. Making use of a variety of 3D simulations of the solar atmosphere, we trained a convolutional neural network, SunnyNet, to learn the translation from LTE to non-LTE atomic populations. Non-LTE populations computed with an existing 3D code were considered as the true values. The network was then used to predict non-LTE populations for other 3D simulations, and synthetic spectra were computed from its predicted non-LTE populations. We used a six-level model atom of hydrogen and H$α$ spectra as test cases. Results. SunnyNet gives reasonable predictions for non-LTE populations with a dramatic speedup of about 10$^5$ times when running on a single GPU and compared to existing codes. When using different snapshots of the same simulation for training and testing, SunnyNet's predictions are within 20-40% of the true values for most points, which results in average differences of a few percent in H$α$ spectra. Predicted H$α$ intensity maps agree very well with existing codes. Most importantly, they show the telltale signs of 3D radiative transfer in the morphology of chromospheric fibrils. The results are not as reliable when the training and testing are done with different families of simulations. SunnyNet is open source and publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Survey of Computational Tools in Solar Physics

The SunPy Project developed a 13-question survey to understand the software and hardware usage of the solar physics community. 364 members of the solar physics community, across 35 countries, responded to our survey. We found that 99$\pm$0.5% of respondents use software in their research and 66% use the Python scientific software stack. Students are twice as likely as faculty, staff scientists, and researchers to use Python rather than Interactive Data Language (IDL). In this respect, the astrophysics and solar physics communities differ widely: 78% of solar physics faculty, staff scientists, and researchers in our sample uses IDL, compared with 44% of astrophysics faculty and scientists sampled by Momcheva and Tollerud (2015). 63$\pm$4% of respondents have not taken any computer-science courses at an undergraduate or graduate level. We also found that most respondents utilize consumer hardware to run software for solar-physics research. Although 82% of respondents work with data from space-based or ground-based missions, some of which (e.g. the Solar Dynamics Observatory and Daniel K. Inouye Solar Telescope) produce terabytes of data a day, 14% use a regional or national cluster, 5% use a commercial cloud provider, and 29% use exclusively a laptop or desktop. Finally, we found that 73$\pm$4% of respondents cite scientific software in their research, although only 42$\pm$3% do so routinely.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Appearance of Spicules in High Resolution Observations of Ca II H and H-alpha

Solar spicules are chromospheric fibrils that appear everywhere on the Sun, yet their origin is not understood. Using high resolution observations of spicules obtained with the Swedish 1-m Solar Telescope, we aim to understand how spicules appear in filtergrams and Dopplergrams, how they compare in Ca II H and H-alpha, and what can make them appear and disappear. We find thatspicules display a rich and detailed spatial structure, and show a distribution of transverse velocities that when aligned with the line of sight can make them appear at different H-alpha wing positions. They become more abundant at positions closer to the line core, reflecting a distribution of Doppler shifts and widths. In H-alpha width maps they stand out as bright features both on disk and off limb, reflecting their large Doppler motions and possibly higher temperatures than in the typical H-alpha formation region. Spicule lifetimes measured from narrowband images at only a few positions will be an underestimate because Doppler shifts can make them disappear prematurely from such images; for such cases, width maps are a more robust tool. In H-alpha and Ca II H filtergrams, off limb spicules essentially have the same properties, appearance, and evolution. We find that the sudden appearance of spicules can be explained by Doppler shifts from their transverse motions, and does not require other convoluted explanations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Internetwork chromospheric bright grains observed with IRIS

The Interface Region Imaging Spectrograph (IRIS) reveals small-scale rapid brightenings in the form of bright grains all over coronal holes and the quiet sun. These bright grains are seen with the IRIS 1330 Å, 1400 Å and 2796 Å slit-jaw filters. We combine coordinated observations with IRIS and from the ground with the Swedish 1-m Solar Telescope (SST) which allows us to have chromospheric (Ca II 8542 Å, Ca II H 3968 Å, Hα, and Mg II k 2796 Å), and transition region (C II 1334 Å, Si IV 1402) spectral imaging, and single-wavelength Stokes maps in Fe I 6302 Åat high spatial (0.33"), temporal and spectral resolution. We conclude that the IRIS slit-jaw grains are the counterpart of so-called acoustic grains, i.e., resulting from chromospheric acoustic waves in a non-magnetic environment. We compare slit-jaw images with spectra from the IRIS spectrograph. We conclude that the grain intensity in the 2796 Å slit-jaw filter comes from both the Mg II k core and wings. The signal in the C II and Si IV lines is too weak to explain the presence of grains in the 1300 and 1400 Å slit-jaw images and we conclude that the grain signal in these passbands comes mostly from the continuum. Even though weak, the characteristic shock signatures of acoustic grains can often be detected in IRIS C II spectra. For some grains, spectral signature can be found in IRIS Si IV. This suggests that upward propagating acoustic waves sometimes reach all the way up to the transition region.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

RH 1.5D: a massively parallel code for multi-level radiative transfer with partial frequency redistribution and Zeeman polarisation

The emergence of three-dimensional magneto-hydrodynamic (MHD) simulations of stellar atmospheres has sparked a need for efficient radiative transfer codes to calculate detailed synthetic spectra. We present RH 1.5D, a massively parallel code based on the RH code and capable of performing Zeeman polarised multi-level non-local thermodynamical equilibrium (NLTE) calculations with partial frequency redistribution for an arbitrary amount of chemical species. The code calculates spectra from 3D, 2D or 1D atmospheric models on a column-by-column basis (or 1.5D). While the 1.5D approximation breaks down in the cores of very strong lines in an inhomogeneous environment, it is nevertheless suitable for a large range of scenarios and allows for faster convergence with finer control over the iteration of each simulation column. The code scales well to at least tens of thousands of CPU cores, and is publicly available. In the present work we briefly describe its inner workings, strategies for convergence optimisation, its parallelism, and some possible applications.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The formation of IRIS diagnostics VIII. IRIS observations in the C II 133.5 nm multiplet

The C II 133.5 nm multiplet has been observed by NASA's Interface Region Imaging Spectrograph (IRIS) in unprecedented spatial resolution. The aims of this work are to characterize these new observations of the C II lines, place them in context with previous work, and to identify any additional value the C II lines bring when compared with other spectral lines. We make use of wide, long exposure IRIS rasters covering the quiet Sun and an active region. Line properties such as velocity shift and width are extracted from individual spectra and analyzed. The lines have a variety of shapes (mostly single-peak or double-peak), are strongest in active regions and weaker in the quiet Sun. The ratio between the 133.4 nm and 133.5 nm components is always less than 1.8, indicating that their radiation is optically thick in all locations. Maps of the C II line widths are a powerful new diagnostic of chromospheric structures, and their line shifts are a robust velocity diagnostic. Compared with earlier quiet Sun observations, we find similar absolute intensities and mean line widths, but smaller red shifts; this difference can perhaps be attributed to differences in spectral resolution and spatial coverage. The C II intensity maps are somewhat similar to those of transition region lines, but also share some features with chromospheric maps such as those from the Mg II k line, indicating that they are formed between the upper chromosphere and transition region. C II intensity, width, and velocity maps can therefore be used to gather additional information about the upper chromosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Formation of IRIS diagnostics. IV. The Mg II triplet lines as a new diagnostic for lower chromospheric heating

A triplet of subordinate lines of Mg II exists in the region around the h&k lines. In solar spectra these lines are seen mostly in absorption, but in some cases can become emission lines. The aim of this work is to study the formation of this triplet, and investigate any diagnostic value they can bring. Using 3D radiative magnetohydrodynamic simulations of quiet Sun and flaring flux emergence, we synthesize spectra and investigate how spectral features respond to the underlying atmosphere. We find that emission in the lines is rare and is typically caused by a steep temperature increase in the lower chromosphere (above 1500 K, with electron densities above 10$^{18}$ m$^{-3}$). In both simulations the lines are sensitive to temperature increases taking place at column masses >= 5e-4 g cm$^{-2}$. Additional information can also be inferred from the peak-to-wing ratio and shape of the line profiles. Using observations from NASA's Interface Region Imaging Spectrograph we find both absorption and emission line profiles with similar shapes to the synthetic spectra, which suggests that these lines represent a useful diagnostic that complements the MgII h&k lines.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Detailed Comparison Between The Observed and Synthesized Properties of a Simulated Type ii Spicule

We performed a 3D radiative MHD simulation of the solar atmosphere. This simulation shows a jet-like feature that shows similarities to the type II spicules observed for the first time with Hinode. Rapid Blueshifted Events (RBEs) on the solar disk are associated with these spicules. Observational results suggest they may contribute significantly in supplying the corona with hot plasma. We perform a detailed comparison of the properties of the simulated jet with those of type II spicules (observed with Hinode) and RBEs (with ground-based instruments). We analyze variety of synthetic emission and absorption lines from the simulations including chromospheric Ca II and Ha to TR and coronal temperatures (10E4 to several 10E6K). We compare their synthetic intensities, line profiles, Doppler shifts, line widths and asymmetries with observations from Hinode/SOT and EIS, SOHO/SUMER, SST and SDO/AIA. Many properties of the synthetic observables resemble the observations, and we describe in detail the physical processes that lead to these observables. Detailed analysis of the synthetic observables provides insight into how observations should be analyzed to derive information about physical variables in such a dynamic event. For example, we find that LOS superposition in the optically thin atmosphere requires the combination of Doppler shifts and spectral line asymmetry to determine the velocity in the jet. Other properties differ from the observations, especially in the chromospheric lines. The mass density of the part of the spicule with a chromospheric temperature is too low to produce significant opacity in chromospheric lines. These and other discrepancies are described in detail, and we discuss which mechanisms and physical processes may need to be included in the MHD simulations to mimic the thermodynamic processes of the chromosphere and corona, in particular to reproduce type II spicules.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

How realistic are solar model atmospheres?

Recently, new solar model atmospheres have been developed to replace classical 1D LTE hydrostatic models and used to for example derive the solar chemical composition. We aim to test various models against key observational constraints. In particular, a 3D model used to derive the solar abundances, a 3D MHD model (with an imposed 10 mT vertical magnetic field), 1D models from the PHOENIX project, the 1D MARCS model, and the 1D semi-empirical model of Holweger & Müller. We confront the models with observational diagnostics of the temperature profile: continuum centre-to-limb variations (CLV), absolute continuum fluxes, and the wings of hydrogen lines. We also test the 3D models for the intensity distribution of the granulation and spectral line shapes. The predictions from the 3D model are in excellent agreement with the continuum CLV observations, performing even better than the Holweger & Müller model (constructed largely to fulfil such observations). The predictions of the 1D theoretical models are worse, given their steeper temperature gradients. For the continuum fluxes, predictions for most models agree well with the observations. No model fits all hydrogen lines perfectly, but again the 3D model comes ahead. The 3D model also reproduces the observed continuum intensity fluctuations and spectral line shapes very well. The excellent agreement of the 3D model with the observables reinforces the view that its temperature structure is realistic. It outperforms the MHD simulation in all diagnostics, implying that recent claims for revised abundances based on MHD modelling are premature. Several weaknesses in the 1D models are exposed. The differences between the PHOENIX LTE and NLTE models are small. We conclude that the 3D hydrodynamical model is superior to any of the tested 1D models, which gives further confidence in the solar abundance analyses based on it.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The formation of IRIS diagnostics. III. Near-ultraviolet Spectra and Images

The Mg II h&k lines are the prime chromospheric diagnostics of NASA's Interface Region Imaging Spectrograph (IRIS). In the previous papers of this series we used a realistic three-dimensional radiative magnetohydrodynamics model to calculate the h&k lines in detail and investigated how their spectral features relate to the underlying atmosphere. In this work, we employ the same approach to investigate how the h&k diagnostics fare when taking into account the finite resolution of IRIS and different noise levels. In addition, we investigate the diagnostic potential of several other photospheric lines and near-continuum regions present in the near-ultraviolet (NUV) window of IRIS and study the formation of the NUV slit-jaw images. We find that the instrumental resolution of IRIS has a small effect on the quality of the h&k diagnostics; the relations between the spectral features and atmospheric properties are mostly unchanged. The peak separation is the most affected diagnostic, but mainly due to limitations of the simulation. The effects of noise start to be noticeable at a signal-to-noise ratio (S/N) of 20, but we show that with noise filtering one can obtain reliable diagnostics at least down to a S/N of 5. The many photospheric lines present in the NUV window provide velocity information for at least eight distinct photospheric heights. Using line-free regions in the h&k far wings we derive good estimates of photospheric temperature for at least three heights. Both of these diagnostics, in particular the latter, can be obtained even at S/Ns as low as 5.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quantifying Spicules

Understanding the dynamic solar chromosphere is fundamental in solar physics. Spicules are an important feature of the chromosphere, connecting the photosphere to the corona, potentially mediating the transfer of energy and mass. The aim of this work is to study the properties of spicules over different regions of the sun. Our goal is to investigate if there is more than one type of spicules, and how spicules behave in the quiet sun, coronal holes, and active regions. We make use of high-cadence and high-spatial resolution Ca II H observations taken by Hinode/SOT. Making use of a semi-automated detection algorithm, we self-consistently track and measure the properties of 519 spicules over different regions. We find clear evidence of two types of spicules. Type I spicules show a rise and fall and have typical lifetimes of 150-400 s and maximum ascending velocities of 15-40 km/s, while type II spicules have shorter lifetimes of 50-150 s, faster velocities of 30-110 km/s, and are not seen to fall down, but rather fade at around their maximum length. Type II spicules are the most common, seen in quiet sun and coronal holes. Type I spicules are seen mostly in active regions. There are regional differences between quiet sun and coronal hole spicules, likely attributable to the different field configurations. The properties of type II spicules are consistent with published results of Rapid Blueshifted Events (RBEs), supporting the hypothesis that RBEs are their disk counterparts. For type I spicules we find the relations between their properties to be consistent with a magnetoacoustic shock wave driver, and with dynamic fibrils as their disk counterpart. The driver of type II spicules remains unclear from limb observations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Oxygen lines in solar granulation. I. Testing 3D models against new observations with high spatial and spectral resolution

Aims: we seek to provide additional tests of the line formation of theoretical 3D solar photosphere models. In particular, we set out to test the spatially-resolved line formation at several viewing angles, from the solar disk-centre to the limb and focusing on atomic oxygen lines. The purpose of these tests is to provide additional information on whether the 3D model is suitable to derive the solar oxygen abundance. We also aim to empirically constrain the NLTE recipes for neutral hydrogen collisions, using the spatially-resolved observations of the OI 777 nm lines. Methods: using the Swedish 1-m Solar Telescope we obtained high-spatial-resolution observations of five atomic oxygen lines (along with lines for other species) for five positions on the solar disk. These observations have a high spatial and spectral resolution, and a continuum intensity contrast up to 9% at 615 nm. The theoretical line profiles were computed using the 3D model, with a full 3D NLTE treatment for oxygen and LTE for the other lines. Results: at disk-centre we find an excellent agreement between predicted and observed line shifts, strengths, FWHM and asymmetries. At other viewing angles the agreement is also good, but the smaller continuum intensity contrast makes a quantitative comparison harder. We use the disk-centre observations we constrain S_H, the scaling factor for the efficiency of collisions with neutral hydrogen. We find that S_H=1 provides the best match to the observations. Conclusions: overall there is a very good agreement between predicted and observed line properties over the solar granulation. This further reinforces the view that the 3D model is realistic and a reliable tool to derive the solar oxygen abundance.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Oxygen lines in solar granulation. II. Centre-to-limb variation, NLTE line formation, blends and the solar oxygen abundance

Context: There is a lively debate about the solar oxygen abundance and the role of 3D models in its recent downward revision. The models have been tested using high resolution solar atlases. Further testing can be done using centre-to-limb variations. Aims: Using high quality observations of oxygen lines across the solar surface we seek to test if the 3D and 1D models reproduce their observed centre-to-limb variation (CLV). In particular we seek to assess whether the 3D model is appropriate to derive the solar oxygen abundance. Methods: We use our recent observations of OI 777 nm, OI 615.81 nm, [OI] 630.03 nm and nine lines of other elements for five viewing angles 0.2<mu<1 of the quiet solar disk. We compare them with the predicted line profiles from the 3D and 1D models computed using the most up-to-date line formation codes, line data and allowing for departures of LTE. The CLV of the OI 777 nm lines is also used to obtain an empirical correction for the poorly known efficiency of the inelastic collisions with H. Results: The 3D model generally reproduces the CLV observations of the lines very well, particularly the oxygen lines. From the OI 777 nm lines we find that the classical Drawin recipe slightly overestimates H collisions. The limb observations of the OI 615.82 nm line allow us to identify a previously unknown contribution of molecules for this line, prevalent at the solar limb. A detailed treatment of the [OI] 630.03 nm line shows that the 3D modeling provides an excellent agreement with the observations. The derived oxygen abundances with the 3D model are 8.68 (OI 777 nm), 8.66 ([OI] 630.03 nm) and 8.62 (OI 615.82 nm). Conclusions: These additional tests have reinforced the trustworthiness of the 3D model and line formation for abundance analyses.

0 citations0 reviewsNo code linkedOpen access