Search

Current and future electron and neutrino scattering experiments will be greatly aided by a better understanding of the role played by short-range correlations in nuclei. Two-body physics, including nucleon-nucleon correlations and two-body electroweak currents, is required to explain the body of experimental data for both static and dynamical nuclear properties. In this work, we focus on examining nucleon-nucleon correlations from a chiral effective field theory perspective and provide a comprehensive set of new variational Monte Carlo calculations of one- and two-body densities and momentum distributions based on the Norfolk many-body nuclear Hamiltonians for A, Comment: 12 pages, 11 figures

We present the equation of state of infinite neutron matter as obtained from highly-realistic Hamiltonians that include nucleon-nucleon and three-nucleon coordinate-space potentials. We benchmark three independent many-body methods: Brueckner-Bethe-Goldstone (BBG), Fermi hypernetted chain/single-operator chain (FHNC/SOC), and auxiliary-field diffusion Monte Carlo (AFDMC). We find them to provide similar equations of state when the Argonne $v_{18}$ and the Argonne $v_{6}^\prime$ nucleon-nucleon potentials are used in combination with the Urbana IX three-body force. Only at densities larger than about 1.5 the nuclear saturation density ($\rho_0 = 0.16\,\rm{fm}^{-3}$) the FHNC/SOC energies are appreciably lower than the other two approaches. The AFDMC calculations carried out with all of the Norfolk potentials fitted to reproduce the experimental trinucleon ground-state energies and $nd$ doublet scattering length yield unphysically bound neutron matter, associated with the formation of neutron droplets. Including tritium $\beta$-decay in the fitting procedure, as in the second family of Norfolk potentials, mitigates but does not completely resolve this problem. An excellent agreement between the BBG and AFDMC results is found for the subset of Norfolk interactions that do not make neutron-matter collapse, while the FHNC/SOC equations of state are moderately softer., Comment: 15 pages, 3 figures. arXiv admin note: text overlap with arXiv:1908.04426

While mean-field approximations, such as the nuclear shell model, provide a good description of many bulk nuclear properties, they fail to capture the important effects of nucleon–nucleon correlations such as the short-distance and high-momentum components of the nuclear many-body wave function1. Here, we study these components using the effective pair-based generalized contact formalism2,3 and ab initio quantum Monte Carlo calculations of nuclei from deuteron to 40Ca (refs. 4–6). We observe a universal factorization of the many-body nuclear wave function at short distance into a strongly interacting pair and a weakly interacting residual system. The residual system distribution is consistent with that of an uncorrelated system, showing that short-distance correlation effects are predominantly embedded in two-body correlations. Spin- and isospin-dependent ‘nuclear contact terms’ are extracted in both coordinate and momentum space for different realistic nuclear potentials. The contact coefficient ratio between two different nuclei shows very little dependence on the nuclear interaction model. These findings thus allow extending the application of mean-field approximations to short-range correlated pair formation by showing that the relative abundance of short-range pairs in the nucleus is a long-range (that is, mean field) quantity that is insensitive to the short-distance nature of the nuclear force. Effects of nucleon–nucleon correlations are studied with the generalized contact formalism and ab initio quantum Monte Carlo calculations. For nuclei from deuteron to 40Ca, the many-body nuclear wave function is shown to factorize at short distances.

Understanding the structure and reactions of nuclei from first principles has been a long-standing goal of nuclear physics. In this respect, few- and many-body systems provide a unique laboratory for studying nuclear interactions. In the past couple of decades, the modeling of nuclear interactions has progressed significantly owing, in particular, to the development of chiral effective field theory ($\chi$EFT), a low-energy effective representation of quantum chromodynamics (QCD). Within $\chi$EFT, many studies have dealt with the construction of both two- and three-nucleon interactions. The aim of the present article is to provide a concise account of chiral interaction models that are local in configuration space, and to report on a selection of recent results for nuclear systems obtained with these interactions., Comment: 14 pages, 2 figures, Contribution to Special Issue of Few-Body Systems: Celebrating 30 years of Steven Weinberg's papers on Nuclear Forces from Chiral Lagrangians

