Showing 1–50 of 229 results for author: Cameron, R

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

The mean solar butterfly diagram and poloidal field generation rate at the surface of the Sun

Authors: Simon Cloutier, Robert H. Cameron, Laurent Gizon

Abstract: The difference between individual solar cycles in the magnetic butterfly diagram can mostly be ascribed to the stochasticity of the emergence process. We aim to obtain the expectation value of the butterfly diagram from observations of four cycles. This allows us to further determine the generation rate of the surface radial magnetic field. We use data from Wilcox Solar Observatory to generate tim… ▽ More The difference between individual solar cycles in the magnetic butterfly diagram can mostly be ascribed to the stochasticity of the emergence process. We aim to obtain the expectation value of the butterfly diagram from observations of four cycles. This allows us to further determine the generation rate of the surface radial magnetic field. We use data from Wilcox Solar Observatory to generate time-latitude diagrams spanning cycles 21 to 24 of the surface radial and toroidal magnetic fields, symmetrize them across the equator and cycle-average them. From the mean butterfly diagram and surface toroidal field we then infer the mean poloidal field generation rate at the surface of the Sun. The averaging procedure removes realization noise from individual cycles. The amount of emerging flux required to account for the evolution of the surface radial field is found to match that provided by the observed surface toroidal field and Joy's law. Cycle-averaging butterfly diagrams removes realization noise and artefacts due to imperfect scale separation, and corresponds to an ensemble average that can be interpreted in the mean-field framework. The result can then be directly compared to $αΩ$-type dynamo models. The Babcock-Leighton $α$-effect is consistent with observations, a result that can be appreciated only if the observational data is averaged in some way. △ Less

Submitted 27 May, 2024; originally announced May 2024.

Comments: 5 pages, 4 figures, submitted to A&A

Non-thermal broadening of coronal lines in a 3D MHD loop model

Authors: C. A. Breu, H. Peter, S. K. Solanki, R. Cameron, I. De Moortel

Abstract: Observed spectral profiles of emission lines from the corona are found to have widths exceeding the thermal line width. To investigate the physical mechanism, we run a 3D MHD model of a single, straightened loop in which we partially resolve turbulent motions that form in response to the driving by self-consistently evolving magneto-convection in the photosphere. The convective motions shear and t… ▽ More Observed spectral profiles of emission lines from the corona are found to have widths exceeding the thermal line width. To investigate the physical mechanism, we run a 3D MHD model of a single, straightened loop in which we partially resolve turbulent motions that form in response to the driving by self-consistently evolving magneto-convection in the photosphere. The convective motions shear and twist the magnetic field lines, leading to heating. From the model we synthesize spectral profiles of emission lines forming at temperatures around and above 1 MK. The coronal heating process generates a range of velocity amplitudes and directions structured on a scale much smaller than the resolving power of current instruments, leading to a broadening of the spectral lines. Our model includes the mass exchange between corona and chromosphere, thus we also capture flows parallel to the loop axis. We find that the spectral lines show a non-thermal line broadening roughly consistent with observations for a viewing angle perpendicular to the axis. The broadening through field-parallel flows is comparable, although slightly smaller. The line broadening is independent of the instrument resolution for a perpendicular line-of-sight. We can connect the non-thermal line broadening to heating events and flows. While small-scale velocities along the line-of-sight are mainly responsible for the broadening observed perpendicular to the loop, chromospheric evaporation is important for the line broadening observed along the loop. The model reproduces observed values for non-thermal line widths. In the model these result from continuous driving by magnetoconvection, without imposing driving motions or starting from an already braided field. △ Less

Submitted 29 March, 2024; originally announced April 2024.

Comments: 18 pages, 18 figures

The Sun's differential rotation is controlled by high-latitude baroclinically unstable inertial modes

Authors: Yuto Bekki, Robert H. Cameron, Laurent Gizon

Abstract: Rapidly rotating fluids have a rotation profile which depends only on the distance from the rotation axis, in accordance with the Taylor-Proudman theorem. Although the Sun was expected to be such a body, helioseismology showed that the rotation rate in the convection zone is closer to constant on radii. It has been postulated that this deviation is due to the poles being warmer than the equator by… ▽ More Rapidly rotating fluids have a rotation profile which depends only on the distance from the rotation axis, in accordance with the Taylor-Proudman theorem. Although the Sun was expected to be such a body, helioseismology showed that the rotation rate in the convection zone is closer to constant on radii. It has been postulated that this deviation is due to the poles being warmer than the equator by a few degrees. Using numerical simulations, we show that the pole-to-equator temperature difference cannot exceed 7 Kelvin as a result of the back-reaction of the high-latitude baroclinically unstable inertial modes. The observed amplitudes of the modes further indicate that this maximum temperature difference is reached in the Sun. We conclude that the Sun's latitudinal differential rotation reaches its maximum allowed value. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Comments: 52 pages, 21 figures, 2 tables. Published in Science Advances on 27 Mar 2024. The authors' version

Journal ref: Science Advances, 10 (13), eadk5643 (2024)

Stellar Surface Magnetic Fields Impact Limb Darkening

Authors: Nadiia M. Kostogryz, Alexander I. Shapiro, Veronika Witzke, Robert H. Cameron, Laurent Gizon, Natalie A. Krivova, Hans-G. Ludwig, Pierre F. L. Maxted, Sara Seager, Sami K. Solanki, Jeff Valenti

Abstract: Stars appear darker at their limbs than at their disk centers because at the limb we are viewing the higher and cooler layers of stellar photospheres. Limb darkening derived from state-of-the-art stellar atmosphere models systematically fails to reproduce recent transiting exoplanet light curves from the Kepler, TESS, and JWST telescopes -- stellar brightness obtained from measurements drops less… ▽ More Stars appear darker at their limbs than at their disk centers because at the limb we are viewing the higher and cooler layers of stellar photospheres. Limb darkening derived from state-of-the-art stellar atmosphere models systematically fails to reproduce recent transiting exoplanet light curves from the Kepler, TESS, and JWST telescopes -- stellar brightness obtained from measurements drops less steeply towards the limb than predicted by models. All previous models assumed atmosphere devoid of magnetic fields. Here we use our new stellar atmosphere models computed with the 3D radiative magneto-hydrodynamic code MURaM to show that small-scale concentration of magnetic fields on the stellar surface affect limb darkening at a level that allows us to explain the observations. Our findings provide a way forward to improve the determination of exoplanet radii and especially the transmission spectroscopy analysis for transiting planets, which relies on a very accurate description of stellar limb darkening from the visible through the infrared. Furthermore, our findings imply that limb darkening allows measuring the small-scale magnetic field on stars with transiting planets. △ Less

Submitted 29 February, 2024; originally announced March 2024.

Comments: 19 pages, 6 figures This is the submission version of manuscript to Nature Astronomy

How well does surface magnetism represent deep Sun-like star dynamo action?

Authors: Adam J. Finley, Sacha A. Brun, Antoine Strugarek, Robert Cameron

Abstract: For Sun-like stars, the generation of toroidal magnetic field from poloidal magnetic field is an essential piece of the dynamo mechanism powering their magnetism. Previous authors have estimated the net toroidal flux generated in each hemisphere of the Sun by exploiting its conservative nature. This only requires observations of the surface magnetic field and differential rotation. We explore this… ▽ More For Sun-like stars, the generation of toroidal magnetic field from poloidal magnetic field is an essential piece of the dynamo mechanism powering their magnetism. Previous authors have estimated the net toroidal flux generated in each hemisphere of the Sun by exploiting its conservative nature. This only requires observations of the surface magnetic field and differential rotation. We explore this approach using a 3D magnetohydrodynamic dynamo simulation of a cool star, for which the magnetic field generation is known throughout the entire star. Changes to the net toroidal flux in each hemisphere were evaluated using a closed line integral bounding the cross-sectional area of each hemisphere, following the application of Stokes-theorem to the induction equation; the individual line segments corresponded to the stellar surface, base, equator, and rotation axis. The influence of the large-scale flows, the fluctuating flows, and magnetic diffusion to each of the line segments was evaluated, along with their depth-dependence. In the simulation, changes to the net toroidal flux via the surface line segment typically dominate the total line integral surrounding each hemisphere, with smaller contributions from the equator and rotation axis. The bulk of the toroidal flux is generated deep inside the convection zone, with the surface observables capturing this due to the conservative nature of the net flux. Surface magnetism and rotation can therefore be used to estimate the net toroidal flux generated in each hemisphere, allowing us to constrain the reservoir of magnetic flux for the next magnetic cycle. However, this methodology cannot identify the physical origin, nor the location, of the toroidal flux generation. In addition, not all dynamo mechanisms depend on the net toroidal field produced in each hemisphere, meaning this method may not be able to characterise every magnetic cycle. △ Less

Submitted 19 January, 2024; originally announced January 2024.

Comments: Accepted to A&A. 8 Pages + Appendix. 6 Figures + 2 Appendix Figures

Small-scale Dynamo in Cool Stars III. Changes in the photospheres of F3V to M0V stars

Authors: Tanayveer S. Bhatia, Robert H. Cameron, Hardi Peter, Sami K. Solanki

Abstract: Some of the quiet solar magnetic flux could be attributed to a small-scale dynamo (SSD) operating in the convection zone. An SSD operating in cool main-sequence stars is expected to affect the atmospheric structure, in particular the convection, and should have observational signatures. We aim to investigate the distribution of these fields as well as their effect on intensity characteristics, vel… ▽ More Some of the quiet solar magnetic flux could be attributed to a small-scale dynamo (SSD) operating in the convection zone. An SSD operating in cool main-sequence stars is expected to affect the atmospheric structure, in particular the convection, and should have observational signatures. We aim to investigate the distribution of these fields as well as their effect on intensity characteristics, velocities and spatial distribution of kinetic (KE) and magnetic energy (ME) in the lower photosphere of spectral types F3V, G2V, K0V and M0V using 3D radiative-MHD simulations. PDFs of field strength at the $τ=1$ surface are quite similar for all cases. The M0V star displays the strongest fields, but relative to the gas pressure, the fields on the F3V star reach the largest values. All stars display an excess of horizontal field relative to vertical field in the middle photosphere, with this excess becoming increasingly prominent towards later spectral types. These fields result in a decrease in upflow velocities, slightly smaller granules as well as the formation of bright points in intergranular lanes. The spatial distribution of KE and ME is also similar for all cases, implying a simple pressure scale height proportionality of important scales. SSD fields have rather similar effects on the photospheres of cool main-sequence stars, namely, a significant reduction in convective velocities as well as a slight reduction in granule size, and concentration of field to kG levels in intergranular lanes associated with the formation of bright points. The distribution of field strengths and energies is also rather similar. △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: 11 pages, 14 figures, accepted for publication in A&A

Testing MURaM and MPS-ATLAS against the quiet solar spectrum

Authors: Veronika Witzke, Alexander I. Shapiro, Nadiia M. Kostogryz, Lucien Mauviard, Tanayveer S. Bhatia, Robert Cameron, Laurent Gizon, Damien Przybylski, Sami K. Solanki, Yvonne C. Unruh, Li Yue

Abstract: Three-dimensional (3D) radiative magnetohydrodynamics (MHD) simulations are the only way to model stellar atmospheres without any ad hoc parameterisations. Several 3D radiative MHD codes have achieved good quantitative agreement with observables for our Sun. We aim to validate the most up-to-date version of the MURaM code against well established quiet Sun measurements, in particular spatially ave… ▽ More Three-dimensional (3D) radiative magnetohydrodynamics (MHD) simulations are the only way to model stellar atmospheres without any ad hoc parameterisations. Several 3D radiative MHD codes have achieved good quantitative agreement with observables for our Sun. We aim to validate the most up-to-date version of the MURaM code against well established quiet Sun measurements, in particular spatially averaged measurements that are relevant for stellar studies. This validation extends the number of solar observables that MURaM can reproduce with high precision. It is also an essential condition for using MURaM to accurately calculate spectra of other cool stars. We simulate the solar upper convection zone and photosphere harbouring a small-scale-dynamo. Using time series of 3D snapshots we calculate the spectral irradiance, limb darkening and selected spectral lines, which we compare to observations. The computed observables agree well with the observations, in particular the limb darkening of the quiet Sun is reproduced remarkably well. △ Less

Submitted 9 October, 2023; originally announced October 2023.

