Showing 1–50 of 57 results for author: Guerrero, G

Simulated Coronal Mass Ejections on a young Solar-Type Star and the Associated Instantaneous Angular Momentum Loss

Authors: Yu Xu, Julián D. Alvarado-Gómez, Hui Tian, Katja Poppenhäger, Gustavo Guerrero, Xianyu Liu

Abstract: Coronal mass ejections (CMEs) on stars can change the stars' magnetic field configurations and mass loss rates during the eruption and propagation and therefore, may affect the stars' rotation properties on long time-scales. The dynamics of stellar CMEs and their influence on the stellar angular momentum loss rate are not yet well understood. In order to start investigating these CME-related aspec… ▽ More Coronal mass ejections (CMEs) on stars can change the stars' magnetic field configurations and mass loss rates during the eruption and propagation and therefore, may affect the stars' rotation properties on long time-scales. The dynamics of stellar CMEs and their influence on the stellar angular momentum loss rate are not yet well understood. In order to start investigating these CME-related aspects on other stars, we conducted a series of magnetohydrodynamic simulations of CMEs on a solar-type star of moderate activity levels. The propagation and evolution of the CMEs were traced in the three-dimensional outputs and the temporal evolution of their dynamic properties (such as masses, velocities, and kinetic energies) were determined. The simulated stellar CMEs are more massive and energetic than their solar analog, which is a result of the stronger magnetic field on the surface of the simulated star than that of the Sun. The simulated CMEs display masses ranging from ~10^16 g to ~10^18 g and kinetic energies from ~10^31 erg to ~10^33 erg. We also investigated the instantaneous influence of the CMEs to the star's angular momentum loss rate. Our results suggest that angular momentum can either be added to or be removed from the star during the evolution of CME events. We found a positive correlation between the amplitude of the angular momentum loss rate variation and the CME's kinetic energy as well as mass, suggesting that more energetic/massive CMEs have higher possibility to add angular momentum to the star. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 16 pages, 10 figures; accepted by ApJ

Dynamics of the tachocline

Authors: Antoine Strugarek, Bernadett Belucz, Allan Sacha Brun, Mausumi Dikpati, Gustavo Guerrero

Abstract: The solar tachocline is an internal region of the Sun possessing strong radial and latitudinal shears straddling the base of the convective envelope. Based on helioseismic inversions, the tachocline is known to be thin (less than 5\% of the solar radius). Since the first theory of the solar tachocline in 1992, this thinness has not ceased to puzzle solar physicists. In this review, we lay out the… ▽ More The solar tachocline is an internal region of the Sun possessing strong radial and latitudinal shears straddling the base of the convective envelope. Based on helioseismic inversions, the tachocline is known to be thin (less than 5\% of the solar radius). Since the first theory of the solar tachocline in 1992, this thinness has not ceased to puzzle solar physicists. In this review, we lay out the grounds of our understanding of this fascinating region of the solar interior. We detail the various physical mechanisms at stake in the solar tachocline, and put a particular focus on the mechanisms that have been proposed to explain its thinness. We also examine the full range of MHD processes including waves and instabilies that are likely to occur in the tachocline, as well as their possible connection with active region patterns observed at the surface. We reflect on the most recent findings for each of them, and highlight the physical understanding that is still missing and that would allow the research community to understand, in a generic sense, how the solar tachocline and stellar tachocline are formed, are sustained, and evolve on secular timescales. △ Less

Submitted 21 November, 2023; originally announced November 2023.

Comments: 46 pages, 10 figures, accepted for publication in Space Science Reviews

Simulations of solar and stellar dynamos and their theoretical interpretation

Authors: Petri J. Käpylä, Matthew K. Browning, Allan Sacha Brun, Gustavo Guerrero, Jörn Warnecke

Abstract: We review the state of the art of three dimensional numerical simulations of solar and stellar dynamos. We summarize fundamental constraints of numerical modelling and the techniques to alleviate these restrictions. Brief summary of the relevant observations that the simulations seek to capture is given. We survey the current progress of simulations of solar convection and the resulting large-scal… ▽ More We review the state of the art of three dimensional numerical simulations of solar and stellar dynamos. We summarize fundamental constraints of numerical modelling and the techniques to alleviate these restrictions. Brief summary of the relevant observations that the simulations seek to capture is given. We survey the current progress of simulations of solar convection and the resulting large-scale dynamo. We continue to studies that model the Sun at different ages and to studies of stars of different masses and evolutionary stages. Both simulations and observations indicate that rotation, measured by the Rossby number which is the ratio of rotation period and convective turnover time, is a key ingredient in setting the overall level and characteristics of magnetic activity. Finally, efforts to understand global 3D simulations in terms of mean-field dynamo theory are discussed. △ Less

Submitted 26 May, 2023; originally announced May 2023.

Comments: 70 pages, 9 figures, submitted to Space Science Reviews

The Puzzling Structure of Solar Convection: Window into the Dynamo

Authors: Nicholas A. Featherstone, Kyle C. Augustson, Jonathan M. Aurnou, Catherine Blume, Benjamin P. Brown, Nicholas Brummell, Keaton J. Burns, Michael A. Calkins, Maria Camisassa, Mausumi Dikpati, Yuhong Fan, J. R. Fuentes, Gustavo Guerrero, Bradley W. Hindman, Keith Julien, Irina N. Kitiashvili, Lydia Korre, Daniel Lecoanet, Bhishek Manek, Loren Matilsky, Mark Miesch, Nicholas J. Nelson, Jeffrey S. Oishi, Whitney T. Powers, Matthias Rempel , et al. (3 additional authors not shown)

Abstract: The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective… ▽ More The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective flows. Those flows are in-turn thought to sustain both the Sun's turbulent EMF and the large-scale flows of differential rotation and meridional circulation suspected of influencing the dynamo's organization and timing. Continued progress in this area is complicated by (i) inconsistencies between helioseismic measurements of convective and meridional flow made with different techniques and instruments, and (ii) a lack of high-latitude data for convection, differential rotation, and meridional flow. We suggest that the path forward to resolving these difficulties is twofold. First, the acquisition of long-term helioseismic and emissivity measurements obtained from a polar vantage point is vital to complete our picture of the Sun's outer convection zone. Second, sustained and expanded investment in theory-oriented and combined theory/observational research initiatives will be crucial to fully exploit these new observations and to resolve inconsistencies between existing measurements. △ Less

Submitted 15 May, 2023; originally announced May 2023.

Comments: White paper submitted to the Decadal Survey for Solar and Space Physics 2024-2033. 7 pages; 4 figures

Global simulations of Tayler instability in stellar interiors: a long-time multi-stage evolution of the magnetic field

Authors: G. Monteiro, G. Guerrero, F. Del Sordo, A. Bonanno, P. K. Smolarkiewicz

Abstract: Magnetic fields have been observed in massive Ap/Bp stars and presumably are also present in the radiative zone of solar-like stars. Yet, to date there is no clear understanding of the dynamics of the magnetic field in stably stratified layers. A purely toroidal magnetic field configuration is known to be unstable, developing mainly non-axisymmetric modes. Rotation and a small poloidal field compo… ▽ More Magnetic fields have been observed in massive Ap/Bp stars and presumably are also present in the radiative zone of solar-like stars. Yet, to date there is no clear understanding of the dynamics of the magnetic field in stably stratified layers. A purely toroidal magnetic field configuration is known to be unstable, developing mainly non-axisymmetric modes. Rotation and a small poloidal field component may lead to a stable configuration. Here we perform global MHD simulations with the EULAG-MHD code to explore the evolution of a toroidal magnetic field located in a layer whose stratification resembles the solar tachocline. Our numerical experiments allow us to explore the initial unstable phase as well as the long-term evolution of the magnetic field. During the first Alfven cycles, we observe the development of the Tayler instability with the prominent longitudinal wavenumber, $m=1$. Rotation decreases the growth rate of the instability, and eventually suppresses it. However, after a stable phase, sudden energy surges lead to the development of higher order modes even for fast rotation. These modes extract energy from the initial toroidal field. Nevertheless, our results show that sufficiently fast rotation leads to a lower saturation energy of the unstable modes, resulting in a magnetic topology with only a small fraction of poloidal field which remains steady for several hundreds of Alfven travel times. At this stage, the system becomes turbulent and the field is prone to turbulent diffusion. The final toroidal-poloidal configuration of the magnetic field may represent an important aspect of the field generation and evolution in stably-stratified layers. △ Less

Submitted 24 March, 2023; v1 submitted 18 November, 2022; originally announced November 2022.

Comments: 15 pages, 16 figures, submitted to MNRAS

Journal ref: Volume 521, Issue 1, May 2023, Pages 1415 1428

Implicit large eddy simulations of global solar convection: effects of numerical resolution in non-rotating and rotating cases

Authors: G. Guerrero, A. M. Stejko, A. G. Kosovichev, P. K Smolarkiewicz, A. Strugarek

Abstract: Simulating deep solar convection and its coupled mean-field motions is a formidable challenge where few observational results constrain models that suffer from the non-physical influence of the grid resolution. We present hydrodynamic global Implicit Large-Eddy simulations (ILES) of deep solar convection performed with the EULAG-MHD code, and explore the effects of grid resolution on the propertie… ▽ More Simulating deep solar convection and its coupled mean-field motions is a formidable challenge where few observational results constrain models that suffer from the non-physical influence of the grid resolution. We present hydrodynamic global Implicit Large-Eddy simulations (ILES) of deep solar convection performed with the EULAG-MHD code, and explore the effects of grid resolution on the properties of rotating and non-rotating convection. The results, based on low-order moments and turbulent spectra reveal that convergence could be achieved in non-rotating simulations provided sufficient resolution in the radial direction. The flow is highly anisotropic, with the energy contained in horizontal divergent motions exceeding by more than three orders of magnitude their radial counterpart. By contrast, in rotating simulations the largest energy is in the toroidal part of the horizontal motions. As the grid resolution increases, the turbulent correlations change in such a way that a solar-like differential rotation, obtained in the simulation with the coarser grid, transitions to the anti-solar differential rotation. The reason for this change is the contribution of the effective viscosity to the balance of the forces driving large-scale flows. As the effective viscosity decreases, the angular momentum balance improves, yet the force balance in the meridional direction lessens, favoring a strong meridional flow that advects angular momentum towards the poles. The results suggest that obtaining the correct distribution of angular momentum may not be a mere issue of numerical resolution. Accounting for additional physics, such as magnetism or the near-surface shear layer, may be necessary in simulating the solar interior. △ Less

Submitted 11 August, 2022; originally announced August 2022.

Comments: 23 pages, 16 figures, submitted to ApJ

Identification of inertial modes in the solar convection zone

Authors: Santiago Andrés Triana, Gustavo Guerrero, Ankit Barik, Jérémy Rekier

Abstract: The observation of global acoustic waves (p modes) in the Sun has been key to unveiling its internal structure and dynamics. A different kind of wave, known as sectoral Rossby modes, have been observed and identified, which potentially opens the door to probing internal processes that are inaccessible through p mode helioseismology. Yet another set of waves, appearing as retrograde-propagating, eq… ▽ More The observation of global acoustic waves (p modes) in the Sun has been key to unveiling its internal structure and dynamics. A different kind of wave, known as sectoral Rossby modes, have been observed and identified, which potentially opens the door to probing internal processes that are inaccessible through p mode helioseismology. Yet another set of waves, appearing as retrograde-propagating, equatorially antisymmetric vorticity waves, have also been observed but their identification remained elusive. Here, through a numerical model implemented as an eigenvalue problem, we provide evidence supporting the identification of those waves as a class of inertial eigenmodes, distinct from the Rossby mode class, with radial velocities comparable to the horizontal ones deep in the convective zone, but still small compared to the horizontal velocities towards the surface. We also suggest that the signature of tesseral-like Rossby modes might be present in the recent observational data. △ Less