Classes of two-nucleon ($2N$) contact interactions are developed in configuration space at leading order (LO), next-to-leading order (NLO), and next-to-next-to-next-to-leading order (N3LO) by fitting the experimental singlet $np$ scattering length and deuteron binding energy at LO, and $np$ and $pp$ scattering data in the laboratory-energy range 0--15 MeV at NLO and 0--25 MeV at N3LO. These interactions are regularized by including two Gaussian cutoffs, one for $T\,$=$\,0$ and the other for $T\,$=$\,1$ channels. The cutoffs are taken to vary in the ranges $R_0\,$=$(1.5$--2.3) fm and $R_1\,$=$(1.5$--3.0) fm. The 780 (1,100) data points up to 15 (25) MeV energy, primarily differential cross sections, are fitted by the NLO (N3LO) models with a $\chi^2$/datum about 1.7 or less (well below 1.5), when harder cutoff values are adopted. As a first application, we report results for the binding energies of nuclei with mass numbers $A\,$=$\,3$--6 and 16 obtained with selected LO and NLO $2N$ models both by themselves as well as in combination with a LO three-nucleon ($3N$) contact interaction. The latter is characterized by a single low-energy constant that is fixed to reproduce the experimental $^3$H binding energy. The inclusion of the $3N$ interaction largely removes the sensitivity to cutoff variations in the few-nucleon systems and leads to predictions for the $^3$He and $^4$He binding energies that cluster around 7.8 MeV and 30 MeV, respectively. However, in $^{16}$O this cutoff sensitivity remains rather strong. Finally, predictions at LO only are also reported for medium-mass nuclei with $A\,$=$\,40$, 48, and 90., Comment: 23 pages, 6 figures

We report Quantum Monte Carlo calculations of weak transitions in $A\leq 10$ nuclei, based on the Norfolk two- and three-nucleon chiral interactions, and associated one- and two-body axial currents. We find that the contribution from two-body currents is at the $2$ - $3\%$ level, with the exception of matrix elements entering the rates of $^8$Li, $^8$B, and $^8$He beta decays. These matrix elements are suppressed in impulse approximation based on the (leading order) Gamow Teller transition operator alone; two-body currents provide a $20$ - $30\%$ correction, which is, however, insufficient to bring theory in agreement with experimental data. For the other transitions, the agreement with the data is satisfactory, and the results exhibit a negligible to mild model dependence when different combinations of Norfolk interactions are utilized to construct the nuclear wave functions. We report a complete study of two-body weak transition densities which reveals the expected universal behavior of two-body currents at short distances throughout the range of $A\,$=$\,3$ to $A\,$=$\,10$ systems considered here., 14 pages, 9 figures; Corrected experimental values for $^8$Li beta decay in Table V to be consistent with Table VI and the text; Corrected an error in row 2 column 4 of Figure 6, conclusions unchanged; Corrected errors in column 3 rows 3 and 4 of Figure 7, results unchanged; Corrected a typo in Equation 19, results in tables unchanged; Corrected $^8$He ground state isospin from $T=1$ to $T=2$

The energy density is calculated in coordinate space for $^{12}$C, $^{40}$Ca, $^{48}$Ca, and $^{208}$Pb using a dispersive optical model constrained by all relevant data including the corresponding energy of the ground state. The energy density of $^{8}$Be is also calculated using the Green's function Monte-Carlo method employing the Argonne/Urbana two and three-body interactions. The nuclear interior minimally contributes to the total binding energy due to the 4$\pi r^2$ phase space factor. Thus, the volume contribution to the energy in the interior is not well constrained. The dispersive-optical-model energy densities are in good agreement with \textit{ab initio} self-consistent Green's function calculations of infinite nuclear matter restricted to treat only short-range and tensor correlations. These results call into question the degree to which the equation of state for nuclear matter is constrained by the empirical mass formula. In particular, the results in this paper indicate that saturated nuclear matter does not require the canonical value of 16 MeV binding per particle but only about 13-14 MeV when the interior of $^{208}$Pb is considered., Comment: 15 pages, 13 figures

A major goal of nuclear theory is to explain the wealth of data and peculiarities exhibited by nuclear systems in a fully microscopic way. In such an approach, which we refer to as the basic model of nuclear theory, the nucleons interact with each other via many-body (primarily, two- and three-body) effective interactions, and with external electroweak probes via effective current operators. These effective interactions and currents are the main inputs to ab-initio methods that are aimed at solving the many-body Schrodinger equation associated with the nuclear system under consideration. In this talk, I will review recent progress in Quantum Monte Carlo calculations of low-lying spectra and electroweak properties of light nuclei as well as nucleonic matter equation of state. Emphasis will be on calculations based on chiral effective field theory approach.