Comments: 9 Page, 7 figures, accepted in A&A

Journal ref: A&A 681, A81 (2024)

A Babcock-Leighton dynamo model of the Sun incorporating toroidal flux loss and the helioseismically-inferred meridional flow

Authors: S. Cloutier, R. H. Cameron, L. Gizon

Abstract: We investigate whether the Babcock-Leighton flux-transport dynamo model remains in agreement with observations if the meridional flow profile is taken from helioseismic inversions. Additionally, we investigate the effect of the loss of toroidal flux through the solar surface. We employ the 2D flux-transport BL dynamo framework. We use the helioseismically-inferred meridional flow profile, and incl… ▽ More We investigate whether the Babcock-Leighton flux-transport dynamo model remains in agreement with observations if the meridional flow profile is taken from helioseismic inversions. Additionally, we investigate the effect of the loss of toroidal flux through the solar surface. We employ the 2D flux-transport BL dynamo framework. We use the helioseismically-inferred meridional flow profile, and include toroidal flux loss in a way that is consistent with the amount of poloidal flux generated by Joy's law. Our model does not impose a preference for emergences at low latitudes, but we do require that the model produces such a preference. We can find solutions in general agreement with observations, including the equatorward drift of the butterfly wings and the cycle's 11 year period. The most important free parameters in the model are the depth to which the radial turbulent pumping extends and the turbulent diffusivity in the lower half of the convection zone. We find that the pumping needs to extend to depths of about $0.80R_{\odot}$ and the bulk turbulent diffusivity needs to be around 10 km$^2$/s or less. We find that the emergences are restricted to low latitudes without the need to impose such a preference. The flux-transport BL model, incorporating the helioseismically inferred meridional flow and toroidal field loss term, is compatible with the properties of the observed butterfly diagram and with the observed toroidal loss rate. Reasonably tight constraints are placed on the remaining free parameters. The pumping needs to be to just below the depth corresponding to the location where the meridional flow changes direction. Our linear model does not however reproduce the observed "rush to the poles" of the diffuse surface radial field resulting from the decay of sunspots -- reproducing this might require the imposition of a preference for flux to emerge near the equator. △ Less

Submitted 29 September, 2023; originally announced September 2023.

Comments: 10 pages, 8 figures, accepted by A&A

Journal ref: A&A 680, A42 (2023)

The Third Fermi Large Area Telescope Catalog of Gamma-ray Pulsars

Authors: David A. Smith, Philippe Bruel, Colin J. Clark, Lucas Guillemot, Matthew T. Kerr, Paul Ray, Soheila Abdollahi, Marco Ajello, Luca Baldini, Jean Ballet, Matthew Baring, Cees Bassa, Josefa Becerra Gonzalez, Ronaldo Bellazzini, Alessandra Berretta, Bhaswati Bhattacharyya, Elisabetta Bissaldi, Raffaella Bonino, Eugenio Bottacini, Johan Bregeon, Marta Burgay, Toby Burnett, Rob Cameron, Fernando Camilo, Regina Caputo , et al. (134 additional authors not shown)

Abstract: We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems co-located with LAT sources also likely harbor gamma-ray M… ▽ More We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems co-located with LAT sources also likely harbor gamma-ray MSPs. This catalog thus reports roughly 340 gamma-ray pulsars and candidates, 10% of all known pulsars, compared to $\leq 11$ known before Fermi. Half of the gamma-ray pulsars are young. Of these, the half that are undetected in radio have a broader Galactic latitude distribution than the young radio-loud pulsars. The others are MSPs, with 6 undetected in radio. Overall, >235 are bright enough above 50 MeV to fit the pulse profile, the energy spectrum, or both. For the common two-peaked profiles, the gamma-ray peak closest to the magnetic pole crossing generally has a softer spectrum. The spectral energy distributions tend to narrow as the spindown power $\dot E$ decreases to its observed minimum near $10^{33}$ erg s$^{-1}$, approaching the shape for synchrotron radiation from monoenergetic electrons. We calculate gamma-ray luminosities when distances are available. Our all-sky gamma-ray sensitivity map is useful for population syntheses. The electronic catalog version provides gamma-ray pulsar ephemerides, properties and fit results to guide and be compared with modeling results. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: 142 pages. Accepted by the Astrophysical Journal Supplement

Spectral variability of photospheric radiation due to faculae II: Facular contrasts for cool main-sequence stars

Authors: Charlotte M. Norris, Yvonne C. Unruh, Veronika Witzke, Sami K. Solanki, Natalie A. Krivova, Alexander I. Shapiro, Robert Cameron, Benjamin Beeck

Abstract: Magnetic features on the surface of stars, such as spots and faculae, cause stellar spectral variability on time-scales of days and longer. For stars other than the Sun, the spectral signatures of faculae are poorly understood, limiting our ability to account for stellar pollution in exoplanet transit observations. Here we present the first facular contrasts derived from magnetoconvection simulati… ▽ More Magnetic features on the surface of stars, such as spots and faculae, cause stellar spectral variability on time-scales of days and longer. For stars other than the Sun, the spectral signatures of faculae are poorly understood, limiting our ability to account for stellar pollution in exoplanet transit observations. Here we present the first facular contrasts derived from magnetoconvection simulations for K0, M0 and M2 main-sequence stars and compare them to previous calculations for G2 main-sequence stars. We simulate photospheres and immediate subsurface layers of main-sequence spectral types between K0 and M2, with different injected vertical magnetic fields (0 G, 100 G, 300 G and 500 G) using MURaM, a 3D radiation-magnetohydrodynamics code. We show synthetic spectra and contrasts from the UV (300 nm) to the IR (10000 nm) calculated using the ATLAS9 radiative transfer code. The calculations are performed for nine viewing angles to characterise the facular radiation across the disc. The brightness contrasts of magnetic regions are found to change significantly across spectral type, wavelength and magnetic field strength, leading to the conclusion that accurate contrasts cannot be found by scaling solar values. This is due to features of different size, apparent structure and spectral brightness emerging in the presence of a given magnetic field for different spectral types. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Comments: 17 pages, 12 figures, to be published in MNRAS

Swirls in the Solar Corona

Authors: C. Breu, H. Peter, R. Cameron, S. K. Solanki

Abstract: Vortex flows have been found in the photosphere, chromosphere and low corona in observations and simulations. It has been suggested that vortices play an important role for channeling energy and plasma into the corona, but the impact of vortex flows on the corona has not directly been studied in a realistic setup. We investigate the role vortices play for coronal heating using high resolution simu… ▽ More Vortex flows have been found in the photosphere, chromosphere and low corona in observations and simulations. It has been suggested that vortices play an important role for channeling energy and plasma into the corona, but the impact of vortex flows on the corona has not directly been studied in a realistic setup. We investigate the role vortices play for coronal heating using high resolution simulations of coronal loops. The vortices are not artificially driven, but arise self-consistently from magnetoconvection. We perform 3D resistive MHD simulations with the MURaM code. Studying an isolated coronal loop in a Cartesian geometry allows us to resolve the structure of the loop interior. We conduct a statistical analysis to determine vortex properties as a function of height from the chromosphere into the corona. We find that the energy injected into the loop is generated by internal coherent motions within strong magnetic elements. A significant part of the resulting Poynting flux is channeled through the chromosphere in vortex tubes forming a magnetic connection between the photosphere and corona. Vortices can form contiguous structures that reach up to coronal heights, but in the corona itself the vortex tubes get deformed and eventually lose their identity with increasing height. Vortices show increased upward directed Poynting flux and heating rate both in the chromosphere and corona, but their effect becomes less pronounced with increasing height. While vortices play an important role for the energy transport and structuring in the chromosphere and low corona, their importance higher up in the atmosphere is less clear since the swirls are less distinguishable from their environment. Vortex tubes reaching the corona show a complex relationship with the coronal emission. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Comments: 18 pages, 12 figures

Journal ref: A&A 675, A94 (2023)

Observationally guided models for the solar dynamo and the role of the surface field

Authors: Robert Cameron, Manfred Schüssler

Abstract: Theoretical models for the solar dynamo range from simple low-dimensional ``toy models'' to complex 3D-MHD simulations. Here we mainly discuss appproaches that are motivated and guided by solar (and stellar) observations. We give a brief overview of the evolution of solar dynamo models since 1950s, focussing upon the development of the Babcock-Leighton approach between its introduction in the 1960… ▽ More Theoretical models for the solar dynamo range from simple low-dimensional ``toy models'' to complex 3D-MHD simulations. Here we mainly discuss appproaches that are motivated and guided by solar (and stellar) observations. We give a brief overview of the evolution of solar dynamo models since 1950s, focussing upon the development of the Babcock-Leighton approach between its introduction in the 1960s and its revival in the 1990s after being long overshadowed by mean-field turbulent dynamo theory. We summarize observations and simple theoretical deliberations that demonstrate the crucial role of the surface fields in the dynamo process and and give quantitative analyses of the generation and loss of toroidal flux in the convection zone as well as of the production of poloidal field resulting from flux emergence at the surface. Furthermore, we discuss possible nonlinearities in the dynamo process suggested by observational results and present models for the long-term variability of solar activity motivated by observations of magnetically active stars and the inherent randomness of the dynamo process. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Comments: submitted to Space Science Reviews

Fermi-GBM Discovery of GRB 221009A: An Extraordinarily Bright GRB from Onset to Afterglow

Authors: S. Lesage, P. Veres, M. S. Briggs, A. Goldstein, D. Kocevski, E. Burns, C. A. Wilson-Hodge, P. N. Bhat, D. Huppenkothen, C. L. Fryer, R. Hamburg, J. Racusin, E. Bissaldi, W. H. Cleveland, S. Dalessi, C. Fletcher, M. M. Giles, B. A. Hristov, C. M. Hui, B. Mailyan, C. Malacaria, S. Poolakkil, O. J. Roberts, A. von Kienlin, J. Wood , et al. (115 additional authors not shown)

Abstract: We report the discovery of GRB 221009A, the highest flux gamma-ray burst ever observed by the Fermi Gamma-ray Burst Monitor (GBM). This GRB has continuous prompt emission lasting more than 600 seconds which smoothly transitions to afterglow visible in the GBM energy range (8 keV--40 MeV), and total energetics higher than any other burst in the GBM sample. By using a variety of new and existing ana… ▽ More We report the discovery of GRB 221009A, the highest flux gamma-ray burst ever observed by the Fermi Gamma-ray Burst Monitor (GBM). This GRB has continuous prompt emission lasting more than 600 seconds which smoothly transitions to afterglow visible in the GBM energy range (8 keV--40 MeV), and total energetics higher than any other burst in the GBM sample. By using a variety of new and existing analysis techniques we probe the spectral and temporal evolution of GRB 221009A. We find no emission prior to the GBM trigger time (t0; 2022 October 9 at 13:16:59.99 UTC), indicating that this is the time of prompt emission onset. The triggering pulse exhibits distinct spectral and temporal properties suggestive of the thermal, photospheric emission of shock-breakout, with significant emission up to $\sim$15 MeV. We characterize the onset of external shock at t0+600 s and find evidence of a plateau region in the early-afterglow phase which transitions to a slope consistent with Swift-XRT afterglow measurements. We place the total energetics of GRB 221009A in context with the rest of the GBM sample and find that this GRB has the highest total isotropic-equivalent energy ($\textrm{E}_{γ,\textrm{iso}}=1.0\times10^{55}$ erg) and second highest isotropic-equivalent luminosity ($\textrm{L}_{γ,\textrm{iso}}=9.9\times10^{53}$ erg/s) based on redshift of z = 0.151. These extreme energetics are what allowed us to observe the continuously emitting central engine of GBM from the beginning of the prompt emission phase through the onset of early afterglow. △ Less

Submitted 12 July, 2023; v1 submitted 24 March, 2023; originally announced March 2023.

Comments: 26 pages 7 figures - accepted for publication in ApJL

First measurement of the nuclear-recoil ionization yield in silicon at 100 eV

Authors: M. F. Albakry, I. Alkhatib, D. Alonso, D. W. P. Amaral, P. An, T. Aralis, T. Aramaki, I. J. Arnquist, I. Ataee Langroudy, E. Azadbakht, S. Banik, P. S. Barbeau, C. Bathurst, R. Bhattacharyya, P. L. Brink, R. Bunker, B. Cabrera, R. Calkins, R. A. Cameron, C. Cartaro, D. G. Cerdeño, Y. -Y. Chang, M. Chaudhuri, R. Chen, N. Chott , et al. (115 additional authors not shown)