Submitted 1 July, 2022; v1 submitted 27 April, 2022; originally announced April 2022.

Comments: 12 pages, 6 figures, accepted for publication in the Astrophysical Journal Letters

Constraining Global Solar Models through Helioseismic Analysis

Authors: Andrey M. Stejko, Alexander G. Kosovichev, Nicholas A. Featherstone, Gustavo Guerrero, Bradley W. Hindman, Loren I. Matilsky, Jörn Warnecke

Abstract: Global hydrodynamic simulations of internal solar dynamics have focused on replicating the conditions for solar-like differential rotation and meridional circulation using the results of helioseismic inversions as a constraint. Inferences of meridional circulation, however, have provided controversial results showing the possibility of one, two, or multiple cells along the radius. To resolve this… ▽ More Global hydrodynamic simulations of internal solar dynamics have focused on replicating the conditions for solar-like differential rotation and meridional circulation using the results of helioseismic inversions as a constraint. Inferences of meridional circulation, however, have provided controversial results showing the possibility of one, two, or multiple cells along the radius. To resolve this controversy and develop a more robust understanding of global flow regimes in the solar interior, we apply a "forward-modeling" approach to the analysis of helioseismic signatures of meridional circulation profiles obtained from numerical simulations. We employ the global acoustic modeling code GALE to simulate the propagation of acoustic waves through regimes of mean mass flows generated by global hydrodynamic and magnetohydrodynamic models: EULAG, the Pencil Code, and the Rayleigh code. These models are used to create synthetic dopplergram data products, used as inputs for local time-distance helioseismology techniques. Helioseismic travel-time signals from solutions obtained through global numerical simulations are compared directly with inferences from solar observations, in order to set additional constraints on global model parameters in a direct way. We show that even though these models are able to replicate solar-like differential rotation, the resulting rotationally-constrained convection develops a multi-cell global meridional circulation profile that is measurably inconsistent with local time-distance inferences of solar observations. However, we find that the development of rotationally-unconstrained convection close to the model surface is able to maintain solar-like differential rotation, while having a significant impact on the helioseismic travel-time signal, replicating solar observations within one standard deviation of the error due to noise. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: 13 pages, 7 figures

Advances and Challenges in Observations and Modeling of the Global-Sun Dynamics and Dynamo

Authors: A. G. Kosovichev, G. Guerrero, A. M. Stejko, V. V. Pipin, A. V. Getling

Abstract: Computational heliophysics has shed light on the fundamental physical processes inside the Sun, such as the differential rotation, meridional circulation, and dynamo-generation of magnetic fields. However, despite the substantial advances, the current results of 3D MHD simulations are still far from reproducing helioseismic inferences and surface observations. The reason is the multi-scale nature… ▽ More Computational heliophysics has shed light on the fundamental physical processes inside the Sun, such as the differential rotation, meridional circulation, and dynamo-generation of magnetic fields. However, despite the substantial advances, the current results of 3D MHD simulations are still far from reproducing helioseismic inferences and surface observations. The reason is the multi-scale nature of the solar dynamics, covering a vast range of scales, which cannot be solved with the current computational resources. In such a situation, significant progress has been achieved by the mean-field approach, based on the separation of small-scale turbulence and large-scale dynamics. The mean-field simulations can reproduce solar observations, qualitatively and quantitatively, and uncover new phenomena. However, they do not reveal the complex physics of large-scale convection, solar magnetic cycles, and the magnetic self-organization that causes sunspots and solar eruptions. Thus, developing a synergy of these approaches seems to be a necessary but very challenging task. △ Less

Submitted 20 March, 2022; originally announced March 2022.

Comments: 20 pages, 11 figures, submitted for publication in Proceedings of IAU Symposium 362 "Predictive Power of Computational Astrophysics as a Discovery Tool"

Leptocline as a Shallow Substructure of Near-Surface Shear Layer in 3D Radiative Hydrodynamic Simulations

Authors: Irina N. Kitiashvili, Alexander G. Kosovichev, Alan A. Wray, Viacheslav M. Sadykov, Gustavo Guerrero

Abstract: Understanding effects driven by rotation in the solar convection zone is essential for many problems related to solar activity, such as the formation of differential rotation, meridional circulation, and others. We analyze realistic 3D radiative hydrodynamics simulations of solar subsurface dynamics in the presence of rotation in a local domain 80 Mm wide and 25 Mm deep, located at 30 degrees lati… ▽ More Understanding effects driven by rotation in the solar convection zone is essential for many problems related to solar activity, such as the formation of differential rotation, meridional circulation, and others. We analyze realistic 3D radiative hydrodynamics simulations of solar subsurface dynamics in the presence of rotation in a local domain 80 Mm wide and 25 Mm deep, located at 30 degrees latitude. The simulation results reveal the development of a shallow 10-Mm deep substructure of the Near-Surface Shear Layer (NSSL), characterized by a strong radial rotational gradient and self-organized meridional flows. This shallow layer ("leptocline") is located in the hydrogen ionization zone associated with enhanced anisotropic overshooting-type flows into a less unstable layer between the H and HeII ionization zones. We discuss current observational evidence of the presence of the leptocline and show that the radial variations of the differential rotation and meridional flow profiles obtained from the simulations in this layer qualitatively agree with helioseismic observations. △ Less

Submitted 27 August, 2022; v1 submitted 2 March, 2022; originally announced March 2022.

Comments: 11 pages, 8 figures, submitted to MNRAS

Numerical convergence of 2D solar convection in implicit large-eddy simulations

Authors: H. D. Nogueira, G. Guerrero, P. K. Smolarkiewicz, A. G. Kosovichev

Abstract: Large-eddy simulations (LES) and implicit LES (ILES) are wise and affordable alternatives to the unfeasible direct numerical simulations (DNS) of turbulent flows at high Reynolds numbers (Re). However, for systems with few observational constraints, it is a formidable challenge to determine if these strategies adequately capture the physics of the system. Here we address this problem by analyzing… ▽ More Large-eddy simulations (LES) and implicit LES (ILES) are wise and affordable alternatives to the unfeasible direct numerical simulations (DNS) of turbulent flows at high Reynolds numbers (Re). However, for systems with few observational constraints, it is a formidable challenge to determine if these strategies adequately capture the physics of the system. Here we address this problem by analyzing numerical convergence of ILES of turbulent convection in 2D, with resolutions between $64^2$ and $2048^2$ grid points, along with the estimation of their effective viscosities, resulting in effective Reynolds numbers between $1$ and $\sim10^4$. The thermodynamic structure of our model resembles the solar interior, including a fraction of the radiative zone and the convection zone. In the convective layer, the ILES solutions converge for the simulations with $\ge 512^2$ grid points, as evidenced by the integral properties of the flow and its power spectra. Most importantly, we found that even a resolution of $128^2$ grid points, Re$\sim10$, is sufficient to capture the dynamics of the large scales accurately. This is a consequence of the ILES method allowing that the energy contained in these scales is the same in simulations with low and high resolution. Special attention is needed in regions with a small density scale height driving the formation of fine structures unresolved by the numerical grid. In the stable layer we found the excitation of internal gravity waves, yet high resolution is needed to capture their development and interaction. △ Less

Submitted 24 February, 2022; v1 submitted 6 February, 2022; originally announced February 2022.

Comments: 15 pages, 9 figures, accepted for publication in The Astrophysical Journal

All-sky Search for Continuous Gravitational Waves from Isolated Neutron Stars in the Early O3 LIGO Data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, A. Adams, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, K. M. Aleman, G. Allen, A. Allocca , et al. (1566 additional authors not shown)

Abstract: We report on an all-sky search for continuous gravitational waves in the frequency band 20-2000\,Hz and with a frequency time derivative in the range of $[-1.0, +0.1]\times10^{-8}$\,Hz/s. Such a signal could be produced by a nearby, spinning and slightly non-axisymmetric isolated neutron star in our galaxy. This search uses the LIGO data from the first six months of Advanced LIGO's and Advanced Vi… ▽ More We report on an all-sky search for continuous gravitational waves in the frequency band 20-2000\,Hz and with a frequency time derivative in the range of $[-1.0, +0.1]\times10^{-8}$\,Hz/s. Such a signal could be produced by a nearby, spinning and slightly non-axisymmetric isolated neutron star in our galaxy. This search uses the LIGO data from the first six months of Advanced LIGO's and Advanced Virgo's third observational run, O3. No periodic gravitational wave signals are observed, and 95\%\ confidence-level (CL) frequentist upper limits are placed on their strengths. The lowest upper limits on worst-case (linearly polarized) strain amplitude $h_0$ are $~1.7\times10^{-25}$ near 200\,Hz. For a circularly polarized source (most favorable orientation), the lowest upper limits are $\sim6.3\times10^{-26}$. These strict frequentist upper limits refer to all sky locations and the entire range of frequency derivative values. For a population-averaged ensemble of sky locations and stellar orientations, the lowest 95\%\ CL upper limits on the strain amplitude are $\sim1.\times10^{-25}$. These upper limits improve upon our previously published all-sky results, with the greatest improvement (factor of $\sim$2) seen at higher frequencies, in part because quantum squeezing has dramatically improved the detector noise level relative to the second observational run, O2. These limits are the most constraining to date over most of the parameter space searched. △ Less

Submitted 8 October, 2021; v1 submitted 1 July, 2021; originally announced July 2021.

Comments: 28 pages, 7 figures

Report number: LIGO-P2000334-v9

Journal ref: Phys. Rev. D 104, 082004 (2021)

Observation of gravitational waves from two neutron star-black hole coalescences

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, A. Adams, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, K. M. Aleman, G. Allen, A. Allocca , et al. (1577 additional authors not shown)

Abstract: We report the observation of gravitational waves from two compact binary coalescences in LIGO's and Virgo's third observing run with properties consistent with neutron star-black hole (NSBH) binaries. The two events are named GW200105_162426 and GW200115_042309, abbreviated as GW200105 and GW200115; the first was observed by LIGO Livingston and Virgo, and the second by all three LIGO-Virgo detecto… ▽ More We report the observation of gravitational waves from two compact binary coalescences in LIGO's and Virgo's third observing run with properties consistent with neutron star-black hole (NSBH) binaries. The two events are named GW200105_162426 and GW200115_042309, abbreviated as GW200105 and GW200115; the first was observed by LIGO Livingston and Virgo, and the second by all three LIGO-Virgo detectors. The source of GW200105 has component masses $8.9^{+1.2}_{-1.5}\,M_\odot$ and $1.9^{+0.3}_{-0.2}\,M_\odot$, whereas the source of GW200115 has component masses $5.7^{+1.8}_{-2.1}\,M_\odot$ and $1.5^{+0.7}_{-0.3}\,M_\odot$ (all measurements quoted at the 90% credible level). The probability that the secondary's mass is below the maximal mass of a neutron star is 89%-96% and 87%-98%, respectively, for GW200105 and GW200115, with the ranges arising from different astrophysical assumptions. The source luminosity distances are $280^{+110}_{-110}$ Mpc and $300^{+150}_{-100}$ Mpc, respectively. The magnitude of the primary spin of GW200105 is less than 0.23 at the 90% credible level, and its orientation is unconstrained. For GW200115, the primary spin has a negative spin projection onto the orbital angular momentum at 88% probability. We are unable to constrain spin or tidal deformation of the secondary component for either event. We infer a NSBH merger rate density of $45^{+75}_{-33}\,\mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$ when assuming GW200105 and GW200115 are representative of the NSBH population, or $130^{+112}_{-69}\,\mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$ under the assumption of a broader distribution of component masses. △ Less

Submitted 29 June, 2021; originally announced June 2021.