To obtain an understanding of the structure and reactions of nuclear systems from first principles has been a long-standing goal of nuclear physics. In this respect, few- and many-body systems provide a unique laboratory for studying nuclear interactions. During the past decades, the development of accurate representations of the nuclear force has undergone substantial progress. Particular emphasis has been devoted to chiral effective field theory (EFT), a low-energy effective representation of quantum chromodynamics (QCD). Within chiral EFT, many studies have been carried out dealing with the construction of both the nucleon-nucleon ($N\!N$) and three-nucleon ($3N$) interactions. The aim of the present article is to give a detailed overview of the chiral interaction models that are local in configuration space, and show recent results for nuclear systems obtained by employing these local chiral forces.

We report benchmark calculations of the energy per particle of pure neutron matter as a function of the baryon density using three independent many-body methods: Brueckner-Bethe-Goldstone, Fermi hypernetted chain/single-operator chain, and auxiliary-field diffusion Monte Carlo. Significant technical improvements are implemented in the latter two methods. The calculations are made for two distinct families of realistic coordinate-space nucleon-nucleon potentials fit to scattering data, including the standard Argonne $v_{18}$ interaction and two of its simplified versions, and four of the new Norfolk $\Delta$-full chiral effective field theory potentials. The results up to twice nuclear matter saturation density show some divergence among the methods, but improved agreement compared to earlier work. We find that the potentials fit to higher-energy nucleon-nucleon scattering data exhibit a much smaller spread of energies., Comment: 16 pages, 7 figures

We present benchmark calculations of Zemach moments and radii of $^{2,3}\mathrm{H}$ and $^{3,4}\mathrm{He}$ using various few-body methods. Zemach moments are required to interpret muonic atom data measured by the CREMA collaboration at the Paul Scherrer Institute. Conversely, radii extracted from spectroscopic measurements can be compared with ab initio computations, posing stringent constraints on the nuclear model. For a given few-body method, different numerical procedures can be applied to compute these quantities. A detailed analysis of the numerical uncertainties entering the total theoretical error is presented. Uncertainties from the few-body method and the calculational procedure are found to be smaller than the dependencies on the dynamical modeling and the single nucleon inputs, which are found to be $\ensuremath{\lesssim}2%$. When relativistic corrections and two-body currents are accounted for, the calculated moments and radii are in very good agreement with the available experimental data.

The process at the heart of neutrinoless double-$\ensuremath{\beta}$ decay, $nn\ensuremath{\rightarrow}pp\phantom{\rule{0.16em}{0ex}}{e}^{\ensuremath{-}}{e}^{\ensuremath{-}}$ induced by a light Majorana neutrino, is investigated in pionless and chiral effective field theory. We show in various regularization schemes the need to introduce a short-range lepton-number-violating operator at leading order, confirming earlier findings. We demonstrate that such a short-range operator is only needed in spin-singlet $S$-wave transitions, while leading-order transitions involving higher partial waves depend solely on long-range currents. Calculations are extended to include next-to-leading-order corrections in perturbation theory, where to this order no additional undetermined parameters appear. We establish a connection based on chiral symmetry between neutrinoless double-$\ensuremath{\beta}$ decay and nuclear charge-independence breaking induced by electromagnetism. Data on the latter confirm the need for a leading-order short-range operator but do not allow for a full determination of the corresponding lepton-number-violating coupling. Using a crude estimate of this coupling, we perform ab initio calculations of the matrix elements for neutrinoless double-$\ensuremath{\beta}$ decay for $^{6}\mathrm{He}$ and $^{12}\mathrm{Be}$. We speculate on the phenomenological impact of the leading short-range operator on the basis of these results.

The magnetic form factors of $^2$H, $^3$H, and $^3$He, deuteron photodisintegration cross sections at low energies, and deuteron threshold electrodisintegration cross sections at backward angles in a wide range of momentum transfers, are calculated with the chiral two-nucleon (and three-nucleon) interactions including $\Delta$ intermediate states that have recently been constructed in configuration space. The $A\,$=$\,$3 wave functions are obtained from hyperspherical-harmonics solutions of the Schr\"odinger equation. The electromagnetic current includes one- and two-body terms, the latter induced by one- and two-pion exchange (OPE and TPE, respectively) mechanisms and contact interactions. The contributions associated with $\Delta$ intermediate states are only retained at the OPE level, and are neglected in TPE loop (tree-level) corrections to two-body (three-body) current operators. Expressions for these currents are derived and regularized in configuration space for consistency with the interactions. The low-energy constants that enter the contact few-nucleon systems. The predicted form factors and deuteron electrodisintegration cross section are in excellent agreement with experiment for momentum transfers up to 2--3 fm$^{-1}$. However, the experimental values for the deuteron photodisintegration cross section are consistently underestimated by theory, unless use is made of the Siegert form of the electric dipole transition operator. A complete analysis of the results is provided, including the clarification of the origin of the aforementioned discrepancy., Comment: 24 pages, 13 figures