Abstract: We measured the nuclear--recoil ionization yield in silicon with a cryogenic phonon-sensitive gram-scale detector. Neutrons from a mono-energetic beam scatter off of the silicon nuclei at angles corresponding to energy depositions from 4\,keV down to 100\,eV, the lowest energy probed so far. The results show no sign of an ionization production threshold above 100\,eV. These results call for furthe… ▽ More We measured the nuclear--recoil ionization yield in silicon with a cryogenic phonon-sensitive gram-scale detector. Neutrons from a mono-energetic beam scatter off of the silicon nuclei at angles corresponding to energy depositions from 4\,keV down to 100\,eV, the lowest energy probed so far. The results show no sign of an ionization production threshold above 100\,eV. These results call for further investigation of the ionization yield theory and a comprehensive determination of the detector response function at energies below the keV scale. △ Less

Submitted 3 March, 2023; originally announced March 2023.

Journal ref: Physical Review Letters 131.9 (2023): 091801

No evidence for synchronization of the solar cycle by a "clock"

Authors: E. Weisshaar, R. H. Cameron, M. Schüssler

Abstract: The length of the solar activity cycle fluctuates considerably. The temporal evolution of the corresponding cycle phase, that is, the deviation of the epochs of activity minima or maxima from strict periodicity, provides relevant information concerning the physical mechanism underlying the cyclic magnetic activity. An underlying strictly periodic process (akin to a perfect "clock"), with the obser… ▽ More The length of the solar activity cycle fluctuates considerably. The temporal evolution of the corresponding cycle phase, that is, the deviation of the epochs of activity minima or maxima from strict periodicity, provides relevant information concerning the physical mechanism underlying the cyclic magnetic activity. An underlying strictly periodic process (akin to a perfect "clock"), with the observer seeing a superposition of the perfect clock and a small random phase perturbation, leads to long-term phase stability in the observations. Such behavior would be expected if cycles were synchronized by tides caused by orbiting planets or by a hypothetical torsional oscillation in the solar radiative interior. Alternatively, in the absence of such synchronization, phase fluctuations accumulate and a random walk of the phase ensues, which is a typical property of randomly perturbed dynamo models. Based on the sunspot record and the reconstruction of solar cycles from cosmogenic C14, we carried out rigorous statistical tests in order to decipher whether there exists phase synchronization or random walk. Synchronization is rejected at significance levels of between 95% (28 cycles from sunspot data) and beyond 99% (84 cycles reconstructed from C14, while the existence of random walk in the phases is consistent with all data sets. This result strongly supports randomly perturbed dynamo models with little inter-cycle memory. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: Astronomy & Astrophysics, in press

Journal ref: A&A 671, A87 (2023)

Exploring the Solar Poles: The Last Great Frontier of the Sun

Authors: Dibyendu Nandy, Dipankar Banerjee, Prantika Bhowmik, Allan Sacha Brun, Robert H. Cameron, S. E. Gibson, Shravan Hanasoge, Louise Harra, Donald M. Hassler, Rekha Jain, Jie Jiang, Laurène Jouve, Duncan H. Mackay, Sushant S. Mahajan, Cristina H. Mandrini, Mathew Owens, Shaonwita Pal, Rui F. Pinto, Chitradeep Saha, Xudong Sun, Durgesh Tripathi, Ilya G. Usoskin

Abstract: Despite investments in multiple space and ground-based solar observatories by the global community, the Sun's polar regions remain unchartered territory - the last great frontier for solar observations. Breaching this frontier is fundamental to understanding the solar cycle - the ultimate driver of short-to-long term solar activity that encompasses space weather and space climate. Magnetohydrodyna… ▽ More Despite investments in multiple space and ground-based solar observatories by the global community, the Sun's polar regions remain unchartered territory - the last great frontier for solar observations. Breaching this frontier is fundamental to understanding the solar cycle - the ultimate driver of short-to-long term solar activity that encompasses space weather and space climate. Magnetohydrodynamic dynamo models and empirically observed relationships have established that the polar field is the primary determinant of the future solar cycle amplitude. Models of solar surface evolution of tilted active regions indicate that the mid to high latitude surges of magnetic flux govern dynamics leading to the reversal and build-up of polar fields. Our theoretical understanding and numerical models of this high latitude magnetic field dynamics and plasma flows - that are a critical component of the sunspot cycle - lack precise observational constraints. This limitation compromises our ability to observe the enigmatic kilo Gauss polar flux patches and constrain the polar field distribution at high latitudes. The lack of these observations handicap our understanding of how high latitude magnetic fields power polar jets, plumes, and the fast solar wind that extend to the boundaries of the heliosphere and modulate solar open flux and cosmic ray flux within the solar system. Accurate observation of the Sun's polar regions, therefore, is the single most outstanding challenge that confronts Heliophysics. This paper argues the scientific case for novel out of ecliptic observations of the Sun's polar regions, in conjunction with existing, or future multi-vantage point heliospheric observatories. Such a mission concept can revolutionize the field of Heliophysics like no other mission concept has - with relevance that transcends spatial regimes from the solar interior to the heliosphere. △ Less

Submitted 30 December, 2022; originally announced January 2023.

Comments: This White Paper was submitted in 2022 to the United States National Academies Solar and Space Physics (Heliophysics) Decadal Survey

Can 1D radiative equilibrium models of faculae be used for calculating contamination of transmission spectra?

Authors: Veronika Witzke, Alexander I. Shapiro, Nadiia M. Kostogryz, Robert Cameron, Benjamin V. Rackham, Sara Seager, Sami K. Solanki, Yvonne C. Unruh

Abstract: The reliable characterization of planetary atmospheres with transmission spectroscopy requires realistic modeling of stellar magnetic features, since features that are attributable to an exoplanet atmosphere could instead stem from the host star's magnetic activity. Current retrieval algorithms for analysing transmission spectra rely on intensity contrasts of magnetic features from 1D radiative-co… ▽ More The reliable characterization of planetary atmospheres with transmission spectroscopy requires realistic modeling of stellar magnetic features, since features that are attributable to an exoplanet atmosphere could instead stem from the host star's magnetic activity. Current retrieval algorithms for analysing transmission spectra rely on intensity contrasts of magnetic features from 1D radiative-convective models. However, magnetic features, especially faculae, are not fully captured by such simplified models. Here we investigate how well such 1D models can reproduce 3D facular contrasts, taking a G2V star as an example. We employ the well established radiative magnetohydrodynamic code MURaM to obtain three-dimensional simulations of the magneto-convection and photosphere harboring a local small-scale-dynamo. Simulations without additional vertical magnetic fields are taken to describe the quiet solar regions, while simulations with initially 100 G, 200 G and 300 G vertical magnetic fields are used to represent different magnetic activity levels. Subsequently, the spectra emergent from the MURaM cubes are calculated with the MPS-ATLAS radiative transfer code. We find that the wavelength dependence of facular contrast from 1D radiative-convective models cannot reproduce facular contrasts obtained from 3D modeling. This has far reaching consequences for exoplanet characterization using transmission spectroscopy, where accurate knowledge of the host star is essential for unbiased inferences of the planetary atmospheric properties. △ Less

Submitted 5 November, 2022; originally announced November 2022.

Comments: 7 pages, 2 figures, submitted to APJL

Small-scale dynamo in cool main sequence stars. II. The effect of metallicity

Authors: V. Witzke, H. B. Duehnen, A. I. Shapiro, D. Przybylski, T. S. Bhatia, R. Cameron, S. K. Solanki

Abstract: All cool main sequence stars including our Sun are thought to have magnetic fields. Observations of the Sun revealed that even in quiet regions small-scale turbulent magnetic fields are present. Simulations further showed that such magnetic fields affect the subsurface and photospheric structure, and thus the radiative transfer and emergent flux. Since small-scale turbulent magnetic fields on othe… ▽ More All cool main sequence stars including our Sun are thought to have magnetic fields. Observations of the Sun revealed that even in quiet regions small-scale turbulent magnetic fields are present. Simulations further showed that such magnetic fields affect the subsurface and photospheric structure, and thus the radiative transfer and emergent flux. Since small-scale turbulent magnetic fields on other stars cannot be directly observed, it is imperative to study their effects on the near surface layers numerically. Until recently comprehensive three-dimensional simulations capturing the effect of small-scale turbulent magnetic fields only exists for the solar case. A series of investigations extending SSD simulations for other stars has been started. Here we aim to examine small-scale turbulent magnetic fields in stars of solar effective temperature but different metallicity. We investigate the properties of three-dimensional simulations of the magneto-convection in boxes covering the upper convection zone and photosphere carried out with the MURaM code for metallicity values of $ \rm M/H = \{-1.0, 0.0, 0.5\}$ with and without a small-scale-dynamo. We find that small-scale turbulent magnetic fields enhanced by a small-scale turbulent dynamo noticeably affect the subsurface dynamics and significantly change the flow velocities in the photosphere. Moreover, significantly stronger magnetic field strengths are present in the convection zone for low metallicity. Whereas, at the optical surface the averaged vertical magnetic field ranges from 64G for M/H = 0.5 to 85G for M/H = -1.0. △ Less

Submitted 4 November, 2022; originally announced November 2022.

Comments: 13 pages, 18 figures, submitted to A&A

Journal ref: A&A 669, A157 (2023)

Toroidal flux loss due to flux emergence explains why solar cycles rise differently but decay in a similar way

Authors: Akash Biswas, Bidya Binay Karak, Robert Cameron

Abstract: A striking feature of the solar cycle is that at the beginning, sunspots appear around mid-latitudes, and over time the latitudes of emergences migrate towards the equator.The maximum level of activity (e.g., sunspot number) varies from cycle to cycle.For strong cycles, the activity begins early and at higher latitudes with wider sunspot distributions than for weak cycles. The activity and the wid… ▽ More A striking feature of the solar cycle is that at the beginning, sunspots appear around mid-latitudes, and over time the latitudes of emergences migrate towards the equator.The maximum level of activity (e.g., sunspot number) varies from cycle to cycle.For strong cycles, the activity begins early and at higher latitudes with wider sunspot distributions than for weak cycles. The activity and the width of sunspot belts increase rapidly and begin to decline when the belts are still at high latitudes. Surprisingly, it has been reported that in the late stages of the cycle the level of activity (sunspot number) as well as the widths and centers of the butterfly wings all have the same statistical properties independent of how strong the cycle was during its rise and maximum phases.We have modeled these features using a Babcock--Leighton type dynamo model and show that the flux loss through magnetic buoyancy is an essential nonlinearity in the solar dynamo.Our study shows that the nonlinearity is effective if the flux emergence becomes efficient at the mean-field strength of the order of $10^4$~G in the lower part of the convection zone. △ Less

Submitted 13 October, 2022; originally announced October 2022.

Comments: Accepted in Phys.Rev.Lett. (PRL)

Three-dimensional non-kinematic simulation of post-emergence evolution of bipolar magnetic regions and Babcock-Leighton dynamo of the Sun

Authors: Yuto Bekki, Robert H. Cameron

Abstract: The Babcock-Leighton (BL) flux-transport model is a widely-accepted dynamo model of the Sun. This dynamo model has been extensively studied in a two-dimensional (2D) mean-field framework in both kinematic and non-kinematic regimes. Recent three-dimensional (3D) models have been restricted to the kinematic regime. In these models, the surface poloidal flux is produced by the emergence of bipolar ma… ▽ More The Babcock-Leighton (BL) flux-transport model is a widely-accepted dynamo model of the Sun. This dynamo model has been extensively studied in a two-dimensional (2D) mean-field framework in both kinematic and non-kinematic regimes. Recent three-dimensional (3D) models have been restricted to the kinematic regime. In these models, the surface poloidal flux is produced by the emergence of bipolar magnetic regions (BMRs) that are tilted according to Joy's law. We investigate the prescription for emergence of a BMR in 3D non-kinematic simulations. We also report initial results of cyclic BL dynamo simulation. We extend a conventional 2D mean-field model of the BL flux-transport dynamo into 3D non-kinematic regime. The large-scale mean flows are driven by the parameterized $Λ$-effect in this model. For the induction equation, we use a BL source term by which the surface BMRs are produced in response to the dynamo-generated toroidal field inside the convection zone. We find that, in the 3D non-kinematic regime, the tilt angle of a newly-emerged BMR is very sensitive to the prescription for the subsurface structure of the BMR. Anti-Joy tilt angles are found unless the BMR is deeply embedded in the convection zone. We also find that the leading spot tends to become stronger than the following spot. The anti-Joy's law trend and the morphological asymmetry of the BMRs can be explained by the Coriolis force acting on the Lorentz-force-driven flows. Furthermore, we demonstrate that the solar-like magnetic cycles can be successfully obtained if the Joy's law is explicitly given in the BL $α$-effect. In these cyclic dynamo simulation, a strong Lorentz force feedback leads to cycle modulations in the differential rotation and meridional circulation. The non-axisymmetric components of the flows are found to exist as inertial modes such as the equatorial Rossby modes. △ Less

Submitted 7 December, 2022; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: 18 pages, 19 figures, accepted for publication in A&A

Journal ref: A&A 670, A101 (2023)

Theory of solar oscillations in the inertial frequency range: Amplitudes of equatorial modes from a nonlinear rotating convection simulation

Authors: Yuto Bekki, Robert H. Cameron, Laurent Gizon

Abstract: Several types of inertial modes have been detected on the Sun. Properties of these inertial modes have been studied in the linear regime but have not been studied in nonlinear simulations of solar rotating convection. Comparing the nonlinear simulations, the linear theory, and the solar observations is important to better understand the differences between the models and the real Sun. We wish to d… ▽ More Several types of inertial modes have been detected on the Sun. Properties of these inertial modes have been studied in the linear regime but have not been studied in nonlinear simulations of solar rotating convection. Comparing the nonlinear simulations, the linear theory, and the solar observations is important to better understand the differences between the models and the real Sun. We wish to detect and characterize the modes present in a nonlinear numerical simulation of solar convection, in particular to understand the amplitudes and lifetimes of the modes. We developed a code with a Yin-Yang grid to carry out fully-nonlinear numerical simulations of rotating convection in a spherical shell. The stratification is solar-like up to 0.96R. The simulations cover a duration of about 15 solar years. Various large-scale modes at low frequencies are extracted from the simulation. Their characteristics are compared to those from the linear model and to the observations. Among other modes, both the equatorial Rossby modes and the columnar convective modes are seen in the simulation. The columnar convective modes contain most of the large-scale velocity power outside the tangential cylinder and substantially contribute to the heat and angular momentum transport. Equatorial Rossby modes with no radial node (n=0) are also found: They have the same spatial structures as the linear eigenfunctions. They are stochastically excited by convection and have the amplitudes of a few m/s and mode linewidths of about 20-30 nHz, which are comparable to those observed on the Sun. We also confirm the existence of the mixed modes between the equatorial Rossby modes and the columnar convective modes in our nonlinear simulation, as predicted by the linear eigenmode analysis. We also see the high-latitude mode with m=1 in our nonlinear simulation but its amplitude is much weaker than that observed on the Sun. △ Less

Submitted 23 August, 2022; originally announced August 2022.

Comments: 18 pages, 21 figures, 1 table. Astronomy & Astrophysics

Journal ref: A&A 666, A135 (2022)

Faculae cancel out on the surfaces of active Suns

Authors: N. -E. Nèmec, A. I. Shapiro, E. Işık, K. Sowmya, S. K. Solanki, N. A. Krivova, R. H. Cameron, L. Gizon

Abstract: Surfaces of the Sun and other cool stars are filled with magnetic fields, which are either seen as dark compact spots or more diffuse bright structures like faculae. Both hamper detection and characterisation of exoplanets, affecting stellar brightness and spectra, as well as transmission spectra. However, the expected facular and spot signals in stellar data are quite different, for instance they… ▽ More Surfaces of the Sun and other cool stars are filled with magnetic fields, which are either seen as dark compact spots or more diffuse bright structures like faculae. Both hamper detection and characterisation of exoplanets, affecting stellar brightness and spectra, as well as transmission spectra. However, the expected facular and spot signals in stellar data are quite different, for instance they have distinct temporal and spectral profiles. Consequently, corrections of stellar data for magnetic activity can greatly benefit from the insight on whether the stellar signal is dominated by spots or faculae. Here, we utilise a surface flux transport model (SFTM) to show that more effective cancellation of diffuse magnetic flux associated with faculae leads to spot area coverages increasing faster with stellar magnetic activity than that by faculae. Our calculations explain the observed dependence between solar spot and facular area coverages and allow its extension to stars more active than the Sun. This extension enables anticipating the properties of stellar signal and its more reliable mitigation, leading to a more accurate characterisation of exoplanets and their atmospheres. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 10 pages, 6 Figures, accepted for publication in ApJL

Impact of spatially correlated fluctuations in sunspots on metrics related to magnetic twist

Authors: C. Baumgartner, A. C. Birch, H. Schunker, R. H. Cameron, L. Gizon

Abstract: The twist of the magnetic field above a sunspot is an important quantity in solar physics. For example, magnetic twist plays a role in the initiation of flares and coronal mass ejections (CMEs). Various proxies for the twist above the photosphere have been found using models of uniformly twisted flux tubes, and are routinely computed from single photospheric vector magnetograms. One class of proxi… ▽ More The twist of the magnetic field above a sunspot is an important quantity in solar physics. For example, magnetic twist plays a role in the initiation of flares and coronal mass ejections (CMEs). Various proxies for the twist above the photosphere have been found using models of uniformly twisted flux tubes, and are routinely computed from single photospheric vector magnetograms. One class of proxies is based on $α_z$, the ratio of the vertical current to the vertical magnetic field. Another class of proxies is based on the so-called twist density, $q$, which depends on the ratio of the azimuthal field to the vertical field. However, the sensitivity of these proxies to temporal fluctuations of the magnetic field has not yet been well characterized. We aim to determine the sensitivity of twist proxies to temporal fluctuations in the magnetic field as estimated from time-series of SDO/HMI vector magnetic field maps. To this end, we introduce a model of a sunspot with a peak vertical field of 2370 Gauss at the photosphere and a uniform twist density $q= -0.024$ Mm$^{-1}$. We add realizations of the temporal fluctuations of the magnetic field that are consistent with SDO/HMI observations, including the spatial correlations. Using a Monte-Carlo approach, we determine the robustness of the different proxies to the temporal fluctuations. The temporal fluctuations of the three components of the magnetic field are correlated for spatial separations up to 1.4 Mm (more than expected from the point spread function alone). The Monte-Carlo approach enables us to demonstrate that several proxies for the twist of the magnetic field are not biased in each of the individual magnetograms. The associated random errors on the proxies have standard deviations in the range between $0.002$ and $0.006$ Mm$^{-1}$, which is smaller by approximately one order of magnitude than the mean value of $q$. △ Less

Submitted 5 July, 2022; originally announced July 2022.

Comments: accepted in A&A

Journal ref: A&A 664, A183 (2022)

Small-scale Dynamo in Cool Stars: I. Changes in stratification and near-surface convection for main-sequence spectral types

Authors: Tanayveer S. Bhatia, Robert H. Cameron, Sami K. Solanki, Hardi Peter, Damien Przybylski, Veronika Witzke, Alexander Shapiro

Abstract: Some of the small-scale solar magnetic flux can be attributed to a small-scale dynamo (SSD) operating in the near-surface convection. The SSD fields have consequences for solar granular convection, basal flux, as well as chromospheric heating. A similar SSD mechanism is expected to be active in the near-surface convection of other cool main-sequence stars, but this has never been investigated. We… ▽ More Some of the small-scale solar magnetic flux can be attributed to a small-scale dynamo (SSD) operating in the near-surface convection. The SSD fields have consequences for solar granular convection, basal flux, as well as chromospheric heating. A similar SSD mechanism is expected to be active in the near-surface convection of other cool main-sequence stars, but this has never been investigated. We aim to investigate changes in stratification and convection due to inclusion of SSD fields for F3V, G2V, K0V and M0V spectral types in the near-surface convection. 3D magnetohydrodynamic (MHD) models of the four stellar boxes, covering the subsurface convection zone up to the lower photosphere in a small cartesian box, are studied using the \textit{MURaM} radiative-MHD simulation code. The SSD runs are compared against reference hydrodynamic runs. An SSD is found to efficiently produce magnetic field with energies ranging between 5\% to 80\% of the plasma kinetic energy at different depths. This ratio tends to be larger for larger $T_{\mathrm{eff}}$. The relative change in density and gas pressure stratification for the deeper convective layers due to SSD magnetic fields is negligible, except for the F-star. For the F-star, there is a substantial reduction in convective velocities due to Lorentz force feedback from magnetic fields, which, in turn, reduces the turbulent pressure. SSD in near-surface convection for cool main-sequence stars introduces small but significant changes in thermodynamic stratification (especially for the F-star) due to reduction in convective velocities. △ Less

Submitted 20 July, 2023; v1 submitted 31 May, 2022; originally announced June 2022.

Comments: Update (20.07.2023): an errata (appendix D) has been included to correct a small error in the derivation of turbulent pressure in appendix B. In addition, the text has been updated with the published version. 10 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 663, A166 (2022)

Effective Field Theory Analysis of CDMSlite Run 2 Data

Authors: SuperCDMS Collaboration, M. F. Albakry, I. Alkhatib, D. W. P. Amaral, T. Aralis, T. Aramaki, I. J. Arnquist, I. Ataee Langroudy, E. Azadbakht, S. Banik, C. Bathurst, D. A. Bauer, L. V. S. Bezerra, R. Bhattacharyya, P. L. Brink, R. Bunker, B. Cabrera, R. Calkins, R. A. Cameron, C. Cartaro, D. G. Cerdeño, Y. -Y. Chang, M. Chaudhuri, R. Chen, N. Chott , et al. (105 additional authors not shown)

Abstract: CDMSlite Run 2 was a search for weakly interacting massive particles (WIMPs) with a cryogenic 600 g Ge detector operated in a high-voltage mode to optimize sensitivity to WIMPs of relatively low mass from 2 - 20 GeV/$c^2$. In this article, we present an effective field theory (EFT) analysis of the CDMSlite Run 2 data using an extended energy range and a comprehensive treatment of the expected back… ▽ More CDMSlite Run 2 was a search for weakly interacting massive particles (WIMPs) with a cryogenic 600 g Ge detector operated in a high-voltage mode to optimize sensitivity to WIMPs of relatively low mass from 2 - 20 GeV/$c^2$. In this article, we present an effective field theory (EFT) analysis of the CDMSlite Run 2 data using an extended energy range and a comprehensive treatment of the expected background. A binned likelihood Bayesian analysis was performed on the recoil energy data, taking into account the parameters of the EFT interactions and optimizing the data selection with respect to the dominant background components. Energy regions within 5$σ$ of known activation peaks were removed from the analysis. The Bayesian evidences resulting from the different operator hypotheses show that the CDMSlite Run 2 data are consistent with the background-only models and do not allow for a signal interpretation assuming any additional EFT interaction. Consequently, upper limits on the WIMP mass and coupling-coefficient amplitudes and phases are presented for each EFT operator. These limits improve previous CDMSlite Run 2 bounds for WIMP masses above 5 GeV/$c^2$. △ Less

Submitted 23 May, 2022; originally announced May 2022.

Comments: 16 pages, 8 figures

Search for new cosmic-ray acceleration sites within the 4FGL catalog Galactic plane sources

Authors: Fermi-LAT Collaboration, S. Abdollahi, F. Acero, M. Ackermann, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, B. Berenji, A. Berretta, E. Bissaldi, R. D. Blandford, R. Bonino, P. Bruel, S. Buson, R. A. Cameron, R. Caputo, P. A. Caraveo, D. Castro, G. Chiaro, N. Cibrario, S. Ciprini, J. Coronado-Blázquez, M. Crnogorcevic , et al. (95 additional authors not shown)

Abstract: Cosmic rays are mostly composed of protons accelerated to relativistic speeds. When those protons encounter interstellar material, they produce neutral pions which in turn decay into gamma rays. This offers a compelling way to identify the acceleration sites of protons. A characteristic hadronic spectrum, with a low-energy break around 200 MeV, was detected in the gamma-ray spectra of four Superno… ▽ More Cosmic rays are mostly composed of protons accelerated to relativistic speeds. When those protons encounter interstellar material, they produce neutral pions which in turn decay into gamma rays. This offers a compelling way to identify the acceleration sites of protons. A characteristic hadronic spectrum, with a low-energy break around 200 MeV, was detected in the gamma-ray spectra of four Supernova Remnants (SNRs), IC 443, W44, W49B and W51C, with the Fermi Large Area Telescope. This detection provided direct evidence that cosmic-ray protons are (re-)accelerated in SNRs. Here, we present a comprehensive search for low-energy spectral breaks among 311 4FGL catalog sources located within 5 degrees from the Galactic plane. Using 8 years of data from the Fermi Large Area Telescope between 50 MeV and 1 GeV, we find and present the spectral characteristics of 56 sources with a spectral break confirmed by a thorough study of systematic uncertainty. Our population of sources includes 13 SNRs for which the proton-proton interaction is enhanced by the dense target material; the high-mass gamma-ray binary LS~I +61 303; the colliding wind binary eta Carinae; and the Cygnus star-forming region. This analysis better constrains the origin of the gamma-ray emission and enlarges our view to potential new cosmic-ray acceleration sites. △ Less

Submitted 6 May, 2022; originally announced May 2022.

Comments: Accepted for publication in The Astrophysical Journal

Investigating the sources of low-energy events in a SuperCDMS-HVeV detector

Authors: SuperCDMS Collaboration, M. F. Albakry, I. Alkhatib, D. W. P. Amaral, T. Aralis, T. Aramaki, I. J. Arnquist, I. Ataee Langroudy, E. Azadbakht, S. Banik, C. Bathurst, D. A. Bauer, R. Bhattacharyya, P. L. Brink, R. Bunker, B. Cabrera, R. Calkins, R. A. Cameron, C. Cartaro, D. G. Cerdeño, Y. -Y. Chang, M. Chaudhuri, R. Chen, N. Chott, J. Cooley , et al. (104 additional authors not shown)

Abstract: Recent experiments searching for sub-GeV/$c^2$ dark matter have observed event excesses close to their respective energy thresholds. Although specific to the individual technologies, the measured excess event rates have been consistently reported at or below event energies of a few-hundred eV, or with charges of a few electron-hole pairs. In the present work, we operated a 1-gram silicon SuperCDMS… ▽ More Recent experiments searching for sub-GeV/$c^2$ dark matter have observed event excesses close to their respective energy thresholds. Although specific to the individual technologies, the measured excess event rates have been consistently reported at or below event energies of a few-hundred eV, or with charges of a few electron-hole pairs. In the present work, we operated a 1-gram silicon SuperCDMS-HVeV detector at three voltages across the crystal (0 V, 60 V and 100 V). The 0 V data show an excess of events in the tens of eV region. Despite this event excess, we demonstrate the ability to set a competitive exclusion limit on the spin-independent dark matter--nucleon elastic scattering cross section for dark matter masses of $\mathcal{O}(100)$ MeV/$c^2$, enabled by operation of the detector at 0 V potential and achievement of a very low $\mathcal{O}(10)$ eV threshold for nuclear recoils. Comparing the data acquired at 0 V, 60 V and 100 V potentials across the crystal, we investigated possible sources of the unexpected events observed at low energy. The data indicate that the dominant contribution to the excess is consistent with a hypothesized luminescence from the printed circuit boards used in the detector holder. △ Less

Submitted 11 October, 2022; v1 submitted 17 April, 2022; originally announced April 2022.

A Gamma-ray Pulsar Timing Array Constrains the Nanohertz Gravitational Wave Background

Authors: M. Ajello, W. B. Atwood, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, A. Berretta, B. Bhattacharyya, E. Bissaldi, R. D. Blandford, E. Bloom, R. Bonino, P. Bruel, R. Buehler, E. Burns, S. Buson, R. A. Cameron, P. A. Caraveo, E. Cavazzuti, N. Cibrario, S. Ciprini, C. J. Clark, I. Cognard, J. Coronado-Blázquez , et al. (107 additional authors not shown)

Abstract: After large galaxies merge, their central supermassive black holes are expected to form binary systems whose orbital motion generates a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background utilize pulsar timing arrays, which perform long-term monitoring of millisecond pulsars (MSPs) at radio wavelengths. We use 12.5 years of Fermi Large Area Telescope data to… ▽ More After large galaxies merge, their central supermassive black holes are expected to form binary systems whose orbital motion generates a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background utilize pulsar timing arrays, which perform long-term monitoring of millisecond pulsars (MSPs) at radio wavelengths. We use 12.5 years of Fermi Large Area Telescope data to form a gamma-ray pulsar timing array. Results from 35 bright gamma-ray pulsars place a 95\% credible limit on the GWB characteristic strain of $1.0\times10^{-14}$ at 1 yr$^{-1}$, which scales as the observing time span $t_{\mathrm{obs}}^{-13/6}$. This direct measurement provides an independent probe of the GWB while offering a check on radio noise models. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: 3 figures in the main text. 3 figures and 8 tables are in the supplementary material

Chromospheric Extension of the MURaM Code

Authors: D. Przybylski, R. Cameron, S. K. Solanki, M. Rempel, J. Leenaarts, L. S. Anusha, V. Witzke, A. I. Shapiro

Abstract: Detailed numerical models of chromosphere and corona are required to understand the heating of the solar atmosphere. An accurate treatment of the solar chromosphere is complicated by the effects arising from Non Local Thermodynamic Equilibrium (NLTE) radiative transfer. A small number of strong, highly scattering lines dominate the cooling and heating in the chromosphere. Additionally, the recombi… ▽ More Detailed numerical models of chromosphere and corona are required to understand the heating of the solar atmosphere. An accurate treatment of the solar chromosphere is complicated by the effects arising from Non Local Thermodynamic Equilibrium (NLTE) radiative transfer. A small number of strong, highly scattering lines dominate the cooling and heating in the chromosphere. Additionally, the recombination times of ionised hydrogen are longer than the dynamical timescales, requiring a non-equilibrium (NE) treatment of hydrogen ionisation. The MURaM code is extended to include the physical process required for accurate simulation of the solar chromosphere, as implemented in the Bifrost code. This includes a time-dependent treatment of hydrogen ionisation, a scattering multi-group radiation transfer scheme and approximations for NLTE radiative cooling. The inclusion of NE and NLTE physics has a large impact on the structure of the chromosphere; the NE treatment of hydrogen ionisation leads to a higher ionisation fraction and enhanced populations in the first excited state throughout cold inter-shock regions of the chromosphere. Additionally this prevents hydrogen ioniation from buffering energy fluctuations, leading to hotter shocks and cooler inter-shock regions. The hydrogen populations in the ground and first excited state are enhanced by $10^2-10^3$ in the upper chromosphere and up to $10^9$ near the transition region. Including the necessary NLTE physics leads to significant differences in chromospheric structure and dynamics. The thermodynamics and hydrogen populations calculated using the extended version of the MURaM code are consistent with previous non-equilibrium simulations. The electron number and temperature calculated using the non-equilibrium treatment of the chromosphere are required to accurately synthesise chromospheric spectral lines. △ Less

Submitted 6 April, 2022; originally announced April 2022.

Comments: 24 pages, 12 figures, to be published in A&A

Journal ref: A&A 664, A91 (2022)

A Strategy for Low-Mass Dark Matter Searches with Cryogenic Detectors in the SuperCDMS SNOLAB Facility

Authors: SuperCDMS Collaboration, M. F. Albakry, I. Alkhatib, D. W. P. Amaral, T. Aralis, T. Aramaki, I. J. Arnquist, I. Ataee Langroudy, E. Azadbakht, S. Banik, C. Bathurst, D. A. Bauer, R. Bhattacharyya, P. L. Brink, R. Bunker, B. Cabrera, R. Calkins, R. A. Cameron, C. Cartaro, D. G. Cerdeno, Y. -Y. Chang, M. Chaudhuri, R. Chen, N. Chott, J. Cooley , et al. (103 additional authors not shown)

Abstract: The SuperCDMS Collaboration is currently building SuperCDMS SNOLAB, a dark matter search focused on nucleon-coupled dark matter in the 1-5 GeV/c$^2$ mass range. Looking to the future, the Collaboration has developed a set of experience-based upgrade scenarios, as well as novel directions, to extend the search for dark matter using the SuperCDMS technology in the SNOLAB facility. The experienced-ba… ▽ More The SuperCDMS Collaboration is currently building SuperCDMS SNOLAB, a dark matter search focused on nucleon-coupled dark matter in the 1-5 GeV/c$^2$ mass range. Looking to the future, the Collaboration has developed a set of experience-based upgrade scenarios, as well as novel directions, to extend the search for dark matter using the SuperCDMS technology in the SNOLAB facility. The experienced-based scenarios are forecasted to probe many square decades of unexplored dark matter parameter space below 5 GeV/c$^2$, covering over 6 decades in mass: 1-100 eV/c$^2$ for dark photons and axion-like particles, 1-100 MeV/c$^2$ for dark-photon-coupled light dark matter, and 0.05-5 GeV/c$^2$ for nucleon-coupled dark matter. They will reach the neutrino fog in the 0.5-5 GeV/c$^2$ mass range and test a variety of benchmark models and sharp targets. The novel directions involve greater departures from current SuperCDMS technology but promise even greater reach in the long run, and their development must begin now for them to be available in a timely fashion. The experienced-based upgrade scenarios rely mainly on dramatic improvements in detector performance based on demonstrated scaling laws and reasonable extrapolations of current performance. Importantly, these improvements in detector performance obviate significant reductions in background levels beyond current expectations for the SuperCDMS SNOLAB experiment. Given that the dominant limiting backgrounds for SuperCDMS SNOLAB are cosmogenically created radioisotopes in the detectors, likely amenable only to isotopic purification and an underground detector life-cycle from before crystal growth to detector testing, the potential cost and time savings are enormous and the necessary improvements much easier to prototype. △ Less

Submitted 1 April, 2023; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: contribution to Snowmass 2021; v2 updated (assorted corrections and improvements to forecasts) October 2022; v3 updated (corrected SuperCDMS SNOLAB sensitivity curves in upgrade forecast plots in body of text) April 2023

Theory of solar oscillations in the inertial frequency range: Linear modes of the convection zone

Authors: Yuto Bekki, Robert H. Cameron, Laurent Gizon

Abstract: On the one hand, several types of global-scale inertial modes of oscillation have been observed on the Sun. They include the equatorial Rossby modes, critical-latitude modes, and high-latitude modes. On the other hand, the columnar convective modes (predicted by simulations; also known as banana cells or thermal Rossby waves) remain elusive. We aim to investigate the influence of turbulent diffusi… ▽ More On the one hand, several types of global-scale inertial modes of oscillation have been observed on the Sun. They include the equatorial Rossby modes, critical-latitude modes, and high-latitude modes. On the other hand, the columnar convective modes (predicted by simulations; also known as banana cells or thermal Rossby waves) remain elusive. We aim to investigate the influence of turbulent diffusivities, non-adiabatic stratification, differential rotation, and a latitudinal entropy gradient on the linear global modes of the rotating solar convection zone. We solve numerically for the eigenmodes of a rotating compressible fluid inside a spherical shell. We identify modes in the inertial frequency range including the columnar convective modes, as well as modes of mixed character. The corresponding mode dispersion relations and eigenfunctions are computed for azimuthal orders $m \leq 16$. The three main results are as follows. Firstly, we find that, for $m \gtrsim 5$, the radial dependence of the equatorial Rossby modes with no radial node ($n=0$) is radically changed from the traditional expectation ($r^m$) for turbulent diffusivities $\gtrsim 10^{12}$ cm$^2$ s$^{-1}$. Secondly, we find mixed modes, i.e. modes that share properties of the equatorial Rossby modes with one radial node ($n=1$) and the columnar convective modes. Thirdly, we show that the $m=1$ high-latitude mode in the model is consistent with the solar observations when the latitudinal entropy gradient corresponding to a thermal wind balance is included (baroclinally unstable mode). To our knowledge, this work is the first realistic eigenvalue calculation of the global modes of the rotating solar convection zone. This calculation reveals a rich spectrum of modes in the inertial frequency range, which can be directly compared to the observations. In turn, the observed modes can inform us about the solar convection zone. △ Less

Submitted 8 March, 2022; originally announced March 2022.

Comments: 24 pages, 29 figures, submitted to Astronomy & Astrophysics on 20 Jan 2022

Journal ref: A&A 662, A16 (2022)

Incremental Fermi Large Area Telescope Fourth Source Catalog

Authors: Fermi-LAT collaboration, :, Soheila Abdollahi, Fabio Acero, Luca Baldini, Jean Ballet, Denis Bastieri, Ronaldo Bellazzini, Bijan Berenji, Alessandra Berretta, Elisabetta Bissaldi, Roger D. Blandford, Elliott Bloom, Raffaella Bonino, Ari Brill, Richard J. Britto, Philippe Bruel, Toby H. Burnett, Sara Buson, Rob A. Cameron, Regina Caputo, Patrizia A. Caraveo, Daniel Castro, Sylvain Chaty, Teddy C. Cheung , et al. (116 additional authors not shown)

Abstract: We present an incremental version (4FGL-DR3, for Data Release 3) of the fourth Fermi-LAT catalog of gamma-ray sources. Based on the first twelve years of science data in the energy range from 50 MeV to 1 TeV, it contains 6658 sources. The analysis improves on that used for the 4FGL catalog over eight years of data: more sources are fit with curved spectra, we introduce a more robust spectral param… ▽ More We present an incremental version (4FGL-DR3, for Data Release 3) of the fourth Fermi-LAT catalog of gamma-ray sources. Based on the first twelve years of science data in the energy range from 50 MeV to 1 TeV, it contains 6658 sources. The analysis improves on that used for the 4FGL catalog over eight years of data: more sources are fit with curved spectra, we introduce a more robust spectral parameterization for pulsars, and we extend the spectral points to 1 TeV. The spectral parameters, spectral energy distributions, and associations are updated for all sources. Light curves are rebuilt for all sources with 1 yr intervals (not 2 month intervals). Among the 5064 original 4FGL sources, 16 were deleted, 112 are formally below the detection threshold over 12 yr (but are kept in the list), while 74 are newly associated, 10 have an improved association, and seven associations were withdrawn. Pulsars are split explicitly between young and millisecond pulsars. Pulsars and binaries newly detected in LAT sources, as well as more than 100 newly classified blazars, are reported. We add three extended sources and 1607 new point sources, mostly just above the detection threshold, among which eight are considered identified, and 699 have a plausible counterpart at other wavelengths. We discuss degree-scale residuals to the global sky model and clusters of soft unassociated point sources close to the Galactic plane, which are possibly related to limitations of the interstellar emission model and missing extended sources. △ Less

Submitted 10 May, 2022; v1 submitted 26 January, 2022; originally announced January 2022.

Comments: accepted in ApJS; follow-up paper to 1902.10045

Journal ref: ApJS 260, 53 (2022)

The crucial role of surface magnetic fields for stellar dynamos: Epsilon Eridani, 61 Cygni A, and the Sun

Authors: S. V. Jeffers, R. H. Cameron, S. C. Marsden, S. Boro Saikia, C. P. Folsom, M. M. Jardine, J. Morin, P. Petit, V. See, A. A. Vidotto, U. Wolter, M. Mittag

Abstract: Cool main-sequence stars, such as the Sun, have magnetic fields which are generated by an internal dynamo mechanism. In the Sun, the dynamo mechanism produces a balance between the amounts of magnetic flux generated and lost over the Sun's 11-year activity cycle and it is visible in the Sun's different atmospheric layers using multi-wavelength observations. We used the same observational diagnosti… ▽ More Cool main-sequence stars, such as the Sun, have magnetic fields which are generated by an internal dynamo mechanism. In the Sun, the dynamo mechanism produces a balance between the amounts of magnetic flux generated and lost over the Sun's 11-year activity cycle and it is visible in the Sun's different atmospheric layers using multi-wavelength observations. We used the same observational diagnostics, spanning several decades, to probe the emergence of magnetic flux on the two close by, active- and low-mass K dwarfs: 61 Cygni A and Epsilon Eridani. Our results show that 61 Cygni A follows the Solar dynamo with a regular cycle at all wavelengths, while Epsilon Eridani represents a more extreme level of the Solar dynamo, while also showing strong Solar-like characteristics. For the first time we show magnetic butterfly diagrams for stars other than the Sun. For the two K stars and the Sun, the rate at which the toroidal field is generated from surface poloidal field is similar to the rate at which toroidal flux is lost through flux emergence. This suggests that the surface field plays a crucial role in the dynamos of all three stars. Finally, for Epsilon Eridani, we show that the two chromospheric cycle periods, of ~3 and ~13 years, correspond to two superimposed magnetic cycles. △ Less

Submitted 19 January, 2022; originally announced January 2022.

Comments: 8 pages, 5 figures: Accepted by A&A

A solar coronal loop in a box: Energy generation and heating

Authors: C. Breu, H. Peter, R. Cameron, S. K. Solanki, D. Przybylski, M. Rempel, L. P. Chitta

Abstract: Coronal loops are the basic building block of the upper solar atmosphere. Comprehending how these are energized, structured, and evolve is key to understanding stellar coronae. Here we investigate how the energy to heat the loop is generated by photospheric magneto-convection, transported into the upper atmosphere, and how the internal structure of a coronal loop forms. In a 3D magnetohydrodyn… ▽ More Coronal loops are the basic building block of the upper solar atmosphere. Comprehending how these are energized, structured, and evolve is key to understanding stellar coronae. Here we investigate how the energy to heat the loop is generated by photospheric magneto-convection, transported into the upper atmosphere, and how the internal structure of a coronal loop forms. In a 3D magnetohydrodynamics (MHD) model, we study an isolated coronal loop rooted with both footpoints in a shallow layer within the convection zone using the MURaM code. To resolve its internal structure, we limited the computational domain to a rectangular box containing a single coronal loop as a straightened magnetic flux tube. Field-aligned heat conduction, gray radiative transfer in the photosphere and chromosphere, and optically thin radiative losses in the corona were taken into account. The footpoints were allowed to interact self-consistently with the granulation surrounding them. The loop is heated by a Poynting flux that is self-consistently generated through small-scale motions within individual magnetic concentrations in the photosphere. Turbulence develops in the upper layers of the atmosphere as a response to the footpoint motions. We see little sign of heating by large-scale braiding of magnetic flux tubes from different photospheric concentrations at a given footpoint. The synthesized emission, as it would be observed by the Atmospheric Imaging Assembly or the X-ray Telescope, reveals transient bright strands that form in response to the heating events. Overall, our model roughly reproduces the properties and evolution of the plasma as observed within coronal loops. With this model we can build a coherent picture of how the energy flux to heat the upper atmosphere is generated near the solar surface and how this process drives and governs the heating and dynamics of a coronal loop. △ Less

Submitted 21 December, 2021; originally announced December 2021.

Comments: 19 pages, 16 figures

Testing solar surface flux transport models in the first days after active region emergence

Authors: Nils Gottschling, Hannah Schunker, Aaron C. Birch, Robert Cameron, Laurent Gizon

Abstract: Active regions (ARs) play an important role in the magnetic dynamics of the Sun. Solar surface flux transport models (SFTMs) are used to describe the evolution of the radial magnetic field at the solar surface. There is however uncertainty about using these models in the early stage of AR evolution. We aim to test the applicability of SFTMs in the first days after the emergence of ARs by comparing… ▽ More Active regions (ARs) play an important role in the magnetic dynamics of the Sun. Solar surface flux transport models (SFTMs) are used to describe the evolution of the radial magnetic field at the solar surface. There is however uncertainty about using these models in the early stage of AR evolution. We aim to test the applicability of SFTMs in the first days after the emergence of ARs by comparing them with observations. The models we employ range from passive evolution to models where the inflows around ARs are included. We simulate the evolution of the surface magnetic field of 17 emerging active regions using a local surface flux transport simulation. We selected regions that do not form fully-fledged sunspots that exhibit moat flows. The simulation includes diffusion and advection. We use observed flows from local correlation tracking of solar granulation, as well as parametrizations of the inflows around ARs. To evaluate our simulations, we measure the cross correlation between the observed and the simulated magnetic field, as well as the total unsigned flux of the ARs, over time. We also test the validity of our simulations by varying the starting time relative to the emergence of flux. We find that the simulations using observed surface flows can reproduce the evolution of the observed magnetic flux. The effect of buffeting of the field by supergranulation can be described as a diffusion process. The SFTM is applicable after 90% of the peak total unsigned flux of the AR has emerged. Diffusivities in the range between $D=250$ to $720$ km$^2$/s are consistent with the evolution of the AR flux in the first five days after this time. We find that the converging flows around emerging ARs are not important for the evolution of the total flux of the AR in these first five days; their effect of increasing flux cancellation is balanced by the decrease of flux transport away from the AR. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: Accepted for publication in Astronomy & Astrophysics. 13 pages, 11 figures

Journal ref: A&A 660, A6 (2022)

Prototype Schwarzschild-Couder Telescope for the Cherenkov Telescope Array: Commissioning the Optical System

Authors: C. B. Adams, G. Ambrosi, M. Ambrosio, C. Aramo, P. I. Batista, W. Benbow, B. Bertucci, E. Bissaldi, M. Bitossi, A. Boiano, C. Bonavolontà, R. Bose, A. Brill, J. H. Buckley, R. A. Cameron, R. Canestrari, M. Capasso, M. Caprai, C. E. Covault, D. Depaoli, L. Di Venere, M. Errando, S. Fegan, Q. Feng, E. Fiandrini , et al. (47 additional authors not shown)

Abstract: A prototype Schwarzschild-Couder Telescope (pSCT) has been constructed at the Fred Lawrence Whipple Observatory as a candidate for the medium-sized telescopes of the Cherenkov Telescope Array Observatory (CTAO). CTAO is currently entering early construction phase of the project and once completed it will vastly improve very high energy gamma-ray detection component in multi-wavelength and multi-me… ▽ More A prototype Schwarzschild-Couder Telescope (pSCT) has been constructed at the Fred Lawrence Whipple Observatory as a candidate for the medium-sized telescopes of the Cherenkov Telescope Array Observatory (CTAO). CTAO is currently entering early construction phase of the project and once completed it will vastly improve very high energy gamma-ray detection component in multi-wavelength and multi-messenger observations due to significantly improved sensitivity, angular resolution and field of view comparing to the current generation of the ground-based gamma-ray observatories H.E.S.S., MAGIC and VERITAS. The pSCT uses a dual aspheric mirror design with a $9.7$ m primary mirror and $5.4$ m secondary mirror, both of which are segmented. The Schwarzschild-Couder (SC) optical system (OS) selected for the prototype telescope achieves wide field of view of $8$ degrees and simultaneously reduces the focal plane plate scale allowing an unprecedented compact ($0.78$m diameter) implementation of the high-resolution camera ($6$mm/ $0.067$deg per imaging pixel with $11,328$ pixels) based on the silicon photo-multipliers (SiPMs). The OS of the telescope is designed to eliminate spherical and comatic aberrations and minimize astigmatism to radically improve off-axis imaging and consequently angular resolution across all the field of view with respect to the conventional single-mirror telescopes. Fast and high imaging resolution OS of the pSCT comes with the challenging submillimeter-precision custom alignment system, which was successfully demonstrated with an on-axis point spread function (PSF) of $2.9$ arcmin prior to the first-light detection of the Crab Nebula in 2020. Ongoing and future commissioning activities are reported. △ Less

Submitted 14 October, 2021; originally announced October 2021.

Journal ref: Proceedings of Science, PoS(ICRC2021)717

Detection of the Crab Nebula by the prototype Schwarzschild-Couder Telescope

Authors: C. B. Adams, G. Ambrosi, M. Ambrosio, C. Aramo, P. I. Batista, W. Benbow, B. Bertucci, E. Bissaldi, M. Bitossi, A. Boiano, C. Bonavolontà, R. Bose, A. Brill, A. M. Brown, J. H. Buckley, R. A. Cameron, R. Canestrari, M. Capasso, M. Caprai, C. E. Covault, D. Depaoli, L. Di Venere, M. Errando, S. Fegan, Q. Feng , et al. (49 additional authors not shown)

Abstract: The Schwarzschild-Couder Telescope (SCT) is a medium-sized telescope technology proposed for the Cherenkov Telescope Array. It uses a novel dual-mirror optical design that removes comatic aberrations across its entire field of view. The SCT camera employs high-resolution silicon photomultiplier (SiPM) sensors with a pixel size of 4 arcminutes. A prototype SCT (pSCT) has been constructed at the Fre… ▽ More The Schwarzschild-Couder Telescope (SCT) is a medium-sized telescope technology proposed for the Cherenkov Telescope Array. It uses a novel dual-mirror optical design that removes comatic aberrations across its entire field of view. The SCT camera employs high-resolution silicon photomultiplier (SiPM) sensors with a pixel size of 4 arcminutes. A prototype SCT (pSCT) has been constructed at the Fred Lawrence Whipple Observatory in Arizona, USA. An observing campaign in 2020, with a partial camera of 1600 pixels (2.7 degrees by 2.7 degrees field of view) resulted in detection of the Crab Nebula at 8.6 sigma statistical significance. Work on the pSCT camera and optical system is ongoing to improve performance and prepare for an upcoming camera upgrade. The pSCT camera upgrade will replace the current camera modules with improved SiPMs and readout electronics and will expand the camera to its full design field of view of 8 degrees in diameter (11,328 pixels). The fully upgraded pSCT will enable next-generation very-high-energy gamma-ray astrophysics through excellent background rejection and angular resolution. In this presentation we describe first results from the successful operation of the pSCT and future plans. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 9 pages, 3 figures, 2 tables, contribution to ICRC 2021, similar to 10.1016/j.astropartphys.2021.102562 (arXiv:2012.08448)

Journal ref: Proceedings of Science, PoS(ICRC2021)830

Design and performance of the prototype Schwarzschild-Couder Telescope camera

Authors: C. B. Adams, G. Ambrosi, M. Ambrosio, C. Aramo, P. I. Batista, W. Benbow, B. Bertucci, E. Bissaldi, M. Bitossi, A. Boiano, C. Bonavolonta, R. Bose, A. Brill, A. M. Brown, J. H. Buckley, R. A. Cameron, M. Capasso, M. Caprai, C. E. Covault, D. Depaoli, L. Di Venere, M. Errando, S. Fegan, Q. Feng, E. Fiandrini , et al. (49 additional authors not shown)

Abstract: The Cherenkov Telescope Array (CTA) is the next-generation ground-based observatory for very-high-energy gamma-ray astronomy. An innovative 9.7 m aperture, dual-mirror Schwarzschild-Couder Telescope (SCT) design is a candidate design for CTA Medium-Sized Telescopes. A prototype SCT (pSCT) has been constructed at the Fred Lawrence Whipple Observatory in Arizona, USA. Its camera is currently partial… ▽ More The Cherenkov Telescope Array (CTA) is the next-generation ground-based observatory for very-high-energy gamma-ray astronomy. An innovative 9.7 m aperture, dual-mirror Schwarzschild-Couder Telescope (SCT) design is a candidate design for CTA Medium-Sized Telescopes. A prototype SCT (pSCT) has been constructed at the Fred Lawrence Whipple Observatory in Arizona, USA. Its camera is currently partially instrumented with 1600 pixels covering a field of view of 2.7 degrees square. The small plate scale of the optical system allows densely packed silicon photomultipliers to be used, which combined with high-density trigger and waveform readout electronics enable the high-resolution camera. The camera's electronics are capable of imaging air shower development at a rate of one billion samples per second. We describe the commissioning and performance of the pSCT camera, including trigger and waveform readout performance, calibration, and absolute GPS time stamping. We also present the upgrade to the camera, which is currently underway. The upgrade will fully populate the focal plane, increasing the field of view to 8 degree diameter, and lower the front-end electronics noise, enabling a lower trigger threshold and improved reconstruction and background rejection. △ Less

Submitted 10 September, 2021; originally announced September 2021.

Comments: 8 pages, 5 figures, Proceedings of the 37th International Cosmic Ray Conference (ICRC 2021), Berlin, Germany

Solar inertial modes: Observations, identification, and diagnostic promise

Authors: Laurent Gizon, Robert H. Cameron, Yuto Bekki, Aaron C. Birch, Richard S. Bogart, Allan Sacha Brun, Cilia Damiani, Damien Fournier, Laura Hyest, Kiran Jain, B. Lekshmi, Zhi-Chao Liang, Bastian Proxauf

Abstract: The oscillations of a slowly rotating star have long been classified into spheroidal and toroidal modes. The spheroidal modes include the well-known 5-min acoustic modes used in helioseismology. Here we report observations of the Sun's toroidal modes, for which the restoring force is the Coriolis force and whose periods are on the order of the solar rotation period. By comparing the observations w… ▽ More The oscillations of a slowly rotating star have long been classified into spheroidal and toroidal modes. The spheroidal modes include the well-known 5-min acoustic modes used in helioseismology. Here we report observations of the Sun's toroidal modes, for which the restoring force is the Coriolis force and whose periods are on the order of the solar rotation period. By comparing the observations with the normal modes of a differentially rotating spherical shell, we are able to identify many of the observed modes. These are the high-latitude inertial modes, the critical-latitude inertial modes, and the equatorial Rossby modes. In the model, the high-latitude and critical-latitude modes have maximum kinetic energy density at the base of the convection zone, and the high-latitude modes are baroclinically unstable due to the latitudinal entropy gradient. As a first application of inertial-mode helioseismology, we constrain the superadiabaticity and the turbulent viscosity in the deep convection zone. △ Less

Submitted 20 July, 2021; originally announced July 2021.

Comments: Forthcoming article, Astronomy and Astrophysics

Fermi Large Area Telescope Performance After 10 Years Of Operation

Authors: The Fermi LAT Collaboration, M. Ajello, W. B. Atwood, M. Axelsson, R. Bagagli, M. Bagni, L. Baldini, D. Bastieri, F. Bellardi, R. Bellazzini, E. Bissaldi, E. D. Bloom, R. Bonino, J. Bregeon, A. Brez, P. Bruel, R. Buehler, S. Buson, R. A. Cameron, P. A. Caraveo, E. Cavazzuti, M. Ceccanti, S. Chen, C. C. Cheung, S. Ciprini , et al. (104 additional authors not shown)

Abstract: The Large Area Telescope (LAT), the primary instrument for the Fermi Gamma-ray Space Telescope (Fermi) mission, is an imaging, wide field-of-view, high-energy gamma-ray telescope, covering the energy range from 30 MeV to more than 300 GeV. We describe the performance of the instrument at the 10-year milestone. LAT performance remains well within the specifications defined during the planning phase… ▽ More The Large Area Telescope (LAT), the primary instrument for the Fermi Gamma-ray Space Telescope (Fermi) mission, is an imaging, wide field-of-view, high-energy gamma-ray telescope, covering the energy range from 30 MeV to more than 300 GeV. We describe the performance of the instrument at the 10-year milestone. LAT performance remains well within the specifications defined during the planning phase, validating the design choices and supporting the compelling case to extend the duration of the Fermi mission. The details provided here will be useful when designing the next generation of high-energy gamma-ray observatories. △ Less

Submitted 6 September, 2021; v1 submitted 23 June, 2021; originally announced June 2021.

Comments: 60 pages, 28 figures. Published in ApJS

Journal ref: ApJS 256 (2021) 12

Catalog of Long-Term Transient Sources in the First 10 Years of Fermi-LAT Data

Authors: L. Baldini, J. Ballet, D. Bastieri, J. Becerra Gonzalez, R. Bellazzini, A. Berretta, E. Bissaldi, R. D. Blandford, E. D. Bloom, R. Bonino, E. Bottacini, P. Bruel, S. Buson, R. A. Cameron, P. A. Caraveo, E. Cavazzuti, S. Chen, G. Chiaro, D. Ciangottini, S. Ciprini, P. Cristarella Orestano, M. Crnogorcevic, S. Cutini, F. D'Ammando, P. de la Torre Luque , et al. (90 additional authors not shown)

Abstract: We present the first Fermi Large Area Telescope (LAT) catalog of long-term $γ$-ray transient sources (1FLT). This comprises sources that were detected on monthly time intervals during the first decade of Fermi-LAT operations. The monthly time scale allows us to identify transient and variable sources that were not yet reported in other Fermi-LAT catalogs. The monthly datasets were analyzed using a… ▽ More We present the first Fermi Large Area Telescope (LAT) catalog of long-term $γ$-ray transient sources (1FLT). This comprises sources that were detected on monthly time intervals during the first decade of Fermi-LAT operations. The monthly time scale allows us to identify transient and variable sources that were not yet reported in other Fermi-LAT catalogs. The monthly datasets were analyzed using a wavelet-based source detection algorithm that provided the candidate new transient sources. The search was limited to the extragalactic regions of the sky to avoid the dominance of the Galactic diffuse emission at low Galactic latitudes. The transient candidates were then analyzed using the standard Fermi-LAT Maximum Likelihood analysis method. All sources detected with a statistical significance above 4$σ$ in at least one monthly bin were listed in the final catalog. The 1FLT catalog contains 142 transient $γ$-ray sources that are not included in the 4FGL-DR2 catalog. Many of these sources (102) have been confidently associated with Active Galactic Nuclei (AGN): 24 are associated with Flat Spectrum Radio Quasars; 1 with a BL Lac object; 70 with Blazars of Uncertain Type; 3 with Radio Galaxies; 1 with a Compact Steep Spectrum radio source; 1 with a Steep Spectrum Radio Quasar; 2 with AGN of other types. The remaining 40 sources have no candidate counterparts at other wavelengths. The median $γ$-ray spectral index of the 1FLT-AGN sources is softer than that reported in the latest Fermi-LAT AGN general catalog. This result is consistent with the hypothesis that detection of the softest $γ$-ray emitters is less efficient when the data are integrated over year-long intervals. △ Less

Submitted 31 May, 2021; originally announced June 2021.

Comments: 41 pages, 17 figures, 7 tables; Accepted by ApJS on 24 May 2021; Contact Authors: I. Mereu, S. Cutini, E. Cavazzuti, G. Tosti

Slow magneto-acoustic waves in simulations of a solar plage region carry enough energy to heat the chromosphere

Authors: Nitin Yadav, Robert H. Cameron, Sami K. Solanki

Abstract: We study the properties of slow magneto-acoustic waves that are naturally excited due to turbulent convection and investigate their role in the energy balance of a plage region using three dimensional (3D) radiation-MHD simulations. We calculate the horizontally averaged (over the whole domain) frequency power spectra for both longitudinal and vertical (i.e. the component perpendicular to the surf… ▽ More We study the properties of slow magneto-acoustic waves that are naturally excited due to turbulent convection and investigate their role in the energy balance of a plage region using three dimensional (3D) radiation-MHD simulations. We calculate the horizontally averaged (over the whole domain) frequency power spectra for both longitudinal and vertical (i.e. the component perpendicular to the surface) components of velocity. To compare our results with the observations we degrade the simulation data with Gaussian kernels having FWHM of 100 km and 200 km, and calculate horizontally averaged power spectra for the vertical component of velocity. The power spectra of the longitudinal component of velocity, averaged over field lines in the core of a kG magnetic flux concentration, reveal that the dominant period of oscillations shifts from around 6.5 minutes in the photosphere to around 4 minutes in the chromosphere. At the same time, the velocity power spectra, averaged horizontally over the whole domain, show that low frequency waves (approximately 6.5 minute period) may reach well into the chromosphere. Importantly, waves with frequencies above 5 mHz propagating along different field lines are found to be out of phase with each other even within a single magnetic concentration. The horizontally averaged power spectra of the vertical component of velocity at various effective resolutions show that the observed acoustic wave energy fluxes are underestimated, by a factor of three even if determined from observations carried out at a high spatial resolution of 200 km. Our results show that longitudinal waves carry (just) sufficient energy to heat the chromosphere in solar plage. We conjecture that current observations (with spatial resolution around 200 km) underestimate the energy flux by roughly a factor of three, or more if the observations have lower spatial resolution. △ Less

Submitted 6 May, 2021; originally announced May 2021.

Comments: Accepted for publication in A & A

Journal ref: A&A 652, A43 (2021)

Non-equilibrium Equation of State in stellar atmospheres

Authors: L. S. Anusha, M. van Noort, R. H. Cameron

Abstract: In the stellar chromospheres, radiative energy transport is dominated by only the strongest spectral lines. For these lines, the approximation of local thermodynamic equilibrium (LTE) is known to be very inaccurate, and a state of equilibrium cannot be assumed in general. To calculate the radiative energy transport under these conditions, the population evolution equation must be evaluated explici… ▽ More In the stellar chromospheres, radiative energy transport is dominated by only the strongest spectral lines. For these lines, the approximation of local thermodynamic equilibrium (LTE) is known to be very inaccurate, and a state of equilibrium cannot be assumed in general. To calculate the radiative energy transport under these conditions, the population evolution equation must be evaluated explicitly, including all time-dependent terms. We develop a numerical method to solve the evolution equation for the atomic-level populations in a time-implicit way, keeping all time-dependent terms to first order. We show that the linear approximation of the time dependence of the populations can handle very large time steps without losing the accuracy. We reproduce the benchmark solutions from earlier, well-established works in terms of non-LTE kinetic equilibrium solution and typical ionization/recombination time-scales in the solar chromosphere. △ Less

Submitted 28 April, 2021; originally announced April 2021.

Comments: 42 pages, 8 figures

Journal ref: ApJ, 911, 71 (2021)

A journey of exploration to the polar regions of a star: probing the solar poles and the heliosphere from high helio-latitude

Authors: Louise Harra, Vincenzo Andretta, Thierry Appourchaux, Frédéric Baudin, Luis Bellot-Rubio, Aaron C. Birch, Patrick Boumier, Robert H. Cameron, Matts Carlsson, Thierry Corbard, Jackie Davies, Andrew Fazakerley, Silvano Fineschi, Wolfgang Finsterle, Laurent Gizon, Richard Harrison, Donald M. Hassler, John Leibacher, Paulett Liewer, Malcolm MacDonald, Milan Maksimovic, Neil Murphy, Giampiero Naletto, Giuseppina Nigro, Christopher Owen , et al. (7 additional authors not shown)

Abstract: A mission to view the solar poles from high helio-latitudes (above 60$^\circ$) will build on the experience of Solar Orbiter as well as a long heritage of successful solar missions and instrumentation (e.g. SOHO \cite{SOHO}, STEREO \cite{stereo}, Hinode \cite{Hinode}, SDO \cite{SDO}), but will focus for the first time on the solar poles, enabling scientific investigations that cannot be done by an… ▽ More A mission to view the solar poles from high helio-latitudes (above 60$^\circ$) will build on the experience of Solar Orbiter as well as a long heritage of successful solar missions and instrumentation (e.g. SOHO \cite{SOHO}, STEREO \cite{stereo}, Hinode \cite{Hinode}, SDO \cite{SDO}), but will focus for the first time on the solar poles, enabling scientific investigations that cannot be done by any other mission. One of the major mysteries of the Sun is the solar cycle. The activity cycle of the Sun drives the structure and behaviour of the heliosphere and is, of course, the driver of space weather. In addition, solar activity and variability provides fluctuating input into the Earth climate models, and these same physical processes are applicable to stellar systems hosting exoplanets. One of the main obstructions to understanding the solar cycle, and hence all solar activity, is our current lack of understanding of the polar regions. In this White Paper, submitted to the European Space Agency in response to the Voyage 2050 call, we describe a mission concept that aims to address this fundamental issue. In parallel, we recognise that viewing the Sun from above the polar regions enables further scientific advantages, beyond those related to the solar cycle, such as unique and powerful studies of coronal mass ejection processes, from a global perspective, and studies of coronal structure and activity in polar regions. Not only will these provide important scientific advances for fundamental stellar physics research, they will feed into our understanding of impacts on the Earth and other planets' space environment. △ Less

Submitted 22 April, 2021; originally announced April 2021.

Comments: 29 pages, 9 figures, White paper for ESA Voyage 2050

On the size distribution of spots within sunspot groups

Authors: Sudip Mandal, Natalie A. Krivova, Robert Cameron, Sami K. Solanki

Abstract: Size distribution of sunspots provides key information about the generation and emergence processes of the solar magnetic field. Previous studies on the size distribution have primarily focused on either the whole group or individual spot areas. In this paper, we investigate the organization of spot areas within sunspot groups. In particular, we analyze the ratio, $\rm{R}$, of the area of the bigg… ▽ More Size distribution of sunspots provides key information about the generation and emergence processes of the solar magnetic field. Previous studies on the size distribution have primarily focused on either the whole group or individual spot areas. In this paper, we investigate the organization of spot areas within sunspot groups. In particular, we analyze the ratio, $\rm{R}$, of the area of the biggest spot ($\rm{A_{big\_spot}}$) inside a group, to the total area of that group ($\rm{A_{group}}$). We use sunspot observations from Kislovodsk, Pulkovo and Debrecen observatories, together covering solar cycles 17 to 24. We find that at the time when the group area reaches its maximum, the single biggest spot in a group typically occupies about 60% of the group area. For half of all groups, $\rm R$ lies in the range between roughly 50% and 70%. We also find R to change with the group area, $\rm{A_{group}}$, such that $\rm{R}$ reaches a maximum of about 0.65 for groups with $\rm{A_{group}}\approx 200μ$Hem and then remains at about 0.6 for lager groups. Our findings imply a scale invariant emergence pattern, providing an observational constraint on the emergence process. Furthermore, extrapolation of our results to larger sunspot groups may have a bearing on the giant unresolved starspot features found in Doppler images of highly active sun-like stars. Our results suggest that such giant features are composed of multiple spots, with the largest spot occupying roughly 55--75% of the total group area (i.e. of the area of the giant starspots seen in Doppler images). △ Less

Submitted 8 April, 2021; originally announced April 2021.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 652, A9 (2021)

Modelling the evolution of the Sun's open and total magnetic flux

Authors: N. A. Krivova, S. K. Solanki, B. Hofer, C. -J. Wu, I. G. Usoskin, R. Cameron

Abstract: Solar activity in all its varied manifestations is driven by the magnetic field. Particularly important for many purposes are two global quantities, the Sun's total and open magnetic flux, which can be computed from sunspot number records using models. Such sunspot-driven models, however, do not take into account the presence of magnetic flux during grand minima, such as the Maunder minimum. Here… ▽ More Solar activity in all its varied manifestations is driven by the magnetic field. Particularly important for many purposes are two global quantities, the Sun's total and open magnetic flux, which can be computed from sunspot number records using models. Such sunspot-driven models, however, do not take into account the presence of magnetic flux during grand minima, such as the Maunder minimum. Here we present a major update of a widely used simple model, which now takes into account the observation that the distribution of all magnetic features on the Sun follows a single power law. The exponent of the power law changes over the solar cycle. This allows for the emergence of small-scale magnetic flux even when no sunspots are present for multiple decades and leads to non-zero total and open magnetic flux also in the deepest grand minima, such as the Maunder minimum, thus overcoming a major shortcoming of the earlier models. The results of the updated model compare well with the available observations and reconstructions of the solar total and open magnetic flux. This opens up the possibility of improved reconstructions of sunspot number from time series of cosmogenic isotope production rate. △ Less

Submitted 29 March, 2021; originally announced March 2021.

Journal ref: A&A 650, A70 (2021)

The First Fermi-LAT Solar Flare Catalog

Authors: M. Ajello, L. Baldini, D. Bastieri, R. Bellazzini, A. Berretta, E. Bissaldi, R. D. Blandford, R. Bonino, P. Bruel, S. Buson, R. A. Cameron, R. Caputo, E. Cavazzuti, C. C. Cheung, G. Chiaro, D. Costantin, S. Cutini, F. D'Ammando, F. de Palma, R. Desiante, N. Di Lalla, L. Di Venere, F. Fana Dirirsa, S. J. Fegan, Y. Fukazawa , et al. (60 additional authors not shown)

Abstract: We present the first Fermi - Large Area Telescope (LAT) solar flare catalog covering the 24 th solar cycle. This catalog contains 45 Fermi -LAT solar flares (FLSFs) with emission in the gamma-ray energy band (30 MeV - 10 GeV) detected with a significance greater than 5 sigma over the years 2010-2018. A subsample containing 37 of these flares exhibit delayed emission beyond the prompt-impulsive har… ▽ More We present the first Fermi - Large Area Telescope (LAT) solar flare catalog covering the 24 th solar cycle. This catalog contains 45 Fermi -LAT solar flares (FLSFs) with emission in the gamma-ray energy band (30 MeV - 10 GeV) detected with a significance greater than 5 sigma over the years 2010-2018. A subsample containing 37 of these flares exhibit delayed emission beyond the prompt-impulsive hard X-ray phase with 21 flares showing delayed emission lasting more than two hours. No prompt-impulsive emission is detected in four of these flares. We also present in this catalog the observations of GeV emission from 3 flares originating from Active Regions located behind the limb (BTL) of the visible solar disk. We report the light curves, spectra, best proton index and localization (when possible) for all the FLSFs. The gamma-ray spectra is consistent with the decay of pions produced by >300 MeV protons. This work contains the largest sample of high-energy gamma-ray flares ever reported and provides the unique opportunity to perform population studies on the different phases of the flare and thus allowing to open a new window in solar physics. △ Less

Submitted 25 January, 2021; originally announced January 2021.

Comments: Published in ApJS

Journal ref: ApJS 252 13 (2021)

Detection of the Crab Nebula with the 9.7 m Prototype Schwarzschild-Couder Telescope

Authors: C. B. Adams, R. Alfaro, G. Ambrosi, M. Ambrosio, C. Aramo, T. Arlen, P. I. Batista, W. Benbow, B. Bertucci, E. Bissaldi, J. Biteau, M. Bitossi, A. Boiano, C. Bonavolontà, R. Bose, A. Bouvier, A. Brill, A. M. Brown, J. H. Buckley, K. Byrum, R. A. Cameron, R. Canestrari, M. Capasso, M. Caprai, C. E. Covault , et al. (83 additional authors not shown)

Abstract: The Schwarzschild-Couder Telescope (SCT) is a telescope concept proposed for the Cherenkov Telescope Array. It employs a dual-mirror optical design to remove comatic aberrations over an $8^{\circ}$ field of view, and a high-density silicon photomultiplier camera (with a pixel resolution of 4 arcmin) to record Cherenkov emission from cosmic ray and gamma-ray initiated particle cascades in the atmos… ▽ More The Schwarzschild-Couder Telescope (SCT) is a telescope concept proposed for the Cherenkov Telescope Array. It employs a dual-mirror optical design to remove comatic aberrations over an $8^{\circ}$ field of view, and a high-density silicon photomultiplier camera (with a pixel resolution of 4 arcmin) to record Cherenkov emission from cosmic ray and gamma-ray initiated particle cascades in the atmosphere. The prototype SCT (pSCT), comprising a 9.7 m diameter primary mirror and a partially instrumented camera with 1536 pixels, has been constructed at the Fred Lawrence Whipple Observatory. The telescope was inaugurated in January 2019, with commissioning continuing throughout 2019. We describe the first campaign of observations with the pSCT, conducted in January and February of 2020, and demonstrate the detection of gamma-ray emission from the Crab Nebula with a statistical significance of $8.6σ$. △ Less

Submitted 15 December, 2020; originally announced December 2020.

Comments: 13 pages, 12 figures, 3 tables, submitted to Astroparticle Physics

Sunspot simulations: penumbra formation and the fluting instability

Authors: Mayukh Panja, Robert Cameron, Sami K. Solanki

Abstract: The fluting instability has been suggested as the driver of the subsurface structure of sunspot flux tubes. We conducted a series of numerical experiments where we used flux tubes with different initial curvatures to study the effect of the fluting instability on the subsurface structure of spots. We used the MURaM code, which has previously been used to simulate complete sunspots, to first comput… ▽ More The fluting instability has been suggested as the driver of the subsurface structure of sunspot flux tubes. We conducted a series of numerical experiments where we used flux tubes with different initial curvatures to study the effect of the fluting instability on the subsurface structure of spots. We used the MURaM code, which has previously been used to simulate complete sunspots, to first compute four sunspots in the slab geometry and then two complete circular spots of opposite polarities. We find that the curvature of a flux tube indeed determines the degree of fluting the flux tube will undergo - the more curved a flux tube is, the more fluted it becomes. In addition, sunspots with strong curvature have strong horizontal fields at the surface and therefore readily form penumbral filaments. The fluted sunspots eventually break up from below, with lightbridges appearing at the surface several hours after fluting commences. △ Less

Submitted 23 November, 2020; originally announced November 2020.

Comments: 8 pages, 5 figures, Accepted to ApJ