Report number: LIGO Document P2000357

Journal ref: ApJL, 915, L5 (2021)

Searches for continuous gravitational waves from young supernova remnants in the early third observing run of Advanced LIGO and Virgo

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, A. Adams, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, K. M. Aleman, G. Allen, A. Allocca , et al. (1567 additional authors not shown)

Abstract: We present results of three wide-band directed searches for continuous gravitational waves from 15 young supernova remnants in the first half of the third Advanced LIGO and Virgo observing run. We use three search pipelines with distinct signal models and methods of identifying noise artifacts. Without ephemerides of these sources, the searches are conducted over a frequency band spanning from 10~… ▽ More We present results of three wide-band directed searches for continuous gravitational waves from 15 young supernova remnants in the first half of the third Advanced LIGO and Virgo observing run. We use three search pipelines with distinct signal models and methods of identifying noise artifacts. Without ephemerides of these sources, the searches are conducted over a frequency band spanning from 10~Hz to 2~kHz. We find no evidence of continuous gravitational radiation from these sources. We set upper limits on the intrinsic signal strain at 95\% confidence level in sample sub-bands, estimate the sensitivity in the full band, and derive the corresponding constraints on the fiducial neutron star ellipticity and $r$-mode amplitude. The best 95\% confidence constraints placed on the signal strain are $7.7\times 10^{-26}$ and $7.8\times 10^{-26}$ near 200~Hz for the supernova remnants G39.2--0.3 and G65.7+1.2, respectively. The most stringent constraints on the ellipticity and $r$-mode amplitude reach $\lesssim 10^{-7}$ and $ \lesssim 10^{-5}$, respectively, at frequencies above $\sim 400$~Hz for the closest supernova remnant G266.2--1.2/Vela Jr. △ Less

Submitted 14 July, 2021; v1 submitted 24 May, 2021; originally announced May 2021.

Comments: https://dcc.ligo.org/P2000479

Search for lensing signatures in the gravitational-wave observations from the first half of LIGO-Virgo's third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, A. Adams, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, K. M. Aleman, G. Allen, A. Allocca, P. A. Altin, A. Amato , et al. (1356 additional authors not shown)

Abstract: We search for signatures of gravitational lensing in the gravitational-wave signals from compact binary coalescences detected by Advanced LIGO and Advanced Virgo during O3a, the first half of their third observing run. We study: 1) the expected rate of lensing at current detector sensitivity and the implications of a non-observation of strong lensing or a stochastic gravitational-wave background o… ▽ More We search for signatures of gravitational lensing in the gravitational-wave signals from compact binary coalescences detected by Advanced LIGO and Advanced Virgo during O3a, the first half of their third observing run. We study: 1) the expected rate of lensing at current detector sensitivity and the implications of a non-observation of strong lensing or a stochastic gravitational-wave background on the merger-rate density at high redshift; 2) how the interpretation of individual high-mass events would change if they were found to be lensed; 3) the possibility of multiple images due to strong lensing by galaxies or galaxy clusters; and 4) possible wave-optics effects due to point-mass microlenses. Several pairs of signals in the multiple-image analysis show similar parameters and, in this sense, are nominally consistent with the strong lensing hypothesis. However, taking into account population priors, selection effects, and the prior odds against lensing, these events do not provide sufficient evidence for lensing. Overall, we find no compelling evidence for lensing in the observed gravitational-wave signals from any of these analyses. △ Less

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

Comments: 31 pages and 6 figures. Accepted by the Astrophysical Journal

Report number: LIGO-P2000400

Constraints from LIGO O3 data on gravitational-wave emission due to r-modes in the glitching pulsar PSR J0537-6910

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, A. Adams, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, K. M. Aleman, G. Allen, A. Allocca , et al. (1574 additional authors not shown)

Abstract: We present a search for continuous gravitational-wave emission due to r-modes in the pulsar PSR J0537-6910 using data from the LIGO-Virgo Collaboration observing run O3. PSR J0537-6910 is a young energetic X-ray pulsar and is the most frequent glitcher known. The inter-glitch braking index of the pulsar suggests that gravitational-wave emission due to r-mode oscillations may play an important role… ▽ More We present a search for continuous gravitational-wave emission due to r-modes in the pulsar PSR J0537-6910 using data from the LIGO-Virgo Collaboration observing run O3. PSR J0537-6910 is a young energetic X-ray pulsar and is the most frequent glitcher known. The inter-glitch braking index of the pulsar suggests that gravitational-wave emission due to r-mode oscillations may play an important role in the spin evolution of this pulsar. Theoretical models confirm this possibility and predict emission at a level that can be probed by ground-based detectors. In order to explore this scenario, we search for r-mode emission in the epochs between glitches by using a contemporaneous timing ephemeris obtained from NICER data. We do not detect any signals in the theoretically expected band of 86-97 Hz, and report upper limits on the amplitude of the gravitational waves. Our results improve on previous amplitude upper limits from r-modes in J0537-6910 by a factor of up to 3 and place stringent constraints on theoretical models for r-mode driven spin-down in PSR J0537-6910, especially for higher frequencies at which our results reach below the spin-down limit defined by energy conservation. △ Less

Submitted 7 January, 2022; v1 submitted 29 April, 2021; originally announced April 2021.

Comments: 28 pages, 19 figures, accepted in ApJ

Report number: LIGO-P2100069

Journal ref: ApJ 922 71 (2021)

Constraints on cosmic strings using data from the third Advanced LIGO-Virgo observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, A. Adams, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, K. M. Aleman, G. Allen, A. Allocca , et al. (1565 additional authors not shown)

Abstract: We search for gravitational-wave signals produced by cosmic strings in the Advanced LIGO and Virgo full O3 data set. Search results are presented for gravitational waves produced by cosmic string loop features such as cusps, kinks and, for the first time, kink-kink collisions.cA template-based search for short-duration transient signals does not yield a detection. We also use the stochastic gravit… ▽ More We search for gravitational-wave signals produced by cosmic strings in the Advanced LIGO and Virgo full O3 data set. Search results are presented for gravitational waves produced by cosmic string loop features such as cusps, kinks and, for the first time, kink-kink collisions.cA template-based search for short-duration transient signals does not yield a detection. We also use the stochastic gravitational-wave background energy density upper limits derived from the O3 data to constrain the cosmic string tension, $Gμ$, as a function of the number of kinks, or the number of cusps, for two cosmic string loop distribution models.cAdditionally, we develop and test a third model which interpolates between these two models. Our results improve upon the previous LIGO-Virgo constraints on $Gμ$ by one to two orders of magnitude depending on the model which is tested. In particular, for one loop distribution model, we set the most competitive constraints to date, $Gμ\lesssim 4\times 10^{-15}$. △ Less

Submitted 28 January, 2021; originally announced January 2021.

Comments: 20 pages, 10 figures

Report number: LIGO-P2000506

Journal ref: Phys. Rev. Lett. 126, 241102 (2021)

Upper Limits on the Isotropic Gravitational-Wave Background from Advanced LIGO's and Advanced Virgo's Third Observing Run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, A. Adams, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, T. Akutsu, K. M. Aleman, G. Allen, A. Allocca, P. A. Altin , et al. (1566 additional authors not shown)

Abstract: We report results of a search for an isotropic gravitational-wave background (GWB) using data from Advanced LIGO's and Advanced Virgo's third observing run (O3) combined with upper limits from the earlier O1 and O2 runs. Unlike in previous observing runs in the advanced detector era, we include Virgo in the search for the GWB. The results are consistent with uncorrelated noise, and therefore we pl… ▽ More We report results of a search for an isotropic gravitational-wave background (GWB) using data from Advanced LIGO's and Advanced Virgo's third observing run (O3) combined with upper limits from the earlier O1 and O2 runs. Unlike in previous observing runs in the advanced detector era, we include Virgo in the search for the GWB. The results are consistent with uncorrelated noise, and therefore we place upper limits on the strength of the GWB. We find that the dimensionless energy density $Ω_{\rm GW}\leq 5.8\times 10^{-9}$ at the 95% credible level for a flat (frequency-independent) GWB, using a prior which is uniform in the log of the strength of the GWB, with 99% of the sensitivity coming from the band 20-76.6 Hz; $\leq 3.4 \times 10^{-9}$ at 25 Hz for a power-law GWB with a spectral index of 2/3 (consistent with expectations for compact binary coalescences), in the band 20-90.6 Hz; and $\leq 3.9 \times 10^{-10}$ at 25 Hz for a spectral index of 3, in the band 20-291.6 Hz. These upper limits improve over our previous results by a factor of 6.0 for a flat GWB. We also search for a GWB arising from scalar and vector modes, which are predicted by alternative theories of gravity; we place upper limits on the strength of GWBs with these polarizations. We demonstrate that there is no evidence of correlated noise of magnetic origin by performing a Bayesian analysis that allows for the presence of both a GWB and an effective magnetic background arising from geophysical Schumann resonances. We compare our upper limits to a fiducial model for the GWB from the merger of compact binaries. Finally, we combine our results with observations of individual mergers andshow that, at design sensitivity, this joint approach may yield stronger constraints on the merger rate of binary black holes at $z \lesssim 2$ than can be achieved with individually resolved mergers alone. [abridged] △ Less

Submitted 28 January, 2021; originally announced January 2021.

Comments: 25 pages, 7 figures, Abstract abridged for arxiv submission

Report number: LIGO-DCC-P2000314

Journal ref: Phys. Rev. D 104, 022004 (2021)

Diving below the spin-down limit: Constraints on gravitational waves from the energetic young pulsar PSR J0537-6910

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, A. Adams, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, K. M. Aleman, G. Allen, A. Allocca , et al. (1568 additional authors not shown)

Abstract: We present a search for continuous gravitational-wave signals from the young, energetic X-ray pulsar PSR J0537-6910 using data from the second and third observing runs of LIGO and Virgo. The search is enabled by a contemporaneous timing ephemeris obtained using NICER data. The NICER ephemeris has also been extended through 2020 October and includes three new glitches. PSR J0537-6910 has the larges… ▽ More We present a search for continuous gravitational-wave signals from the young, energetic X-ray pulsar PSR J0537-6910 using data from the second and third observing runs of LIGO and Virgo. The search is enabled by a contemporaneous timing ephemeris obtained using NICER data. The NICER ephemeris has also been extended through 2020 October and includes three new glitches. PSR J0537-6910 has the largest spin-down luminosity of any pulsar and is highly active with regards to glitches. Analyses of its long-term and inter-glitch braking indices provided intriguing evidence that its spin-down energy budget may include gravitational-wave emission from a time-varying mass quadrupole moment. Its 62 Hz rotation frequency also puts its possible gravitational-wave emission in the most sensitive band of LIGO/Virgo detectors. Motivated by these considerations, we search for gravitational-wave emission at both once and twice the rotation frequency. We find no signal, however, and report our upper limits. Assuming a rigidly rotating triaxial star, our constraints reach below the gravitational-wave spin-down limit for this star for the first time by more than a factor of two and limit gravitational waves from the $l=m=2$ mode to account for less than 14% of the spin-down energy budget. The fiducial equatorial ellipticity is limited to less than about 3e-5, which is the third best constraint for any young pulsar. △ Less

Submitted 10 June, 2021; v1 submitted 23 December, 2020; originally announced December 2020.

Comments: 21 pages, 5 figures, published in ApJL

Report number: LIGO-P2000407

All-sky search in early O3 LIGO data for continuous gravitational-wave signals from unknown neutron stars in binary systems

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, A. Adams, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, K. M. Aleman, G. Allen, A. Allocca, P. A. Altin, A. Amato , et al. (1347 additional authors not shown)