A major goal of nuclear theory is to explain the spectra and stability of nuclei in terms of effective many-body interactions amongst the nucleus' constituents-the nucleons, i.e., protons and neutrons. Such an approach, referred to below as the basic model of nuclear theory, is formulated in terms of point-like nucleons, which emerge as effective degrees of freedom, at sufficiently low energy, as a result of a decimation process, starting from the fundamental quarks and gluons, described by Quantum Chromodynamics (QCD). A systematic way to account for the constraints imposed by the symmetries of QCD, in particular chiral symmetry, is provided by chiral effective field theory, in the framework of a low-energy expansion. Here we show, in quantum Monte Carlo calculations accurate to $\leq\!2\%$ of the binding energy, that two- and three-body chiral interactions fitted {\sl only} to bound- and scattering-state observables in, respectively, the two- and three-nucleon sectors, lead to predictions for the energy levels and level ordering of nuclei in the mass range $A\,$=$\,$4-12 in very satisfactory agreement with experimental data. Our findings provide strong support for the fundamental assumptions of the basic model, and pave the way to its systematic application to the electroweak structure and response of these systems as well as to more complex nuclei., Comment: 21 pages, 3 figures

We review our recent work on the derivation of the nuclear electromagnetic charge and current operators in chiral perturbation theory, based on time-ordered perturbation theory. We then discuss the strategies for fixing the relevant low energy constants, and compare the resulting predictions for the electric and magnetic form factors of the deuteron and trinucleons with experimental data, using as input accurate nuclear wave functions derived with realistic potentials.

We present fully local versions of the minimally nonlocal nucleon-nucleon potentials constructed in a previous paper [Piarulli et al., Phys. Rev. C 91, 024003 (2015)], and use them in hypersperical harmonics and quantum Monte Carlo calculations of ground and excited states of $^{3}\mathrm{H}$, $^{3}\mathrm{He}$, $^{4}\mathrm{He}$, $^{6}\mathrm{He}$, and $^{6}\mathrm{Li}$ nuclei. The long-range part of these local potentials includes one- and two-pion exchange contributions without and with $\mathrm{\ensuremath{\Delta}}$ isobars in the intermediate states up to order ${Q}^{3}$ ($Q$ denotes generically the low momentum scale) in the chiral expansion, while the short-range part consists of contact interactions up to order ${Q}^{4}$. The low-energy constants multiplying these contact interactions are fitted to the 2013 Granada database in two different ranges of laboratory energies, either 0--125 MeV or 0--200 MeV, and to the deuteron binding energy and $nn$ singlet scattering length. Fits to these data are performed for three models characterized by long- and short-range cutoffs, ${R}_{\mathrm{L}}$ and ${R}_{\mathrm{S}}$, respectively, ranging from $({R}_{\mathrm{L}},{R}_{\mathrm{S}})=(1.2,0.8)$ fm down to $(0.8,0.6)$ fm. The long-range (short-range) cutoff regularizes the one- and two-pion exchange (contact) part of the potential.

The objectives of the present work are twofold. The first is to address and resolve some of the differences present in independent, chiral-effective-field-theory (\chiEFT) derivations up to one loop, recently appeared in the literature, of the nuclear charge and current operators. The second objective is to provide a complete set of \chiEFT predictions for the structure functions and tensor polarization of the deuteron, for the charge and magnetic form factors of 3He and 3H, and for the charge and magnetic radii of these few-nucleon systems. The calculations use wave functions derived from high-order chiral two- and three-nucleon potentials and Monte Carlo methods to evaluate the relevant matrix elements. Predictions based on conventional potentials in combination with \chiEFT charge and current operators are also presented. There is excellent agreement between theory and experiment for all these observables for momentum transfers up to q< 2.0-2.5 (1/fm); for a subset of them, this agreement extends to momentum transfers as high as q~5-6 (1/fm). A complete analysis of the results is provided., 34 pages, Revtex

