-
Small-scale and large-scale dynamos in global convection simulations of solar-like stars
Authors:
Jörn Warnecke,
Maarit J. Korpi-Lagg,
Matthias Rheinhard,
Mariangela Viviani,
Ameya Prabhu
Abstract:
It has been recently shown that a small-scale dynamo (SSD) instability could be possible in solar-like low magnetic Prandtl number Pm plasmas. It has been proposed that the presence of SSD can potentially have a significant impact on the dynamics of the large-scale dynamo (LSD) in the stellar convection zones. Studying these two dynamos, SSD and LSD, together in a global magnetoconvection model re…
▽ More
It has been recently shown that a small-scale dynamo (SSD) instability could be possible in solar-like low magnetic Prandtl number Pm plasmas. It has been proposed that the presence of SSD can potentially have a significant impact on the dynamics of the large-scale dynamo (LSD) in the stellar convection zones. Studying these two dynamos, SSD and LSD, together in a global magnetoconvection model requires high-resolution simulations and large amounts of computational resources. Starting from a well-studied global convective dynamo model that produces cyclic magnetic fields, we systematically increased the resolution and lowered the diffusivities to enter the regime of Reynolds numbers that allow for the excitation of SSD on top of the LSD. We studied how the properties of convection, generated differential rotation profiles, and LSD solutions change with the presence of SSD. We performed convective dynamo simulations in a spherical wedge with the Pencil Code. The resolutions of the models were increased in 4 steps by a total factor of 16 to achieve maximal fluid and magnetic Reynolds numbers of over 500. We found that the differential rotation is strongly quenched by the presence of the LSD and SSD. Even though the small-scale magnetic field only mildly decreases increasing Re, the large-scale field strength decreases significantly. We do not find the SSD dynamo significantly quenching the convective flows as claimed recently by other authors; in contrast, the convective flows first grow and then saturate for increasing Re. Furthermore, the angular momentum transport is highly affected by the presence of small-scale magnetic fields, which are mostly generated by LSD. These fields not only change the Reynolds stresses, but also generate dynamically important Maxwell stresses. The LSD evolution in terms of its pattern and field distribution is rather independent of the increase in Rm.
△ Less
Submitted 13 June, 2024;
originally announced June 2024.
-
From convective stellar dynamo simulations to Zeeman-Doppler images
Authors:
Thomas Hackman,
Oleg Kochukhov,
Mariangela Viviani,
Jörn Warnecke,
Maarit J. Korpi-Lagg,
Jyri J. Lehtinen
Abstract:
Zeeman-Doppler imaging (ZDI) is used to reconstruct the surface magnetic field of late-type stars from high resolution spectropolarimetric observations. The results are usually described in terms of characteristics of the field topology, i.e. poloidality vs. toroidality and axi-symmetry vs. non-axisymmetry in addition to the field strength. We want to test how well these characteristics are preser…
▽ More
Zeeman-Doppler imaging (ZDI) is used to reconstruct the surface magnetic field of late-type stars from high resolution spectropolarimetric observations. The results are usually described in terms of characteristics of the field topology, i.e. poloidality vs. toroidality and axi-symmetry vs. non-axisymmetry in addition to the field strength. We want to test how well these characteristics are preserved when applying the ZDI method on simulated data, i.e. how accurately the field topology is preserved and to what extent stellar parameters influence the reconstruction. We use published magnetic field data from direct numerical MHD simulations. These have variable rotation rates, and hence represent different levels of activity, of an otherwise Sun-like setup. Our ZDI reconstruction is based on spherical harmonics expansion. By comparing the original values to those of the reconstructed images, we study the ability to reconstruct the surface magnetic field in terms of various characteristics of the field. The main large-scale features are reasonably well recovered, but the strength of the recovered magnetic field is just a fraction of the original input. The quality of the reconstruction shows clear correlations with the data quality. Furthermore, there are some spurious dependencies between stellar parameters and the characteristics of the field. Our study uncovers some limits of ZDI. Firstly, the recovered field strength will generally be lower than the "real" value as smaller structures with opposite polarities will be blurred in the inversion. Secondly, the axisymmetry is overestimated. The poloidality vs. toroidality is better recovered. The reconstruction works better for a stronger field and faster rotation velocity. Still, the ZDI method works surprisingly well even for a weaker field and slow rotation, provided the data has a high signal-to-noise and good rotation phase coverage.
△ Less
Submitted 22 January, 2024; v1 submitted 13 June, 2023;
originally announced June 2023.
-
Turbulent magneto-genesis in a collisionless plasma
Authors:
F. Pucci,
M. Viviani,
F. Valentini,
G. Lapenta,
W. H. Matthaeus,
S. Servidio
Abstract:
We investigate an efficient mechanism for generating magnetic fields in turbulent, collisionless plasmas. By using fully kinetic, particle-in-cell simulations of an initially non-magnetized plasma, we inspect the genesis of magnetization, in a nonlinear regime. The complex motion is initiated via a Taylor-Green vortex, and the plasma locally develops strong electron temperature anisotropy, due to…
▽ More
We investigate an efficient mechanism for generating magnetic fields in turbulent, collisionless plasmas. By using fully kinetic, particle-in-cell simulations of an initially non-magnetized plasma, we inspect the genesis of magnetization, in a nonlinear regime. The complex motion is initiated via a Taylor-Green vortex, and the plasma locally develops strong electron temperature anisotropy, due to the strain tensor of the turbulent flow. Subsequently, in a domino effect, the anisotropy triggers a Weibel instability, localized in space. In such active wave-particle interaction regions, the magnetic field seed grows exponentially and spreads to larger scales due to the interaction with the underlying stirring motion. Such a self-feeding process might explain magneto-genesis in a variety of astrophysical plasmas, wherever turbulence is present.
△ Less
Submitted 15 September, 2021;
originally announced September 2021.
-
Implementation and Optimization of the PTOLEMY Transverse Drift Electromagnetic Filter
Authors:
A. Apponi,
M. G. Betti,
M. Borghesi,
A. Boscá,
F. Calle,
N. Canci,
G. Cavoto,
C. Chang,
W. Chung,
A. G. Cocco,
A. P. Colijn,
N. D'Ambrosio,
N. de Groot,
M. Faverzani,
A. Ferella,
E. Ferri,
L. Ficcadenti,
P. Garcia-Abia,
G. Garcia Gomez-Tejedor,
S. Gariazzo,
F. Gatti,
C. Gentile,
A. Giachero,
Y. Hochberg,
Y. Kahn
, et al. (31 additional authors not shown)
Abstract:
The PTOLEMY transverse drift filter is a new concept to enable precision analysis of the energy spectrum of electrons near the tritium beta-decay endpoint. This paper details the implementation and optimization methods for successful operation of the filter. We present the first demonstrator that produces the required magnetic field properties with an iron return-flux magnet. Two methods for the s…
▽ More
The PTOLEMY transverse drift filter is a new concept to enable precision analysis of the energy spectrum of electrons near the tritium beta-decay endpoint. This paper details the implementation and optimization methods for successful operation of the filter. We present the first demonstrator that produces the required magnetic field properties with an iron return-flux magnet. Two methods for the setting of filter electrode voltages are detailed. The challenges of low-energy electron transport in cases of low field are discussed, such as the growth of the cyclotron radius with decreasing magnetic field, which puts a ceiling on filter performance relative to fixed filter dimensions. Additionally, low pitch angle trajectories are dominated by motion parallel to the magnetic field lines and introduce non-adiabatic conditions and curvature drift. To minimize these effects and maximize electron acceptance into the filter, we present a three-potential-well design to simultaneously drain the parallel and transverse kinetic energies throughout the length of the filter. These optimizations are shown, in simulation, to achieve low-energy electron transport from a 1 T iron core (or 3 T superconducting) starting field with initial kinetic energy of 18.6 keV drained to <10 eV (<1 eV) in about 80 cm. This result for low field operation paves the way for the first demonstrator of the PTOLEMY spectrometer for measurement of electrons near the tritium endpoint to be constructed at the Gran Sasso National Laboratary (LNGS) in Italy.
△ Less
Submitted 24 January, 2022; v1 submitted 23 August, 2021;
originally announced August 2021.
-
Investigating global convective dynamos with mean-field models: full spectrum of turbulent effects required
Authors:
Jörn Warnecke,
Matthias Rheinhardt,
Mariangela Viviani,
Frederick Gent,
Simo Tuomisto,
Maarit J. Käpylä
Abstract:
The role of turbulent effects for dynamos in the Sun and stars continues to be debated. Mean-field (MF) theory provides a broadly used framework to connect these effects to fundamental magnetohydrodynamics. While inaccessible observationally, turbulent effects can be directly studied using global convective dynamo (GCD) simulations. We measure the turbulent effects in terms of turbulent transport…
▽ More
The role of turbulent effects for dynamos in the Sun and stars continues to be debated. Mean-field (MF) theory provides a broadly used framework to connect these effects to fundamental magnetohydrodynamics. While inaccessible observationally, turbulent effects can be directly studied using global convective dynamo (GCD) simulations. We measure the turbulent effects in terms of turbulent transport coefficients, based on the MF framework, from an exemplary GCD simulation using the test-field method. These coefficients are then used as an input into an MF model. We find a good agreement between the MF and GCD solutions, which validates our theoretical approach. This agreement requires all turbulent effects to be included, even those which have been regarded as unimportant so far. Our results suggest that simple dynamo models, as are commonly used in the solar and stellar community, relying on very few, precisely fine-tuned turbulent effects, may not be representative of the full dynamics of dynamos in global convective simulations and astronomical objects.
△ Less
Submitted 5 October, 2021; v1 submitted 17 May, 2021;
originally announced May 2021.
-
Hunting down the cause of solar magnetism
Authors:
M. Viviani,
A. Prabhu,
J. Warnecke,
L. Duarte,
J. Pekkilä,
M. Rheinhardt,
M. J. Käpylä
Abstract:
To understand solar and stellar dynamos combining local and global numerical modelling with long-term observations is a challenging task: even with state of the art computational methods and resources, the stellar parameter regime remains unattainable. Our goal is to relax some approximations, in order to simulate more realistic systems, and try to connect the results with theoretical predictions…
▽ More
To understand solar and stellar dynamos combining local and global numerical modelling with long-term observations is a challenging task: even with state of the art computational methods and resources, the stellar parameter regime remains unattainable. Our goal is to relax some approximations, in order to simulate more realistic systems, and try to connect the results with theoretical predictions and state-of-the-art observations. We present here the first test-field measurements from our higher-resolution runs with improved heat conduction description. They indicate significant changes in the profiles of the most crucial inductive effect related to solar and stellar dynamo mechanisms. Higher resolution runs, currently undertaken, will bring us into an even more turbulent regime, in which we will be able to study, for the first time, the interaction of small- and large-scale dynamos in a quantitative way.
△ Less
Submitted 4 February, 2021;
originally announced February 2021.
-
Physically motivated heat conduction treatment in simulations of solar-like stars: effects on dynamo transitions
Authors:
M. Viviani,
M. J. Käpylä
Abstract:
Context. Results from global magnetoconvection simulations of solar-like stars are at odds with observations in many respects: They show a surplus of energy in the kinetic power spectrum at large scales, anti-solar differential rotation profiles, with accelerated poles and a slow equator, for the solar rotation rate, and a transition from axi- to non-axisymmetric dynamos at a much lower rotation r…
▽ More
Context. Results from global magnetoconvection simulations of solar-like stars are at odds with observations in many respects: They show a surplus of energy in the kinetic power spectrum at large scales, anti-solar differential rotation profiles, with accelerated poles and a slow equator, for the solar rotation rate, and a transition from axi- to non-axisymmetric dynamos at a much lower rotation rate than what is observed. Even though the simulations reproduce the observed active longitudes in fast rotators, their motion in the rotational frame (the so-called azimuthal dynamo wave, ADW) is retrograde, in contrast to the prevalent prograde motion in observations. Aims. We study the effect of a more realistic treatment of heat conductivity in alleviating the discrepancies between observations and simulations. Methods. We use physically-motivated heat conduction, by applying Kramers opacity law, on a semi-global spherical setup describing convective envelopes of solar-like stars, instead of a prescribed heat conduction profile from mixing-length arguments. Results. We find that some aspects of the results now better correspond to observations: The axi- to non-axisymmetric transition point is shifted towards higher rotation rates. We also find a change in the propagation direction of ADWs so that also prograde waves are now found. The transition from anti-solar to solar-like rotation profile, however, is also shifted towards higher rotation rates, leaving the models into an even more unrealistic regime. Conclusions. Although a Kramers-based heat conduction does not help in reproducing the solar rotation profile, it does help in the faster rotation regime, where the dynamo solutions now match better with observations.
△ Less
Submitted 11 June, 2020; v1 submitted 8 June, 2020;
originally announced June 2020.
-
Shapes of stellar activity cycles
Authors:
T. Willamo,
T. Hackman,
J. J. Lehtinen,
M. J. Käpylä,
N. Olspert,
M. Viviani,
J. Warnecke
Abstract:
Magnetic activity cycles are an important phenomenon in both the Sun and other stars. The shape of the solar cycle is commonly characterised by a fast rise and a slower decline, but not much attention has been paid to the shape of cycles in other stars. Our aim is to study whether the asymmetric shape of the solar cycle is common in other stars as well, and compare the cycle asymmetry to other ste…
▽ More
Magnetic activity cycles are an important phenomenon in both the Sun and other stars. The shape of the solar cycle is commonly characterised by a fast rise and a slower decline, but not much attention has been paid to the shape of cycles in other stars. Our aim is to study whether the asymmetric shape of the solar cycle is common in other stars as well, and compare the cycle asymmetry to other stellar parameters. We also study the differences in the shape of the solar cycle, depending on what activity indicator is used. The observations are also compared to simulated activity cycles. We use the chromospheric Ca II H&K data from the Mount Wilson Observatory HK Project. From this data set we identify 47 individual cycles from 18 stars. We use the statistical skewness of a cycle as a measure of its asymmetry, and compare this to other stellar parameters. A similar analysis has been done to magnetic cycles extracted from direct numerical magnetohydrodynamic simulations of solar-type convection zones. The shape of the solar cycle (fast rise and slower decline) is common in other stars as well, although the Sun has particularly asymmetric cycles. Cycle-to-cycle variations are large, but the average shape of a cycle is still fairly well represented by a sinusoid. We find only slight correlations between the cycle asymmetry and other stellar parameters. There are large differences in the shape of the solar cycle, depending on what activity indicator is used. In the simulated cycles, there is a difference in the symmetry of global simulations covering the full longitudinal range, hence capable of exciting non-axisymmetric dynamo modes, versus wedge simulations covering a partial extent in longitude, where only axisymmetric modes are possible. The former produce preferentially positive skewness, while the latter a negative one.
△ Less
Submitted 29 April, 2020; v1 submitted 11 February, 2020;
originally announced February 2020.
-
Neutrino physics with the PTOLEMY project: active neutrino properties and the light sterile case
Authors:
PTOLEMY collaboration,
M. G. Betti,
M. Biasotti,
A. Boscá,
F. Calle,
G. Cavoto,
C. Chang,
A. G. Cocco,
A. P. Colijn,
J. Conrad,
N. D'Ambrosio,
N. De Groot,
P. F. de Salas,
M. Faverzani,
A. Ferella,
E. Ferri,
P. Garcia-Abia,
I. García-Cortés,
G. Garcia Gomez-Tejedor,
S. Gariazzo,
F. Gatti,
C. Gentile,
A. Giachero,
J. E. Gudmundsson,
Y. Hochberg
, et al. (32 additional authors not shown)
Abstract:
The PTOLEMY project aims to develop a scalable design for a Cosmic Neutrino Background (CNB) detector, the first of its kind and the only one conceived that can look directly at the image of the Universe encoded in neutrino background produced in the first second after the Big Bang. The scope of the work for the next three years is to complete the conceptual design of this detector and to validate…
▽ More
The PTOLEMY project aims to develop a scalable design for a Cosmic Neutrino Background (CNB) detector, the first of its kind and the only one conceived that can look directly at the image of the Universe encoded in neutrino background produced in the first second after the Big Bang. The scope of the work for the next three years is to complete the conceptual design of this detector and to validate with direct measurements that the non-neutrino backgrounds are below the expected cosmological signal. In this paper we discuss in details the theoretical aspects of the experiment and its physics goals. In particular, we mainly address three issues. First we discuss the sensitivity of PTOLEMY to the standard neutrino mass scale. We then study the perspectives of the experiment to detect the CNB via neutrino capture on tritium as a function of the neutrino mass scale and the energy resolution of the apparatus. Finally, we consider an extra sterile neutrino with mass in the eV range, coupled to the active states via oscillations, which has been advocated in view of neutrino oscillation anomalies. This extra state would contribute to the tritium decay spectrum, and its properties, mass and mixing angle, could be studied by analyzing the features in the beta decay electron spectrum.
△ Less
Submitted 30 September, 2019; v1 submitted 14 February, 2019;
originally announced February 2019.
-
Stellar dynamos in the transition regime: multiple dynamo modes and anti-solar differential rotation
Authors:
M. Viviani,
M. J. Käpylä,
J. Warnecke,
P. J. Käpylä,
M. Rheinhardt
Abstract:
Global and semi-global convective dynamo simulations of solar-like stars are known to show a transition from an anti-solar (fast poles, slow equator) to solar-like (fast equator, slow poles) differential rotation (DR) for increasing rotation rate. The dynamo solutions in the latter regime can exhibit regular cyclic modes, whereas in the former one, only stationary or temporally irregular solutions…
▽ More
Global and semi-global convective dynamo simulations of solar-like stars are known to show a transition from an anti-solar (fast poles, slow equator) to solar-like (fast equator, slow poles) differential rotation (DR) for increasing rotation rate. The dynamo solutions in the latter regime can exhibit regular cyclic modes, whereas in the former one, only stationary or temporally irregular solutions have been obtained so far. In this paper we present a semi-global dynamo simulation in the transition region, exhibiting two coexisting dynamo modes, a cyclic and a stationary one, both being dynamically significant. We seek to understand how such a dynamo is driven by analyzing the large-scale flow properties (DR and meridional circulation) together with the turbulent transport coefficients obtained with the test-field method. Neither an $αΩ$ dynamo wave nor an advection-dominated dynamo are able to explain the cycle period and the propagation direction of the mean magnetic field. Furthermore, we find that the $α$ effect is comparable or even larger than the $Ω$ effect in generating the toroidal magnetic field, and therefore, the dynamo seems to be $α^2Ω$ or $α^2$ type. We further find that the effective large-scale flows are significantly altered by turbulent pumping.
△ Less
Submitted 22 November, 2019; v1 submitted 11 February, 2019;
originally announced February 2019.
-
Effects of a subadiabatic layer on convection and dynamos in spherical wedge simulations
Authors:
Petri J. Käpylä,
Mariangela Viviani,
Maarit J. Käpylä,
Axel Brandenburg,
Federico Spada
Abstract:
We consider the effect of a subadiabatic layer at the base of the convection zone on convection itself and the associated large-scale dynamos in spherical wedge geometry. We use a heat conduction prescription based on the Kramers opacity law which allows the depth of the convection zone to dynamically adapt to changes in the physical characteristics such as rotation rate and magnetic fields. We fi…
▽ More
We consider the effect of a subadiabatic layer at the base of the convection zone on convection itself and the associated large-scale dynamos in spherical wedge geometry. We use a heat conduction prescription based on the Kramers opacity law which allows the depth of the convection zone to dynamically adapt to changes in the physical characteristics such as rotation rate and magnetic fields. We find that the convective heat transport is strongly concentrated toward the equatorial and polar regions in the cases without a substantial radiative layer below the convection zone. The presence of a stable layer below the convection zone significantly reduces the anisotropy of radial enthalpy transport. Furthermore, the dynamo solutions are sensitive to subtle changes in the convection zone structure. We find that the kinetic helicity changes sign in the deeper parts of the convection zone at high latitudes in all runs. This region expands progressively toward the equator in runs with a thicker stably stratified layer.
△ Less
Submitted 2 January, 2019; v1 submitted 15 March, 2018;
originally announced March 2018.
-
Transition from axi- to nonaxisymmetric dynamo modes in spherical convection models of solar-like stars
Authors:
Mariangela Viviani,
Jörn Warnecke,
Maarit J. Käpylä,
Petri J. Käpylä,
Nigul Olspert,
Elizabeth M. Cole-Kodikara,
Jyri J. Lehtinen,
Axel Brandenburg
Abstract:
We seek to understand the transition from nearly axisymmetric configurations at solar rotation rates to nonaxisymmetric configurations for rapid rotation using 3D numerical simulations of turbulent convection and considering rotation rates between 1 and 30 times the solar value. We find a transition from axi- to nonaxisymmetric solutions at around 1.8 times the solar rotation rate. This transition…
▽ More
We seek to understand the transition from nearly axisymmetric configurations at solar rotation rates to nonaxisymmetric configurations for rapid rotation using 3D numerical simulations of turbulent convection and considering rotation rates between 1 and 30 times the solar value. We find a transition from axi- to nonaxisymmetric solutions at around 1.8 times the solar rotation rate. This transition coincides with a change in the rotation profile from antisolar- to solar-like differential rotation with a faster equator and slow poles. In the solar-like rotation regime, the field configuration consists of an axisymmetric oscillatory field accompanied by an m=1 azimuthal mode (two active longitudes), which also shows temporal variability. At slow (rapid) rotation, the axisymmetric (nonaxisymmetric) mode dominates. The axisymmetric mode produces latitudinal dynamo waves with polarity reversals, while the nonaxisymmetric mode often exhibits a drift in the rotating reference frame and the strength of the active longitudes changes cyclically over time between the different hemispheres. Most of the obtained dynamo solutions exhibit cyclic variability either caused by latitudinal or azimuthal dynamo waves. In an activity-period diagram, the cycle lengths normalized by the rotation period form two different populations as a function of rotation rate or magnetic activity level. The slowly rotating axisymmetric population lies close to what is called the inactive branch in observations, while the rapidly rotating models are close to the superactive branch with a declining cycle to rotation frequency ratio with increasing rotation rate. We can successfully reproduce the transition from axi- to nonaxisymmetric dynamo solutions for high rotation rates, but high-resolution simulations are required to limit the effect of rotational quenching of convection at rotation rates above 20 times the solar value.
△ Less
Submitted 7 May, 2018; v1 submitted 27 October, 2017;
originally announced October 2017.
-
Implication of the proton-deuteron radiative capture for Big Bang Nucleosynthesis
Authors:
L. E. Marcucci,
G. Mangano,
A. Kievsky,
M. Viviani
Abstract:
The astrophysical $S$-factor for the radiative capture $d(p,γ)^3$He in the energy-range of interest for Big Bang Nucleosynthesis (BBN) is calculated using an {\it ab-initio} approach. The nuclear Hamiltonian retains both two- and three-nucleon interactions - the Argonne $v_{18}$ and the Urbana IX, respectively. Both one- and many-body contributions to the nuclear current operator are included. The…
▽ More
The astrophysical $S$-factor for the radiative capture $d(p,γ)^3$He in the energy-range of interest for Big Bang Nucleosynthesis (BBN) is calculated using an {\it ab-initio} approach. The nuclear Hamiltonian retains both two- and three-nucleon interactions - the Argonne $v_{18}$ and the Urbana IX, respectively. Both one- and many-body contributions to the nuclear current operator are included. The former retain for the first time, besides the $1/m$ leading order contribution ($m$ is the nucleon mass), also the next-to-leading order term, proportional to $1/m^3$. The many-body currents are constructed in order to satisfy the current conservation relation with the adopted Hamiltonian model. The hyperspherical harmonics technique is applied to solve the $A=3$ bound and scattering states. A particular attention is used in this second case in order to obtain, in the energy range of BBN, an uncertainty on the astrophysical $S$-factor of the order or below $\sim$1 %. Then, in this energy range, the $S$-factor is found to be $\sim$10 % larger than the currently adopted values.Part of this increase (1-3 %) is due to the $1/m^3$ one-body operator, while the remaining is due to the new more accurate scattering wave functions. We have studied the implication of this new determination for the $d(p,γ)^3$He $S$-factor on deuterium primordial abundance. We find that the predicted theoretical value for $^2$H/H is in excellent agreement with its experimental determination, using the most recent determination of baryon density of Planck experiment, and with a standard number of relativistic degrees of freedom $N_{\rm eff}=3.046$ during primordial nucleosynthesis.
△ Less
Submitted 27 October, 2015;
originally announced October 2015.
-
Parameter-free effective field theory calculation for the solar proton-fusion and hep processes
Authors:
T. -S. Park,
L. E. Marcucci,
R. Schiavilla,
M. Viviani,
A. Kievsky,
S. Rosati,
K. Kubodera,
D. -P. Min,
M. Rho
Abstract:
Spurred by the recent complete determination of the weak currents in two-nucleon systems up to ${\cal O}(Q^3)$ in heavy-baryon chiral perturbation theory, we carry out a parameter-free calculation of the threshold $S$-factors for the solar $pp$ (proton-fusion) and $hep$ processes in an effective field theory that {\it combines} the merits of the standard nuclear physics method and systematic chi…
▽ More
Spurred by the recent complete determination of the weak currents in two-nucleon systems up to ${\cal O}(Q^3)$ in heavy-baryon chiral perturbation theory, we carry out a parameter-free calculation of the threshold $S$-factors for the solar $pp$ (proton-fusion) and $hep$ processes in an effective field theory that {\it combines} the merits of the standard nuclear physics method and systematic chiral expansion. The power of the EFT adopted here is that one can correlate in a unified formalism the weak-current matrix elements of two-, three- and four-nucleon systems. Using the tritium $β$-decay rate as an input to fix the only unknown parameter in the theory, we can evaluate the threshold $S$ factors with drastically improved precision; the results are $S_{pp}(0) = 3.94\times(1 \pm 0.004) \times 10^{-25} {MeV-b}$ and $S_{hep}(0) = (8.6\pm 1.3)\times 10^{-20} {keV-b}$. The dependence of the calculated $S$-factors on the momentum cutoff parameter $Λ$ has been examined for a physically reasonable range of $Λ$. This dependence is found to be extremely small for the $pp$ process, and to be within acceptable levels for the $hep$ process, substantiating the consistency of our calculational scheme.
△ Less
Submitted 26 August, 2002;
originally announced August 2002.
-
The Solar hep Process in Effective Field Theory
Authors:
T. -S. Park,
L. E. Marcucci,
R. Schiavilla,
M. Viviani,
A. Kievsky,
S. Rosati,
K. Kubodera,
D. -P. Min,
M. Rho
Abstract:
Using effective field theory, we calculate the S-factor for the hep process in a totally parameter-free formulation. The transition operators are organized according to chiral counting, and their matrix elements are evaluated using the realistic nuclear wave functions obtained in the correlated-hyperspherical-harmonics method. Terms of up to next-to-next-to-next-to-leading order in heavy-baryon…
▽ More
Using effective field theory, we calculate the S-factor for the hep process in a totally parameter-free formulation. The transition operators are organized according to chiral counting, and their matrix elements are evaluated using the realistic nuclear wave functions obtained in the correlated-hyperspherical-harmonics method. Terms of up to next-to-next-to-next-to-leading order in heavy-baryon chiral perturbation theory are considered. Fixing the only parameter in the theory by fitting the tritium β-decay rate, we predict the hep S-factor with accuracy better than \sim 20 %.
△ Less
Submitted 3 August, 2001; v1 submitted 4 July, 2001;
originally announced July 2001.
-
Parameter-Free Calculation of the Solar Proton Fusion Rate in Effective Field Theory
Authors:
T. -S. Park,
L. E. Marcucci,
R. Schiavilla,
M. Viviani,
A. Kievsky,
S. Rosati,
K. Kubodera,
D. -P. Min,
M. Rho
Abstract:
Spurred by the recent complete determination of the weak currents in two-nucleon systems up to ${\cal O}(Q^3)$ in heavy-baryon chiral perturbation theory, we carry out a parameter-free calculation of the solar proton fusion rate in an effective field theory that combines the merits of the standard nuclear physics method and systematic chiral expansion. Using the tritium beta-decay rate as an inp…
▽ More
Spurred by the recent complete determination of the weak currents in two-nucleon systems up to ${\cal O}(Q^3)$ in heavy-baryon chiral perturbation theory, we carry out a parameter-free calculation of the solar proton fusion rate in an effective field theory that combines the merits of the standard nuclear physics method and systematic chiral expansion. Using the tritium beta-decay rate as an input to fix the only unknown parameter in the effective Lagrangian, we can evaluate with drastically improved precision the ratio of the two-body contribution to the well established one-body contribution; the ratio is determined to be (0.86\pm 0.05) %. This result is essentially independent of the cutoff parameter for a wide range of its variation (500 MeV \le Λ\le 800 MeV), a feature that substantiates the consistency of the calculation.
△ Less
Submitted 16 July, 2001; v1 submitted 12 June, 2001;
originally announced June 2001.
-
Weak proton capture on 3He
Authors:
L. E. Marcucci,
R. Schiavilla,
M. Viviani,
A. Kievsky,
S. Rosati,
J. F. Beacom
Abstract:
The astrophysical S-factor for the proton weak capture on 3He is calculated with correlated-hyperspherical-harmonics bound and continuum wave functions corresponding to realistic Hamiltonians consisting of the Argonne v14 or Argonne v18 two-nucleon and Urbana-VIII or Urbana-IX three-nucleon interactions. The nuclear weak charge and current operators have vector and axial-vector components, that…
▽ More
The astrophysical S-factor for the proton weak capture on 3He is calculated with correlated-hyperspherical-harmonics bound and continuum wave functions corresponding to realistic Hamiltonians consisting of the Argonne v14 or Argonne v18 two-nucleon and Urbana-VIII or Urbana-IX three-nucleon interactions. The nuclear weak charge and current operators have vector and axial-vector components, that include one- and many-body terms. All possible multipole transitions connecting any of the p 3He S- and P-wave channels to the 4He bound state are considered. The S-factor at a p 3He center-of-mass energy of 10 keV, close to the Gamow-peak energy, is predicted to be 10.1 10^{-20} keV b with the AV18/UIX Hamiltonian, a factor of about 4.5 larger than the value adopted in the standard solar model. The P-wave transitions are found to be important, contributing about 40 % of the calculated S-factor. The energy dependence is rather weak: the AV18/UIX zero-energy S-factor is 9.64 10^{-20} keV b, only 5 % smaller than the 10 keV result quoted above. The model dependence is also found to be weak: the zero-energy S-factor is calculated to be 10.2 10^{-20} keV b with the older AV14/UVIII model, only 6 % larger than the AV18/UIX result. Our best estimate for the S-factor at 10 keV is therefore (10.1 \pm 0.6) 10^{-20} keV b, when the theoretical uncertainty due to the model dependence is included. This value for the calculated S-factor is not as large as determined in fits to the Super-Kamiokande data in which the hep flux normalization is free. However, the precise calculation of the S-factor and the consequent absolute prediction for the hep neutrino flux will allow much greater discrimination among proposed solar neutrino oscillation solutions.
△ Less
Submitted 5 June, 2000;
originally announced June 2000.
-
Realistic Calculation of the hep Astrophysical Factor
Authors:
L. E. Marcucci,
R. Schiavilla,
M. Viviani,
A. Kievsky,
S. Rosati
Abstract:
The astrophysical factor for the proton weak capture on 3He is calculated with correlated-hyperspherical-harmonics bound and continuum wave functions corresponding to a realistic Hamiltonian consisting of the Argonne v18 two-nucleon and Urbana-IX three-nucleon interactions. The nuclear weak charge and current operators have vector and axial-vector components, that include one- and many-body term…
▽ More
The astrophysical factor for the proton weak capture on 3He is calculated with correlated-hyperspherical-harmonics bound and continuum wave functions corresponding to a realistic Hamiltonian consisting of the Argonne v18 two-nucleon and Urbana-IX three-nucleon interactions. The nuclear weak charge and current operators have vector and axial-vector components, that include one- and many-body terms. All possible multipole transitions connecting any of the p-3He S- and P-wave channels to the 4He bound state are considered. The S-factor at a p-3He center-of-mass energy of 10 keV, close to the Gamow-peak energy, is predicted to be 10.1 10^{-20} keV b, a factor of five larger than the standard-solar-model value. The P-wave transitions are found to be important, contributing about 40 % of the calculated S-factor.
△ Less
Submitted 29 March, 2000;
originally announced March 2000.