Abstract: Rapidly spinning neutron stars are promising sources of persistent, continuous gravitational waves. Detecting such a signal would allow probing of the physical properties of matter under extreme conditions. A significant fraction of the known pulsar population belongs to binary systems. Searching for unknown neutron stars in binary systems requires specialized algorithms to address unknown orbital… ▽ More Rapidly spinning neutron stars are promising sources of persistent, continuous gravitational waves. Detecting such a signal would allow probing of the physical properties of matter under extreme conditions. A significant fraction of the known pulsar population belongs to binary systems. Searching for unknown neutron stars in binary systems requires specialized algorithms to address unknown orbital frequency modulations. We present a search for continuous gravitational waves emitted by neutron stars in binary systems in early data from the third observing run of the Advanced LIGO and Advanced Virgo detectors using the semicoherent, GPU-accelerated, BinarySkyHough pipeline. The search analyzes the most sensitive frequency band of the LIGO detectors, 50 - 300 Hz. Binary orbital parameters are split into four regions, comprising orbital periods of 3 - 45 days and projected semimajor axes of 2 - 40 light-seconds. No detections are reported. We estimate the sensitivity of the search using simulated continuous wave signals, achieving the most sensitive results to date across the analyzed parameter space. △ Less

Submitted 19 March, 2021; v1 submitted 22 December, 2020; originally announced December 2020.

Comments: 23 pages, 12 figures, 7 tables

Report number: LIGO-P2000298

Journal ref: Phys. Rev. D 103, 064017 (2021)

Population Properties of Compact Objects from the Second LIGO-Virgo Gravitational-Wave Transient Catalog

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, A. Adams, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, P. A. Altin, A. Amato, S. Anand, A. Ananyeva , et al. (1316 additional authors not shown)

Abstract: We report on the population of the 47 compact binary mergers detected with a false-alarm rate 1/yr in the second LIGO--Virgo Gravitational-Wave Transient Catalog, GWTC-2. We observe several characteristics of the merging binary black hole (BBH) population not discernible until now. First, we find that the primary mass spectrum contains structure beyond a power-law with a sharp high-mass cut-off; i… ▽ More We report on the population of the 47 compact binary mergers detected with a false-alarm rate 1/yr in the second LIGO--Virgo Gravitational-Wave Transient Catalog, GWTC-2. We observe several characteristics of the merging binary black hole (BBH) population not discernible until now. First, we find that the primary mass spectrum contains structure beyond a power-law with a sharp high-mass cut-off; it is more consistent with a broken power law with a break at $39.7^{+20.3}_{-9.1}\,M_\odot$, or a power law with a Gaussian feature peaking at $33.1^{+4.0}_{-5.6}\,M_\odot$ (90\% credible interval). While the primary mass distribution must extend to $\sim65\,M_\odot$ or beyond, only $2.9^{+3.5}_{1.7}\%$ of systems have primary masses greater than $45\,M_\odot$. Second, we find that a fraction of BBH systems have component spins misaligned with the orbital angular momentum, giving rise to precession of the orbital plane. Moreover, 12% to 44% of BBH systems have spins tilted by more than $90^\circ$, giving rise to a negative effective inspiral spin parameter $χ_\mathrm{eff}$. Under the assumption that such systems can only be formed by dynamical interactions, we infer that between 25% and 93% of BBH with non-vanishing $|χ_\mathrm{eff}| > 0.01$ are dynamically assembled. Third, we estimate merger rates, finding $\mathcal{R}_\text{BBH} = 23.9^{+14.3}_{8.6}$ Gpc$^{-3}$ yr$^{-1}$ for BBH and $\mathcal{R}_\text{BNS}= 320^{+490}_{-240}$ Gpc$^{-3}$ yr$^{-1}$ for binary neutron stars. We find that the BBH rate likely increases with redshift ($85\%$ credibility), but not faster than the star-formation rate ($86\%$ credibility). Additionally, we examine recent exceptional events in the context of our population models, finding that the asymmetric masses of GW190412 and the high component masses of GW190521 are consistent with our models, but the low secondary mass of GW190814 makes it an outlier. △ Less

Submitted 25 February, 2021; v1 submitted 27 October, 2020; originally announced October 2020.

Comments: 53 pages, including 24 pages main text, 18 pages appendix, 30 figures

Report number: LIGO-P2000077

Tests of General Relativity with Binary Black Holes from the second LIGO-Virgo Gravitational-Wave Transient Catalog

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, A. Adams, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, P. A. Altin, A. Amato, S. Anand, A. Ananyeva , et al. (1322 additional authors not shown)

Abstract: Gravitational waves enable tests of general relativity in the highly dynamical and strong-field regime. Using events detected by LIGO-Virgo up to 1 October 2019, we evaluate the consistency of the data with predictions from the theory. We first establish that residuals from the best-fit waveform are consistent with detector noise, and that the low- and high-frequency parts of the signals are in ag… ▽ More Gravitational waves enable tests of general relativity in the highly dynamical and strong-field regime. Using events detected by LIGO-Virgo up to 1 October 2019, we evaluate the consistency of the data with predictions from the theory. We first establish that residuals from the best-fit waveform are consistent with detector noise, and that the low- and high-frequency parts of the signals are in agreement. We then consider parametrized modifications to the waveform by varying post-Newtonian and phenomenological coefficients, improving past constraints by factors of ${\sim}2$; we also find consistency with Kerr black holes when we specifically target signatures of the spin-induced quadrupole moment. Looking for gravitational-wave dispersion, we tighten constraints on Lorentz-violating coefficients by a factor of ${\sim}2.6$ and bound the mass of the graviton to $m_g \leq 1.76 \times 10^{-23} \mathrm{eV}/c^2$ with 90% credibility. We also analyze the properties of the merger remnants by measuring ringdown frequencies and damping times, constraining fractional deviations away from the Kerr frequency to $δ\hat{f}_{220} = 0.03^{+0.38}_{-0.35}$ for the fundamental quadrupolar mode, and $δ\hat{f}_{221} = 0.04^{+0.27}_{-0.32}$ for the first overtone; additionally, we find no evidence for postmerger echoes. Finally, we determine that our data are consistent with tensorial polarizations through a template-independent method. When possible, we assess the validity of general relativity based on collections of events analyzed jointly. We find no evidence for new physics beyond general relativity, for black hole mimickers, or for any unaccounted systematics. △ Less

Submitted 16 June, 2021; v1 submitted 27 October, 2020; originally announced October 2020.

Comments: 24 pages + appendices, 19 figures; journal version

Report number: LIGO-P2000091

Journal ref: Phys. Rev. D 103, 122002 (2021)

GWTC-2: Compact Binary Coalescences Observed by LIGO and Virgo During the First Half of the Third Observing Run

Authors: R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, A. Adams, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, S. Akcay, G. Allen, A. Allocca, P. A. Altin, A. Amato, S. Anand, A. Ananyeva, S. B. Anderson , et al. (1327 additional authors not shown)

Abstract: We report on gravitational wave discoveries from compact binary coalescences detected by Advanced LIGO and Advanced Virgo in the first half of the third observing run (O3a) between 1 April 2019 15:00 UTC and 1 October 2019 15:00. By imposing a false-alarm-rate threshold of two per year in each of the four search pipelines that constitute our search, we present 39 candidate gravitational wave event… ▽ More We report on gravitational wave discoveries from compact binary coalescences detected by Advanced LIGO and Advanced Virgo in the first half of the third observing run (O3a) between 1 April 2019 15:00 UTC and 1 October 2019 15:00. By imposing a false-alarm-rate threshold of two per year in each of the four search pipelines that constitute our search, we present 39 candidate gravitational wave events. At this threshold, we expect a contamination fraction of less than 10%. Of these, 26 candidate events were reported previously in near real-time through GCN Notices and Circulars; 13 are reported here for the first time. The catalog contains events whose sources are black hole binary mergers up to a redshift of ~0.8, as well as events whose components could not be unambiguously identified as black holes or neutron stars. For the latter group, we are unable to determine the nature based on estimates of the component masses and spins from gravitational wave data alone. The range of candidate events which are unambiguously identified as binary black holes (both objects $\geq 3~M_\odot$) is increased compared to GWTC-1, with total masses from $\sim 14~M_\odot$ for GW190924_021846 to $\sim 150~M_\odot$ for GW190521. For the first time, this catalog includes binary systems with significantly asymmetric mass ratios, which had not been observed in data taken before April 2019. We also find that 11 of the 39 events detected since April 2019 have positive effective inspiral spins under our default prior (at 90% credibility), while none exhibit negative effective inspiral spin. Given the increased sensitivity of Advanced LIGO and Advanced Virgo, the detection of 39 candidate events in ~26 weeks of data (~1.5 per week) is consistent with GWTC-1. △ Less

Submitted 8 March, 2021; v1 submitted 27 October, 2020; originally announced October 2020.

Comments: This version updates with minor revisions to typographical errors. We would also like to call attention to the updated parameter estimation samples data release here: https://dcc.ligo.org/LIGO-P2000223/public

Report number: P2000061

Journal ref: Phys. Rev. X 11, 021053 (2021)

The Pencil Code, a modular MPI code for partial differential equations and particles: multipurpose and multiuser-maintained

Authors: A. Brandenburg, A. Johansen, P. A. Bourdin, W. Dobler, W. Lyra, M. Rheinhardt, S. Bingert, N. E. L. Haugen, A. Mee, F. Gent, N. Babkovskaia, C. -C. Yang, T. Heinemann, B. Dintrans, D. Mitra, S. Candelaresi, J. Warnecke, P. J. Käpylä, A. Schreiber, P. Chatterjee, M. J. Käpylä, X. -Y. Li, J. Krüger, J. R. Aarnes, G. R. Sarson , et al. (12 additional authors not shown)

Abstract: The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the… ▽ More The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the framework is general and can also be applied to situations not related to hydrodynamics or even PDEs, for example when just the message passing interface or input/output strategies of the code are to be used. The code can also evolve Lagrangian (inertial and noninertial) particles, their coagulation and condensation, as well as their interaction with the fluid. △ Less

Submitted 17 September, 2020; originally announced September 2020.

Comments: 7 pages, submitted to the Journal for Open Source Software (JOSS)

Report number: NORDITA-2020-087

Journal ref: Journal of Open Source Software 6, 2807 (2021)

Gravitational-wave constraints on the equatorial ellipticity of millisecond pulsars

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, A. Adams, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, P. A. Altin, A. Amato, S. Anand, A. Ananyeva , et al. (1311 additional authors not shown)

Abstract: We present a search for continuous gravitational waves from five radio pulsars, comprising three recycled pulsars (PSR J0437-4715, PSR J0711-6830, and PSR J0737-3039A) and two young pulsars: the Crab pulsar (J0534+2200) and the Vela pulsar (J0835-4510). We use data from the third observing run of Advanced LIGO and Virgo combined with data from their first and second observing runs. For the first t… ▽ More We present a search for continuous gravitational waves from five radio pulsars, comprising three recycled pulsars (PSR J0437-4715, PSR J0711-6830, and PSR J0737-3039A) and two young pulsars: the Crab pulsar (J0534+2200) and the Vela pulsar (J0835-4510). We use data from the third observing run of Advanced LIGO and Virgo combined with data from their first and second observing runs. For the first time we are able to match (for PSR J0437-4715) or surpass (for PSR J0711-6830) the indirect limits on gravitational-wave emission from recycled pulsars inferred from their observed spin-downs, and constrain their equatorial ellipticities to be less than $10^{-8}$. For each of the five pulsars, we perform targeted searches that assume a tight coupling between the gravitational-wave and electromagnetic signal phase evolution. We also present constraints on PSR J0711-6830, the Crab pulsar and the Vela pulsar from a search that relaxes this assumption, allowing the gravitational-wave signal to vary from the electromagnetic expectation within a narrow band of frequencies and frequency derivatives. △ Less