The muon-capture reactions 2H(μ-,νμ)nn and 3He(μ-,νμ)3H are studied with conventional or chiral realistic potentials and consistent weak currents. The initial and final A=2 and A=3 nuclear wave functions are obtained from the Argonne v18 or chiral next-to-next-to-next-to leading order (N3LO) two-nucleon potential, in combination with, respectively, the Urbana IX or chiral next-to-next-to leading order (N2LO) three-nucleon potential in the case of A=3. The weak current consists of polar- and axial-vector components. The former are related to the isovector piece of the electromagnetic current via the conserved-vector-current hypothesis. These and the axial currents are derived either in a meson-exchange or in a chiral effective field theory (χEFT) framework. There is one parameter (either the N-to-Δ axial coupling constant in the meson-exchange model, or the strength of a contact term in the χEFT model) that is fixed by reproducing the Gamow-Teller matrix element in tritium β decay. The model dependence relative to the adopted interactions and currents (and cutoff sensitivity in the χEFT currents) is weak, resulting in total rates of 392.0±2.3 s-1 for A=2, and 1484±13 s-1 for A=3, where the spread accounts for this model dependence.

We construct a coordinate-space chiral potential, including $\Delta$-isobar intermediate states in its two-pion-exchange component. The contact interactions entering at next-to-leading and next-to-next-to-next-to-leading orders ($Q^2$ and $Q^4$, respectively, $Q$ denoting generically the low momentum scale) are rearranged by Fierz transformations to yield terms at most quadratic in the relative momentum operator of the two nucleons. The low-energy constants multiplying these contact interactions are fitted to the 2013 Granada database, consisting of 2309 $pp$ and 2982 $np$ data (including, respectively, 148 and 218 normalizations) in the laboratory-energy range 0--300 MeV. For the total 5291 $pp$ and $np$ data in this range, we obtain a $\chi^2$/datum of roughly 1.3 for a set of three models characterized by long- and short-range cutoffs, $R_{\rm L}$ and $R_{\rm S}$ respectively, ranging from $(R_{\rm L},R_{\rm S})=(1.2,0.8)$ fm down to $(0.8,0.6)$ fm. The long-range (short-range) cutoff regularizes the one- and two-pion exchange (contact) part of the potential., Comment: 32 pages, 19 figures, accepted for publication as a Regular Article in Physical Review C

Background: The $A=2$ and 3 form factors are among the observables of choice for testing models of nuclear interactions and associated electromagnetic charge and current operators. Here we investigate the validity of the chiral-effective-field-theory ($\ensuremath{\chi}$EFT) approach to describe the strong-interaction dynamics in these few-nucleon systems and their response to electromagnetic probes.Purpose: The objectives of the present work are twofold. The first is to address and resolve some of the differences present in independent, $\ensuremath{\chi}$EFT derivations up to one loop, recently appearing in the literature, of the nuclear charge and current operators. The second objective is to provide a complete set of $\ensuremath{\chi}$EFT and hybrid predictions for the structure functions and tensor polarization of the deuteron, for the charge and magnetic form factors of ${}^{3}$He and ${}^{3}$H, and for the charge and magnetic radii of these few-nucleon systems.Methods: The calculations use wave functions derived from either chiral or conventional two- and three-nucleon potentials and Monte Carlo methods to evaluate the relevant matrix elements.Results: In reference to the two objectives mentioned earlier, we find that (i) there are no differences between the $\ensuremath{\chi}$EFT magnetic dipole operator at one loop derived in our formalism and that obtained by K\"olling et al. [Phys. Rev. C 80, 045502 (2009)] with the unitary transformation method and (ii) there is excellent agreement between theory and experiment for the static properties and elastic form factors of these $A=2$ and 3 nuclei up to momentum transfers $q\ensuremath{\lesssim}2.0$--2.5 fm${}^{\ensuremath{-}1}$. A complete analysis of the results is provided.Conclusions: Nuclear $\ensuremath{\chi}$EFT provides a very satisfactory description of the isoscalar and isovector charge and magnetic structure of the $A=2$ and 3 nuclei at low momentum transfers $q\ensuremath{\lesssim}3{m}_{\ensuremath{\pi}}$. In particular, contributions from two-body charge and current operators are crucial for bringing theory into close agreement with experiment. At higher $q$ values the present $\ensuremath{\chi}$EFT predictions are similar to those obtained in the hybrid approach, as well as in older studies based on the conventional meson-exchange picture, and fail to reproduce the measured $A=2$ and 3 form factors in the diffraction region.