Submitted 13 October, 2020; v1 submitted 28 July, 2020; originally announced July 2020.

Comments: 23 pages, 6 figures, 3 tables

Report number: LIGO-P2000029

Journal ref: 2020 ApJL 902 L21

Diffusion of large-scale magnetic fields by reconnection in MHD turbulence

Authors: R. Santos-Lima, G. Guerrero, E. M. de Gouveia Dal Pino, A. Lazarian

Abstract: The rate of magnetic field diffusion plays an essential role in several astrophysical plasma processes. It has been demonstrated that the omnipresent turbulence in astrophysical media induces fast magnetic reconnection, which consequently leads to large-scale magnetic flux diffusion at a rate independent of the plasma microphysics. This process is called ``reconnection diffusion'' (RD) and allows… ▽ More The rate of magnetic field diffusion plays an essential role in several astrophysical plasma processes. It has been demonstrated that the omnipresent turbulence in astrophysical media induces fast magnetic reconnection, which consequently leads to large-scale magnetic flux diffusion at a rate independent of the plasma microphysics. This process is called ``reconnection diffusion'' (RD) and allows for the diffusion of fields which are dynamically important. The current theory describing RD is based on incompressible magnetohydrodynamic (MHD) turbulence. In this work, we have tested quantitatively the predictions of the RD theory when magnetic forces are dominant in the turbulence dynamics (Alfvénic Mach number $M_A < 1$). We employed the \textsc{Pencil Code} to perform numerical simulations of forced MHD turbulence, extracting the values of the diffusion coefficient $η_{RD}$ using the Test-Field method. Our results are consistent with the RD theory ($η_{RD} \sim M_A^{3}$ for $M_A < 1$) when turbulence approaches the incompressible limit (sonic Mach number $M_S \lesssim 0.02$), while for larger $M_S$ the diffusion is faster ($η_{RD} \sim M_A^{2}$). This work shows for the first time simulations of compressible MHD turbulence with the suppression of the cascade in the direction parallel to the mean magnetic field, which is consistent with incompressible weak turbulence theory. We also verified that in our simulations the energy cascading time does not follow the scaling with $M_A$ predicted for the weak regime, in contradiction with the RD theory assumption. Our results generally support and expand the RD theory predictions. △ Less

Submitted 15 February, 2021; v1 submitted 15 May, 2020; originally announced May 2020.

Comments: 22 pages, 13 figures. Accepted for publication in MNRAS

Is the primary CoRoT target HD43587 under a Maunder minimum phase?

Authors: Rafael Ramon Ferreira, Rafaella Barbosa, Matthieu Castro, Gustavo Guerrero, Leandro de Almeida, Patrick Boumier, José Dias do Nascimento Jr

Abstract: One of the most enigmatic phenomena related to solar activity is the so-called Maunder minimum phase. It consists of the lowest sunspot's counting ever registered for the Sun, and never so far confirmed for other stars. Since the spectroscopic observations of stellar activity at the Mount Wilson Observatory, the solar analog HD43587 has shown a very low and apparently invariant activity level, whi… ▽ More One of the most enigmatic phenomena related to solar activity is the so-called Maunder minimum phase. It consists of the lowest sunspot's counting ever registered for the Sun, and never so far confirmed for other stars. Since the spectroscopic observations of stellar activity at the Mount Wilson Observatory, the solar analog HD43587 has shown a very low and apparently invariant activity level, which makes it a Maunder minimum candidate. We analyze the chromospheric activity evolution of HD43587 and its evolutive status, intending to unravel the reasons for this low and flat activity. We use an activity measurements dataset available in the literature, and compute the activity S-index from HARPS and NARVAL spectra, to infer a cycle period. Besides, we analyze the CoRoT light-curve of HD43587, and apply gyrochronology and activity calibrations, to determine its rotation period. Finally, based on an evolutionary model and the inferred rotation period, we use the EULAG-MHD code to perform global MHD simulations of HD43587 with the aim of getting some insights about its dynamo process. We confirm the almost flat activity profile, with a cycle period $P_{\mathrm{cyc}} = 10.44\pm3.03$ yrs deduced from the S-index time series, and a long-term trend that might be a period larger than 50 yrs. It was impossible to define a rotation period from the light-curve, however gyrochronology and activity calibrations allow to infer an indirect estimate of $\overline{P}_{\mathrm{rot}} = 22.6 \pm 1.9$ d. Furthermore, the MHD simulations confirm an oscillatory dynamo with a cycle period in good agreement with the observations and a low level of surface magnetic activity. We conclude that this object might be experiencing a "natural" decrease in magnetic activity as a consequence of its age. Nevertheless, the possibility that HD43587 is in a Maunder minimum phase cannot be ruled out. △ Less

Submitted 13 August, 2020; v1 submitted 11 May, 2020; originally announced May 2020.

Comments: 12 pages, 10 figures, 2 tables, has been published in Astronomy & Astrophysics, Volume 640

Journal ref: A&A, 640, A46 (2020)

Global simulations of stellar dynamos

Authors: G Guerrero

Abstract: The dynamo mechanism, responsible for the solar magnetic activity, is still an open problem in astrophysics. Different theories proposed to explain such phenomena have failed in reproducing the observational properties of the solar magnetism. Thus, ab-initio computational modeling of the convective dynamo in a spherical shell turns out as the best alternative to tackle this problem. In this work w… ▽ More The dynamo mechanism, responsible for the solar magnetic activity, is still an open problem in astrophysics. Different theories proposed to explain such phenomena have failed in reproducing the observational properties of the solar magnetism. Thus, ab-initio computational modeling of the convective dynamo in a spherical shell turns out as the best alternative to tackle this problem. In this work we review the efforts performed in global simulations over the past decades. Regarding the development and sustain of mean-flows, as well as mean magnetic field, we discuss the points of agreement and divergence between the different modeling strategies. Special attention is given to the implicit large-eddy simulations performed with the EULAG-MHD code. △ Less

Submitted 28 January, 2020; originally announced January 2020.

Comments: 21 pages, 9 figures, proceedings of IAU symposium 354, Solar and Stellar Magnetic Fields: Origins and Manifestations

3D MHD Modeling of the Impact of Subsurface Stratification on the Solar Dynamo

Authors: Andrey M. Stejko, Gustavo Guerrero, Alexander G. Kosovichev, Piotr K. Smolarkiewicz

Abstract: Various models of solar subsurface stratification are tested in the global EULAG-MHD solver to simulate diverse regimes of near-surface convective transport. Sub- and superadiabacity are altered at the surface of the model ($ r > 0.95~R_{\odot}$) to either suppress or enhance convective flow speeds in an effort to investigate the impact of the near-surface layer on global dynamics. A major consequ… ▽ More Various models of solar subsurface stratification are tested in the global EULAG-MHD solver to simulate diverse regimes of near-surface convective transport. Sub- and superadiabacity are altered at the surface of the model ($ r > 0.95~R_{\odot}$) to either suppress or enhance convective flow speeds in an effort to investigate the impact of the near-surface layer on global dynamics. A major consequence of increasing surface convection rates appears to be a significant alteration of the distribution of angular momentum, especially below the tachocline where the rotational frequency predominantly increases at higher latitudes. These hydrodynamic changes correspond to large shifts in the development of the current helicity in this stable layer ($r<0.72R_{\odot}$), significantly altering its impact on the generation of poloidal and toroidal fields at the tachocline and below, acting as a major contributor towards transitions in the dynamo cycle. The enhanced near-surface flow speed manifests in a global shift of the toroidal field ($B_φ$) in the butterfly diagram - from a North-South symmetric pattern to a staggered anti-symmetric emergence. △ Less

Submitted 21 November, 2019; originally announced November 2019.

Global simulations of Tayler instability in stellar interiors: The stabilizing effect of gravity

Authors: G. Guerrero, F. Del Sordo, A. Bonanno, P. K. Smolarkiewicz

Abstract: Unveiling the evolution of toroidal field instability, known as Tayler instability, is essential to understand the strength and topology of the magnetic fields observed in early-type stars, in the core of the red giants, or in any stellar radiative zone. We want to study the non-linear evolution of the instability of a toroidal field stored in a stably stratified layer, in spherical symmetry and i… ▽ More Unveiling the evolution of toroidal field instability, known as Tayler instability, is essential to understand the strength and topology of the magnetic fields observed in early-type stars, in the core of the red giants, or in any stellar radiative zone. We want to study the non-linear evolution of the instability of a toroidal field stored in a stably stratified layer, in spherical symmetry and in the absence of rotation. In particular, we intend to quantify the suppression of the instability as a function of the Brunt-Väisäla ($ω_{\rm BV}$) and the Alfvén ($ω_{\rm A}$) frequencies. We use the MHD equations as implemented in the anelastic approximation in the EULAG-MHD code and perform a large series of numerical simulations of the instability exploring the parameter space for the $ω_{\rm BV}$ and $ω_{\rm A}$. We show that beyond a critical value gravity strongly suppress the instability, in agreement with the linear analysis. The intensity of the initial field also plays an important role: weaker fields show much slower growth rates. Moreover, in the case of very low gravity, the fastest growing modes have a large characteristic radial scale, at variance with the case of strong gravity, where the instability is characterized by horizontal displacements. Our results illustrate that the anelastic approximation can efficiently describe the evolution of toroidal field instability in stellar interiors. The suppression of the instability as a consequence of increasing values of $ω_{\rm BV}$ might play a role to explain the magnetic desert in Ap/Bp stars since weak fields are only marginally unstable in the case of strong gravity. △ Less

Submitted 7 October, 2019; v1 submitted 6 September, 2019; originally announced September 2019.

Comments: 11 pages, 11 figures, Accepted for publication in MNRAS

What sets the magnetic field strength and cycle period in solar-type stars?

Authors: G. Guerrero, B. Zaire, P. K. Smolarkiewicz, E. M. de Gouveia Dal Pino, A. G. Kosovichev, N. N. Mansour

Abstract: Two fundamental properties of stellar magnetic fields have been determined by observations for solar-like stars with different Rossby numbers (Ro), namely, the magnetic field strength and the magnetic cycle period. The field strength exhibits two regimes: 1) for fast rotation it is independent of Ro, 2) for slow rotation it decays with Ro following a power law. For the magnetic cycle period two re… ▽ More Two fundamental properties of stellar magnetic fields have been determined by observations for solar-like stars with different Rossby numbers (Ro), namely, the magnetic field strength and the magnetic cycle period. The field strength exhibits two regimes: 1) for fast rotation it is independent of Ro, 2) for slow rotation it decays with Ro following a power law. For the magnetic cycle period two regimes of activity, the active and inactive branches, also have been identified. For both of them, the longer the rotation period, the longer the activity cycle. Using global dynamo simulations of solar like stars with Rossby numbers between ~0.4 and ~2, this paper explores the relevance of rotational shear layers in determining these observational properties. Our results, consistent with non-linear alpha^2-Omega dynamos, show that the total magnetic field strength is independent of the rotation period. Yet at surface levels, the origin of the magnetic field is determined by Ro. While for Ro<1 it is generated in the convection zone, for Ro>1 strong toroidal fields are generated at the tachocline and rapidly emerge towards the surface. In agreement with the observations, the magnetic cycle period increases with the rotational period. However, a bifurcation is observed for Ro~1, separating a regime where oscillatory dynamos operate mainly in the convection zone, from the regime where the tachocline has a predominant role. In the latter the cycles are believed to result from the periodic energy exchange between the dynamo and the magneto-shear instabilities developing in the tachocline and the radiative interior. △ Less

Submitted 13 May, 2019; v1 submitted 18 October, 2018; originally announced October 2018.

Comments: 43 pages, 14 figures, accepted for publication in The Astrophysical Journal

Magnetic field generation in PMS stars with and without radiative core

Authors: Bonnie Zaire, Gustavo Guerrero, Alexander Kosovichev, Piotr Smolarkiewicz, Natalia Landin

Abstract: Recent observations of the magnetic field in pre-main sequence stars suggest that the magnetic field topology changes as a function of age. The presence of a tachocline could be an important factor in the development of magnetic field with higher multipolar modes. In this work we performed MHD simulations using the EULAG-MHD code to study the magnetic field generation and evolution in models that… ▽ More Recent observations of the magnetic field in pre-main sequence stars suggest that the magnetic field topology changes as a function of age. The presence of a tachocline could be an important factor in the development of magnetic field with higher multipolar modes. In this work we performed MHD simulations using the EULAG-MHD code to study the magnetic field generation and evolution in models that mimic stars at two evolutionary stages. The stratification for both stellar phases was computed by fitting stellar structure profiles obtained with the ATON stellar evolution code. The first stage is at $1.1 \mathrm{Myr}$, when the star is completely convective. The second stage is at $14 \mathrm{Myrs}$, when the star is partly convective, with a radiative core developed up to $30\%$ of the stellar radius. In this proceedings we present a preliminary analysis of the resulting mean-flows and magnetic field. The mean-flow analysis shown that the star rotate almost rigidly on the fully convective phase, whereas at the partially convective phase there is differential rotation with conical contours of iso-rotation. As for the mean magnetic field both simulations show similarities with respect to the field evolution. However, the topology of the magnetic field is different. △ Less

Submitted 6 November, 2017; originally announced November 2017.

Comments: 8 pages, 10 figures, Symposium

Journal ref: Proceedings of the International Astronomical Union, 12(S328), 30-37 (2016)

Meridional Circulation Dynamics in a Cyclic Convective Dynamo

Authors: D. Passos, M. Miesch, G. Guerrero, P. Charbonneau

Abstract: Surface observations indicate that the speed of the solar meridional circulation in the photosphere varies in anti-phase with the solar cycle. The current explanation for the source of this variation is that inflows into active regions alter the global surface pattern of the meridional circulation. When these localized inflows are integrated over a full hemisphere, they contribute to the slow down… ▽ More Surface observations indicate that the speed of the solar meridional circulation in the photosphere varies in anti-phase with the solar cycle. The current explanation for the source of this variation is that inflows into active regions alter the global surface pattern of the meridional circulation. When these localized inflows are integrated over a full hemisphere, they contribute to the slow down of the axisymmetric poleward horizontal component. The behavior of this large scale flow deep inside the convection zone remains largely unknown. Present helioseismic techniques are not sensitive enough to capture the dynamics of this weak large scale flow. Moreover, the large time of integration needed to map the meridional circulation inside the convection zone, also masks some of the possible dynamics on shorter timescales. In this work we examine the dynamics of the meridional circulation that emerges from a 3D MHD global simulation of the solar convection zone. Our aim is to assess and quantify the behavior of meridional circulation deep inside the convection zone, where the cyclic large-scale magnetic field can reach considerable strength. Our analyses indicate that the meridional circulation morphology and amplitude are both highly influenced by the magnetic field, via the impact of magnetic torques on the global angular momentum distribution. A dynamic feature induced by these magnetic torques is the development of a prominent upward flow at mid latitudes in the lower convection zone that occurs near the equatorward edge of the toroidal bands and that peaks during cycle maximum. Globally, the dynamo-generated large-scale magnetic field drives variations in the meridional flow, in stark contrast to the conventional kinematic flux transport view of the magnetic field being advected passively by the flow. △ Less

Submitted 18 July, 2017; v1 submitted 8 February, 2017; originally announced February 2017.

Comments: 20 pages, 13 figures, Accepted in A&A

Journal ref: A&A 607, A120 (2017)

3D Global MHD Modeling of Solar Subsurface Convection

Authors: Andrey M. Stejko, Gustavo Guerrero, Alexander G. Kosovichev, Piotr K. Smolarkiewicz

Abstract: Various models of solar sub-surface stratification are tested in the global EULAG-MHD solver to simulate diverse regimes of near-surface convective transport. Sub- and superadiabacity are altered at the surface of the model ($ r > 0.95~R_{\odot}$) to either suppress or enhance convective flow speeds in an effort to investigate the impact of the near-surface layer on global dynamics. A major conseq… ▽ More Various models of solar sub-surface stratification are tested in the global EULAG-MHD solver to simulate diverse regimes of near-surface convective transport. Sub- and superadiabacity are altered at the surface of the model ($ r > 0.95~R_{\odot}$) to either suppress or enhance convective flow speeds in an effort to investigate the impact of the near-surface layer on global dynamics. A major consequence of increasing surface convection rates appears to be a significant alteration of the distribution of angular momentum, especially below the tachocline where the rotational frequency predominantly increases at higher latitudes. These hydrodynamic changes correspond to large shifts in the development of the current helicity in this stable layer ($r<0.72R_{\odot}$), significantly altering its impact on the generation of poloidal and toroidal fields at the tachocline and below, acting as a major contributor towards transitions in the dynamo cycle. The enhanced near-surface flow speed manifests in a global shift of the toroidal field ($B_φ$) in the butterfly diagram -- from a North-South symmetric pattern to a staggered anti-symmetric emergence. △ Less

Submitted 18 November, 2019; v1 submitted 29 January, 2017; originally announced January 2017.

Comments: The wrong article was accidentally submitted, it will be resubmitted in its full form along with a link to the published article in ApJ, so that it can be more easily linked to NASA ADS

Preliminary design of the INPE's Solar Vector Magnetograph

Authors: L. E. A. Vieira, A. L. Clúa de Gonzalez, A. Dal Lago, C. Wrasse, E. Echer, F. L. Guarnieri, F. Reis Cardoso, G. Guerrero, J. Rezende Costa, J. Palacios, L. Balmaceda, L. Ribeiro Alves, L. da Silva, L. L. Costa, M. Sampaio, M. C. Rabello Soares, M. Barbosa, M. Domingues, N. Rigozo, O. Mendes Jr., P. Jauer, R. Dallaqua, R. H. Branco, T. Stekel, W. Gonzalez , et al. (1 additional authors not shown)

Abstract: We describe the preliminary design of a magnetograph and visible-light imager instrument to study the solar dynamo processes through observations of the solar surface magnetic field distribution. The instrument will provide measurements of the vector magnetic field and of the line-of-sight velocity in the solar photosphere. As the magnetic field anchored at the solar surface produces most of the s… ▽ More We describe the preliminary design of a magnetograph and visible-light imager instrument to study the solar dynamo processes through observations of the solar surface magnetic field distribution. The instrument will provide measurements of the vector magnetic field and of the line-of-sight velocity in the solar photosphere. As the magnetic field anchored at the solar surface produces most of the structures and energetic events in the upper solar atmosphere and significantly influences the heliosphere, the development of this instrument plays an important role in reaching the scientific goals of The Atmospheric and Space Science Coordination (CEA) at the Brazilian National Institute for Space Research (INPE). In particular, the CEA's space weather program will benefit most from the development of this technology. We expect that this project will be the starting point to establish a strong research program on Solar Physics in Brazil. Our main aim is acquiring progressively the know-how to build state-of-art solar vector magnetograph and visible-light imagers for space-based platforms to contribute to the efforts of the solar-terrestrial physics community to address the main unanswered questions on how our nearby Star works. △ Less

Submitted 4 September, 2016; originally announced September 2016.

Comments: 5 pages, 1 figure

Journal ref: Proceedings IAU Symposium No. 305, 2014

Understanding solar torsional oscillations from global dynamo models

Authors: G. Guerrero, P. K. Smolarkiewicz, E. M. de Gouveia Dal Pino, A. G. Kosovichev, N. N. Mansour

Abstract: The phenomenon of solar "torsional oscillations" (TO) represents migratory zonal flows associated with the solar cycle. These flows are observed on the solar surface and, according to helioseismology, extend through the convection zone. We study the origin of the TO using results from a global MHD simulation of the solar interior that reproduces several of the observed characteristics of the mean-… ▽ More The phenomenon of solar "torsional oscillations" (TO) represents migratory zonal flows associated with the solar cycle. These flows are observed on the solar surface and, according to helioseismology, extend through the convection zone. We study the origin of the TO using results from a global MHD simulation of the solar interior that reproduces several of the observed characteristics of the mean-flows and magnetic fields. Our results indicate that the magnetic tension (MT) in the tachocline region is a key factor for the periodic changes in the angular momentum transport that causes the TO. The torque induced by the MT at the base of the convection zone is positive at the poles and negative at the equator. A rising MT torque at higher latitudes causes the poles to speed-up, whereas a declining negative MT torque at the lower latitudes causes the equator to slow-down. These changes in the zonal flows propagate through the convection zone up to the surface. Additionally, our results suggest that it is the magnetic field at the tachocline that modulates the amplitude of the surface meridional flow rather than the opposite as assumed by flux-transport dynamo models of the solar cycle. △ Less

Submitted 26 August, 2016; v1 submitted 7 August, 2016; originally announced August 2016.

Comments: 16 pages, 4 figures, accepted for publication in The Astrophysical Journal Letters

Journal ref: ApJL, 828, L3, 2016

On the role of tachoclines in solar and stellar dynamos

Authors: G. Guerrero, P. K. Smolarkiewicz, E. M. de Gouveia Dal Pino, A. G. Kosovichev, N. N. Mansour

Abstract: Rotational shear layers at the boundary between radiative and convective zones, tachoclines, play a key role in the process of magnetic field generation in solar-like stars. We present two sets of global simulations of rotating turbulent convection and dynamo. The first set considers a stellar convective envelope only; the second one, aiming at the formation of a tachocline, considers also the upp… ▽ More Rotational shear layers at the boundary between radiative and convective zones, tachoclines, play a key role in the process of magnetic field generation in solar-like stars. We present two sets of global simulations of rotating turbulent convection and dynamo. The first set considers a stellar convective envelope only; the second one, aiming at the formation of a tachocline, considers also the upper part of the radiative zone. Our results indicate that the resulting mean-flows and dynamo properties like the growth rate, saturation energy and mode depend on the Rossby (Ro) number. For the first set of models either oscillatory (with ~2 yr period) or steady dynamo solutions are obtained. The models in the second set naturally develop a tachocline which, in turn, leads to the generation of strong mean magnetic field. Since the field is also deposited into the stable deeper layer, its evolutionary time-scale is much longer than in the models without a tachocline. Surprisingly, the magnetic field in the upper turbulent convection zone evolves in the same time scale as the deep field. These models result in either an oscillatory dynamo with ~30 yr period or in a steady dynamo depending on Ro. In terms of the mean-field dynamo coefficients computed using FOSA, the field evolution in the oscillatory models without a tachocline seems to be consistent with dynamo waves propagating according to the Parker-Yoshimura sign rule. In the models with tachoclines the dynamics is more complex involving other transport mechanisms as well as tachocline instabilities. △ Less

Submitted 18 December, 2015; v1 submitted 15 July, 2015; originally announced July 2015.

Comments: 42 pages, 9 figures. Accepted for publication in ApJ

Differential rotation in solar-like stars from global simulations

Authors: G. Guerrero, P. K. Smolarkiewicz, A. G. Kosovichev, N. N. Mansour