This work investigates the muon capture reactions 2H(\mu^-,\nu_\mu)nn and 3He(\mu^-,\nu_\mu)3H and the contribution to their total capture rates arising from the axial two-body currents obtained imposing the partially-conserved-axial-current (PCAC) hypothesis. The initial and final A=2 and 3 nuclear wave functions are obtained from the Argonne v_{18} two-nucleon potential, in combination with the Urbana IX three-nucleon potential in the case of A=3. The weak current consists of vector and axial components derived in chiral effective field theory. The low-energy constant entering the vector (axial) component is determined by reproducting the isovector combination of the trinucleon magnetic moment (Gamow-Teller matrix element of tritium beta-decay). The total capture rates are 393.1(8) s^{-1} for A=2 and 1488(9) s^{-1} for A=3, where the uncertainties arise from the adopted fitting procedure., Comment: 6 pages, submitted to Few-Body Syst

The aim of this study was to verify the existence of synergistic antibacterial effect between four essential oils (Aniba rosaeodora, Melaleuca alternifolia, Origanum vulgare, and Pelargonium graveolens) individually combined with the antibacterial drug Gentamicin. We investigated the effectiveness in vitro of the association of essential oil/Gentamicin, against fifteen different strains of Gram positive and Gram negative bacteria. The antibacterial effects of these oils in combination with Gentamicin were evaluated by using the MHB microdilution method, while gas chromatography (GC) and GC/Mass spectrometry were used to analyze the chemical composition of the oils. A synergistic interaction was observed against all tested strains with the associations between the essential oils Aniba rosaeodora/Gentamicin and Pelargonium graveolens/Gentamicin. In particular a very strong synergistic interaction was observed against Acinetobacter baumannii ATCC 19606 (FIC index = 0.11). In contrast, the essential oils Origanum vulgare and Melaleuca alternifolia in association with Gentamicin were less effective on bacterial species growth. In vitro interaction can improve the antimicrobial effectiveness of the Gentamicin and may contribute to reduce its dose correlated to side effects.

Experimental form factors of the hydrogen and helium isotopes, extracted from an up-to-date global analysis of cross sections and polarization observables measured in elastic electron scattering from these systems, are compared to predictions obtained in three different theoretical approaches: the first is based on realistic interactions and currents, including relativistic corrections (labeled as the conventional approach); the second relies on a chiral effective field theory description of the strong and electromagnetic interactions in nuclei (labeled $\chi$EFT); the third utilizes a fully relativistic treatment of nuclear dynamics as implemented in the covariant spectator theory (labeled CST). For momentum transfers below $Q \lesssim 5$ fm$^{-1}$ there is satisfactory agreement between experimental data and theoretical results in all three approaches. However, at $Q \gtrsim 5$ fm$^{-1}$, particularly in the case of the deuteron, a relativistic treatment of the dynamics, as is done in the CST, is necessary. The experimental data on the deuteron $A$ structure function extend to $Q \simeq 12$ fm$^{-1}$, and the close agreement between these data and the CST results suggests that, even in this extreme kinematical regime, there is no evidence for new effects coming from quark and gluon degrees of freedom at short distances., Comment: 66 pages, 23 figures, topical review submitted to J. Phys. G: Nucl. Part. Phys

A coordinate-space nucleon-nucleon potential is constructed in chiral effective field theory (χ EFT) retaining pions, nucleons and Δ-isobars as explicit degrees of freedom. The calculation of the potential is carried out by including one- and two-pionexchange contributions up to next-to-next-to-leading order (N2LO) and contact interactions up to next-to-next-to-next-to-leading order (N3LO). The low-energy constants multiplying these contact interactions are fitted to the 2013 Granada database in the laboratory-energy range 0–300 MeV.

We describe our method for deriving the nuclear electromagnetic charge and current operators in chiral perturbation theory, based on time-ordered perturbation theory. We then discuss possible strategies for fixing the relevant low-energy constants, from the magnetic moments of the deuteron and of the trinucleons, and from the radiative np capture cross sections, and identify a scheme which, partly relying on {Delta} resonance saturation, leads to a reasonable pattern of convergence of the chiral expansion.