Abstract: To explore the physics of large-scale flows in solar-like stars, we perform 3D anelastic simulations of rotating convection for global models with stratification resembling the solar interior. The numerical method is based on an implicit large-eddy simulation approach designed to capture effects from non-resolved small scales. We obtain two regimes of differential rotation, with equatorial zonal f… ▽ More To explore the physics of large-scale flows in solar-like stars, we perform 3D anelastic simulations of rotating convection for global models with stratification resembling the solar interior. The numerical method is based on an implicit large-eddy simulation approach designed to capture effects from non-resolved small scales. We obtain two regimes of differential rotation, with equatorial zonal flows accelerated either in the direction of rotation (solar-like) or in the opposite direction (anti-solar). While the models with the solar-like differential rotation tend to produce multiple cells of meridional circulation, the models with anti-solar differential rotation result in only one or two meridional cells. Our simulations indicate that the rotation and large-scale flow patterns critically depend on the ratio between buoyancy and Coriolis forces. By including a subadiabatic layer at the bottom of the domain, corresponding to the stratification of a radiative zone, we reproduce a layer of strong radial shear similar to the solar tachocline. Similarly, enhanced superadiabaticity at the top results in a near-surface shear layer located mainly at lower latitudes. The models reveal a latitudinal entropy gradient localized at the base of the convection zone and in the stable region, which however does not propagate across the convection zone. In consequence, baroclinicity effects remain small and the rotation iso-contours align in cylinders along the rotation axis. Our results confirm the alignment of large convective cells along the rotation axis in the deep convection zone, and suggest that such "banana-cell" pattern can be hidden beneath the supergranulation layer. △ Less

Submitted 30 October, 2013; originally announced October 2013.

Comments: 46 pages, 13 figures, Accepted for publication in ApJ

Solar differential rotation: hints to reproduce a near-surface shear layer in global simulations

Authors: G. Guerrero, P. K. Smolarkiewicz, A. Kosovichev, N. Mansour

Abstract: Convective turbulent motions in the solar interior, as well as the mean flows resulting from them, determine the evolution of the solar magnetic field. With the aim to get a better understanding of these flows we study anelastic rotating convection in a spherical shell whose stratification resembles that of the solar interior. This study is done through numerical simulations performed with the EUL… ▽ More Convective turbulent motions in the solar interior, as well as the mean flows resulting from them, determine the evolution of the solar magnetic field. With the aim to get a better understanding of these flows we study anelastic rotating convection in a spherical shell whose stratification resembles that of the solar interior. This study is done through numerical simulations performed with the EULAG code. Due to the numerical formulation, these simulations are known as implicit large eddy simulations (ILES), since they intrinsically capture the contribution of, non-resolved, small scales at the same time maximizing the effective Reynolds number. We reproduce some previous results and find a transition between buoyancy and rotation dominated regimes which results in anti-solar or solar like rotation patterns. Even thought the rotation profiles are dominated by Taylor-Proudman columnar rotation, we are able to reproduce the tachocline and a low latitude near-surface shear layer. We find that simulations results depend on the grid resolution as a consequence of a different sub-grid scale contribution. △ Less

Submitted 8 January, 2013; v1 submitted 7 January, 2013; originally announced January 2013.

Comments: Proceedings of IAUS 294 "Solar and Astrophysical Dynamos and Magnetic Activity" Editors A.G. Kosovichev, E.M. de Gouveia Dal Pino and Y.Yan, Cambridge University Press, to appear (2013). 9 pages, 4 figures; corrected typos in Metadata

Turbulent dynamo with advective magnetic helicity flux

Authors: Fabio Del Sordo, Gustavo Guerrero, Axel Brandenburg

Abstract: Many astrophysical bodies harbor magnetic fields that are thought to be sustained by a dynamo process. However, it has been argued that the production of large-scale magnetic fields by mean-field dynamo action is strongly suppressed at large magnetic Reynolds numbers owing to the conservation of magnetic helicity. This phenomenon is known as {\it catastrophic quenching}. Advection of magnetic fiel… ▽ More Many astrophysical bodies harbor magnetic fields that are thought to be sustained by a dynamo process. However, it has been argued that the production of large-scale magnetic fields by mean-field dynamo action is strongly suppressed at large magnetic Reynolds numbers owing to the conservation of magnetic helicity. This phenomenon is known as {\it catastrophic quenching}. Advection of magnetic fields by stellar and galactic winds toward the outer boundaries and away from the dynamo is expected to alleviate such quenching. Here we explore the relative roles played by advective and turbulent--diffusive fluxes of magnetic helicity in the dynamo. In particular, we study how the dynamo is affected by advection. We do this by performing direct numerical simulations of a turbulent dynamo of $α^2$ type driven by forced turbulence in a Cartesian domain in the presence of a flow away from the equator where helicity changes sign. Our results indicate that in the presence of advection, the dynamo, otherwise stationary, becomes oscillatory. We confirm an earlier result for turbulent--diffusive magnetic helicity fluxes that for small magnetic Reynolds numbers ($\Rm\lesssim 100...200$, based on the wavenumber of the energy-carrying eddies) the magnetic helicity flux scales less strongly with magnetic Reynolds number ($\Rm^{-1/2}$) than the term describing magnetic helicity destruction by resistivity ($\Rm^{-1}$). Our new results now suggest that for larger $\Rm$ the former becomes approximately independent of $\Rm$, while the latter falls off more slowly. We show for the first time that both for weak and stronger winds, the magnetic helicity flux term becomes comparable to the resistive term for $\Rm\gtrsim 1000$, which is necessary for alleviating catastrophic quenching. △ Less

Submitted 13 November, 2012; v1 submitted 15 May, 2012; originally announced May 2012.

Comments: 9 pages, 9 figures, accepted for publication in MNRAS

Report number: NORDITA 2012 - 38

Journal ref: Mon. Not. Roy. Astron. Soc. 429, 1686-1694 (2013)

Magnetic Flux Transport by turbulent reconnection in astrophysical flows

Authors: Elisabete M. de Gouveia Dal Pino, Márcia R. M. Leão, Reinaldo Santos-Lima, Gustavo Guerrero, Grzegorz Kowal, Alex Lazarian

Abstract: The role of MHD turbulence in astrophysical environments is still highly debated. An important question that permeates this debate is the transport of magnetic flux. This is particularly important, for instance, in the context of star formation. When clouds collapse gravitationally to form stars, there must be some magnetic flux transport. otherwise the new born stars would have magnetic fields se… ▽ More The role of MHD turbulence in astrophysical environments is still highly debated. An important question that permeates this debate is the transport of magnetic flux. This is particularly important, for instance, in the context of star formation. When clouds collapse gravitationally to form stars, there must be some magnetic flux transport. otherwise the new born stars would have magnetic fields several orders of magnitude larger than the observed ones. Also, the magnetic flux that is dragged in the late stages of the formation of a star can remove all the rotational support from the accretion disk that grows around the protostar. The efficiency of the mechanism which is often invoked to allow the transport of magnetic fields in the different stages of star formation, namely, the ambipolar diffusion, has been lately put in check. We here discuss an alternative mechanism for magnetic flux transport which is based on turbulent fast magnetic reconnection. We review recent results obtained from 3D MHD numerical simulations that indicate that this mechanism is very efficient for decoupling and transport magnetic flux from the inner denser regions to the outskirts of collapsing clouds in the different stages of star formation. We also discuss this mechanism in the context of dynamo processes and speculate that it can play a role both in the solar dynamo and in accretion disk dynamo processes. △ Less

Submitted 20 December, 2011; originally announced December 2011.

Comments: 10 pages, 2 figures, review submitted to Physica Scripta

Magnetic feature tracking, what determines the speed?

Authors: G. Guerrero, M. Rheinhardt, M. Dikpati

Abstract: Recent observations revealed that small magnetic elements abundant at the solar surface move poleward with a velocity which seems to be lower than the plasma velocity $U_θ$. Guerrero et al. (2011) explained this discrepancy as a consequence of diffusive spreading of the magnetic elements due to a positive radial gradient of $|U_θ|$. As the gradient's sign (inferred by local helioseismology) is sti… ▽ More Recent observations revealed that small magnetic elements abundant at the solar surface move poleward with a velocity which seems to be lower than the plasma velocity $U_θ$. Guerrero et al. (2011) explained this discrepancy as a consequence of diffusive spreading of the magnetic elements due to a positive radial gradient of $|U_θ|$. As the gradient's sign (inferred by local helioseismology) is still unclear, cases with a negative gradient are studied in this paper. Under this condition, the velocity of the magnetic tracers turns out to be larger than the plasma velocity, in disagreement with the observations. Alternative mechanisms for explaining them independently are proposed. For the turbulent magnetic pumping it is shown that it has to be unrealistically strong to reconcile the model with the observations. △ Less

Submitted 15 November, 2011; originally announced November 2011.

Comments: 4 pages, 2 figures, to appear in Proceedings of IAU Symp 286, Comparative Magnetic Minima, C. H. Mandrini, eds

Report number: Nordita-2011-098

Cycles and cycle modulations

Authors: Axel Brandenburg, Gustavo Guerrero

Abstract: Some selected concepts for the solar activity cycle are briefly reviewed. Cycle modulations through a stochastic alpha effect are being identified with limited scale separation ratios. Three-dimensional turbulence simulations with helicity and shear are compared at two different scale separation ratios. In both cases the level of fluctuations shows relatively little variation with the dynamo cycle… ▽ More Some selected concepts for the solar activity cycle are briefly reviewed. Cycle modulations through a stochastic alpha effect are being identified with limited scale separation ratios. Three-dimensional turbulence simulations with helicity and shear are compared at two different scale separation ratios. In both cases the level of fluctuations shows relatively little variation with the dynamo cycle. Prospects for a shallow origin of sunspots are discussed in terms of the negative effective magnetic pressure instability. Tilt angles of bipolar active regions are discussed as a consequence of shear rather than the Coriolis force. △ Less

Submitted 14 November, 2011; originally announced November 2011.

Comments: 12 pages, 5 figures, to appear in Proceedings of IAU Symp 286, Comparative Magnetic Minima, C. H. Mandrini, eds

Report number: NORDITA-2011-103

Journal ref: in Comparative Magnetic Minima: Characterizing quiet times in the Sun and stars, ed. C. H. Mandrini & D. F. Webb, Proc. IAU Symp., Vol. 286, pp. 37-48 (2012)

Plasma flow vs. magnetic feature-tracking speeds in the Sun

Authors: G. Guerrero, M. Rheinhardt, A. Brandenburg, M. Dikpati

Abstract: We simulate the magnetic feature tracking (MFT) speed using advective-diffusive transport models in both one and two dimensions. By depositing magnetic bipolar regions at different latitudes at the Sun's surface and following their evolution for a prescribed meridional circulation and magnetic diffusivity profiles, we derive the MFT speed as a function of latitude. We find that in a one dimensiona… ▽ More We simulate the magnetic feature tracking (MFT) speed using advective-diffusive transport models in both one and two dimensions. By depositing magnetic bipolar regions at different latitudes at the Sun's surface and following their evolution for a prescribed meridional circulation and magnetic diffusivity profiles, we derive the MFT speed as a function of latitude. We find that in a one dimensional surface-transport model the simulated MFT speed at the surface is always the same as the meridional flow-speed used as input to the model, but is different in a two-dimensional transport model in the meridional (r,theta) plane. The difference depends on the value of the magnetic diffusivity and on the radial gradient of the latitudinal velocity. We have confirmed our results with two different codes in spherical and Cartesian coordinates. △ Less

Submitted 24 July, 2011; originally announced July 2011.

Comments: 5 pages, 7 figures. Submitted to MNRAS: Letters

Report number: NORDITA-2011-61

Journal ref: Monthly Notices Roy. Astron. Soc. 420, L1-L5 (2012)

Dynamo action and magnetic buoyancy in convection simulations with vertical shear

Authors: G. Guerrero, P. Käpylä

Abstract: A hypothesis for sunspot formation is the buoyant emergence of magnetic flux tubes created by the strong radial shear at the tachocline. In this scenario, the magnetic field has to exceed a threshold value before it becomes buoyant and emerges through the whole convection zone. We follow the evolution of a random seed magnetic field with the aim of study under what conditions it is possible to exc… ▽ More A hypothesis for sunspot formation is the buoyant emergence of magnetic flux tubes created by the strong radial shear at the tachocline. In this scenario, the magnetic field has to exceed a threshold value before it becomes buoyant and emerges through the whole convection zone. We follow the evolution of a random seed magnetic field with the aim of study under what conditions it is possible to excite the dynamo instability and whether the dynamo generated magnetic field becomes buoyantly unstable and emerges to the surface as expected in the flux-tube context. We perform numerical simulations of compressible turbulent convection that include a vertical shear layer. Like the solar tachocline, the shear is located at the interface between convective and stable layers. We find that shear and convection are able to amplify the initial magnetic field and form large-scale elongated magnetic structures. The magnetic field strength depends on several parameters such as the shear amplitude, the thickness and location of the shear layer, and the magnetic Reynolds number ($\Rm$). Whenever the toroidal magnetic field reaches amplitudes greater a threshold value which is close to the equipartition value, it becomes buoyant and rises into the convection zone where it expands and forms mushroom shape structures. Some events of emergence, i.e. those with the largest amplitudes of the initial field, are able to reach the very uppermost layers of the domain. These episodes are able to modify the convective pattern forming either broader convection cells or convective eddies elongated in the direction of the field. However, in none of these events the field preserves its initial structure. △ Less

Submitted 16 February, 2011; originally announced February 2011.

Comments: 13 pages, 10 figures, Submitted to A&A

Reynolds stress and heat flux in spherical shell convection

Authors: P. J. Käpylä, M. J. Mantere, G. Guerrero, A. Brandenburg, P. Chatterjee

Abstract: (abridged) Context. Turbulent fluxes of angular momentum and heat due to rotationally affected convection play a key role in determining differential rotation of stars. Here we perform a systematic comparison between Cartesian and spherical geometries as a function of the rotation rate. Aims. We extend the earlier studies by using spherical wedges to obtain turbulent angular momentum and heat tran… ▽ More (abridged) Context. Turbulent fluxes of angular momentum and heat due to rotationally affected convection play a key role in determining differential rotation of stars. Here we perform a systematic comparison between Cartesian and spherical geometries as a function of the rotation rate. Aims. We extend the earlier studies by using spherical wedges to obtain turbulent angular momentum and heat transport as functions of the rotation rate from stratified convection. We compare results from spherical and Cartesian models in the same parameter regime. In particular, we want to clarify whether the sharp equatorial profile of the horizontal Reynolds stress found in earlier Cartesian models is reproduced in spherical models. Methods. We employ direct numerical simulations of turbulent convection. In order to reach as high spatial resolution as possible in the spherical runs, we model only parts of the latitude and longitude. The rotational influence, measured by the Coriolis number, is varied from zero to roughly seven, which is the regime that is likely to be realised in the solar convection zone. Cartesian simulations are performed in overlapping parameter regimes. Results. For slow rotation we find that the radial and latitudinal turbulent angular momentum fluxes are directed inward and equatorward, respectively. In the rapid rotation regime the radial flux changes sign in contradiction with theory. The latitudinal flux remains mostly equatorward and develops a maximum close to the equator. In Cartesian simulations this peak can be explained by the strong `banana cells'. The latitudinal heat flux is mostly equatorward for slow rotation but changes sign for rapid rotation. The rotation profiles vary from anti-solar (slow equator) for slow and intermediate rotation to solar-like (fast equator) for rapid rotation. The solar-like profiles are dominated by the Taylor--Proudman balance. △ Less

Submitted 6 July, 2011; v1 submitted 6 October, 2010; originally announced October 2010.

Comments: 17 pages, 19 figures, updated to match the published version. Version with higher resolution figures can be found at http://www.helsinki.fi/~kapyla/publ.html

Report number: NORDITA-2010-85

Journal ref: Astron. Astrophys., 531, A162 (2011)

Can catastrophic quenching be alleviated by separating shear and alpha effect?

Authors: Piyali Chatterjee, Axel Brandenburg, Gustavo Guerrero

Abstract: The small-scale magnetic helicity produced as a by-product of the large-scale dynamo is believed to play a major role in dynamo saturation. In a mean-field model the generation of small-scale magnetic helicity can be modelled by using the dynamical quenching formalism. Catastrophic quenching refers to a decrease of the saturation field strength with increasing Reynolds number. It has been suggeste… ▽ More The small-scale magnetic helicity produced as a by-product of the large-scale dynamo is believed to play a major role in dynamo saturation. In a mean-field model the generation of small-scale magnetic helicity can be modelled by using the dynamical quenching formalism. Catastrophic quenching refers to a decrease of the saturation field strength with increasing Reynolds number. It has been suggested that catastrophic quenching only affects the region of non-zero helical turbulence (i.e. where the kinematic alpha operates) and that it is possible to alleviate catastrophic quenching by separating the region of strong shear from the alpha layer. We perform a systematic study of a simple axisymmetric two-layer alpha-omega dynamo in a spherical shell for Reynolds numbers in the range 1 < Rm < 10^5. In the framework of dynamical quenching we show that this may not be the case, suggesting that magnetic helicity fluxes would be necessary. △ Less

Submitted 31 May, 2010; originally announced May 2010.

Comments: 8 pages, 5 figures (Accepted in Geophysical and Astrophysical Fluid Dynamics)

Report number: NORDITA-2010-034

Journal ref: Geophys. Astrophys. Fluid Dyn. 104, 591-599 (2010)

Magnetic helicity fluxes in interface and flux transport dynamos

Authors: Piyali Chatterjee, Gustavo Guerrero, Axel Brandenburg

Abstract: Dynamos in the Sun and other bodies tend to produce magnetic fields that possess magnetic helicity of opposite sign at large and small scales, respectively. The build-up of magnetic helicity at small scales provides an important saturation mechanism. In order to understand the nature of the solar dynamo we need to understand the details of the saturation mechanism in spherical geometry. In particu… ▽ More Dynamos in the Sun and other bodies tend to produce magnetic fields that possess magnetic helicity of opposite sign at large and small scales, respectively. The build-up of magnetic helicity at small scales provides an important saturation mechanism. In order to understand the nature of the solar dynamo we need to understand the details of the saturation mechanism in spherical geometry. In particular, we want to understand the effects of magnetic helicity fluxes from turbulence and meridional circulation. We consider a model with just radial shear confined to a thin layer (tachocline) at the bottom of the convection zone. The kinetic alpha owing to helical turbulence is assumed to be localized in a region above the convection zone. The dynamical quenching formalism is used to describe the build-up of mean magnetic helicity in the model, which results in a magnetic alpha effect that feeds back on the kinetic alpha effect. In some cases we compare with results obtained using a simple algebraic alpha quenching formula. In agreement with earlier findings, the magnetic alpha effect in the dynamical alpha quenching formalism has the opposite sign compared with the kinetic alpha effect and leads to a catastrophic decrease of the saturation field strength with increasing magnetic Reynolds numbers. However, at high latitudes this quenching effect can lead to secondary dynamo waves that propagate poleward due to the opposite sign of alpha. Magnetic helicity fluxes both from turbulent mixing and from meridional circulation alleviate catastrophic quenching. △ Less

Submitted 28 May, 2010; originally announced May 2010.

Comments: 9 pages, 14 figures, submitted to A & A

Report number: NORDITA-2010-35

Journal ref: Astron. Astrophys. 525, A5 (2011)

Shear-driven and diffusive helicity fluxes in alpha-Omega dynamos

Authors: Gustavo Guerrero, Piyali Chatterjee, Axel Brandenburg

Abstract: We present nonlinear mean-field alpha-Omega dynamo simulations in spherical geometry with simplified profiles of kinematic alpha effect and shear. We take magnetic helicity evolution into account by solving a dynamical equation for the magnetic alpha effect. This gives a consistent description of the quenching mechanism in mean-field dynamo models. The main goal of this work is to explore the effe… ▽ More We present nonlinear mean-field alpha-Omega dynamo simulations in spherical geometry with simplified profiles of kinematic alpha effect and shear. We take magnetic helicity evolution into account by solving a dynamical equation for the magnetic alpha effect. This gives a consistent description of the quenching mechanism in mean-field dynamo models. The main goal of this work is to explore the effects of this quenching mechanism in solar-like geometry, and in particular to investigate the role of magnetic helicity fluxes, specifically diffusive and Vishniac-Cho (VC) fluxes, at large magnetic Reynolds numbers (Rm). For models with negative radial shear or positive latitudinal shear, the magnetic alpha effect has predominantly negative (positive) sign in the northern (southern) hemisphere. In the absence of fluxes, we find that the magnetic energy follows an Rm^-1 dependence, as found in previous works. This catastrophic quenching is alleviated in models with diffusive magnetic helicity fluxes resulting in magnetic fields comparable to the equipartition value even for Rm=10^7. On the other hand, models with a shear-driven Vishniac-Cho flux show an increase of the amplitude of the magnetic field with respect to models without fluxes, but only for Rm<10^4. This is mainly a consequence of assuming a vacuum outside the Sun which cannot support a significant VC flux across the boundary. However, in contrast with the diffusive flux, the VC flux modifies the distribution of the magnetic field. In addition, if an ill-determined scaling factor in the expression for the VC flux is large enough, subcritical dynamo action is possible that is driven by the action of shear and the divergence of current helicity flux. △ Less

Submitted 26 May, 2010; originally announced May 2010.

Comments: 12 pages, 12 figures. Submitted to MNRAS

Report number: NORDITA-2010-36

Journal ref: Monthly Notices Roy. Astron. Soc. 409, 1619-1630 (2010)

The role of diffusivity quenching in flux-transport dynamo models

Authors: G. A. Guerrero, M. Dikpati, E. M. de Gouveia Dal Pino

Abstract: In the non-linear phase of a dynamo process, the back-reaction of the magnetic field upon the turbulent motion results in a decrease of the turbulence level and therefore in a suppression of both the magnetic field amplification (the alpha-quenching effect) and the turbulent magnetic diffusivity (the eta-quenching effect). While the former has been widely explored, the effects of eta-quenching i… ▽ More In the non-linear phase of a dynamo process, the back-reaction of the magnetic field upon the turbulent motion results in a decrease of the turbulence level and therefore in a suppression of both the magnetic field amplification (the alpha-quenching effect) and the turbulent magnetic diffusivity (the eta-quenching effect). While the former has been widely explored, the effects of eta-quenching in the magnetic field evolution have rarely been considered. In this work we investigate the role of the suppression of diffusivity in a flux-transport solar dynamo model that also includes a non-linear alpha-quenching term. Our results indicate that, although for alpha-quenching the dependence of the magnetic field amplification with the quenching factor is nearly linear, the magnetic field response to eta-quenching is non-linear and spatially non-uniform. We have found that the magnetic field can be locally amplified in this case, forming long-lived structures whose maximum amplitude can be up to ~2.5 times larger at the tachocline and up to ~2 times larger at the center of the convection zone than in models without quenching. However, this amplification leads to unobservable effects and to a worse distribution of the magnetic field in the butterfly diagram. Since the dynamo cycle period increases when the efficiency of the quenching increases, we have also explored whether the eta-quenching can cause a diffusion-dominated model to drift into an advection-dominated regime. We have found that models undergoing a large suppression of the diffusivity produce a strong segregation of magnetic fields that may lead to unsteady dynamo-oscillations. On the other hand, an initially diffusion-dominated model undergoing a small suppression in eta remains in the diffusion-dominated regime. △ Less

Submitted 19 June, 2009; originally announced June 2009.

Comments: 36 pages, 12 figures, accepted for publication in ApJ

Journal ref: Astrophys.J.701:725-736,2009